BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDsToKSKmPipPip Class Reference

#include <EvtDsToKSKmPipPip.hh>

Inheritance diagram for EvtDsToKSKmPipPip:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtDsToKSKmPipPip()

EvtDsToKSKmPipPip::EvtDsToKSKmPipPip ( )
inline

Definition at line 11 of file EvtDsToKSKmPipPip.hh.

11{}

Referenced by clone().

◆ ~EvtDsToKSKmPipPip()

EvtDsToKSKmPipPip::~EvtDsToKSKmPipPip ( )
virtual

Definition at line 33 of file EvtDsToKSKmPipPip.cc.

33{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDsToKSKmPipPip::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 37 of file EvtDsToKSKmPipPip.cc.

37{ return new EvtDsToKSKmPipPip; }

◆ decay()

void EvtDsToKSKmPipPip::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 162 of file EvtDsToKSKmPipPip.cc.

162 {
163 /*
164 double maxprob=0;
165 for(int ir=0;ir<=60000000;ir++){
166 p->initializePhaseSpace(getNDaug(),getDaugs());
167 EvtVector4R _ks = p->getDaug(0)->getP4();
168 EvtVector4R _km = p->getDaug(1)->getP4();
169 EvtVector4R _pip = p->getDaug(2)->getP4();
170 EvtVector4R _pip2 = p->getDaug(3)->getP4();
171
172 double _Pip[4],_Pip2[4],_Ks[4],_Km[4];
173 _Ks[0] = _ks.get(0); _Km[0] = _km.get(0); _Pip[0] = _pip.get(0); _Pip2[0] =
174 _pip2.get(0); _Ks[1] = _ks.get(1); _Km[1] = _km.get(1); _Pip[1] = _pip.get(1); _Pip2[1] =
175 _pip2.get(1); _Ks[2] = _ks.get(2); _Km[2] = _km.get(2); _Pip[2] = _pip.get(2); _Pip2[2] =
176 _pip2.get(2); _Ks[3] = _ks.get(3); _Km[3] = _km.get(3); _Pip[3] = _pip.get(3); _Pip2[3] =
177 _pip2.get(3);
178
179 double _prob;
180 int g0[11]={1,1,1,1,1,1,1,1,1,1,1};
181 int nstates=11;
182
183 calEvaMy(_Ks, _Km, _Pip, _Pip2, mass, width, rho, phi, g0, modetype, nstates, _prob);
184 if(_prob>maxprob) {
185 maxprob=_prob;
186 printf("ir = %d,maxprob = %.10f,prob = %.10f\n\n",ir,maxprob,_prob);
187 }
188 }
189 printf("maxprob = %.10f\n", maxprob);*/
190
191 p->initializePhaseSpace( getNDaug(), getDaugs() );
192 EvtVector4R ks = p->getDaug( 0 )->getP4();
193 EvtVector4R km = p->getDaug( 1 )->getP4();
194 EvtVector4R pip = p->getDaug( 2 )->getP4();
195 EvtVector4R pip2 = p->getDaug( 3 )->getP4();
196
197 double Pip[4], Pip2[4], Km[4], Ks[4];
198 Ks[0] = ks.get( 0 );
199 Km[0] = km.get( 0 );
200 Pip[0] = pip.get( 0 );
201 Pip2[0] = pip2.get( 0 );
202 Ks[1] = ks.get( 1 );
203 Km[1] = km.get( 1 );
204 Pip[1] = pip.get( 1 );
205 Pip2[1] = pip2.get( 1 );
206 Ks[2] = ks.get( 2 );
207 Km[2] = km.get( 2 );
208 Pip[2] = pip.get( 2 );
209 Pip2[2] = pip2.get( 2 );
210 Ks[3] = ks.get( 3 );
211 Km[3] = km.get( 3 );
212 Pip[3] = pip.get( 3 );
213 Pip2[3] = pip2.get( 3 );
214
215 // Ks[0] = 0.6096125554; Ks[1] = 0.2036145752; Ks[2] = -0.2256821642; Ks[3] =
216 // 0.1884652565 ; Km[0] = 0.7106732604; Km[1] = -0.2599714599; Km[2] =
217 // 0.3643877399; Km[3] = 0.2469330230 ; Pip[0] = 0.4177584141; Pip[1] =
218 // -0.0562004660; Pip[2] = 0.2821155804; Pip[3] =-0.0020089709 ; Pip2[0] =
219 // 0.3945120988; Pip2[1] = -0.1225632042; Pip2[2] = -0.3645928242; Pip2[3] = 0.0521817527
220 // ;
221
222 // Ks[0] = 0.6439438863; Ks[1] = 0.3368405824; Ks[2] = 0.2308441601; Ks[3] = 0.0171298406;
223 // Km[0] = 0.7146391007; Km[1] = 0.4061832461; Km[2] = -0.1841669581; Km[3] = 0.2609401579;
224 // Pip[0] = 0.2817222681; Pip[1] = -0.1005161765;Pip[2] = -0.1172908044;Pip[3] =
225 // -0.1898077097; Pip2[0] = 0.3638598699;Pip2[1] = -0.3106229764;Pip2[2] =
226 // 0.0905020525;Pip2[3] = 0.0907574503;
227
228 double value;
229 int g0[11] = { 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 };
230 int nstates = 11;
231
232 calEvaMy( Ks, Km, Pip, Pip2, mass, width, rho, phi, g0, modetype, nstates, value );
233 setProb( value );
234
235 return;
236}
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition EvtDecayBase.hh:68
EvtDecayProb::setProb
void setProb(double prob)
Definition EvtDecayProb.hh:40
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition EvtParticle.cc:108
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition EvtParticle.cc:80
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition EvtParticle.cc:1042
EvtVector4R::get
double get(int i) const
Definition EvtVector4R.hh:171

◆ getName()

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

Implements EvtDecayBase.

Definition at line 35 of file EvtDsToKSKmPipPip.cc.

35{ model_name = "DsToKSKmPipPip"; }

◆ init()

void EvtDsToKSKmPipPip::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 39 of file EvtDsToKSKmPipPip.cc.

39 {
40 // check that there are 0 arguments
41 checkNArg( 0 );
42 checkNDaug( 4 );
43
44 checkSpinParent( EvtSpinType::SCALAR );
45 checkSpinDaughter( 0, EvtSpinType::SCALAR );
46 checkSpinDaughter( 1, EvtSpinType::SCALAR );
47 checkSpinDaughter( 2, EvtSpinType::SCALAR );
48 checkSpinDaughter( 3, EvtSpinType::SCALAR );
49
50 int mode = 0;
51 phi[1] = -1.5871e+00;
52 phi[2] = -6.4308e+00;
53 phi[3] = -4.4006e+00;
54 phi[4] = 4.7540e+00;
55 phi[5] = 4.2785e+00;
56 phi[6] = 6.3523e+00;
57 phi[7] = 6.3523e+00;
58 phi[8] = 5.4902e+00;
59 phi[9] = 5.2024e+00;
60 phi[10] = -3.9776e+00;
61
62 rho[1] = 3.6108e-01;
63 rho[2] = 4.5990e-01;
64 rho[3] = 1.7722e+00;
65 rho[4] = 2.1434e+00;
66 rho[5] = 1.2504e+00;
67 rho[6] = 1.0069e+00;
68 rho[7] = 1.0069e+00;
69 rho[8] = 1.7351e+00;
70 rho[9] = 4.1005e+00;
71 rho[10] = 3.7821e+00;
72
73 modetype[0] = 1;
74 modetype[1] = 2;
75 modetype[2] = 3;
76 modetype[3] = 21;
77 modetype[4] = 22;
78 modetype[5] = 8;
79 modetype[6] = 4;
80 modetype[7] = 14;
81 modetype[8] = 17;
82 modetype[9] = 5;
83 modetype[10] = 20;
84
85 phi[0] = 0;
86 rho[0] = 1;
87
88 // cout << "DsToKSKmPipPip :" << endl;
89 // for (int i=0; i<11; i++) {
90 // cout << i << " rho= " << rho[i] << " phi= " << phi[i] << endl;
91 // }
92
93 G1510 = 1.0000e+00;
94 GKst0 = 4.7300e-02;
95 GKstp = 5.0300e-02;
96 Ga0_980 = 5.0000e-02;
97 Geta1475 = 9.0000e-02;
98 Gf1285 = 2.2700e-02;
99 m1510 = 1.6200e+00;
100 mKst0 = 8.9555e-01;
101 mKstp = 8.9176e-01;
102 ma0_980 = 9.9000e-01;
103 meta1475 = 1.4750e+00;
104 mf1285 = 1.2819e+00;
105
106 mD = 1.86486;
107 metap = 0.95778;
108 mk0 = 0.497614;
109 mass_Kaon = 0.49368;
110 mass_Eta = 0.547862;
111 mass_Pion = 0.13957;
112 math_pi = 3.1415926;
113 mass_Pion2 = 0.0194797849;
114 mass_2Pion = 0.27914;
115 math_2pi = 6.2831852;
116 rD2 = 25.0; // 5*5
117 rRes2 = 9.0; // 3*3
118 g1 = 0.5468;
119 g2 = 0.23;
120 GS1 = 0.636619783;
121 GS2 = 0.01860182466;
122 GS3 = 0.1591549458; // 1/(2*math_2pi)
123 GS4 = 0.00620060822; // mass_Pion2/math_pi
124
125 int GG[4][4] = { { 1, 0, 0, 0 }, { 0, -1, 0, 0 }, { 0, 0, -1, 0 }, { 0, 0, 0, -1 } };
126 int EE[4][4][4][4] = {
127 { { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
128 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 1 }, { 0, 0, -1, 0 } },
129 { { 0, 0, 0, 0 }, { 0, 0, 0, -1 }, { 0, 0, 0, 0 }, { 0, 1, 0, 0 } },
130 { { 0, 0, 0, 0 }, { 0, 0, 1, 0 }, { 0, -1, 0, 0 }, { 0, 0, 0, 0 } } },
131 { { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, -1 }, { 0, 0, 1, 0 } },
132 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
133 { { 0, 0, 0, 1 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { -1, 0, 0, 0 } },
134 { { 0, 0, -1, 0 }, { 0, 0, 0, 0 }, { 1, 0, 0, 0 }, { 0, 0, 0, 0 } } },
135 { { { 0, 0, 0, 0 }, { 0, 0, 0, 1 }, { 0, 0, 0, 0 }, { 0, -1, 0, 0 } },
136 { { 0, 0, 0, -1 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 1, 0, 0, 0 } },
137 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
138 { { 0, 1, 0, 0 }, { -1, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } } },
139 { { { 0, 0, 0, 0 }, { 0, 0, -1, 0 }, { 0, 1, 0, 0 }, { 0, 0, 0, 0 } },
140 { { 0, 0, 1, 0 }, { 0, 0, 0, 0 }, { -1, 0, 0, 0 }, { 0, 0, 0, 0 } },
141 { { 0, -1, 0, 0 }, { 1, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
142 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } } } };
143 for ( int i = 0; i < 4; i++ )
144 {
145 for ( int j = 0; j < 4; j++ )
146 {
147 G[i][j] = GG[i][j];
148 for ( int k = 0; k < 4; k++ )
149 {
150 for ( int l = 0; l < 4; l++ ) { E[i][j][k][l] = EE[i][j][k][l]; }
151 }
152 }
153 }
154}
EvtDecayBase::checkSpinDaughter
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:492
EvtDecayBase::checkSpinParent
void checkSpinParent(EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:479
EvtDecayBase::checkNDaug
void checkNDaug(int d1, int d2=-1)
Definition EvtDecayBase.cc:465
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition EvtDecayBase.cc:449

◆ initProbMax()

void EvtDsToKSKmPipPip::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 156 of file EvtDsToKSKmPipPip.cc.

156 {
157 // setProbMax(75000);
158 setProbMax( 95050 );
159 // setProbMax(1.0);
160}
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition EvtDecayBase.cc:261
The documentation for this class was generated from the following files: