BOSS 7.0.9
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtD0ToKpipi0pi0 Class Reference

#include <EvtD0ToKpipi0pi0.hh>

+ Inheritance diagram for EvtD0ToKpipi0pi0:

Public Member Functions

 EvtD0ToKpipi0pi0 ()
 
virtual ~EvtD0ToKpipi0pi0 ()
 
void getName (std::string &name)
 
EvtDecayBaseclone ()
 
void init ()
 
void initProbMax ()
 
void decay (EvtParticle *p)
 
 EvtD0ToKpipi0pi0 ()
 
virtual ~EvtD0ToKpipi0pi0 ()
 
void getName (std::string &name)
 
EvtDecayBaseclone ()
 
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 ()
 
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 EvtDecayBaseclone ()=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 ()
 
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 ()
 
virtual void getName (std::string &name)=0
 
virtual void decay (EvtParticle *p)=0
 
virtual void makeDecay (EvtParticle *p)=0
 
virtual EvtDecayBaseclone ()=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 ()
 
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)
 
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 ()
 
bool daugsDecayedByParentModel ()
 
- Protected Attributes inherited from EvtDecayBase
bool _daugsDecayedByParentModel
 

Detailed Description

Constructor & Destructor Documentation

◆ EvtD0ToKpipi0pi0() [1/2]

EvtD0ToKpipi0pi0::EvtD0ToKpipi0pi0 ( )
inline

Definition at line 12 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtD0ToKpipi0pi0.hh.

12{}

Referenced by clone().

◆ ~EvtD0ToKpipi0pi0() [1/2]

EvtD0ToKpipi0pi0::~EvtD0ToKpipi0pi0 ( )
virtual

◆ EvtD0ToKpipi0pi0() [2/2]

EvtD0ToKpipi0pi0::EvtD0ToKpipi0pi0 ( )
inline

◆ ~EvtD0ToKpipi0pi0() [2/2]

virtual EvtD0ToKpipi0pi0::~EvtD0ToKpipi0pi0 ( )
virtual

Member Function Documentation

◆ clone() [1/2]

◆ clone() [2/2]

EvtDecayBase * EvtD0ToKpipi0pi0::clone ( )
virtual

Implements EvtDecayBase.

◆ decay() [1/2]

void EvtD0ToKpipi0pi0::decay ( EvtParticle p)
virtual

Implements EvtDecayBase.

Definition at line 181 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtD0ToKpipi0pi0.cc.

181 {
182/*
183 double maxprob = 0.0;
184 for(int ir=0;ir<=60000000;ir++){
185 p->initializePhaseSpace(getNDaug(),getDaugs());
186 EvtVector4R Km0 = p->getDaug(0)->getP4();
187 EvtVector4R pi1 = p->getDaug(1)->getP4();
188 EvtVector4R pi2 = p->getDaug(2)->getP4();
189 EvtVector4R pi3 = p->getDaug(3)->getP4();
190 double Km[4],Pip[4],Pi01[4],Pi02[4];
191 Km[0] = Km0.get(0); Pip[0] = pi1.get(0); Pi01[0] = pi2.get(0); Pi02[0] = pi3.get(0);
192 Km[1] = Km0.get(1); Pip[1] = pi1.get(1); Pi01[1] = pi2.get(1); Pi02[1] = pi3.get(1);
193 Km[2] = Km0.get(2); Pip[2] = pi1.get(2); Pi01[2] = pi2.get(2); Pi02[2] = pi3.get(2);
194 Km[3] = Km0.get(3); Pip[3] = pi1.get(3); Pi01[3] = pi2.get(3); Pi02[3] = pi3.get(3);
195 double Prob = calPDF(Km, Pip, Pi01, Pi02);
196 if(Prob>maxprob) {
197 maxprob=Prob;
198 std::cout << "Max PDF = " << ir << " prob= " << Prob << std::endl;
199 }
200 }
201 std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
202*/
204 EvtVector4R Km0 = p->getDaug(0)->getP4();
205 EvtVector4R pi1 = p->getDaug(1)->getP4();
206 EvtVector4R pi2 = p->getDaug(2)->getP4();
207 EvtVector4R pi3 = p->getDaug(3)->getP4();
208
209 double Km[4],Pip[4],Pi01[4],Pi02[4];
210 Km[0] = Km0.get(0); Pip[0] = pi1.get(0); Pi01[0] = pi2.get(0); Pi02[0] = pi3.get(0);
211 Km[1] = Km0.get(1); Pip[1] = pi1.get(1); Pi01[1] = pi2.get(1); Pi02[1] = pi3.get(1);
212 Km[2] = Km0.get(2); Pip[2] = pi1.get(2); Pi01[2] = pi2.get(2); Pi02[2] = pi3.get(2);
213 Km[3] = Km0.get(3); Pip[3] = pi1.get(3); Pi01[3] = pi2.get(3); Pi02[3] = pi3.get(3);
214 double prob = calPDF(Km, Pip, Pi01, Pi02);
215 setProb(prob);
216 return;
217}
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)

◆ decay() [2/2]

void EvtD0ToKpipi0pi0::decay ( EvtParticle p)
virtual

Implements EvtDecayBase.

◆ getName() [1/2]

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

Implements EvtDecayBase.

Definition at line 32 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtD0ToKpipi0pi0.cc.

32 {
33 model_name="D0ToKpipi0pi0";
34}

◆ getName() [2/2]

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

Implements EvtDecayBase.

◆ init() [1/2]

void EvtD0ToKpipi0pi0::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 40 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtD0ToKpipi0pi0.cc.

40 {
41 checkNArg(0);
42 checkNDaug(4);
44/*
45 checkSpinDaughter(0,EvtSpinType::SCALAR);
46 checkSpinDaughter(1,EvtSpinType::SCALAR);
47 checkSpinDaughter(3,EvtSpinType::SCALAR);
48 checkSpinDaughter(4,EvtSpinType::SCALAR);
49*/
50 std::cout << "EvtD0ToKpipi0pi0 ==> Initialization !" << std::endl;
51
52 mod[0]= 1; rho[0]= 2.02; phi[0]= -0.75;
53 mod[1]= 1; rho[1]= 1.66; phi[1]= -2.90;
54 mod[2]= 0; rho[2]= 0; phi[2]= 0;
55 mod[3]= 0; rho[3]= 0; phi[3]= 0;
56 mod[4]= 0; rho[4]= 0; phi[4]= 0;
57 mod[5]= 0; rho[5]= 0; phi[5]= 0;
58 mod[6]= 0; rho[6]= 0; phi[6]= 0;
59 mod[7]= 1; rho[7]= 1; phi[7]= 0;
60 mod[8]= 1; rho[8]= 0.842; phi[8]= -2.05;
61 mod[9]= 1; rho[9]= 0.0218; phi[9]= 1.84;
62 mod[10]=0; rho[10]= 0; phi[10]= 0;
63 mod[11]=0; rho[11]= 0; phi[11]= 0;
64 mod[12]=0; rho[12]= 0; phi[12]= 0;
65 mod[13]=1; rho[13]= 0.0336; phi[13]= -1.55;
66 mod[14]=1; rho[14]= 0.109; phi[14]= -1.35;
67 mod[15]=1; rho[15]= 0.196; phi[15]= -2.07;
68 mod[16]=0; rho[16]= 0; phi[16]= 0;
69 mod[17]=0; rho[17]= 0; phi[17]= 0;
70 mod[18]=0; rho[18]= 0; phi[18]= 0;
71 mod[19]=1; rho[19]= 0.363; phi[19]= 1.93;
72 mod[20]=0; rho[20]= 0; phi[20]= 0;
73 mod[21]=0; rho[21]= 0; phi[21]= 0;
74 mod[22]=0; rho[22]= 0; phi[22]= 0;
75 mod[23]=1; rho[23]= 0.555; phi[23]= 0.44;
76 mod[24]=1; rho[24]= 0.526; phi[24]= -1.84;
77 mod[25]=0; rho[25]= 0; phi[25]= 0;
78 mod[26]=1; rho[26]= 1; phi[26]= 0.64;
79 mod[27]=0; rho[27]= 0; phi[27]= 0;
80 mod[28]=0; rho[28]= 0; phi[28]= 0;
81 mod[29]=1; rho[29]= 3.34; phi[29]= -0.02;
82 mod[30]=0; rho[30]= 0; phi[30]= 0;
83 mod[31]=0; rho[31]= 0; phi[31]= 0;
84 mod[32]=0; rho[32]= 0; phi[32]= 0;
85 mod[33]=0; rho[33]= 0; phi[33]= 0;
86 mod[34]=1; rho[34]= 1.76; phi[34]= -2.39;
87 mod[35]=1; rho[35]= 0.175; phi[35]= 1.59;
88 mod[36]=1; rho[36]= 0.397; phi[36]= 1.45;
89 mod[37]=0; rho[37]= 0; phi[37]= 0;
90 mod[38]=0; rho[38]= 0; phi[38]= 0;
91 mod[39]=0; rho[39]= 0; phi[39]= 0;
92 mod[40]=0; rho[40]= 0; phi[40]= 0;
93 mod[41]=1; rho[41]= 1.02; phi[41]= 0.52;
94 mod[42]=0; rho[42]= 0; phi[42]= 0;
95 mod[43]=0; rho[43]= 0; phi[43]= 0;
96 mod[44]=0; rho[44]= 0; phi[44]= 0;
97 mod[45]=0; rho[45]= 0; phi[45]= 0;
98 mod[46]=1; rho[46]= 0.146; phi[46]= 1.24;
99 mod[47]=1; rho[47]= 0.0978; phi[47]= -2.89;
100 mod[48]=1; rho[48]= 0.233; phi[48]= 2.41;
101 mod[49]=0; rho[49]= 0; phi[49]= 0;
102 mod[50]=1; rho[50]= 0.424; phi[50]= -0.94;
103 mod[51]=1; rho[51]= 1.03; phi[51]= -1.93;
104 mod[52]=0; rho[52]= 0; phi[52]= 0;
105 mod[53]=0; rho[53]= 0; phi[53]= 0;
106 mod[54]=1; rho[54]= 0.474; phi[54]= -1.17;
107 mod[55]=0; rho[55]= 0; phi[55]= 0;
108 mod[56]=0; rho[56]= 0; phi[56]= 0;
109 mod[57]=0; rho[57]= 0; phi[57]= 0;
110 mod[58]=0; rho[58]= 0; phi[58]= 0;
111 mod[59]=0; rho[59]= 0; phi[59]= 0;
112 mod[60]=0; rho[60]= 0; phi[60]= 0;
113 mod[61]=1; rho[61]= 6.74; phi[61]= -1.74;
114 mod[62]=0; rho[62]= 0; phi[62]= 0;
115 mod[63]=0; rho[63]= 0; phi[63]= 0;
116 mod[64]=0; rho[64]= 0; phi[64]= 0;
117 mod[65]=0; rho[65]= 0; phi[65]= 0;
118 mod[66]=1; rho[66]= 1.54; phi[66]= -2.93;
119 mod[67]=1; rho[67]= 1.36; phi[67]= 2.23;
120
121 mass[0]= 1.23; width[0]= 0.50204;
122 mass[1]= 1.2723; width[1]= 0.09;
123 mass[2]= 0.89166; width[2]= 0.0508;
124 mass[3]= 0.89581; width[3]= 0.0474;
125 mass[4]= 0.77511; width[4]= 0.1491;
126
127 for (int i=0; i<69; i++) {
128 cout << i << "rho,phi = " << rho[i] << ", "<< phi[i] << endl;
129 }
130 mD = 1.86486;
131 rRes = 3.0;
132 rD = 5.0;
133 metap = 0.95778;
134 mk0 = 0.497614;
135 mass_Kaon = 0.49368;
136 mass_Pion = 0.13957;
137 math_pi = 3.1415926;
138 mass_Pi0 = 0.1349766;
139 mkstrm = 0.89594;
140 mkstr0 = 0.89594;
141
142 pi = 3.1415926;
143 mpi = 0.13957;
144 g1 = 0.5468;
145 g2 = 0.23;
146
147 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
148 int EE[4][4][4][4] =
149 { { {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
150 {{0,0,0,0}, {0,0,0,0}, {0,0,0,1}, {0,0,-1,0}},
151 {{0,0,0,0}, {0,0,0,-1}, {0,0,0,0}, {0,1,0,0} },
152 {{0,0,0,0}, {0,0,1,0}, {0,-1,0,0}, {0,0,0,0} } },
153 { {{0,0,0,0}, {0,0,0,0}, {0,0,0,-1}, {0,0,1,0} },
154 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
155 {{0,0,0,1}, {0,0,0,0}, {0,0,0,0}, {-1,0,0,0}},
156 {{0,0,-1,0}, {0,0,0,0}, {1,0,0,0}, {0,0,0,0} } },
157 { {{0,0,0,0}, {0,0,0,1}, {0,0,0,0}, {0,-1,0,0}},
158 {{0,0,0,-1}, {0,0,0,0}, {0,0,0,0}, {1,0,0,0} },
159 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
160 {{0,1,0,0}, {-1,0,0,0}, {0,0,0,0}, {0,0,0,0} } },
161 { {{0,0,0,0}, {0,0,-1,0}, {0,1,0,0}, {0,0,0,0} },
162 {{0,0,1,0}, {0,0,0,0}, {-1,0,0,0}, {0,0,0,0} },
163 {{0,-1,0,0}, {1,0,0,0}, {0,0,0,0}, {0,0,0,0} },
164 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} } } };
165 for (int i=0; i<4; i++) {
166 for (int j=0; j<4; j++) {
167 G[i][j] = GG[i][j];
168 for (int k=0; k<4; k++) {
169 for (int l=0; l<4; l++) {
170 E[i][j][k][l] = EE[i][j][k][l];
171 }
172 }
173 }
174 }
175}
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)

◆ init() [2/2]

void EvtD0ToKpipi0pi0::init ( )
virtual

Reimplemented from EvtDecayBase.

◆ initProbMax() [1/2]

void EvtD0ToKpipi0pi0::initProbMax ( )
virtual

◆ initProbMax() [2/2]

void EvtD0ToKpipi0pi0::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.


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