BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtChi2BB2 Class Reference

#include <EvtChi2BB2.hh>

Inheritance diagram for EvtChi2BB2:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtChi2BB2()

EvtChi2BB2::EvtChi2BB2 ( )
inline

Definition at line 31 of file EvtChi2BB2.hh.

31{}

Referenced by clone().

◆ ~EvtChi2BB2()

EvtChi2BB2::~EvtChi2BB2 ( )
virtual

Definition at line 36 of file EvtChi2BB2.cc.

36{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtChi2BB2::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 40 of file EvtChi2BB2.cc.

40{ return new EvtChi2BB2; }

◆ decay()

void EvtChi2BB2::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 59 of file EvtChi2BB2.cc.

59 {
60
61 p->initializePhaseSpace( getNDaug(), getDaugs() );
62
63 EvtParticle *v, *s1;
64 EvtVector4R pv, ps, ppr;
65
66 v = p->getDaug( 0 );
67 s1 = p->getDaug( 1 );
68 pv = v->getP4();
69 ps = s1->getP4();
70 ppr = p->getP4();
71
72 // Put phase space results into the daughters.
73 EvtHelSys angles( ppr, pv ); // using helicity sys. angles
74 double theta = angles.getHelAng( 1 );
75 double phi = angles.getHelAng( 2 );
76 double gamma = 0;
77 double F00, F01, F02, F03, F10, F11, F13, F20, F22, F23, F30, F31, F32, F33;
78 F00 = getArg( 0 );
79 F01 = getArg( 1 );
80 F02 = getArg( 2 );
81 F03 = getArg( 3 );
82 F11 = getArg( 4 );
83
84 F33 = F00;
85 F32 = F01;
86 F31 = F02;
87 F30 = F03;
88 F22 = F11;
89 F10 = F01;
90 F20 = F02;
91 F23 = F10;
92 F13 = F20;
93 // cout<<"********F00,F01,F02,F03,F10,F11,F13,F20,F22,F23,F30,F31,F32,F33="<<F00<<F01<<F02<<F03<<F10<<F11<<F13<<F20<<F22<<F23<<F30<<F31<<F32<<F33<<endl;
94
95 vertex( 0, 0, 0, Djmn( 2, 0, 0, phi, theta, gamma ) * F00 ); // Chi2 helicity =0
96 // corresponding index 0;
97 vertex( 0, 0, 1, Djmn( 2, 0, -1, phi, theta, gamma ) * F01 );
98 vertex( 0, 0, 2, Djmn( 2, 0, 2, phi, theta, gamma ) * F02 );
99 vertex( 0, 0, 3, Djmn( 2, 0, 1, phi, theta, gamma ) * F03 );
100 vertex( 0, 1, 0, Djmn( 2, 0, 1, phi, theta, gamma ) * F10 );
101 vertex( 0, 1, 1, Djmn( 2, 0, 0, phi, theta, gamma ) * F11 );
102 vertex( 0, 1, 2, 0.0 );
103 vertex( 0, 1, 3, Djmn( 2, 0, 2, phi, theta, gamma ) * F13 );
104 vertex( 0, 2, 0, Djmn( 2, 0, -2, phi, theta, gamma ) * F20 );
105 vertex( 0, 2, 1, 0.0 );
106 vertex( 0, 2, 2, Djmn( 2, 0, 0, phi, theta, gamma ) * F22 );
107 vertex( 0, 2, 3, Djmn( 2, 0, -1, phi, theta, gamma ) * F23 );
108 vertex( 0, 3, 0, Djmn( 2, 0, -1, phi, theta, gamma ) * F30 );
109 vertex( 0, 3, 1, Djmn( 2, 0, -2, phi, theta, gamma ) * F31 );
110 vertex( 0, 3, 2, Djmn( 2, 0, 1, phi, theta, gamma ) * F32 );
111 vertex( 0, 3, 3, Djmn( 2, 0, 0, phi, theta, gamma ) * F33 );
112
113 vertex( 1, 0, 0, Djmn( 2, 2, 0, phi, theta, gamma ) * F00 ); // Chi2 helicity =2
114 // corresponding index 1;
115 vertex( 1, 0, 1, Djmn( 2, 2, -1, phi, theta, gamma ) * F01 );
116 vertex( 1, 0, 2, Djmn( 2, 2, 2, phi, theta, gamma ) * F02 );
117 vertex( 1, 0, 3, Djmn( 2, 2, 1, phi, theta, gamma ) * F03 );
118 vertex( 1, 1, 0, Djmn( 2, 2, 1, phi, theta, gamma ) * F10 );
119 vertex( 1, 1, 1, Djmn( 2, 2, 0, phi, theta, gamma ) * F11 );
120 vertex( 1, 1, 2, 0.0 );
121 vertex( 1, 1, 3, Djmn( 2, 2, 2, phi, theta, gamma ) * F13 );
122 vertex( 1, 2, 0, Djmn( 2, 2, -2, phi, theta, gamma ) * F20 );
123 vertex( 1, 2, 1, 0.0 );
124 vertex( 1, 2, 2, Djmn( 2, 2, 0, phi, theta, gamma ) * F22 );
125 vertex( 1, 2, 3, Djmn( 2, 2, -1, phi, theta, gamma ) * F23 );
126 vertex( 1, 3, 0, Djmn( 2, 2, -1, phi, theta, gamma ) * F30 );
127 vertex( 1, 3, 1, Djmn( 2, 2, -2, phi, theta, gamma ) * F31 );
128 vertex( 1, 3, 2, Djmn( 2, 2, 1, phi, theta, gamma ) * F32 );
129 vertex( 1, 3, 3, Djmn( 2, 2, 0, phi, theta, gamma ) * F33 );
130
131 vertex( 2, 0, 0, Djmn( 2, -2, 0, phi, theta, gamma ) * F00 ); // Chi2 helicity =-2
132 // corresponding index 2;
133 vertex( 2, 0, 1, Djmn( 2, -2, -1, phi, theta, gamma ) * F01 );
134 vertex( 2, 0, 2, Djmn( 2, -2, 2, phi, theta, gamma ) * F02 );
135 vertex( 2, 0, 3, Djmn( 2, -2, 1, phi, theta, gamma ) * F03 );
136 vertex( 2, 1, 0, Djmn( 2, -2, 1, phi, theta, gamma ) * F10 );
137 vertex( 2, 1, 1, Djmn( 2, -2, 0, phi, theta, gamma ) * F11 );
138 vertex( 2, 1, 2, 0.0 );
139 vertex( 2, 1, 3, Djmn( 2, -2, 2, phi, theta, gamma ) * F13 );
140 vertex( 2, 2, 0, Djmn( 2, -2, -2, phi, theta, gamma ) * F20 );
141 vertex( 2, 2, 1, 0.0 );
142 vertex( 2, 2, 2, Djmn( 2, -2, 0, phi, theta, gamma ) * F22 );
143 vertex( 2, 2, 3, Djmn( 2, -2, -1, phi, theta, gamma ) * F23 );
144 vertex( 2, 3, 0, Djmn( 2, -2, -1, phi, theta, gamma ) * F30 );
145 vertex( 2, 3, 1, Djmn( 2, -2, -2, phi, theta, gamma ) * F31 );
146 vertex( 2, 3, 2, Djmn( 2, -2, 1, phi, theta, gamma ) * F32 );
147 vertex( 2, 3, 3, Djmn( 2, -2, 0, phi, theta, gamma ) * F33 );
148
149 vertex( 3, 0, 0, Djmn( 2, 1, 0, phi, theta, gamma ) * F00 ); // Chi2 helicity =1
150 // corresponding index 3;
151 vertex( 3, 0, 1, Djmn( 2, 1, -1, phi, theta, gamma ) * F01 );
152 vertex( 3, 0, 2, Djmn( 2, 1, 2, phi, theta, gamma ) * F02 );
153 vertex( 3, 0, 3, Djmn( 2, 1, 1, phi, theta, gamma ) * F03 );
154 vertex( 3, 1, 0, Djmn( 2, 1, 1, phi, theta, gamma ) * F10 );
155 vertex( 3, 1, 1, Djmn( 2, 1, 0, phi, theta, gamma ) * F11 );
156 vertex( 3, 1, 2, 0.0 );
157 vertex( 3, 1, 3, Djmn( 2, 1, 2, phi, theta, gamma ) * F13 );
158 vertex( 3, 2, 0, Djmn( 2, 1, -2, phi, theta, gamma ) * F20 );
159 vertex( 3, 2, 1, 0.0 );
160 vertex( 3, 2, 2, Djmn( 2, 1, 0, phi, theta, gamma ) * F22 );
161 vertex( 3, 2, 3, Djmn( 2, 1, -1, phi, theta, gamma ) * F23 );
162 vertex( 3, 3, 0, Djmn( 2, 1, -1, phi, theta, gamma ) * F30 );
163 vertex( 3, 3, 1, Djmn( 2, 1, -2, phi, theta, gamma ) * F31 );
164 vertex( 3, 3, 2, Djmn( 2, 1, 1, phi, theta, gamma ) * F32 );
165 vertex( 3, 3, 3, Djmn( 2, 1, 0, phi, theta, gamma ) * F33 );
166
167 vertex( 4, 0, 0, Djmn( 2, -1, 0, phi, theta, gamma ) * F00 ); // Chi2 helicity =-1
168 // corresponding index 4;
169 vertex( 4, 0, 1, Djmn( 2, -1, -1, phi, theta, gamma ) * F01 );
170 vertex( 4, 0, 2, Djmn( 2, -1, 2, phi, theta, gamma ) * F02 );
171 vertex( 4, 0, 3, Djmn( 2, -1, 1, phi, theta, gamma ) * F03 );
172 vertex( 4, 1, 0, Djmn( 2, -1, 1, phi, theta, gamma ) * F10 );
173 vertex( 4, 1, 1, Djmn( 2, -1, 0, phi, theta, gamma ) * F11 );
174 vertex( 4, 1, 2, 0.0 );
175 vertex( 4, 1, 3, Djmn( 2, -1, 2, phi, theta, gamma ) * F13 );
176 vertex( 4, 2, 0, Djmn( 2, -1, -2, phi, theta, gamma ) * F20 );
177 vertex( 4, 2, 1, 0.0 );
178 vertex( 4, 2, 2, Djmn( 2, -1, 0, phi, theta, gamma ) * F22 );
179 vertex( 4, 2, 3, Djmn( 2, -1, -1, phi, theta, gamma ) * F23 );
180 vertex( 4, 3, 0, Djmn( 2, -1, -1, phi, theta, gamma ) * F30 );
181 vertex( 4, 3, 1, Djmn( 2, -1, -2, phi, theta, gamma ) * F31 );
182 vertex( 4, 3, 2, Djmn( 2, -1, 1, phi, theta, gamma ) * F32 );
183 vertex( 4, 3, 3, Djmn( 2, -1, 0, phi, theta, gamma ) * F33 );
184
185 return;
186}
Djmn
EvtComplex Djmn(int j, int m, int n, double phi, double theta, double gamma)
Definition EvtHelSys.cc:165
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:36
EvtDecayBase::getArg
double getArg(int j)
Definition EvtDecayBase.cc:520
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition EvtDecayBase.hh:68
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition EvtParticle.cc:108
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition EvtParticle.cc:80
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition EvtParticle.cc:1042

◆ getName()

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

Implements EvtDecayBase.

Definition at line 38 of file EvtChi2BB2.cc.

38{ model_name = "Chi2BB2"; }

◆ init()

void EvtChi2BB2::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 42 of file EvtChi2BB2.cc.

42 {
43 // check that there are 0 arguments
44 checkNArg( 5 );
45 checkNDaug( 2 );
46 checkSpinParent( EvtSpinType::TENSOR );
47 checkSpinDaughter( 0, EvtSpinType::RARITASCHWINGER );
48 checkSpinDaughter( 1, EvtSpinType::RARITASCHWINGER );
49}
EvtDecayBase::checkSpinDaughter
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:492
EvtDecayBase::checkSpinParent
void checkSpinParent(EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:479
EvtDecayBase::checkNDaug
void checkNDaug(int d1, int d2=-1)
Definition EvtDecayBase.cc:465
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition EvtDecayBase.cc:449
The documentation for this class was generated from the following files: