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

#include <EvtDsToKSKpPipPim.hh>

+ Inheritance diagram for EvtDsToKSKpPipPim:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtDsToKSKpPipPim()

EvtDsToKSKpPipPim::EvtDsToKSKpPipPim ( )
inline

Definition at line 12 of file EvtDsToKSKpPipPim.hh.

12{}

Referenced by clone().

◆ ~EvtDsToKSKpPipPim()

EvtDsToKSKpPipPim::~EvtDsToKSKpPipPim ( )
virtual

Definition at line 33 of file EvtDsToKSKpPipPim.cc.

33{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDsToKSKpPipPim::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 39 of file EvtDsToKSKpPipPim.cc.

39 {
40 return new EvtDsToKSKpPipPim;
41}

◆ decay()

void EvtDsToKSKpPipPim::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 156 of file EvtDsToKSKpPipPim.cc.

156 {
157/*
158 double maxprob=0;
159 for(int ir=0;ir<=60000000;ir++){
160 p->initializePhaseSpace(getNDaug(),getDaugs());
161 EvtVector4R _ks = p->getDaug(0)->getP4();
162 EvtVector4R _kp = p->getDaug(1)->getP4();
163 EvtVector4R _pip = p->getDaug(2)->getP4();
164 EvtVector4R _pim = p->getDaug(3)->getP4();
165
166 double _Pip[4],_Pim[4],_Ks[4],_Kp[4];
167 _Ks[0] = _ks.get(0); _Kp[0] = _kp.get(0); _Pip[0] = _pip.get(0); _Pim[0] = _pim.get(0);
168 _Ks[1] = _ks.get(1); _Kp[1] = _kp.get(1); _Pip[1] = _pip.get(1); _Pim[1] = _pim.get(1);
169 _Ks[2] = _ks.get(2); _Kp[2] = _kp.get(2); _Pip[2] = _pip.get(2); _Pim[2] = _pim.get(2);
170 _Ks[3] = _ks.get(3); _Kp[3] = _kp.get(3); _Pip[3] = _pip.get(3); _Pim[3] = _pim.get(3);
171
172 double _prob;
173 int g0[8]={1,1,1,1,1,1,1,1};
174 int nstates=10;
175
176 calEvaMy(_Ks, _Kp, _Pip, _Pim, mass, width, rho, phi, g0, modetype, nstates, _prob);
177 if(_prob>maxprob) {
178 maxprob=_prob;
179 printf("ir = %d,maxprob = %.10f,prob = %.10f\n\n",ir,maxprob,_prob);
180 }
181 }
182 printf("maxprob = %.10f\n", maxprob);
183*/
184 p->initializePhaseSpace(getNDaug(),getDaugs());
185 EvtVector4R ks = p->getDaug(0)->getP4();
186 EvtVector4R kp = p->getDaug(1)->getP4();
187 EvtVector4R pip = p->getDaug(2)->getP4();
188 EvtVector4R pim = p->getDaug(3)->getP4();
189
190 double Pip[4],Pim[4],Kp[4],Ks[4];
191 Ks[0] = ks.get(0); Kp[0] = kp.get(0); Pip[0] = pip.get(0); Pim[0] = pim.get(0);
192 Ks[1] = ks.get(1); Kp[1] = kp.get(1); Pip[1] = pip.get(1); Pim[1] = pim.get(1);
193 Ks[2] = ks.get(2); Kp[2] = kp.get(2); Pip[2] = pip.get(2); Pim[2] = pim.get(2);
194 Ks[3] = ks.get(3); Kp[3] = kp.get(3); Pip[3] = pip.get(3); Pim[3] = pim.get(3);
195
196 //Ks[0] = 0.6096125554; Ks[1] = 0.2036145752; Ks[2] = -0.2256821642; Ks[3] = 0.1884652565;
197 //Kp[0] = 0.7106732604; Kp[1] = -0.2599714599; Kp[2] = 0.3643877399; Kp[3] = 0.2469330230;
198 //Pip[0] = 0.4177584141; Pip[1] = -0.0562004660; Pip[2] = 0.2821155804; Pip[3] =-0.0020089709;
199 //Pim[0] = 0.3945120988; Pim[1] = -0.1225632042; Pim[2] = -0.3645928242; Pim[3] = 0.0521817527;
200
201 //Ks[0] = 0.6439438863; Ks[1] = 0.3368405824; Ks[2] = 0.2308441601; Ks[3] = 0.0171298406;
202 //Kp[0] = 0.7146391007; Kp[1] = 0.4061832461; Kp[2] = -0.1841669581; Kp[3] = 0.2609401579;
203 //Pip[0] = 0.2817222681; Pip[1] = -0.1005161765;Pip[2] = -0.1172908044;Pip[3] = -0.1898077097;
204 //Pim[0] = 0.3638598699;Pim[1] = -0.3106229764;Pim[2] = 0.0905020525;Pim[3] = 0.0907574503;
205
206 double value;
207 int g0[8]={1,1,1,1,1,1,1,1};
208 int nstates=8;
209
210 calEvaMy(Ks,Kp,Pip,Pim,mass,width,rho,phi,g0,modetype,nstates,value);
211 setProb(value);
212
213 return;
214}
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 EvtDsToKSKpPipPim::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 35 of file EvtDsToKSKpPipPim.cc.

35 {
36 model_name="DsToKSKpPipPim";
37}

◆ init()

void EvtDsToKSKpPipPim::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 43 of file EvtDsToKSKpPipPim.cc.

43 {
44 // check that there are 0 arguments
45 checkNArg(0);
46 checkNDaug(4);
47
48 checkSpinParent(EvtSpinType::SCALAR);
49 checkSpinDaughter(0,EvtSpinType::SCALAR);
50 checkSpinDaughter(1,EvtSpinType::SCALAR);
51 checkSpinDaughter(2,EvtSpinType::SCALAR);
52 checkSpinDaughter(3,EvtSpinType::SCALAR);
53
54 int mode = 0;
55 phi[1] = -3.1173e-03;
56 phi[2] = 5.6119e+00;
57 phi[3] = -7.0591e-01;
58 phi[4] = -3.1220e+00;
59 phi[5] = 3.5914e+00;
60 phi[6] = 3.8510e+00;
61 phi[7] = -1.0830e+00;
62
63 rho[1] = 1.8516e+00;
64 rho[2] = 2.8585e+00;
65 rho[3] = 1.1065e+00;
66 rho[4] = 5.4084e-01;
67 rho[5] = 1.3816e+00;
68 rho[6] = 1.9177e+00;
69 rho[7] = 2.2119e+00;
70 modetype[0]= 7;
71 modetype[1]= 4;
72 modetype[2]= 3;
73 modetype[3]= 15;
74 modetype[4]= 8;
75 modetype[5]= 11;
76 modetype[6]= 12;
77 modetype[7]= 17;
78 phi[0] = 0;
79 rho[0] = 1;
80
81 //cout << "DsToKSKpPipPim :" << endl;
82 //for (int i=0; i<8; i++) {
83 // cout << i << " rho= " << rho[i] << " phi= " << phi[i] << endl;
84 //}
85
86 mD = 1.86486;
87 metap = 0.95778;
88 mk0 = 0.497614;
89 mass_Kaon = 0.49368;
90 mass_Eta = 0.547862;
91 mass_Pion = 0.13957;
92 math_pi = 3.1415926;
93 mass_Pion2 = 0.0194797849;
94 mass_2Pion = 0.27914;
95 math_2pi = 6.2831852;
96
97 rD2 = 25.0;
98 rRes2 = 9.0;
99 GK1270 = 8.7000e-02;
100 GK1400 = 1.7400e-01;
101 GKst0 = 4.7300e-02;
102 GKstp = 5.0300e-02;
103 GP1680 = 1.5000e-01;
104 GRho = 1.4780e-01;
105 Ga0_980 = 5.0000e-02;
106 Geta1475= 9.0000e-02;
107 mK1270 = 1.2720e+00;
108 mK1400 = 1.4030e+00;
109 mKst0 = 8.9555e-01;
110 mKstp = 8.9176e-01;
111 mP1680 = 1.6800e+00;
112 mRho = 7.7526e-01;
113 ma0_980 = 9.9000e-01;
114 meta1475= 1.4750e+00;
115
116 GS1 = 0.636619783;
117 GS2 = 0.01860182466;
118 GS3 = 0.1591549458; // 1/(2*math_2pi)
119 GS4 = 0.00620060822; // mass_Pion2/math_pi
120
121 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
122 int EE[4][4][4][4] =
123 { { {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
124 {{0,0,0,0}, {0,0,0,0}, {0,0,0,1}, {0,0,-1,0}},
125 {{0,0,0,0}, {0,0,0,-1}, {0,0,0,0}, {0,1,0,0} },
126 {{0,0,0,0}, {0,0,1,0}, {0,-1,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,0}, {0,0,0,0}, {0,0,0,0} },
129 {{0,0,0,1}, {0,0,0,0}, {0,0,0,0}, {-1,0,0,0}},
130 {{0,0,-1,0}, {0,0,0,0}, {1,0,0,0}, {0,0,0,0} } },
131 { {{0,0,0,0}, {0,0,0,1}, {0,0,0,0}, {0,-1,0,0}},
132 {{0,0,0,-1}, {0,0,0,0}, {0,0,0,0}, {1,0,0,0} },
133 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
134 {{0,1,0,0}, {-1,0,0,0}, {0,0,0,0}, {0,0,0,0} } },
135 { {{0,0,0,0}, {0,0,-1,0}, {0,1,0,0}, {0,0,0,0} },
136 {{0,0,1,0}, {0,0,0,0}, {-1,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,0,0}, {0,0,0,0}, {0,0,0,0} } } };
139 for (int i=0; i<4; i++) {
140 for (int j=0; j<4; j++) {
141 G[i][j] = GG[i][j];
142 for (int k=0; k<4; k++) {
143 for (int l=0; l<4; l++) {
144 E[i][j][k][l] = EE[i][j][k][l];
145 }
146 }
147 }
148 }
149
150}
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 EvtDsToKSKpPipPim::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 152 of file EvtDsToKSKpPipPim.cc.

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