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

#include <EvtDToKSpipipi.hh>

+ Inheritance diagram for EvtDToKSpipipi:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtDToKSpipipi()

EvtDToKSpipipi::EvtDToKSpipipi ( )
inline

Definition at line 11 of file EvtDToKSpipipi.hh.

11{}

Referenced by clone().

◆ ~EvtDToKSpipipi()

EvtDToKSpipipi::~EvtDToKSpipipi ( )
virtual

Definition at line 31 of file EvtDToKSpipipi.cc.

31{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDToKSpipipi::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 37 of file EvtDToKSpipipi.cc.

37 {
38 return new EvtDToKSpipipi;
39}

◆ decay()

void EvtDToKSpipipi::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 159 of file EvtDToKSpipipi.cc.

159 {
160/*
161 double maxprob = 0.0;
162 for(int ir=0;ir<=60000000;ir++){
163 p->initializePhaseSpace(getNDaug(),getDaugs());
164 EvtVector4R Ks0 = p->getDaug(0)->getP4();
165 EvtVector4R pi1 = p->getDaug(1)->getP4();
166 EvtVector4R pi2 = p->getDaug(2)->getP4();
167 EvtVector4R pi3 = p->getDaug(3)->getP4();
168 double Ks[4],Pip1[4],Pip2[4],Pim[4];
169 Ks[0] = Ks0.get(0); Pip1[0] = pi1.get(0); Pip2[0] = pi2.get(0); Pim[0] = pi3.get(0);
170 Ks[1] = Ks0.get(1); Pip1[1] = pi1.get(1); Pip2[1] = pi2.get(1); Pim[1] = pi3.get(1);
171 Ks[2] = Ks0.get(2); Pip1[2] = pi1.get(2); Pip2[2] = pi2.get(2); Pim[2] = pi3.get(2);
172 Ks[3] = Ks0.get(3); Pip1[3] = pi1.get(3); Pip2[3] = pi2.get(3); Pim[3] = pi3.get(3);
173 double Prob = calPDF(Ks, Pip1, Pip2, Pim);
174 if(Prob>maxprob) {
175 maxprob=Prob;
176 std::cout << "Max PDF = " << ir << " prob= " << Prob << std::endl;
177 }
178 }
179 std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
180*/
181 p->initializePhaseSpace(getNDaug(),getDaugs());
182 EvtVector4R Ks0 = p->getDaug(0)->getP4();
183 EvtVector4R pi1 = p->getDaug(1)->getP4();
184 EvtVector4R pi2 = p->getDaug(2)->getP4();
185 EvtVector4R pi3 = p->getDaug(3)->getP4();
186
187 double Ks[4],Pip1[4],Pip2[4],Pim[4];
188 Ks[0] = Ks0.get(0); Pip1[0] = pi1.get(0); Pip2[0] = pi2.get(0); Pim[0] = pi3.get(0);
189 Ks[1] = Ks0.get(1); Pip1[1] = pi1.get(1); Pip2[1] = pi2.get(1); Pim[1] = pi3.get(1);
190 Ks[2] = Ks0.get(2); Pip1[2] = pi1.get(2); Pip2[2] = pi2.get(2); Pim[2] = pi3.get(2);
191 Ks[3] = Ks0.get(3); Pip1[3] = pi1.get(3); Pip2[3] = pi2.get(3); Pim[3] = pi3.get(3);
192 double prob = calPDF(Ks, Pip1, Pip2, Pim);
193 setProb(prob);
194 return;
195}
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 EvtDToKSpipipi::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 33 of file EvtDToKSpipipi.cc.

33 {
34 model_name="DToKSpipipi";
35}

◆ init()

void EvtDToKSpipipi::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 41 of file EvtDToKSpipipi.cc.

41 {
42 checkNArg(0);
43 checkNDaug(4);
44 checkSpinParent(EvtSpinType::SCALAR);
45/*
46 checkSpinDaughter(0,EvtSpinType::SCALAR);
47 checkSpinDaughter(1,EvtSpinType::SCALAR);
48 checkSpinDaughter(3,EvtSpinType::SCALAR);
49 checkSpinDaughter(4,EvtSpinType::SCALAR);
50*/
51 //std::cout << "Initializing EvtDToKSpipipi" << std::endl;
52
53 mK1270 = 1.272; mK1400 = 1.403;
54 GK1270 = 0.09; GK1400 = 0.174;
55 mKstr = 0.89166; mrho = 0.77549;
56 GKstr = 0.0462; Grho = 0.1491;
57 msigma = 0.472;
58 Gsigma = 0.542;
59 phi_omega = -0.02;
60 mK1650 = 1.65;
61 GK1650 = 0.158;
62 rho[0] = 1.0;
63 phi[0] = 0.0;
64
65 ma1 = 1.22;
66 Ga1 = 0.4282;
67 mK1460 = 1.4152;
68 GK1460 = 0.2485;
69 rho_omega = 0.00294;
70
71 phi[1] = -1.55;
72 rho[1] = 0.5843;
73
74 phi[2] = -1.8223;
75 rho[2] = 2.0974;
76
77 phi[3] = -2.6751;
78 rho[3] = 0.46642;
79
80 phi[4] = -2.2429;
81 rho[4] = 0.33334;
82
83 phi[5] = -0.55888;
84 rho[5] = 0.15549;
85
86 phi[6] = -1.8778;
87 rho[6] = 0.94452;
88
89 phi[7] = 2.7724;
90 rho[7] = 0.99411;
91
92 phi[8] = -2.6461;
93 rho[8] = 0.21231;
94
95 phi[9] = -0.95137;
96 rho[9] = 0.29895;
97
98 phi[10] = -3.0843;
99 rho[10] = 3.6361;
100
101 phi[11] = 2.0954;
102 rho[11] = 0.96472;
103
104 phi[12] = -2.4965;
105 rho[12] = 0.48470;
106
107 //for (int i=0; i<13; i++) {
108 // cout << i << " rho,phi = " << rho[i] << ", "<< phi[i] << endl;
109 //}
110
111 mD = 1.86486;
112 rD = 5;
113 metap = 0.95778;
114 mkstr = 0.89594;
115 mk0 = 0.497614;
116 mass_Kaon = 0.49368;
117 mass_Pion = 0.13957;
118 math_pi = 3.1415926;
119
120 pi = 3.1415926;
121 mpi = 0.13957;
122 g1 = 0.5468;
123 g2 = 0.23;
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 for (int j=0; j<4; j++) {
145 G[i][j] = GG[i][j];
146 for (int k=0; k<4; k++) {
147 for (int l=0; l<4; l++) {
148 E[i][j][k][l] = EE[i][j][k][l];
149 }
150 }
151 }
152 }
153}
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 EvtDToKSpipipi::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 155 of file EvtDToKSpipipi.cc.

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