EvtCPUtil.hh

//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of the EvtGen package developed jointly
//      for the BaBar and CLEO collaborations.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/COPYRIGHT
//      Copyright (C) 1998      Caltech, UCSB
//
// Module: EvtGen/EvtCPUtil.hh
//
// Description:Class to hold CP physics utilities.
//
// Modification history:
//
//    RYD     March 24, 1998         Module created
//
//------------------------------------------------------------------------

#ifndef EVTCPUTIL_HH
#define EVTCPUTIL_HH

#include "EvtGenBase/EvtComplex.hh"
#include "EvtGenBase/EvtPatches.hh"
class EvtParticle;
class EvtId;

class EvtCPUtil {

public:

  EvtCPUtil(int mixingType);
  ~EvtCPUtil();

  enum MixingType {Coherent = 0, Incoherent = 1};

  static EvtCPUtil* getInstance();

  void setMixingType(int mixingType) {_mixingType = mixingType;}
  int getMixingType() {return _mixingType;}

  void fractB0CP(EvtComplex Af, EvtComplex Abarf, 
                        double deltam, double beta, double &fract);

  void fractB0nonCP(EvtComplex Af, EvtComplex Abarf, 
                    EvtComplex Afbar, EvtComplex Abarfbar, 
                    double deltam, double beta, int flip, 
                    double &fract);

  // Mark Whitehead 7/12/2009
  // Add required lines from EvtIncoherentMixing.hh to fix CPV

  // Functions to check if a B has mixed (comes from a B)
  bool isB0Mixed( EvtParticle * ) ;
  bool isBsMixed( EvtParticle * ) ;

  bool flipIsEnabled() ;
  void enableFlip() ;
  void disableFlip() ;

  void OtherB(EvtParticle *p, double &t, EvtId &otherb);

  void OtherCoherentB(EvtParticle *p, double &t, EvtId &otherb, double probB0);
  void OtherIncoherentB(EvtParticle *p, double &t, EvtId &otherb, double probB0);

  void OtherB(EvtParticle *p, double &t, EvtId &otherb, double probB0);

  //id is the produced particle
  //t returns the lifetime of the particle
  //and mix will be 1 if it mixed otherwise 0
  void incoherentMix(const EvtId id, double &t, int &mix);

  double getDeltaGamma(const EvtId id);
  double getDeltaM(const EvtId id);

private:

  bool _enableFlip;
  int _mixingType;

};


#endif