BOSS 7.1.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDToKSpipipi Class Reference

#include <EvtDToKSpipipi.hh>

+ Inheritance diagram for EvtDToKSpipipi:

Public Member Functions

 EvtDToKSpipipi ()
 
virtual ~EvtDToKSpipipi ()
 
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 void getName (std::string &name)=0
 
virtual void decay (EvtParticle *p)=0
 
virtual void makeDecay (EvtParticle *p)=0
 
virtual EvtDecayBaseclone ()=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 ()
 
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 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*/
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}
EvtId * getDaugs()
Definition: EvtDecayBase.hh:65
void setProb(double prob)
Definition: EvtDecayProb.hh:34
const EvtVector4R & getP4() const
Definition: EvtParticle.cc:120
EvtParticle * getDaug(int i)
Definition: EvtParticle.cc:84
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
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);
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}
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 EvtDToKSpipipi::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 155 of file EvtDToKSpipipi.cc.

155 {
156 setProbMax(3700.0);
157}
void setProbMax(double prbmx)

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