CGEM BOSS 6.6.5.h
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)
 
EvtDecayBaseclone ()
 
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 ()
 
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 28 of file EvtChi2BB2.hh.

Constructor & Destructor Documentation

◆ EvtChi2BB2()

EvtChi2BB2::EvtChi2BB2 ( )
inline

Definition at line 32 of file EvtChi2BB2.hh.

32{}

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 45 of file EvtChi2BB2.cc.

45 {
46
47 return new EvtChi2BB2;
48
49}

◆ decay()

void EvtChi2BB2::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 69 of file EvtChi2BB2.cc.

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

◆ getName()

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

Implements EvtDecayBase.

Definition at line 38 of file EvtChi2BB2.cc.

38 {
39
40 model_name="Chi2BB2";
41
42}

◆ init()

void EvtChi2BB2::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 51 of file EvtChi2BB2.cc.

51 {
52// check that there are 0 arguments
53 checkNArg(5);
54 checkNDaug(2);
58
59 }
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
void checkSpinParent(EvtSpinType::spintype sp)
void checkNDaug(int d1, int d2=-1)
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)

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