1 //--------------------------------------------------------------------------
4 // This software is part of the EvtGen package developed jointly
5 // for the BaBar and CLEO collaborations. If you use all or part
6 // of it, please give an appropriate acknowledgement.
8 // Copyright Information: See EvtGen/COPYRIGHT
9 // Copyright (C) 1998 Caltech, UCSB
11 // Module: EvtSSSCPT.cc
13 // Description: Routine to decay scalar -> 2 scalars (CPT)
15 // Modification history:
17 // SHY April 28, 1997 Module created
19 //------------------------------------------------------------------------
21 #include "EvtGenBase/EvtPatches.hh"
23 #include "EvtGenBase/EvtParticle.hh"
24 #include "EvtGenBase/EvtGenKine.hh"
25 #include "EvtGenBase/EvtCPUtil.hh"
26 #include "EvtGenBase/EvtPDL.hh"
27 #include "EvtGenBase/EvtReport.hh"
28 #include "EvtGenModels/EvtSSSCPT.hh"
29 #include "EvtGenBase/EvtId.hh"
31 #include "EvtGenBase/EvtConst.hh"
33 EvtSSSCPT::~EvtSSSCPT() {}
35 std::string EvtSSSCPT::getName(){
42 EvtDecayBase* EvtSSSCPT::clone(){
48 void EvtSSSCPT::init(){
50 // check that there are 8 arguments
57 void EvtSSSCPT::decay( EvtParticle *p){
59 //added by Lange Jan4,2000
60 static EvtId B0=EvtPDL::getId("B0");
61 static EvtId B0B=EvtPDL::getId("anti-B0");
67 EvtCPUtil::OtherB(p,t,other_b);
69 p->initializePhaseSpace(getNDaug(),getDaugs());
76 P=EvtComplex(cos(-getArg(0)),sin(-getArg(0)));
77 Q=EvtComplex(cos(getArg(0)),sin(getArg(0)));
78 D=EvtComplex(getArg(6)*cos(getArg(7)),getArg(6)*sin(getArg(7)));
79 Im=EvtComplex(0.0,1.0);
81 A=EvtComplex(getArg(2)*cos(getArg(3)),getArg(2)*sin(getArg(3)));
82 Abar=EvtComplex(getArg(4)*cos(getArg(5)),getArg(4)*sin(getArg(5)));
85 amp=A*cos(getArg(1)*t/(2*EvtConst::c))+
86 Im*sin(getArg(1)*t/(2*EvtConst::c))*
90 amp=Abar*cos(getArg(1)*t/(2*EvtConst::c))+
91 Im*sin(getArg(1)*t/(2*EvtConst::c))*