BOSS 7.1.0
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | List of all members
EvtD0ToKSpipipi0pi0 Class Reference

#include <EvtD0ToKSpipipi0pi0.hh>

+ Inheritance diagram for EvtD0ToKSpipipi0pi0:

Public Member Functions

 EvtD0ToKSpipipi0pi0 ()
 
virtual ~EvtD0ToKSpipipi0pi0 ()
 
void getName (std::string &name)
 
EvtDecayBase * clone ()
 
void init ()
 
void initProbMax ()
 
void decay (EvtParticle *p)
 
- Public Member Functions inherited from EvtDecayProb
void makeDecay (EvtParticle *p)
 
void setProb (double prob)
 
double getProb ()
 
void setWeight (double weight)
 
virtual ~EvtDecayProb ()
 
- Public Member Functions inherited from EvtDecayBase
virtual void getName (std::string &name)=0
 
virtual void decay (EvtParticle *p)=0
 
virtual void makeDecay (EvtParticle *p)=0
 
virtual EvtDecayBase * clone ()=0
 
virtual void init ()
 
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 EvtDecayBase
bool _daugsDecayedByParentModel
 

Detailed Description

Definition at line 9 of file EvtD0ToKSpipipi0pi0.hh.

Constructor & Destructor Documentation

◆ EvtD0ToKSpipipi0pi0()

EvtD0ToKSpipipi0pi0::EvtD0ToKSpipipi0pi0 ( )
inline

Definition at line 12 of file EvtD0ToKSpipipi0pi0.hh.

12{}

Referenced by clone().

◆ ~EvtD0ToKSpipipi0pi0()

EvtD0ToKSpipipi0pi0::~EvtD0ToKSpipipi0pi0 ( )
virtual

Definition at line 35 of file EvtD0ToKSpipipi0pi0.cc.

35{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtD0ToKSpipipi0pi0::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 41 of file EvtD0ToKSpipipi0pi0.cc.

41 {
42 return new EvtD0ToKSpipipi0pi0;
43 }

◆ decay()

void EvtD0ToKSpipipi0pi0::decay ( EvtParticle *  p)
virtual

Implements EvtDecayBase.

Definition at line 184 of file EvtD0ToKSpipipi0pi0.cc.

184 {
185/*
186 double maxprob = 0.0;
187 for(int ir=0;ir<=60000000;ir++){
188
189 p->initializePhaseSpace(getNDaug(),getDaugs());
190 EvtVector4R D1 = p->getDaug(0)->getP4();
191 EvtVector4R D2 = p->getDaug(1)->getP4();
192 EvtVector4R D3 = p->getDaug(2)->getP4();
193 EvtVector4R D4 = p->getDaug(3)->getP4();
194 EvtVector4R D5 = p->getDaug(4)->getP4();
195
196 double P1[4], P2[4], P3[4], P4[4], P5[4];
197 P1[0] = D1.get(0); P1[1] = D1.get(1); P1[2] = D1.get(2); P1[3] = D1.get(3);
198 P2[0] = D2.get(0); P2[1] = D2.get(1); P2[2] = D2.get(2); P2[3] = D2.get(3);
199 P3[0] = D3.get(0); P3[1] = D3.get(1); P3[2] = D3.get(2); P3[3] = D3.get(3);
200 P4[0] = D4.get(0); P4[1] = D4.get(1); P4[2] = D4.get(2); P4[3] = D4.get(3);
201 P5[0] = D5.get(0); P5[1] = D5.get(1); P5[2] = D5.get(2); P5[3] = D5.get(3);
202
203 double value;
204 int g0[8]={1,1,1,0,1,1,1,1};
205 int g1[8]={1,1,1,0,1,1,1,1};
206 double r0[8]={1,1,1,0,1,1,1,1};
207 double r1[8]={1,1,1,0,1,1,1,1};
208 int spin[8]={0,1,2,0,1020,1020,0,1};
209
210 int nstates=8;
211 calEva(P1, P4, P2, P3, P5, mass1, mass2, width1, width2, rho, phi, r0, r1, g0, g1, spin, modetype, nstates, value);
212
213 if (value<0) continue;
214 if(value>maxprob) {
215 maxprob=value;
216 cout << "ir= " << ir << endl;
217 cout << "double P1[4] = {" << P1[0] <<","<< P1[1] <<","<< P1[2] <<","<< P1[3] <<"};"<< endl;
218 cout << "double P2[4] = {" << P2[0] <<","<< P2[1] <<","<< P2[2] <<","<< P2[3] <<"};"<< endl;
219 cout << "double P3[4] = {" << P3[0] <<","<< P3[1] <<","<< P3[2] <<","<< P3[3] <<"};"<< endl;
220 cout << "double P4[4] = {" << P4[0] <<","<< P4[1] <<","<< P4[2] <<","<< P4[3] <<"};"<< endl;
221 cout << "double P5[4] = {" << P5[0] <<","<< P5[1] <<","<< P5[2] <<","<< P5[3] <<"};"<< endl;
222 cout << "MAX====> " << maxprob << endl;
223
224 }
225 }
226 printf("MAXprob = %.10f\n",maxprob);
227*/
228 p->initializePhaseSpace(getNDaug(),getDaugs());
229 EvtVector4R D1 = p->getDaug(0)->getP4();
230 EvtVector4R D2 = p->getDaug(1)->getP4();
231 EvtVector4R D3 = p->getDaug(2)->getP4();
232 EvtVector4R D4 = p->getDaug(3)->getP4();
233 EvtVector4R D5 = p->getDaug(4)->getP4();
234
235 double P1[4], P2[4], P3[4], P4[4], P5[4];
236 P1[0] = D1.get(0); P1[1] = D1.get(1); P1[2] = D1.get(2); P1[3] = D1.get(3);
237 P2[0] = D2.get(0); P2[1] = D2.get(1); P2[2] = D2.get(2); P2[3] = D2.get(3);
238 P3[0] = D3.get(0); P3[1] = D3.get(1); P3[2] = D3.get(2); P3[3] = D3.get(3);
239 P4[0] = D4.get(0); P4[1] = D4.get(1); P4[2] = D4.get(2); P4[3] = D4.get(3);
240 P5[0] = D5.get(0); P5[1] = D5.get(1); P5[2] = D5.get(2); P5[3] = D5.get(3);
241
242// P1[0] = 0.699073767764678; P1[1] = 0.116528634106469; P1[2] = -0.416631548481574; P1[3] = -0.232030327801847;
243// P2[0] = 0.431675469128871; P2[1] = 0.090021693909732; P2[2] = 0.149419303854623; P2[3] =0.371064840093330;
244// P3[0] = 0.188455625514340; P3[1] = -0.100081780055657;P3[2] = -0.052767368589371; P3[3] = 0.056876884593353;
245// P4[0] = 0.311123417560542; P4[1] = 0.231980078879333; P4[2] = 0.149145989336585; P4[3] = -0.035478344416011;
246// P5[0] = 0.257852771533528; P5[1] = -0.185420043160991;P5[2] = -0.073435020040710; P5[3] = -0.092139879037438;
247
248 double value;
249 int g0[8]={1,1,1,0,1,1,1,1};
250 int g1[8]={1,1,1,0,1,1,1,1};
251 double r0[8]={1,1,1,0,1,1,1,1};
252 double r1[8]={1,1,1,0,1,1,1,1};
253 int spin[8]={0,1,2,0,1020,1020,0,0};
254
255 int nstates=8;
256 calEva(P1, P4, P2, P3, P5, mass1, mass2, width1, width2, rho, phi, r0, r1, g0, g1, spin, modetype, nstates, value);
257 setProb(value);
258
259 return ;
260}
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition: EvtDecayBase.hh:65
EvtDecayProb::setProb
void setProb(double prob)
Definition: EvtDecayProb.hh:34
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition: EvtParticle.cc:120
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition: EvtParticle.cc:84
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition: EvtParticle.cc:1070
EvtVector4R::get
double get(int i) const
Definition: EvtVector4R.hh:179

◆ getName()

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

Implements EvtDecayBase.

Definition at line 37 of file EvtD0ToKSpipipi0pi0.cc.

37 {
38 model_name="D0ToKSpipipi0pi0";
39 }

◆ init()

void EvtD0ToKSpipipi0pi0::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 44 of file EvtD0ToKSpipipi0pi0.cc.

44 {
45 checkNArg(0);
46 checkNDaug(5);
47 checkSpinParent(EvtSpinType::SCALAR);
48 checkSpinDaughter(0,EvtSpinType::SCALAR);
49 checkSpinDaughter(1,EvtSpinType::SCALAR);
50 checkSpinDaughter(2,EvtSpinType::SCALAR);
51 checkSpinDaughter(3,EvtSpinType::SCALAR);
52 checkSpinDaughter(4,EvtSpinType::SCALAR);
53
54 string name;
55 Double_t mean_v;
56
57 phi[0] = 0.0;
58 phi[1] = 1.808202793;
59 phi[2] = 1.353154265;
60 phi[3] = 0.02534778524;
61 phi[4] = 2.572695445;
62 phi[5] = 2.572695445;
63 phi[6] = -2.683237531;
64 phi[7] = 0.1209896124;
65
66 rho[0] = 1.0;
67 rho[1] = 0.5632373979;
68 rho[2] = 1.541400902;
69 rho[3] = 163.725261;
70 rho[4] = 3.207856388;
71 rho[5] = 3.207856388;
72 rho[6] = 2.943199635;
73 rho[7] = 3.023783561;
74
75 modetype[0]= 1;
76 modetype[1]= 1;
77 modetype[2]= 1;
78 modetype[3]= 41;
79 modetype[4]= 20;
80 modetype[5]= 33;
81 modetype[6]= 4;
82 modetype[7]= 40;
83
84 //std::cout << "Initializing EvtD0ToKSpipipi0pi0" << std::endl;
85 //for (int i=0; i<8; i++) {
86 // cout << i << " rho= " << rho[i] << " phi= " << phi[i] << endl;
87 //}
88 width1[0] = 0.00849;
89 width1[1] = 0.00849;
90 width1[2] = 0.00849;
91 width1[3] = 0.00849;//ks4pi
92 width1[4] = 0.420;
93 width1[5] = 0.420;
94 width1[6] = 0.142;
95 width1[7] = 0.00849;
96
97 width2[0] = 0.0473;
98 width2[1] = 0.0473;
99 width2[2] = 0.0473;
100 width2[3] = 0.0473;
101 width2[4] = 0.0473;
102 width2[5] = 0.0473;
103 width2[6] = 0.00849;
104 width2[7] = 0.0473;
105
106 mass1[0] = 0.78265;
107 mass1[1] = 0.78265;
108 mass1[2] = 0.78265;
109 mass1[3] = 0.78265;
110 mass1[4] = 1.230;
111 mass1[5] = 1.230;
112 mass1[6] = 1.2295;
113 mass1[7] = 0.78265;
114
115 mass2[0] = 0.89555;
116 mass2[1] = 0.89555;
117 mass2[2] = 0.89555;
118 mass2[3] = 0.89555;
119 mass2[4] = 0.89555;
120 mass2[5] = 0.89555;
121 mass2[6] = 0.78265;
122 mass2[7] = 0.89555;
123
124 mD0M = 1.86486;
125 mD = 1.86486;
126 metap = 0.95778;
127 mkstr = 0.89594;
128 mk0 = 0.497614;
129 mass_Kaon = 0.49368;
130 mass_Pion = 0.13957;
131 mass_Pion2 = 0.0194797849;
132 mass_2Pion = 0.27914;
133 math_2pi = 6.2831852;
134 rD2 = 25.0; // 5*5
135 rDs2 = 25.0; // 5*5
136 rRes2 = 9.0; // 3*3
137 gg1 = 0.5468;
138 gg2 = 0.23; // K*0(1430)
139 GS1 = 0.636619783;
140 GS2 = 0.01860182466;
141 GS3 = 0.1591549458; // 1/(2*math_2pi)
142 GS4 = 0.00620060822;
143
144 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
145 for (int i=0; i<4; i++) {
146 for (int j=0; j<4; j++) {
147 G[i][j] = GG[i][j];
148 }
149 }
150 double EE[4][4][4][4] =
151{ { {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
152 {{0,0,0,0}, {0,0,0,0}, {0,0,0,1}, {0,0,-1,0}},
153 {{0,0,0,0}, {0,0,0,-1}, {0,0,0,0}, {0,1,0,0} },
154 {{0,0,0,0}, {0,0,1,0}, {0,-1,0,0}, {0,0,0,0} } },
155 { {{0,0,0,0}, {0,0,0,0}, {0,0,0,-1}, {0,0,1,0} },
156 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
157 {{0,0,0,1}, {0,0,0,0}, {0,0,0,0}, {-1,0,0,0}},
158 {{0,0,-1,0}, {0,0,0,0}, {1,0,0,0}, {0,0,0,0} } },
159 { {{0,0,0,0}, {0,0,0,1}, {0,0,0,0}, {0,-1,0,0}},
160 {{0,0,0,-1}, {0,0,0,0}, {0,0,0,0}, {1,0,0,0} },
161 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
162 {{0,1,0,0}, {-1,0,0,0}, {0,0,0,0}, {0,0,0,0} } },
163 { {{0,0,0,0}, {0,0,-1,0}, {0,1,0,0}, {0,0,0,0} },
164 {{0,0,1,0}, {0,0,0,0}, {-1,0,0,0}, {0,0,0,0} },
165 {{0,-1,0,0}, {1,0,0,0}, {0,0,0,0}, {0,0,0,0} },
166 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} } } };
167
168 for (int i=0; i<4; i++) {
169 for (int j=0; j<4; j++) {
170 for (int k=0; k<4; k++) {
171 for (int l=0; l<4; l++) {
172 E[i][j][k][l] = EE[i][j][k][l];
173 }
174 }
175 }
176 }
177
178}
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 EvtD0ToKSpipipi0pi0::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 180 of file EvtD0ToKSpipipi0pi0.cc.

180 {
181 setProbMax(52000000);
182}
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition: EvtDecayBase.cc:297
The documentation for this class was generated from the following files: