EvtCPUtil.hh

Go to the documentation of this file.

 
/***********************************************************************
* Copyright 1998-2020 CERN for the benefit of the EvtGen authors       *
*                                                                      *
* This file is part of EvtGen.                                         *
*                                                                      *
* EvtGen is free software: you can redistribute it and/or modify       *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation, either version 3 of the License, or    *
* (at your option) any later version.                                  *
*                                                                      *
* EvtGen is distributed in the hope that it will be useful,            *
* but WITHOUT ANY WARRANTY; without even the implied warranty of       *
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the        *
* GNU General Public License for more details.                         *
*                                                                      *
* You should have received a copy of the GNU General Public License    *
* along with EvtGen.  If not, see <https://www.gnu.org/licenses/>.     *
***********************************************************************/
 
#ifndef EVTCPUTIL_HH
#define EVTCPUTIL_HH
 
#include "EvtGenBase/EvtComplex.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 ) { m_mixingType = mixingType; }
    int getMixingType() { return m_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 m_enableFlip;
    int m_mixingType;
};
 
#endif