Page MenuHomeHEPForge
Files F19252186

D137.1759122900.diff
No OneTemporary
Actions

Size
15 KB
Referenced Files
None
Subscribers
None

D137.1759122900.diff
View Options

diff --git a/EvtGenBase/EvtCPUtil.hh b/EvtGenBase/EvtCPUtil.hh
--- a/EvtGenBase/EvtCPUtil.hh
+++ b/EvtGenBase/EvtCPUtil.hh
@@ -48,9 +48,6 @@
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* );
diff --git a/EvtGenBase/EvtIncoherentMixing.hh b/EvtGenBase/EvtIncoherentMixing.hh
deleted file mode 100644
--- a/EvtGenBase/EvtIncoherentMixing.hh
+++ /dev/null
@@ -1,93 +0,0 @@
-
-/***********************************************************************
-* 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 EVTINCOHERENTMIXING_HH
-#define EVTINCOHERENTMIXING_HH 1
-
-// Include files
-
-#include "EvtGenBase/EvtParticle.hh"
-
-/** @class EvtIncoherentMixing EvtIncoherentMixing.hh EvtGenBase/EvtIncoherentMixing.hh
- * Class to contain the parameters of the incoherent B0 and B0S mixing
- *
- * @author Patrick Robbe
- * @date 2003-10-09
- */
-class EvtIncoherentMixing {
- public:
- /// Standard constructor
- EvtIncoherentMixing();
-
- ~EvtIncoherentMixing(); ///< Destructor
-
- // activate or desactivate the Bs mixing
- static void setB0Mixing();
- static void unsetB0Mixing();
-
- // activate or desactivate the B0 mixing
- static void setBsMixing();
- static void unsetBsMixing();
-
- // is mixing activated ?
- static bool doB0Mixing();
- static bool doBsMixing();
-
- // set values for the mixing
- static void setdGammad( double value );
- static void setdeltamd( double value );
- static void setdGammas( double value );
- static void setdeltams( double value );
-
- // get parameters for mixing
- static double getdGammad();
- static double getdeltamd();
- static double getdGammas();
- static double getdeltams();
-
- // Functions to obtain time and type of the mixing ( 1 mix, 0 unmix )
- static void incoherentB0Mix( const EvtId id, double& t, int& mix );
- static void incoherentBsMix( const EvtId id, double& t, int& mix );
-
- // Functions to check if a B has mixed (comes from a B)
- static bool isB0Mixed( EvtParticle* );
- static bool isBsMixed( EvtParticle* );
-
- // Functions for CP models
- // returns for particle p the decay time t , the flavour of the tag
- // side (B0, anti-B0, B_s0, anti-B_s0). The flavour of p can flip
- // with probability probB
- static void OtherB( EvtParticle* p, double& t, EvtId& otherb, double probB );
- static void OtherB( EvtParticle* p, double& t, EvtId& otherb );
-
- static bool flipIsEnabled();
- static void enableFlip();
- static void disableFlip();
-
- private:
- static bool m_doB0Mixing;
- static bool m_doBsMixing;
- static double m_dGammad;
- static double m_deltamd;
- static double m_dGammas;
- static double m_deltams;
- static bool m_enableFlip;
-};
-#endif // EVTGENBASE_EVTINCOHERENTMIXING_HH
diff --git a/src/EvtGenBase/EvtIncoherentMixing.cpp b/src/EvtGenBase/EvtIncoherentMixing.cpp
deleted file mode 100644
--- a/src/EvtGenBase/EvtIncoherentMixing.cpp
+++ /dev/null
@@ -1,326 +0,0 @@
-
-/***********************************************************************
-* 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/>. *
-***********************************************************************/
-
-#include "EvtGenBase/EvtIncoherentMixing.hh"
-
-#include "EvtGenBase/EvtId.hh"
-#include "EvtGenBase/EvtPDL.hh"
-#include "EvtGenBase/EvtRandom.hh"
-
-#include <cmath>
-#include <cstdlib>
-#include <iostream>
-
-//-----------------------------------------------------------------------------
-// Implementation file for class : EvtIncoherentMixing
-//
-// 2003-10-09 : Patrick Robbe
-//-----------------------------------------------------------------------------
-
-bool EvtIncoherentMixing::m_doB0Mixing = false;
-bool EvtIncoherentMixing::m_doBsMixing = false;
-bool EvtIncoherentMixing::m_enableFlip = false;
-double EvtIncoherentMixing::m_dGammad = 0.;
-double EvtIncoherentMixing::m_deltamd = 0.502e12;
-// dGamma_s corresponds to DeltaGamma / Gamma = 10 %
-double EvtIncoherentMixing::m_dGammas = 6.852e10;
-double EvtIncoherentMixing::m_deltams = 20.e12;
-
-//=============================================================================
-// Standard constructor, initializes variables
-//=============================================================================
-EvtIncoherentMixing::EvtIncoherentMixing()
-{
- m_doB0Mixing = false;
- m_doBsMixing = false;
- m_dGammad = 0.;
- // dGammas corresponds to DeltaGamma / Gamma = 10 %
- m_dGammas = 6.852e10;
- m_deltamd = 0.502e12;
- m_deltams = 20.e12;
- m_enableFlip = false;
-}
-//=============================================================================
-void EvtIncoherentMixing::incoherentB0Mix( const EvtId id, double& t, int& mix )
-{
- static const EvtId B0 = EvtPDL::getId( "B0" );
- static const EvtId B0B = EvtPDL::getId( "anti-B0" );
-
- if ( ( B0 != id ) && ( B0B != id ) ) {
- EvtGenReport( EVTGEN_ERROR, "EvtGen" )
- << "Bad configuration in incoherentB0Mix" << std::endl;
- ::abort();
- }
-
- double x = getdeltamd() * EvtPDL::getctau( B0 ) / EvtConst::c;
-
- double y = getdGammad() * ( EvtPDL::getctau( B0 ) / EvtConst::c ) / 2.;
-
- double fac = 1.; // No CP violation
-
- double mixprob = ( x * x + y * y ) /
- ( x * x + y * y + ( 1. / fac ) * ( 2. + x * x - y * y ) );
-
- int mixsign;
-
- // decide if state is mixed
- mixsign = ( mixprob > EvtRandom::Flat( 0., 1. ) ) ? -1 : 1;
-
- double prob;
-
- do {
- t = -log( EvtRandom::Flat() ) * EvtPDL::getctau( B0 ) / ( 1. - y );
- prob = ( 1. + exp( -2. * y * t / EvtPDL::getctau( B0 ) ) +
- mixsign * 2. * exp( -y * t / EvtPDL::getctau( B0 ) ) *
- cos( getdeltamd() * t / EvtConst::c ) ) /
- 2.;
- } while ( prob < 2. * EvtRandom::Flat() );
-
- mix = 0;
- if ( mixsign == -1 )
- mix = 1;
-
- return;
-}
-// ============================================================================
-void EvtIncoherentMixing::incoherentBsMix( const EvtId id, double& t, int& mix )
-{
- static const EvtId BS = EvtPDL::getId( "B_s0" );
- static const EvtId BSB = EvtPDL::getId( "anti-B_s0" );
-
- if ( ( BS != id ) && ( BSB != id ) ) {
- EvtGenReport( EVTGEN_ERROR, "EvtGen" )
- << "Bad configuration in incoherentBsMix" << std::endl;
- ::abort();
- }
-
- double x = getdeltams() * EvtPDL::getctau( BS ) / EvtConst::c;
-
- double y = getdGammas() * ( EvtPDL::getctau( BS ) / EvtConst::c ) / 2.;
-
- double fac = 1.; // No CP violation
-
- double mixprob = ( x * x + y * y ) /
- ( x * x + y * y + ( 1. / fac ) * ( 2. + x * x - y * y ) );
-
- int mixsign;
-
- // decide if state is mixed
- mixsign = ( mixprob > EvtRandom::Flat( 0., 1. ) ) ? -1 : 1;
-
- double prob;
-
- do {
- t = -log( EvtRandom::Flat() ) * EvtPDL::getctau( BS ) / ( 1. - y );
- prob = ( 1. + exp( -2. * y * t / EvtPDL::getctau( BS ) ) +
- mixsign * 2. * exp( -y * t / EvtPDL::getctau( BS ) ) *
- cos( getdeltams() * t / EvtConst::c ) ) /
- 2.;
- } while ( prob < 2. * EvtRandom::Flat() );
-
- mix = 0;
- if ( mixsign == -1 )
- mix = 1;
-
- return;
-}
-
-// ========================================================================
-bool EvtIncoherentMixing::isBsMixed( EvtParticle* p )
-{
- if ( !( p->getParent() ) )
- return false;
-
- static const EvtId BS0 = EvtPDL::getId( "B_s0" );
- static const EvtId BSB = EvtPDL::getId( "anti-B_s0" );
-
- if ( ( p->getId() != BS0 ) && ( p->getId() != BSB ) )
- return false;
-
- if ( ( p->getParent()->getId() == BS0 ) || ( p->getParent()->getId() == BSB ) )
- return true;
-
- return false;
-}
-
-// ========================================================================
-bool EvtIncoherentMixing::isB0Mixed( EvtParticle* p )
-{
- if ( !( p->getParent() ) )
- return false;
-
- static const EvtId B0 = EvtPDL::getId( "B0" );
- static const EvtId B0B = EvtPDL::getId( "anti-B0" );
-
- if ( ( p->getId() != B0 ) && ( p->getId() != B0B ) )
- return false;
-
- if ( ( p->getParent()->getId() == B0 ) || ( p->getParent()->getId() == B0B ) )
- return true;
-
- return false;
-}
-//============================================================================
-// Return the tag of the event (ie the anti-flavour of the produced
-// B meson). Flip the flavour of the event with probB probability
-//============================================================================
-void EvtIncoherentMixing::OtherB( EvtParticle* p, double& t, EvtId& otherb,
- double probB )
-{
- //if(p->getId() == B0 || p->getId() == B0B)
- //added by liming Zhang
- enableFlip();
- if ( ( isB0Mixed( p ) ) || ( isBsMixed( p ) ) ) {
- p->getParent()->setLifetime();
- t = p->getParent()->getLifetime();
- } else {
- p->setLifetime();
- t = p->getLifetime();
- }
-
- if ( flipIsEnabled() ) {
- //std::cout << " liming << flipIsEnabled " << std::endl;
- // Flip the flavour of the particle with probability probB
- bool isFlipped = ( EvtRandom::Flat( 0., 1. ) < probB );
-
- if ( isFlipped ) {
- if ( ( isB0Mixed( p ) ) || ( isBsMixed( p ) ) ) {
- p->getParent()->setId(
- EvtPDL::chargeConj( p->getParent()->getId() ) );
- p->setId( EvtPDL::chargeConj( p->getId() ) );
- } else {
- p->setId( EvtPDL::chargeConj( p->getId() ) );
- }
- }
- }
-
- if ( ( isB0Mixed( p ) ) || ( isBsMixed( p ) ) ) {
- // if B has mixed, tag flavour is charge conjugate of parent of B-meson
- otherb = EvtPDL::chargeConj( p->getParent()->getId() );
- } else {
- // else it is opposite flavour than this B hadron
- otherb = EvtPDL::chargeConj( p->getId() );
- }
-
- return;
-}
-//============================================================================
-// Return the tag of the event (ie the anti-flavour of the produced
-// B meson). No flip
-//============================================================================
-void EvtIncoherentMixing::OtherB( EvtParticle* p, double& t, EvtId& otherb )
-{
- if ( ( isB0Mixed( p ) ) || ( isBsMixed( p ) ) ) {
- p->getParent()->setLifetime();
- t = p->getParent()->getLifetime();
- } else {
- p->setLifetime();
- t = p->getLifetime();
- }
-
- if ( ( isB0Mixed( p ) ) || ( isBsMixed( p ) ) ) {
- // if B has mixed, tag flavour is charge conjugate of parent of B-meson
- otherb = EvtPDL::chargeConj( p->getParent()->getId() );
- } else {
- // else it is opposite flavour than this B hadron
- otherb = EvtPDL::chargeConj( p->getId() );
- }
-
- return;
-}
-
-// activate or desactivate the Bs mixing
-void EvtIncoherentMixing::setB0Mixing()
-{
- m_doB0Mixing = true;
-}
-void EvtIncoherentMixing::unsetB0Mixing()
-{
- m_doB0Mixing = false;
-}
-
-// activate or desactivate the B0 mixing
-void EvtIncoherentMixing::setBsMixing()
-{
- m_doBsMixing = true;
-}
-void EvtIncoherentMixing::unsetBsMixing()
-{
- m_doBsMixing = false;
-}
-
-// is mixing activated ?
-bool EvtIncoherentMixing::doB0Mixing()
-{
- return m_doB0Mixing;
-}
-bool EvtIncoherentMixing::doBsMixing()
-{
- return m_doBsMixing;
-}
-
-// set values for the mixing
-void EvtIncoherentMixing::setdGammad( double value )
-{
- m_dGammad = value;
-}
-void EvtIncoherentMixing::setdeltamd( double value )
-{
- m_deltamd = value;
-}
-void EvtIncoherentMixing::setdGammas( double value )
-{
- m_dGammas = value;
-}
-void EvtIncoherentMixing::setdeltams( double value )
-{
- m_deltams = value;
-}
-
-// get parameters for mixing
-double EvtIncoherentMixing::getdGammad()
-{
- return m_dGammad;
-}
-double EvtIncoherentMixing::getdeltamd()
-{
- return m_deltamd;
-}
-double EvtIncoherentMixing::getdGammas()
-{
- return m_dGammas;
-}
-double EvtIncoherentMixing::getdeltams()
-{
- return m_deltams;
-}
-
-bool EvtIncoherentMixing::flipIsEnabled()
-{
- return m_enableFlip;
-}
-void EvtIncoherentMixing::enableFlip()
-{
- m_enableFlip = true;
-}
-void EvtIncoherentMixing::disableFlip()
-{
- m_enableFlip = false;
-}
diff --git a/src/EvtGenModels/EvtSSDCP.cpp b/src/EvtGenModels/EvtSSDCP.cpp
--- a/src/EvtGenModels/EvtSSDCP.cpp
+++ b/src/EvtGenModels/EvtSSDCP.cpp
@@ -23,7 +23,6 @@
#include "EvtGenBase/EvtCPUtil.hh"
#include "EvtGenBase/EvtConst.hh"
#include "EvtGenBase/EvtGenKine.hh"
-#include "EvtGenBase/EvtIncoherentMixing.hh"
#include "EvtGenBase/EvtPDL.hh"
#include "EvtGenBase/EvtParticle.hh"
#include "EvtGenBase/EvtRandom.hh"
@@ -208,7 +207,6 @@
EvtComplex amp;
EvtCPUtil::getInstance()->OtherB( p, t, other_b, 0.5 ); // t is c*Dt (mm)
- // EvtIncoherentMixing::OtherB( p , t , other_b , 0.5 ) ;
//if (flip) t=-t;

File Metadata

Mime Type
text/plain
Expires
Mon, Sep 29, 6:15 AM (17 h, 36 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
6531231
Default Alt Text
D137.1759122900.diff (15 KB)

Event Timeline

HEPForge.org · .