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

#include <EvtD0ToKSpi0pi0pi0.hh>

+ Inheritance diagram for EvtD0ToKSpi0pi0pi0:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtD0ToKSpi0pi0pi0()

EvtD0ToKSpi0pi0pi0::EvtD0ToKSpi0pi0pi0 ( )
inline

Definition at line 11 of file EvtD0ToKSpi0pi0pi0.hh.

11{}

Referenced by clone().

◆ ~EvtD0ToKSpi0pi0pi0()

EvtD0ToKSpi0pi0pi0::~EvtD0ToKSpi0pi0pi0 ( )
virtual

Definition at line 33 of file EvtD0ToKSpi0pi0pi0.cc.

33{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtD0ToKSpi0pi0pi0::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 39 of file EvtD0ToKSpi0pi0pi0.cc.

39 {
40 return new EvtD0ToKSpi0pi0pi0;
41}

◆ decay()

void EvtD0ToKSpi0pi0pi0::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 158 of file EvtD0ToKSpi0pi0pi0.cc.

158 {
159 //-----------for max value------------------
160/* double maxprob = 0.0;
161 for(int ir=0;ir<=60000000;ir++){
162 p->initializePhaseSpace(getNDaug(),getDaugs());
163 EvtVector4R ks = p->getDaug(0)->getP4();
164 EvtVector4R kp = p->getDaug(1)->getP4();
165 EvtVector4R pim = p->getDaug(2)->getP4();
166 EvtVector4R pi0 = p->getDaug(3)->getP4();
167 mother_c=EvtPDL::getStdHep(p->getId());
168 int cc;
169 if(mother_c==421){
170 cc=1;
171 } else if(mother_c==-421){
172 cc=-1;
173 }
174 double value;
175 double KS[4],Pi01[4],Pi02[4],Pi03[4];
176 KS[0] = ks.get(0); Pi01[0] = kp.get(0); Pi02[0] = pim.get(0); Pi03[0] = pi0.get(0);
177 KS[1] = cc*ks.get(1); Pi01[1] = cc*kp.get(1); Pi02[1] = cc*pim.get(1); Pi03[1] = cc*pi0.get(1);
178 KS[2] = cc*ks.get(2); Pi01[2] = cc*kp.get(2); Pi02[2] = cc*pim.get(2); Pi03[2] = cc*pi0.get(2);
179 KS[3] = cc*ks.get(3); Pi01[3] = cc*kp.get(3); Pi02[3] = cc*pim.get(3); Pi03[3] = cc*pi0.get(3);
180 calPDF(KS, Pi01, Pi02, Pi03, value);
181 if(value>maxprob) {
182 maxprob=value;
183 std::cout << "Max PDF = " << ir << " prob= " << value << std::endl;
184 }
185 }
186 std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
187 return;*/
188 //-----------------------------------------------
189 p->initializePhaseSpace(getNDaug(),getDaugs());
190 EvtVector4R ks = p->getDaug(0)->getP4();
191 EvtVector4R pi01 = p->getDaug(1)->getP4();
192 EvtVector4R pi02 = p->getDaug(2)->getP4();
193 EvtVector4R pi03 = p->getDaug(3)->getP4();
194
195 mother_c=EvtPDL::getStdHep(p->getId());
196 int cc;
197 if(mother_c==421){
198 cc=1;
199 } else if(mother_c==-421){
200 cc=-1;
201 }
202
203 double KS[4],Pi01[4],Pi02[4],Pi03[4];
204 KS[0] = ks.get(0); Pi01[0] = pi01.get(0); Pi02[0] = pi02.get(0); Pi03[0] = pi03.get(0);
205 KS[1] = cc*ks.get(1); Pi01[1] = cc*pi01.get(1); Pi02[1] = cc*pi02.get(1); Pi03[1] = cc*pi03.get(1);
206 KS[2] = cc*ks.get(2); Pi01[2] = cc*pi01.get(2); Pi02[2] = cc*pi02.get(2); Pi03[2] = cc*pi03.get(2);
207 KS[3] = cc*ks.get(3); Pi01[3] = cc*pi01.get(3); Pi02[3] = cc*pi02.get(3); Pi03[3] = cc*pi03.get(3);
208
209 double value;
210 calPDF(KS, Pi01, Pi02, Pi03, value);
211 setProb(value);
212 return;
213
214}
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition EvtDecayBase.hh:65
EvtDecayProb::setProb
void setProb(double prob)
Definition EvtDecayProb.hh:34
EvtPDL::getStdHep
static int getStdHep(EvtId id)
Definition EvtPDL.hh:56
EvtParticle::getId
EvtId getId() const
Definition EvtParticle.cc:112
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 EvtD0ToKSpi0pi0pi0::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 35 of file EvtD0ToKSpi0pi0pi0.cc.

35 {
36 model_name="D0ToKSpi0pi0pi0";
37}

◆ init()

void EvtD0ToKSpi0pi0pi0::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 43 of file EvtD0ToKSpi0pi0pi0.cc.

43 {
44 checkNArg(0);
45 checkNDaug(4);
46 checkSpinParent(EvtSpinType::SCALAR);
47 //cout << "Initializing EvtD0ToKSpi0pi0pi0" << endl;
48
49 mD = 1.86484; mPi = 0.13957; mKa = 0.493677;
50 math_pi = 3.1415926;
51 mKstr0 = 0.89555; mrhop = 0.77511;
52 GKstr0 = 0.0473; Grhop = 0.1491;
53
54 mass_Pion = 0.13957;
55 mass_Eta = 0.547862;
56 rD2 = 25.0; // 5*5
57 rRes2 = 9.0; // 3*3
58
59 ma1_1260 = 1.195; mK1_1270 = 1.2898;
60 Ga1_1260 = 0.422; GK1_1270 = 0.1161;
61 mpi_1300 = 1.371; mK1_1400 = 1.403;
62 Gpi_1300 = 0.314; GK1_1400 = 0.174;
63 mf0_500 = 0.5; Gf0_500 = 0.335;
64 mK1460 = 1.4824; GK1460 = 0.335;
65//--------fixed---------------//
66 phi01 = 0.;
67 phi04 = 0.;
68 phi07 = 0.;
69 phi08 = 0.;
70 phi10 = 0.;
71 phi126 = 0.;
72
73 rho01 = 1.0;
74 rho04 = 1.0;
75 rho07 = 1.0;
76 rho08 = 1.0;
77 rho10 = 1.0;
78 rho126 = 1.0;
79//---------------------------//
80
81
82 K1460b1 = 1.7104e+00;
83 K1460b2 = 7.4066e+00;
84 K1460b6 = -1.4986e+00;
85 K1460b7 = -8.0174e-01;
86 K1460f2 = 1.7438e+01;
87 K1460f7 = -4.2315e+00;
88 Kstb2 = 1.2922e+00;
89 Kstb7 = 4.9340e+00;
90 Kstf1 = 2.7165e-01;
91 Kstf6 = 9.3143e-01;
92 Phspb2 = 2.1656e+00;
93 Phspb7 = 4.6264e+00;
94 Phspf1 = 1.2052e-01;
95 Phspf6 = 1.2506e+00;
96 a1260b1 = 2.3244e+00;
97 a1260b2 = 1.4414e+01;
98 a1260b6 = 4.3493e+00;
99 a1260b7 = 3.1104e+00;
100 a1260f1 = 3.3290e+00;
101 a1260f6 = 4.7215e+00;
102 a321 = 1.8517e+01;
103 a322 = 6.3000e+00;
104 ake1 = 2.0000e+02;
105 ake2 = 3.1394e+00;
106 phi02 = 9.5172e+00;
107 phi03 = 7.7804e+00;
108 phi06 = 5.6253e+00;
109 rho02 = 2.0961e-01;
110 rho03 = 7.5251e-01;
111 rho06 = 1.9247e+00;
112
113
114
115 modetype[0]=0;
116 modetype[1]=1;
117 modetype[2]=1;
118 modetype[3]=2;
119 modetype[4]=2;
120 modetype[5]=1;
121 modetype[6]=33;
122 modetype[7]=11;
123
124 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
125 int EE[4][4][4][4] =
126 { { {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
127 {{0,0,0,0}, {0,0,0,0}, {0,0,0,1}, {0,0,-1,0}},
128 {{0,0,0,0}, {0,0,0,-1}, {0,0,0,0}, {0,1,0,0} },
129 {{0,0,0,0}, {0,0,1,0}, {0,-1,0,0}, {0,0,0,0} } },
130 { {{0,0,0,0}, {0,0,0,0}, {0,0,0,-1}, {0,0,1,0} },
131 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
132 {{0,0,0,1}, {0,0,0,0}, {0,0,0,0}, {-1,0,0,0}},
133 {{0,0,-1,0}, {0,0,0,0}, {1,0,0,0}, {0,0,0,0} } },
134 { {{0,0,0,0}, {0,0,0,1}, {0,0,0,0}, {0,-1,0,0}},
135 {{0,0,0,-1}, {0,0,0,0}, {0,0,0,0}, {1,0,0,0} },
136 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
137 {{0,1,0,0}, {-1,0,0,0}, {0,0,0,0}, {0,0,0,0} } },
138 { {{0,0,0,0}, {0,0,-1,0}, {0,1,0,0}, {0,0,0,0} },
139 {{0,0,1,0}, {0,0,0,0}, {-1,0,0,0}, {0,0,0,0} },
140 {{0,-1,0,0}, {1,0,0,0}, {0,0,0,0}, {0,0,0,0} },
141 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} } } };
142 for (int i=0; i<4; i++) {
143 for (int j=0; j<4; j++) {
144 G[i][j] = GG[i][j];
145 for (int k=0; k<4; k++) {
146 for (int l=0; l<4; l++) {
147 E[i][j][k][l] = EE[i][j][k][l];
148 }
149 }
150 }
151 }
152}
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 EvtD0ToKSpi0pi0pi0::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 154 of file EvtD0ToKSpi0pi0pi0.cc.

154 {
155 setProbMax(77195.0);
156}
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition EvtDecayBase.cc:297
The documentation for this class was generated from the following files: