CGEM BOSS 6.6.5.h
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | List of all members
EvtKstarstargamma Class Reference

#include <EvtKstarstargamma.hh>

+ Inheritance diagram for EvtKstarstargamma:

Public Member Functions

 EvtKstarstargamma ()
 
virtual ~EvtKstarstargamma ()
 
void getName (std::string &name)
 
EvtDecayBase * clone ()
 
void decay (EvtParticle *p)
 
void init ()
 
void initProbMax ()
 
- 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 ()
 
EvtId * getDaugs ()
 
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 28 of file EvtKstarstargamma.hh.

Constructor & Destructor Documentation

◆ EvtKstarstargamma()

EvtKstarstargamma::EvtKstarstargamma ( )
inline

Definition at line 32 of file EvtKstarstargamma.hh.

32{}

Referenced by clone().

◆ ~EvtKstarstargamma()

EvtKstarstargamma::~EvtKstarstargamma ( )
virtual

Definition at line 39 of file EvtKstarstargamma.cc.

39{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtKstarstargamma::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 48 of file EvtKstarstargamma.cc.

48 {
49
50 return new EvtKstarstargamma;
51
52}

◆ decay()

void EvtKstarstargamma::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 75 of file EvtKstarstargamma.cc.

75 {
76
77 p->initializePhaseSpace(getNDaug(),getDaugs());
78
79 EvtParticle* kaon = p->getDaug(0);
80 EvtParticle* pion = p->getDaug(1);
81 EvtParticle* photon = p->getDaug(2);
82
83
84 EvtComplexPtrPtr Hd1=new EvtComplexPtr[5];
85 Hd1[0]=new EvtComplex[2];
86 Hd1[1]=new EvtComplex[2];
87 Hd1[2]=new EvtComplex[2];
88 Hd1[3]=new EvtComplex[2];
89 Hd1[4]=new EvtComplex[2];
90
91 Hd1[0][0]=0.0;
92 Hd1[0][1]=0.0;
93 Hd1[1][0]=0.0;
94 Hd1[1][1]=0.0;
95 Hd1[2][0]=0.0;
96 Hd1[2][1]=0.0;
97 Hd1[3][0]=0.0;
98 Hd1[3][1]=1.0;
99 Hd1[4][0]=0.0;
100 Hd1[4][1]=0.0;
101
102 EvtEvalHelAmp d1(EvtSpinType::SCALAR,EvtSpinType::TENSOR,
103 EvtSpinType::PHOTON,Hd1);
104
105 EvtScalarParticle theB;
106
107 theB.init(p->getId(),p->getP4Restframe());
108
109 EvtVector4R theKstarP4=kaon->getP4()+pion->getP4();
110
111 EvtTensorParticle theKstar;
112 theKstar.init(EvtPDL::getId(std::string("K_2*0")),theKstarP4);
113
114 EvtPhotonParticle thePhoton;
115 thePhoton.init(EvtPDL::getId(std::string("K_2*0")),photon->getP4());
116
117 theKstar.addDaug(&theB);
118 thePhoton.addDaug(&theB);
119
120 EvtAmp amp1;
121
122 d1.evalAmp(&theB,amp1);
123
124 EvtComplexPtrPtr Hd2=new EvtComplexPtr[1];
125 Hd2[0]=new EvtComplex[1];
126
127 Hd2[0][0]=1.0;
128
129
130 EvtEvalHelAmp d2(EvtSpinType::TENSOR,EvtSpinType::SCALAR,
131 EvtSpinType::SCALAR,Hd2);
132
133
134 EvtVector4R theKstarP4boost(theKstarP4.get(0),-theKstarP4.get(1),-theKstarP4.get(2),-theKstarP4.get(3));
135
136 EvtScalarParticle theKaon;
137 theKaon.init(EvtPDL::getId(std::string("K+")),boostTo(kaon->getP4(),theKstarP4boost));
138
139 EvtScalarParticle thePion;
140 thePion.init(EvtPDL::getId(std::string("pi+")),boostTo(pion->getP4(),theKstarP4boost));
141
142 theKaon.addDaug(&theKstar);
143 thePion.addDaug(&theKstar);
144
145 // Calculate the propagator
146
147 double m = theKstarP4.mass();
148 EvtTwoBodyVertex v(0.5,0.14,1.4,2);
149 EvtTwoBodyKine v1(0.5,0.14,m);
150 EvtPropBreitWignerRel prop(1.4,0.2);
151
152 // Mass-dependent width correction and amplitude calculation
153
154 double width = prop.g0() * v.widthFactor(v1);
155 prop.set_g0(width);
156 EvtComplex bwamp = prop.evaluate(m);
157
158
159 EvtAmp amp2;
160
161 d2.evalAmp(&theKstar,amp2);
162
163 vertex(0,bwamp*(amp1._amp[0]*amp2._amp[0]+
164 amp1._amp[1]*amp2._amp[1]+
165 amp1._amp[2]*amp2._amp[2]+
166 amp1._amp[3]*amp2._amp[3]+
167 amp1._amp[4]*amp2._amp[4]));
168
169 vertex(1,bwamp*(amp1._amp[5]*amp2._amp[0]+
170 amp1._amp[6]*amp2._amp[1]+
171 amp1._amp[7]*amp2._amp[2]+
172 amp1._amp[8]*amp2._amp[3]+
173 amp1._amp[9]*amp2._amp[4]));
174
175 return;
176}
pion
@ pion
Definition DstMdcDedx.h:9
kaon
@ kaon
Definition DstMdcDedx.h:9
boostTo
EvtDiracSpinor boostTo(const EvtDiracSpinor &sp, const EvtVector4R p4)
Definition EvtDiracSpinor.cc:75
v
**********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
EvtDecayAmp::vertex
void vertex(const EvtComplex &amp)
Definition EvtDecayAmp.hh:37
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition EvtDecayBase.hh:65
EvtPDL::getId
static EvtId getId(const std::string &name)
Definition EvtPDL.cc:287
EvtParticle::getId
EvtId getId() const
Definition EvtParticle.cc:113
EvtParticle::getP4Restframe
EvtVector4R getP4Restframe()
Definition EvtParticle.cc:698
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition EvtParticle.cc:121
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition EvtParticle.cc:85
EvtParticle::addDaug
void addDaug(EvtParticle *node)
Definition EvtParticle.cc:104
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition EvtParticle.cc:1069
EvtPhotonParticle::init
void init(EvtId part_n, double e, double px, double py, double pz)
Definition EvtPhotonParticle.cc:39
EvtScalarParticle::init
void init(EvtId part_n, double e, double px, double py, double pz)
Definition EvtScalarParticle.cc:29
EvtTensorParticle::init
void init(EvtId part_n, double e, double px, double py, double pz)
Definition EvtTensorParticle.cc:44
EvtVector4R::mass
double mass() const
Definition EvtVector4R.cc:39
EvtVector4R::get
double get(int i) const
Definition EvtVector4R.hh:179

◆ getName()

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

Implements EvtDecayBase.

Definition at line 41 of file EvtKstarstargamma.cc.

41 {
42
43 model_name="KSTARSTARGAMMA";
44
45}

◆ init()

void EvtKstarstargamma::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 54 of file EvtKstarstargamma.cc.

54 {
55
56 // check that there are 0 arguments
57 checkNArg(0);
58
59 // check that there are 3 daughters
60 checkNDaug(3);
61
62 // check the parent and daughter spins
63 checkSpinParent(EvtSpinType::SCALAR);
64 checkSpinDaughter(0,EvtSpinType::SCALAR);
65 checkSpinDaughter(1,EvtSpinType::SCALAR);
66 checkSpinDaughter(2,EvtSpinType::PHOTON);
67}
EvtDecayBase::checkSpinDaughter
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:533
EvtDecayBase::checkSpinParent
void checkSpinParent(EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:520
EvtDecayBase::checkNDaug
void checkNDaug(int d1, int d2=-1)
Definition EvtDecayBase.cc:504
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition EvtDecayBase.cc:482

◆ initProbMax()

void EvtKstarstargamma::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 69 of file EvtKstarstargamma.cc.

69 {
70
71 //setProbMax(1.0);
72
73}
The documentation for this class was generated from the following files: