BOSS 7.0.6
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 30 of file EvtDToKSpipipi.cc.

30{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDToKSpipipi::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 36 of file EvtDToKSpipipi.cc.

36 {
37 return new EvtDToKSpipipi;
38}

◆ decay()

void EvtDToKSpipipi::decay ( EvtParticle p)
virtual

Implements EvtDecayBase.

Definition at line 158 of file EvtDToKSpipipi.cc.

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

32 {
33 model_name="DToKSpipipi";
34}

◆ init()

void EvtDToKSpipipi::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 40 of file EvtDToKSpipipi.cc.

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

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

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