BOSS 7.1.2
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDsToKKpipipi Class Reference

#include <EvtDsToKKpipipi.hh>

+ Inheritance diagram for EvtDsToKKpipipi:

Public Member Functions

 EvtDsToKKpipipi ()
 
virtual ~EvtDsToKKpipipi ()
 
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 ()
 
- 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 ()
 
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)
 
- Protected Member Functions inherited from EvtDecayBase
bool daugsDecayedByParentModel ()
 
- Protected Attributes inherited from EvtDecayBase
bool _daugsDecayedByParentModel
 

Detailed Description

Definition at line 8 of file EvtDsToKKpipipi.hh.

Constructor & Destructor Documentation

◆ EvtDsToKKpipipi()

EvtDsToKKpipipi::EvtDsToKKpipipi ( )
inline

Definition at line 12 of file EvtDsToKKpipipi.hh.

12{}

Referenced by clone().

◆ ~EvtDsToKKpipipi()

EvtDsToKKpipipi::~EvtDsToKKpipipi ( )
virtual

Definition at line 34 of file EvtDsToKKpipipi.cc.

34{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDsToKKpipipi::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 40 of file EvtDsToKKpipipi.cc.

40 {
41 return new EvtDsToKKpipipi;
42}

◆ decay()

void EvtDsToKKpipipi::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 141 of file EvtDsToKKpipipi.cc.

141 {
142/*
143 double maxprob = 0.0;
144 for(int ir=0;ir<=60000000;ir++){
145 p->initializePhaseSpace(getNDaug(),getDaugs());
146 EvtVector4R _km = p->getDaug(0)->getP4();
147 EvtVector4R _kp = p->getDaug(1)->getP4();
148 EvtVector4R _pip1 = p->getDaug(2)->getP4();
149 EvtVector4R _pip2 = p->getDaug(3)->getP4();
150 EvtVector4R _pim = p->getDaug(4)->getP4();
151
152 double _Km[4], _Kp[4], _Pip1[4], _Pip2[4], _Pim[4];
153 _Km[0] = _km.get(0); _Kp[0] = _kp.get(0); _Pip1[0] = _pip1.get(0); _Pip2[0] = _pip2.get(0); _Pim[0] = _pim.get(0);
154 _Km[1] = _km.get(1); _Kp[1] = _kp.get(1); _Pip1[1] = _pip1.get(1); _Pip2[1] = _pip2.get(1); _Pim[1] = _pim.get(1);
155 _Km[2] = _km.get(2); _Kp[2] = _kp.get(2); _Pip1[2] = _pip1.get(2); _Pip2[2] = _pip2.get(2); _Pim[2] = _pim.get(2);
156 _Km[3] = _km.get(3); _Kp[3] = _kp.get(3); _Pip1[3] = _pip1.get(3); _Pip2[3] = _pip2.get(3); _Pim[3] = _pim.get(3);
157
158 double _prob;
159 int nstates=3;
160 int g0[3]={1,1,1};
161 int g1[3]={1,1,1};
162 int g2[3]={0,0,0};
163 int g3[3]={0,1,0};
164 int g4[3]={0,0,0};
165 calEvaMy(_Km,_Kp,_Pip1,_Pip2,_Pim,mass1,mass2,mass3,width1,width2,width3,rho,phi,modetype,nstates,_prob);
166 if(_prob>maxprob) {
167 maxprob=_prob;
168 cout << "Max PDF = " << ir << " " << _prob << endl;
169 }
170 }
171 cout << "Max!!!!!!!!!!! " << maxprob<< endl;
172*/
174 EvtVector4R km = p->getDaug(0)->getP4();
175 EvtVector4R kp = p->getDaug(1)->getP4();
176 EvtVector4R pip1 = p->getDaug(2)->getP4();
177 EvtVector4R pip2 = p->getDaug(3)->getP4();
178 EvtVector4R pim = p->getDaug(4)->getP4();
179
180 double Km[4],Kp[4], Pip1[4],Pip2[4],Pim[4];
181 Km[0] = km.get(0); Kp[0] = kp.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0); Pim[0] = pim.get(0);
182 Km[1] = km.get(1); Kp[1] = kp.get(1); Pip1[1] = pip1.get(1); Pip2[1] = pip2.get(1); Pim[1] = pim.get(1);
183 Km[2] = km.get(2); Kp[2] = kp.get(2); Pip1[2] = pip1.get(2); Pip2[2] = pip2.get(2); Pim[2] = pim.get(2);
184 Km[3] = km.get(3); Kp[3] = kp.get(3); Pip1[3] = pip1.get(3); Pip2[3] = pip2.get(3); Pim[3] = pim.get(3);
185
186 double prob;
187 int nstates=3;
188 calEvaMy(Km,Kp,Pip1,Pip2,Pim,mass1,mass2,mass3,width1,width2,width3,rho,phi,modetype,nstates,prob);
189
190 setProb(prob);
191
192 return;
193}
EvtId * getDaugs()
void setProb(double prob)
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
double get(int i) const

◆ getName()

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

Implements EvtDecayBase.

Definition at line 36 of file EvtDsToKKpipipi.cc.

36 {
37 model_name="DsToKKpipipi";
38}

◆ init()

void EvtDsToKKpipipi::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 44 of file EvtDsToKKpipipi.cc.

44 {
45 // check that there are 0 arguments
46 checkNArg(0);
47 checkNDaug(5);
48
55
56 int mode = 0;
57//--------------------amp---------------
58 phi[1] = 1.4732;
59 phi[2] = -4.2901;
60 rho[1] = 0.23396;
61 rho[2] = 7.4744;
62//--------------------mass---------------------
63 mass1[0] = 1.019461;
64 mass1[1] = 1.019461;
65 mass1[2] = 1.019461;
66 mass2[0] = 0.77526;
67 mass2[1] = 0.77526;
68 mass2[2] = 0.77526;
69 mass3[0] = 1.23;
70 mass3[1] = 1.23;
71 mass3[2] = 1.23;
72
73 width1[0] = 0.004249;
74 width1[1] = 0.004249;
75 width1[2] = 0.004249;
76 width2[0] = 0.1478;
77 width2[1] = 0.1478;
78 width2[2] = 0.1478;
79 width3[0] = 0.42;
80 width3[1] = 0.42;
81 width3[2] = 0.42;
82
83 modetype[0]= 1;
84 modetype[1]= 3;
85 modetype[2]= 13;
86
87 phi[0] = 0;
88 rho[0] = 1;
89
90 //cout << "DsToKKpipipi :" << endl;
91 //for (int i=0; i<3; i++) {
92 // cout << i << " rho= " << rho[i] << " phi= " << phi[i] << endl;
93 //}
94
95 mass_Pion = 0.13957;
96 mass_Pion_N = 0.134977;
97 mass_Eta = 0.547862;
98 math_pi = 3.1415926;
99 rD2 = 25.0; // 5*5
100 rRes2 = 9.0; // 3*3
101 GS1 = 0.636619783;
102 GS2 = 0.01860182466;
103 GS3 = 0.1591549458; // 1/(2*math_2pi)
104 GS4 = 0.00620060822; // mass_Pion2/math_pi
105
106 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
107 int EE[4][4][4][4] =
108 { { {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
109 {{0,0,0,0}, {0,0,0,0}, {0,0,0,1}, {0,0,-1,0}},
110 {{0,0,0,0}, {0,0,0,-1}, {0,0,0,0}, {0,1,0,0} },
111 {{0,0,0,0}, {0,0,1,0}, {0,-1,0,0}, {0,0,0,0} } },
112 { {{0,0,0,0}, {0,0,0,0}, {0,0,0,-1}, {0,0,1,0} },
113 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
114 {{0,0,0,1}, {0,0,0,0}, {0,0,0,0}, {-1,0,0,0}},
115 {{0,0,-1,0}, {0,0,0,0}, {1,0,0,0}, {0,0,0,0} } },
116 { {{0,0,0,0}, {0,0,0,1}, {0,0,0,0}, {0,-1,0,0}},
117 {{0,0,0,-1}, {0,0,0,0}, {0,0,0,0}, {1,0,0,0} },
118 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
119 {{0,1,0,0}, {-1,0,0,0}, {0,0,0,0}, {0,0,0,0} } },
120 { {{0,0,0,0}, {0,0,-1,0}, {0,1,0,0}, {0,0,0,0} },
121 {{0,0,1,0}, {0,0,0,0}, {-1,0,0,0}, {0,0,0,0} },
122 {{0,-1,0,0}, {1,0,0,0}, {0,0,0,0}, {0,0,0,0} },
123 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} } } };
124 for (int i=0; i<4; i++) {
125 for (int j=0; j<4; j++) {
126 G[i][j] = GG[i][j];
127 for (int k=0; k<4; k++) {
128 for (int l=0; l<4; l++) {
129 E[i][j][k][l] = EE[i][j][k][l];
130 }
131 }
132 }
133 }
134
135}
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
void checkSpinParent(EvtSpinType::spintype sp)
void checkNDaug(int d1, int d2=-1)
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)

◆ initProbMax()

void EvtDsToKKpipipi::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 137 of file EvtDsToKKpipipi.cc.

137 {
138 setProbMax(585000.0);
139}
void setProbMax(double prbmx)

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