Page MenuHomeHEPForge
Files F19244894

EvtBToDDalitzCPK.cpp
No OneTemporary
Actions

Size
4 KB
Referenced Files
None
Subscribers
None

EvtBToDDalitzCPK.cpp
View Options

// $Id: EvtBToDDalitzCPK.cpp,v 1.2 2009-03-16 16:28:42 robbep Exp $
// Include files
#include "EvtGenBase/EvtPDL.hh"
#include "EvtGenBase/EvtParticle.hh"
#include "EvtGenBase/EvtVector4C.hh"
// local
#include "EvtGenModels/EvtBToDDalitzCPK.hh"
//-----------------------------------------------------------------------------
// Implementation file for class : EvtBToDDalitzCPK
// Decay Model for B->D0K with D0->Ks pi+ pi-
// it is just a way to get the arguments...
// Works also for other B->D0K decays...
// 2003-12-08 : Patrick Robbe
//-----------------------------------------------------------------------------
//=============================================================================
// Standard destructor
//=============================================================================
EvtBToDDalitzCPK::~EvtBToDDalitzCPK ( ) { }
//=============================================================================
// Name of the model
//=============================================================================
std::string EvtBToDDalitzCPK::getName( )
{
return "BTODDALITZCPK" ;
}
//=============================================================================
// Clone method
//=============================================================================
EvtDecayBase * EvtBToDDalitzCPK::clone ( )
{
return new EvtBToDDalitzCPK ;
}
//=============================================================================
// Initialisation method
//=============================================================================
void EvtBToDDalitzCPK::init ( )
{
// Check that there are 3 arguments
checkNArg( 3 ) ;
// Check that there are 2 daughters
checkNDaug( 2 ) ;
// Check that the particles of the decay are :
// B+/- -> D0/bar K+/-
// B+/- -> K+/- D0/bar
// B0/bar -> K*0/bar D0/bar
// and nothing else ...
static EvtId BP = EvtPDL::getId( "B+" ) ;
static EvtId BM = EvtPDL::getId( "B-" ) ;
static EvtId B0 = EvtPDL::getId( "B0" ) ;
static EvtId B0B = EvtPDL::getId( "anti-B0" ) ;
static EvtId KP = EvtPDL::getId( "K+" ) ;
static EvtId KM = EvtPDL::getId( "K-" ) ;
static EvtId KS = EvtPDL::getId( "K*0" ) ;
static EvtId KSB = EvtPDL::getId( "anti-K*0" ) ;
static EvtId D0 = EvtPDL::getId( "D0" ) ;
static EvtId D0B = EvtPDL::getId( "anti-D0" ) ;
_flag = 0 ;
EvtId parent = getParentId() ;
EvtId d1 = getDaug( 0 ) ;
EvtId d2 = getDaug( 1 ) ;
if ( ( ( parent == BP ) || ( parent == BM ) ) &&
( ( d1 == D0 ) || ( d1 == D0B ) ) &&
( ( d2 == KP ) || ( d2 == KM ) ) ) {
_flag = 1 ;
// PHSP Decay
}
else if ( ( ( parent == BP ) || ( parent == BM ) ) &&
( ( d1 == KP ) || ( d1 == KM ) ) &&
( ( d2 == D0 ) || ( d2 == D0B ) ) ) {
_flag = 1 ;
// also PHSP decay
}
else if ( ( ( parent == B0 ) || ( parent == B0B ) ) &&
( ( d1 == KS ) || ( d1 == KSB ) ) &&
( ( d2 == D0 ) || ( d2 == D0B ) ) ) {
_flag = 2 ;
// SVS Decay
}
if ( _flag == 0 ) {
report( ERROR , "EvtGen" ) << "EvtBToDDalitzCPK : Invalid mode."
<< std::endl ;
assert( 0 ) ;
}
}
//=============================================================================
// Set prob max
//=============================================================================
void EvtBToDDalitzCPK::initProbMax( )
{
if ( _flag == 1 ) {
// PHSP
setProbMax ( 0. ) ;
}
else if ( _flag == 2 ) {
// SVS
setProbMax ( 1.0 ) ;
}
}
//=============================================================================
// decay particle
//=============================================================================
void EvtBToDDalitzCPK::decay( EvtParticle * p )
{
if ( _flag == 1 ) {
// PHSP
p -> initializePhaseSpace( getNDaug() , getDaugs() ) ;
vertex ( 0. ) ;
}
else if ( _flag == 2 ) {
// SVS
p->initializePhaseSpace(getNDaug(),getDaugs());
EvtParticle *v;
v = p->getDaug(0);
double massv = v->mass();
EvtVector4R momv = v->getP4();
EvtVector4R moms = p->getDaug(1)->getP4();
double m_parent = p->mass();
EvtVector4R p4_parent = momv+moms;
double norm=massv/(momv.d3mag()*m_parent);
p4_parent = norm*p4_parent;
vertex(0,p4_parent*(v->epsParent(0)));
vertex(1,p4_parent*(v->epsParent(1)));
vertex(2,p4_parent*(v->epsParent(2)));
}
}

File Metadata

Mime Type
text/x-c
Expires
Tue, Sep 30, 4:46 AM (11 h, 56 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
6494660
Default Alt Text
EvtBToDDalitzCPK.cpp (4 KB)

Event Timeline

HEPForge.org · .