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

#include <EvtDToKmPipPipPi0.hh>

+ Inheritance diagram for EvtDToKmPipPipPi0:

Public Member Functions

 EvtDToKmPipPipPi0 ()
 
virtual ~EvtDToKmPipPipPi0 ()
 
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 EvtDToKmPipPipPi0.hh.

Constructor & Destructor Documentation

◆ EvtDToKmPipPipPi0()

EvtDToKmPipPipPi0::EvtDToKmPipPipPi0 ( )
inline

Definition at line 13 of file EvtDToKmPipPipPi0.hh.

13{}

Referenced by clone().

◆ ~EvtDToKmPipPipPi0()

EvtDToKmPipPipPi0::~EvtDToKmPipPipPi0 ( )
virtual

Definition at line 37 of file EvtDToKmPipPipPi0.cc.

37{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDToKmPipPipPi0::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 43 of file EvtDToKmPipPipPi0.cc.

43 {
44 return new EvtDToKmPipPipPi0;
45}

◆ decay()

void EvtDToKmPipPipPi0::decay ( EvtParticle *  p)
virtual

Implements EvtDecayBase.

Definition at line 148 of file EvtDToKmPipPipPi0.cc.

148 {
149//--------------max
150/*
151 double maxprob=0,maxq1270;
152 for(int ir=0;ir<=20000000;ir++){
153 p->initializePhaseSpace(getNDaug(),getDaugs());
154 EvtVector4R km = p->getDaug(0)->getP4();
155 EvtVector4R pip1 = p->getDaug(1)->getP4();
156 EvtVector4R pip2 = p->getDaug(2)->getP4();
157 EvtVector4R pi0 = p->getDaug(3)->getP4();
158
159 double Pip1[4],Pip2[4],Km[4],Pi0[4];
160 Km[0] = km.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0);Pi0[0] = pi0.get(0);
161 Km[1] = km.get(1); Pip1[1] = pip1.get(1); Pip2[1] = pip2.get(1);Pi0[1] = pi0.get(1);
162 Km[2] = km.get(2); Pip1[2] = pip1.get(2); Pip2[2] = pip2.get(2);Pi0[2] = pi0.get(2);
163 Km[3] = km.get(3); Pip1[3] = pip1.get(3); Pip2[3] = pip2.get(3);Pi0[3] = pi0.get(3);
164
165
166 double _prob;
167 double value;
168 int nstates=14;
169 int modetype[14]= {1,1,8,4,12,13,12,17,8,4,13,9,4,11};
170 int g0[14] = {1,1,1,1,1,0,0,1,1,1,0,1,1,1};
171 int g1[14] = {1,1,1,1,1,1,1,1,0,0,0,0,1,1};
172 int g2[14] = {0,1,0,0,1,1,1,0,2,0,1,0,2,1};
173 double mass1[14] = {mKst0,mKst0,mrho, mKst0, mKst0, mrho, mKst0, mKst0, mrho, mKst0, mrho, mKst0, mKst0, mKst0};
174 double mass2[14] = {mrho, mrho, mK1270,mK1400,mK1460,mK1460,mK1460,mrho, mK1270,mK1400,mK1460,mK1270,mK1400, mK1680};
175 double width1[14] = {GKst0,GKst0,Grho, GKst0, GKst0, Grho, GKst0, GKst0, Grho, GKst0, Grho, GKst0, GKst0, GKst0};
176 double width2[14] = {Grho, Grho, GK1270,GK1400,GK1460,GK1460,GK1460,Grho, GK1270,GK1400,GK1460,GK1270,GK1400, GK1680};
177
178 calEvaMy(Km,Pip1,Pip2,Pi0,mass1,mass2,width1,width2,rho,phi,g0,g1,g2,modetype,nstates,_prob);
179 if(_prob>maxprob) {
180 maxprob=_prob;
181 }
182 }
183 printf("maxprob = %.10f\n", maxprob);
184*/
185//--------------------main--------------
186
187 p->initializePhaseSpace(getNDaug(),getDaugs());
188 EvtVector4R km = p->getDaug(0)->getP4();
189 EvtVector4R pip1 = p->getDaug(1)->getP4();
190 EvtVector4R pip2 = p->getDaug(2)->getP4();
191 EvtVector4R pi0 = p->getDaug(3)->getP4();
192 int num = EvtPDL::getStdHep(p->getId());
193
194 double Pip1[4],Pip2[4],Km[4],Pi0[4];
195 //----------------------------------------------------------------------------------------
196 if(num>0){
197 Km[0] = km.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0);Pi0[0] = pi0.get(0);
198 Km[1] = km.get(1); Pip1[1] = pip1.get(1); Pip2[1] = pip2.get(1);Pi0[1] = pi0.get(1);
199 Km[2] = km.get(2); Pip1[2] = pip1.get(2); Pip2[2] = pip2.get(2);Pi0[2] = pi0.get(2);
200 Km[3] = km.get(3); Pip1[3] = pip1.get(3); Pip2[3] = pip2.get(3);Pi0[3] = pi0.get(3);
201}
202 if(num<0){
203 Km[0] = km.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0);Pi0[0] = pi0.get(0);
204 Km[1] = -km.get(1); Pip1[1] = -pip1.get(1); Pip2[1] = -pip2.get(1);Pi0[1] = -pi0.get(1);
205 Km[2] = -km.get(2); Pip1[2] = -pip1.get(2); Pip2[2] = -pip2.get(2);Pi0[2] = -pi0.get(2);
206 Km[3] = -km.get(3); Pip1[3] = -pip1.get(3); Pip2[3] = -pip2.get(3);Pi0[3] = -pi0.get(3);
207}
208 double value;
209 int nstates=14;
210 int modetype[14]= {1,1,8,4,12,13,12,17,8,4,13,9,4,11};
211 int g0[14] = {1,1,1,1,1,0,0,1,1,1,0,1,1,1};
212 int g1[14] = {1,1,1,1,1,1,1,1,0,0,0,0,1,1};
213 int g2[14] = {0,1,0,0,1,1,1,0,2,0,1,0,2,1};
214 double mass1[14] = {mKst0,mKst0,mrho, mKst0, mKst0, mrho, mKst0, mKst0, mrho, mKst0, mrho, mKst0, mKst0, mKst0};
215 double mass2[14] = {mrho, mrho, mK1270,mK1400,mK1460,mK1460,mK1460,mrho, mK1270,mK1400,mK1460,mK1270,mK1400, mK1680};
216 double width1[14] = {GKst0,GKst0,Grho, GKst0, GKst0, Grho, GKst0, GKst0, Grho, GKst0, Grho, GKst0, GKst0, GKst0};
217 double width2[14] = {Grho, Grho, GK1270,GK1400,GK1460,GK1460,GK1460,Grho, GK1270,GK1400,GK1460,GK1270,GK1400, GK1680};
218
219 calEvaMy(Km,Pip1,Pip2,Pi0,mass1,mass2,width1,width2,rho,phi,g0,g1,g2,modetype,nstates,value);
220 setProb(value);
221
222 return;
223}
g1
TF1 * g1
Definition: DataBase/tau_mode.c:66
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
rb::Grho
const double Grho
Definition: TConstant.h:12

◆ getName()

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

Implements EvtDecayBase.

Definition at line 39 of file EvtDToKmPipPipPi0.cc.

39 {
40 model_name="DToKmPipPipPi0";
41}

◆ init()

void EvtDToKmPipPipPi0::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 47 of file EvtDToKmPipPipPi0.cc.

47 {
48 // check that there are 0 arguments
49 checkNArg(0);
50 checkNDaug(4);
51
52 checkSpinParent(EvtSpinType::SCALAR);
53 checkSpinDaughter(0,EvtSpinType::SCALAR);
54 checkSpinDaughter(1,EvtSpinType::SCALAR);
55 checkSpinDaughter(2,EvtSpinType::SCALAR);
56 checkSpinDaughter(3,EvtSpinType::SCALAR);
57
58//------------------parameters-----------------
59 rho[0] = 1.00; phi[0] = 0.00;
60 rho[1] = 1.7545e-01; phi[1] = 1.4102e+00;
61 rho[2] = -5.1947e-01; phi[2] = 3.1246e+00;
62 rho[3] = 9.9483e-01; phi[3] = 4.6896e-01;
63 rho[4] = 7.1106e-01; phi[4] = 4.3110e-01;
64 rho[5] = 6.3661e+00; phi[5] = 4.9169e+00;
65 rho[6] = 1.4064e+01; phi[6] = -2.3110e+00;
66 rho[7] = 3.0275e+00; phi[7] = 2.9577e+00;
67 rho[8] = -6.9322e+00; phi[8] = 4.6902e+00;
68 rho[9] = 1.1069e+00; phi[9] = 3.4885e-01;
69 rho[10]= 5.5398e+00; phi[10]= -8.4782e-01;
70 rho[11]= -7.5436e-01; phi[11]= 1.0600e+00;
71 rho[12]= 7.1533e-01; phi[12]= -5.3175e+00;
72 rho[13]= 1.7440e+00; phi[13]= -5.1549e+00;
73 mD = 1.86486;
74
75 //cout << "Initializing DToKpipipi0" << endl;
76 //for (int i=0; i<14; i++) {
77 // cout << i << " rho,phi = " << rho[i] << ", "<< phi[i] << endl;
78 //}
79
80//------------------------new----------------------
81 mKst0 = 0.89555;
82 mrho = 0.77511;
83 mrho1450 = 1.465;
84 mK1270 = 1.289;
85 mK1400 = 1.403;
86 mK1460 = 1.482;
87 mK1650 = 1.672;
88 mK1680 = 1.718;
89
90 GKst0 = 0.0473;
91 Grho = 0.1491;
92 Grho1450 = 0.400;
93 GK1270 = 0.116;
94 GK1400 = 0.174;
95 GK1460 = 0.3356;
96 GK1650 = 0.158;
97 GK1680 = 0.322;
98
99 mass_Pion = 0.13957;
100 mass_Pion_N = 0.134977;
101 mass_Eta = 0.547862;
102 mass_Kaon = 0.493677;
103 math_pi = 3.1415926;
104
105 rD2 = 5.0; // 5*5
106 rRes1 = 3.0; // 3*3
107 rRes2 = 3.0; // 3*3
108
109 GS1 = 0.636619783;
110 GS2 = 0.01860182466;
111 GS3 = 0.1591549458; // 1/(2*math_2pi)
112 GS4 = 0.00620060822; // mass_Pion2/math_pi
113
114 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
115 int EE[4][4][4][4] =
116 { { {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
117 {{0,0,0,0}, {0,0,0,0}, {0,0,0,1}, {0,0,-1,0}},
118 {{0,0,0,0}, {0,0,0,-1}, {0,0,0,0}, {0,1,0,0} },
119 {{0,0,0,0}, {0,0,1,0}, {0,-1,0,0}, {0,0,0,0} } },
120 { {{0,0,0,0}, {0,0,0,0}, {0,0,0,-1}, {0,0,1,0} },
121 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
122 {{0,0,0,1}, {0,0,0,0}, {0,0,0,0}, {-1,0,0,0}},
123 {{0,0,-1,0}, {0,0,0,0}, {1,0,0,0}, {0,0,0,0} } },
124 { {{0,0,0,0}, {0,0,0,1}, {0,0,0,0}, {0,-1,0,0}},
125 {{0,0,0,-1}, {0,0,0,0}, {0,0,0,0}, {1,0,0,0} },
126 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
127 {{0,1,0,0}, {-1,0,0,0}, {0,0,0,0}, {0,0,0,0} } },
128 { {{0,0,0,0}, {0,0,-1,0}, {0,1,0,0}, {0,0,0,0} },
129 {{0,0,1,0}, {0,0,0,0}, {-1,0,0,0}, {0,0,0,0} },
130 {{0,-1,0,0}, {1,0,0,0}, {0,0,0,0}, {0,0,0,0} },
131 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} } } };
132 for (int i=0; i<4; i++) {
133 for (int j=0; j<4; j++) {
134 G[i][j] = GG[i][j];
135 for (int k=0; k<4; k++) {
136 for (int l=0; l<4; l++) {
137 E[i][j][k][l] = EE[i][j][k][l];
138 }
139 }
140 }
141 }
142}
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 EvtDToKmPipPipPi0::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 144 of file EvtDToKmPipPipPi0.cc.

144 {
145 setProbMax(100000);
146}
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition: EvtDecayBase.cc:297
The documentation for this class was generated from the following files: