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

#include <EvtDsToKKpipipi.hh>

+ Inheritance diagram for EvtDsToKKpipipi:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtDsToKKpipipi()

EvtDsToKKpipipi::EvtDsToKKpipipi ( )
inline

Definition at line 12 of file EvtDsToKKpipipi.hh.

12{}

Referenced by clone().

◆ ~EvtDsToKKpipipi()

EvtDsToKKpipipi::~EvtDsToKKpipipi ( )
virtual

Definition at line 34 of file EvtDsToKKpipipi.cc.

34{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDsToKKpipipi::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 40 of file EvtDsToKKpipipi.cc.

40 {
41 return new EvtDsToKKpipipi;
42}

◆ decay()

void EvtDsToKKpipipi::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 141 of file EvtDsToKKpipipi.cc.

141 {
142/*
143 double maxprob = 0.0;
144 for(int ir=0;ir<=60000000;ir++){
145 p->initializePhaseSpace(getNDaug(),getDaugs());
146 EvtVector4R _km = p->getDaug(0)->getP4();
147 EvtVector4R _kp = p->getDaug(1)->getP4();
148 EvtVector4R _pip1 = p->getDaug(2)->getP4();
149 EvtVector4R _pip2 = p->getDaug(3)->getP4();
150 EvtVector4R _pim = p->getDaug(4)->getP4();
151
152 double _Km[4], _Kp[4], _Pip1[4], _Pip2[4], _Pim[4];
153 _Km[0] = _km.get(0); _Kp[0] = _kp.get(0); _Pip1[0] = _pip1.get(0); _Pip2[0] = _pip2.get(0); _Pim[0] = _pim.get(0);
154 _Km[1] = _km.get(1); _Kp[1] = _kp.get(1); _Pip1[1] = _pip1.get(1); _Pip2[1] = _pip2.get(1); _Pim[1] = _pim.get(1);
155 _Km[2] = _km.get(2); _Kp[2] = _kp.get(2); _Pip1[2] = _pip1.get(2); _Pip2[2] = _pip2.get(2); _Pim[2] = _pim.get(2);
156 _Km[3] = _km.get(3); _Kp[3] = _kp.get(3); _Pip1[3] = _pip1.get(3); _Pip2[3] = _pip2.get(3); _Pim[3] = _pim.get(3);
157
158 double _prob;
159 int nstates=3;
160 int g0[3]={1,1,1};
161 int g1[3]={1,1,1};
162 int g2[3]={0,0,0};
163 int g3[3]={0,1,0};
164 int g4[3]={0,0,0};
165 calEvaMy(_Km,_Kp,_Pip1,_Pip2,_Pim,mass1,mass2,mass3,width1,width2,width3,rho,phi,modetype,nstates,_prob);
166 if(_prob>maxprob) {
167 maxprob=_prob;
168 cout << "Max PDF = " << ir << " " << _prob << endl;
169 }
170 }
171 cout << "Max!!!!!!!!!!! " << maxprob<< endl;
172*/
173 p->initializePhaseSpace(getNDaug(),getDaugs());
174 EvtVector4R km = p->getDaug(0)->getP4();
175 EvtVector4R kp = p->getDaug(1)->getP4();
176 EvtVector4R pip1 = p->getDaug(2)->getP4();
177 EvtVector4R pip2 = p->getDaug(3)->getP4();
178 EvtVector4R pim = p->getDaug(4)->getP4();
179
180 double Km[4],Kp[4], Pip1[4],Pip2[4],Pim[4];
181 Km[0] = km.get(0); Kp[0] = kp.get(0); Pip1[0] = pip1.get(0); Pip2[0] = pip2.get(0); Pim[0] = pim.get(0);
182 Km[1] = km.get(1); Kp[1] = kp.get(1); Pip1[1] = pip1.get(1); Pip2[1] = pip2.get(1); Pim[1] = pim.get(1);
183 Km[2] = km.get(2); Kp[2] = kp.get(2); Pip1[2] = pip1.get(2); Pip2[2] = pip2.get(2); Pim[2] = pim.get(2);
184 Km[3] = km.get(3); Kp[3] = kp.get(3); Pip1[3] = pip1.get(3); Pip2[3] = pip2.get(3); Pim[3] = pim.get(3);
185
186 double prob;
187 int nstates=3;
188 calEvaMy(Km,Kp,Pip1,Pip2,Pim,mass1,mass2,mass3,width1,width2,width3,rho,phi,modetype,nstates,prob);
189
190 setProb(prob);
191
192 return;
193}
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 EvtDsToKKpipipi::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 36 of file EvtDsToKKpipipi.cc.

36 {
37 model_name="DsToKKpipipi";
38}

◆ init()

void EvtDsToKKpipipi::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 44 of file EvtDsToKKpipipi.cc.

44 {
45 // check that there are 0 arguments
46 checkNArg(0);
47 checkNDaug(5);
48
49 checkSpinParent(EvtSpinType::SCALAR);
50 checkSpinDaughter(0,EvtSpinType::SCALAR);
51 checkSpinDaughter(1,EvtSpinType::SCALAR);
52 checkSpinDaughter(2,EvtSpinType::SCALAR);
53 checkSpinDaughter(3,EvtSpinType::SCALAR);
54 checkSpinDaughter(4,EvtSpinType::SCALAR);
55
56 int mode = 0;
57//--------------------amp---------------
58 phi[1] = 1.4732;
59 phi[2] = -4.2901;
60 rho[1] = 0.23396;
61 rho[2] = 7.4744;
62//--------------------mass---------------------
63 mass1[0] = 1.019461;
64 mass1[1] = 1.019461;
65 mass1[2] = 1.019461;
66 mass2[0] = 0.77526;
67 mass2[1] = 0.77526;
68 mass2[2] = 0.77526;
69 mass3[0] = 1.23;
70 mass3[1] = 1.23;
71 mass3[2] = 1.23;
72
73 width1[0] = 0.004249;
74 width1[1] = 0.004249;
75 width1[2] = 0.004249;
76 width2[0] = 0.1478;
77 width2[1] = 0.1478;
78 width2[2] = 0.1478;
79 width3[0] = 0.42;
80 width3[1] = 0.42;
81 width3[2] = 0.42;
82
83 modetype[0]= 1;
84 modetype[1]= 3;
85 modetype[2]= 13;
86
87 phi[0] = 0;
88 rho[0] = 1;
89
90 //cout << "DsToKKpipipi :" << endl;
91 //for (int i=0; i<3; i++) {
92 // cout << i << " rho= " << rho[i] << " phi= " << phi[i] << endl;
93 //}
94
95 mass_Pion = 0.13957;
96 mass_Pion_N = 0.134977;
97 mass_Eta = 0.547862;
98 math_pi = 3.1415926;
99 rD2 = 25.0; // 5*5
100 rRes2 = 9.0; // 3*3
101 GS1 = 0.636619783;
102 GS2 = 0.01860182466;
103 GS3 = 0.1591549458; // 1/(2*math_2pi)
104 GS4 = 0.00620060822; // mass_Pion2/math_pi
105
106 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
107 int EE[4][4][4][4] =
108 { { {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
109 {{0,0,0,0}, {0,0,0,0}, {0,0,0,1}, {0,0,-1,0}},
110 {{0,0,0,0}, {0,0,0,-1}, {0,0,0,0}, {0,1,0,0} },
111 {{0,0,0,0}, {0,0,1,0}, {0,-1,0,0}, {0,0,0,0} } },
112 { {{0,0,0,0}, {0,0,0,0}, {0,0,0,-1}, {0,0,1,0} },
113 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
114 {{0,0,0,1}, {0,0,0,0}, {0,0,0,0}, {-1,0,0,0}},
115 {{0,0,-1,0}, {0,0,0,0}, {1,0,0,0}, {0,0,0,0} } },
116 { {{0,0,0,0}, {0,0,0,1}, {0,0,0,0}, {0,-1,0,0}},
117 {{0,0,0,-1}, {0,0,0,0}, {0,0,0,0}, {1,0,0,0} },
118 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
119 {{0,1,0,0}, {-1,0,0,0}, {0,0,0,0}, {0,0,0,0} } },
120 { {{0,0,0,0}, {0,0,-1,0}, {0,1,0,0}, {0,0,0,0} },
121 {{0,0,1,0}, {0,0,0,0}, {-1,0,0,0}, {0,0,0,0} },
122 {{0,-1,0,0}, {1,0,0,0}, {0,0,0,0}, {0,0,0,0} },
123 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} } } };
124 for (int i=0; i<4; i++) {
125 for (int j=0; j<4; j++) {
126 G[i][j] = GG[i][j];
127 for (int k=0; k<4; k++) {
128 for (int l=0; l<4; l++) {
129 E[i][j][k][l] = EE[i][j][k][l];
130 }
131 }
132 }
133 }
134
135}
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 EvtDsToKKpipipi::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 137 of file EvtDsToKKpipipi.cc.

137 {
138 setProbMax(585000.0);
139}
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition EvtDecayBase.cc:297
The documentation for this class was generated from the following files: