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

#include <EvtD0ToKpipi0pi0.hh>

+ Inheritance diagram for EvtD0ToKpipi0pi0:

Public Member Functions

 EvtD0ToKpipi0pi0 ()
 
virtual ~EvtD0ToKpipi0pi0 ()
 
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 8 of file EvtD0ToKpipi0pi0.hh.

Constructor & Destructor Documentation

◆ EvtD0ToKpipi0pi0()

EvtD0ToKpipi0pi0::EvtD0ToKpipi0pi0 ( )
inline

Definition at line 12 of file EvtD0ToKpipi0pi0.hh.

12{}

Referenced by clone().

◆ ~EvtD0ToKpipi0pi0()

EvtD0ToKpipi0pi0::~EvtD0ToKpipi0pi0 ( )
virtual

Definition at line 31 of file EvtD0ToKpipi0pi0.cc.

31{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtD0ToKpipi0pi0::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 37 of file EvtD0ToKpipi0pi0.cc.

37 {
38 return new EvtD0ToKpipi0pi0;
39}

◆ decay()

void EvtD0ToKpipi0pi0::decay ( EvtParticle *  p)
virtual

Implements EvtDecayBase.

Definition at line 182 of file EvtD0ToKpipi0pi0.cc.

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

Implements EvtDecayBase.

Definition at line 33 of file EvtD0ToKpipi0pi0.cc.

33 {
34 model_name="D0ToKpipi0pi0";
35}

◆ init()

void EvtD0ToKpipi0pi0::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 41 of file EvtD0ToKpipi0pi0.cc.

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

Reimplemented from EvtDecayBase.

Definition at line 178 of file EvtD0ToKpipi0pi0.cc.

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