CGEM BOSS 6.6.5.h
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtSVSCPLH Class Reference

#include <EvtSVSCPLH.hh>

+ Inheritance diagram for EvtSVSCPLH:

Public Member Functions

 EvtSVSCPLH ()
 
virtual ~EvtSVSCPLH ()
 
void getName (std::string &name)
 
EvtDecayBaseclone ()
 
void initProbMax ()
 
void init ()
 
void decay (EvtParticle *p)
 
- Public Member Functions inherited from EvtDecayAmp
void makeDecay (EvtParticle *p)
 
void setWeight (double weight)
 
void vertex (const EvtComplex &amp)
 
void vertex (int i1, const EvtComplex &amp)
 
void vertex (int i1, int i2, const EvtComplex &amp)
 
void vertex (int i1, int i2, int i3, const EvtComplex &amp)
 
void vertex (int *i1, const EvtComplex &amp)
 
virtual ~EvtDecayAmp ()
 
- Public Member Functions inherited from EvtDecayBase
virtual std::string commandName ()
 
virtual void command (std::string cmd)
 
double getProbMax (double prob)
 
double resetProbMax (double prob)
 
 EvtDecayBase ()
 
virtual ~EvtDecayBase ()
 
virtual bool matchingDecay (const EvtDecayBase &other) const
 
EvtId getParentId ()
 
double getBranchingFraction ()
 
void disableCheckQ ()
 
void checkQ ()
 
int getNDaug ()
 
EvtIdgetDaugs ()
 
EvtId getDaug (int i)
 
int getNArg ()
 
int getPHOTOS ()
 
void setPHOTOS ()
 
void setVerbose ()
 
void setSummary ()
 
double * getArgs ()
 
std::string * getArgsStr ()
 
double getArg (int j)
 
std::string getArgStr (int j)
 
std::string getModelName ()
 
int getDSum ()
 
int summary ()
 
int verbose ()
 
void saveDecayInfo (EvtId ipar, int ndaug, EvtId *daug, int narg, std::vector< std::string > &args, std::string name, double brfr)
 
void printSummary ()
 
void setProbMax (double prbmx)
 
void noProbMax ()
 
void checkNArg (int a1, int a2=-1, int a3=-1, int a4=-1)
 
void checkNDaug (int d1, int d2=-1)
 
void checkSpinParent (EvtSpinType::spintype sp)
 
void checkSpinDaughter (int d1, EvtSpinType::spintype sp)
 
virtual int nRealDaughters ()
 

Additional Inherited Members

- Static Public Member Functions inherited from EvtDecayBase
static void findMasses (EvtParticle *p, int ndaugs, EvtId daugs[10], double masses[10])
 
static void findMass (EvtParticle *p)
 
static double findMaxMass (EvtParticle *p)
 
- Protected Member Functions inherited from EvtDecayBase
bool daugsDecayedByParentModel ()
 
- Protected Attributes inherited from EvtDecayAmp
EvtAmp _amp2
 
- Protected Attributes inherited from EvtDecayBase
bool _daugsDecayedByParentModel
 

Detailed Description

Definition at line 29 of file EvtSVSCPLH.hh.

Constructor & Destructor Documentation

◆ EvtSVSCPLH()

EvtSVSCPLH::EvtSVSCPLH ( )
inline

Definition at line 33 of file EvtSVSCPLH.hh.

33{}

Referenced by clone().

◆ ~EvtSVSCPLH()

EvtSVSCPLH::~EvtSVSCPLH ( )
virtual

Definition at line 37 of file EvtSVSCPLH.cc.

37{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtSVSCPLH::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 46 of file EvtSVSCPLH.cc.

46 {
47
48 return new EvtSVSCPLH;
49
50}

◆ decay()

void EvtSVSCPLH::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 96 of file EvtSVSCPLH.cc.

96 {
97
99
100 static EvtId B0=EvtPDL::getId("B0");
101 static EvtId B0B=EvtPDL::getId("anti-B0");
102
103 double t;
104 EvtId other_b;
105
106 EvtCPUtil::OtherB(p,t,other_b,0.5);
107
108 //convert time from mm to seconds
109 t/=EvtConst::c;
110
111 //sign convention is dm=Mheavy-Mlight
112 // dGamma=Gammalight-Gammaheavy
113 //such that in the standard model both of these are positive.
114 EvtComplex gp=0.5*(exp(EvtComplex(0.25*t*_dgamma,-0.5*t*_dm))+exp(EvtComplex(-0.25*t*_dgamma,0.5*t*_dm)));
115 EvtComplex gm=0.5*(exp(EvtComplex(0.25*t*_dgamma,-0.5*t*_dm))-exp(EvtComplex(-0.25*t*_dgamma,0.5*t*_dm)));
116
117 EvtComplex amp;
118
119 if (other_b==B0B){
120 amp=gp*_Af+_qop*gm*_Abarf;
121 }
122 else if (other_b==B0){
123 amp=gp*_Abarf+_poq*gm*_Af;
124 }
125 else{
126 report(ERROR,"EvtGen") << "other_b was not B0 or B0B!"<<endl;
127 ::abort();
128 }
129
130 EvtVector4R p4_parent=p->getP4Restframe();;
131
132 double norm=p->getDaug(0)->mass()/(p->getDaug(0)->getP4().d3mag()*p4_parent.mass());
133
134 EvtParticle* v=p->getDaug(0);
135
136 vertex(0,amp*norm*(p4_parent*(v->epsParent(0))));
137 vertex(1,amp*norm*(p4_parent*(v->epsParent(1))));
138 vertex(2,amp*norm*(p4_parent*(v->epsParent(2))));
139
140
141 return ;
142}
EvtComplex exp(const EvtComplex &c)
ostream & report(Severity severity, const char *facility)
Definition EvtReport.cc:36
@ ERROR
Definition EvtReport.hh:49
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyWtm INTEGER m_idyfs DOUBLE PRECISION m_zini DOUBLE PRECISION m_q2 DOUBLE PRECISION m_Wt_KF DOUBLE PRECISION m_WtCut INTEGER m_KFfin *COMMON c_KarLud $ !Input CMS energy[GeV] $ !CMS energy after beam spread beam strahlung[GeV] $ !Beam energy spread[GeV] $ !z boost due to beam spread $ !electron beam mass *ff pair spectrum $ !minimum v
Definition KarLud.h:35
HepLorentzVector getP4(RecEmcShower *gtrk)
Definition MyUtil.cxx:16
static void OtherB(EvtParticle *p, double &t, EvtId &otherb)
Definition EvtCPUtil.cc:229
static const double c
Definition EvtConst.hh:32
void vertex(const EvtComplex &amp)
EvtId * getDaugs()
Definition EvtId.hh:27
static EvtId getId(const std::string &name)
Definition EvtPDL.cc:287
EvtVector4R getP4Restframe()
EvtParticle * getDaug(int i)
double mass() const
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
double mass() const
int t()
Definition t.c:1

◆ getName()

void EvtSVSCPLH::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 39 of file EvtSVSCPLH.cc.

39 {
40
41 model_name="SVS_CPLH";
42
43}

◆ init()

void EvtSVSCPLH::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 52 of file EvtSVSCPLH.cc.

52 {
53
54 // check that there are 8 arguments
55 checkNArg(8);
56 checkNDaug(2);
57
59
62
63 static double ctau=EvtPDL::getctau(EvtPDL::getId("B0"));
64
65 // hbar/s
66 _dm=getArg(0);
67 _dgamma=EvtConst::c*getArg(1)/ctau;
68
69 _qop=getArg(2)*exp(EvtComplex(0.0,getArg(3)));
70
71 _poq=1.0/_qop;
72
73 _Af=getArg(4)*exp(EvtComplex(0.0,getArg(5)));
74 _Abarf=getArg(6)*exp(EvtComplex(0.0,getArg(7)));
75
76 if (verbose()){
77 report(INFO,"EvtGen")<<":EvtSVSCPLH:dm="<<_dm<<endl;
78 report(INFO,"EvtGen")<<":EvtSVSCPLH:dGamma="<<_dgamma<<endl;
79 report(INFO,"EvtGen")<<":EvtSVSCPLH:q/p="<<_qop<<endl;
80 report(INFO,"EvtGen")<<":EvtSVSCPLH:Af="<<_Af<<endl;
81 report(INFO,"EvtGen")<<":EvtSVSCPLH:Abarf="<<_Abarf<<endl;
82 }
83
84
85}
@ INFO
Definition EvtReport.hh:52
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
void checkSpinParent(EvtSpinType::spintype sp)
double getArg(int j)
void checkNDaug(int d1, int d2=-1)
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
static double getctau(EvtId i)
Definition EvtPDL.hh:55

◆ initProbMax()

void EvtSVSCPLH::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 87 of file EvtSVSCPLH.cc.

87 {
88
89 //This is probably not quite right, but it should do as a start...
90 //Anders
91
92 setProbMax(4.0*(getArg(4)*getArg(4)+getArg(6)*getArg(6)));
93
94}
void setProbMax(double prbmx)

The documentation for this class was generated from the following files: