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

#include <EvtDTopipi0pi0.hh>

+ Inheritance diagram for EvtDTopipi0pi0:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtDTopipi0pi0()

EvtDTopipi0pi0::EvtDTopipi0pi0 ( )
inline

Definition at line 11 of file EvtDTopipi0pi0.hh.

11{}

Referenced by clone().

◆ ~EvtDTopipi0pi0()

EvtDTopipi0pi0::~EvtDTopipi0pi0 ( )
virtual

Definition at line 39 of file EvtDTopipi0pi0.cc.

39{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDTopipi0pi0::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 45 of file EvtDTopipi0pi0.cc.

45 {
46 return new EvtDTopipi0pi0;
47}

◆ decay()

void EvtDTopipi0pi0::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 127 of file EvtDTopipi0pi0.cc.

127 {
128/*
129 double maxprob = 0.0;
130 for(int ir=0;ir<=60000000;ir++){
131 p->initializePhaseSpace(getNDaug(),getDaugs());
132 EvtVector4R D1 = p->getDaug(0)->getP4();
133 EvtVector4R D2 = p->getDaug(1)->getP4();
134 EvtVector4R D3 = p->getDaug(2)->getP4();
135
136 double P1[4], P2[4], P3[4];
137 P1[0] = D1.get(0); P1[1] = D1.get(1); P1[2] = D1.get(2); P1[3] = D1.get(3);
138 P2[0] = D2.get(0); P2[1] = D2.get(1); P2[2] = D2.get(2); P2[3] = D2.get(3);
139 P3[0] = D3.get(0); P3[1] = D3.get(1); P3[2] = D3.get(2); P3[3] = D3.get(3);
140
141 double value;
142 int g0[6]={2,3,1,1,2,2};
143 int spin[6]={0,0,0,2,1,1};
144 int nstates=6;
145 double r0[6]={3,3,3,3,3,3};
146 double r1[6]={5,5,5,5,5,5};
147 calEva(P1, P2, P3, mass, width, rho, phi, g0, spin, modetype, nstates, value,r0,r1);
148
149 if (value<0) continue;
150 if(value>maxprob) {
151 maxprob=value;
152 cout << "ir= " << ir << endl;
153 // cout << "double P1[4] = {" << P1[0] <<","<< P1[1] <<","<< P1[2] <<","<< P1[3] <<"};"<< endl;
154 // cout << "double P2[4] = {" << P2[0] <<","<< P2[1] <<","<< P2[2] <<","<< P2[3] <<"};"<< endl;
155 // cout << "double P3[4] = {" << P3[0] <<","<< P3[1] <<","<< P3[2] <<","<< P3[3] <<"};"<< endl;
156 cout << "MAX====> " << maxprob << endl;
157 }
158 }
159 printf("MAXprob = %.10f\n",maxprob);
160*/
161
162 p->initializePhaseSpace(getNDaug(),getDaugs());
163 EvtVector4R D1 = p->getDaug(0)->getP4();
164 EvtVector4R D2 = p->getDaug(1)->getP4();
165 EvtVector4R D3 = p->getDaug(2)->getP4();
166
167 double P1[4], P2[4], P3[4];
168 P1[0] = D1.get(0); P1[1] = D1.get(1); P1[2] = D1.get(2); P1[3] = D1.get(3);
169 P2[0] = D2.get(0); P2[1] = D2.get(1); P2[2] = D2.get(2); P2[3] = D2.get(3);
170 P3[0] = D3.get(0); P3[1] = D3.get(1); P3[2] = D3.get(2); P3[3] = D3.get(3);
171
172 double value;
173 int g0[4]={1,2,2,6};
174 int spin[4]={2,1,1,0};
175 int nstates=4;
176 double r0[4]={3,3,3,3};
177 double r1[4]={5,5,5,5};
178
179 calEva(P1, P2, P3, mass, width, rho, phi, g0, spin, modetype, nstates, value,r0,r1);
180
181 setProb(value);
182 return ;
183
184}
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 EvtDTopipi0pi0::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 41 of file EvtDTopipi0pi0.cc.

41 {
42 model_name="DTopipi0pi0";
43}

◆ init()

void EvtDTopipi0pi0::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 49 of file EvtDTopipi0pi0.cc.

49 {
50 // check that there are 0 arguments
51 checkNArg(0);
52 checkNDaug(3);
53 checkSpinParent(EvtSpinType::SCALAR);
54 checkSpinDaughter(0,EvtSpinType::SCALAR);
55 checkSpinDaughter(1,EvtSpinType::SCALAR);
56 checkSpinDaughter(2,EvtSpinType::SCALAR);
57
58 phi[0] =-3.2360;
59 phi[1] = 0.0;
60 phi[2] = 2.6124;
61 phi[3] = 1.2792;
62
63 rho[0] = 1.2606;
64 rho[1] = 1.0;
65 rho[2] = 0.98286;
66 rho[3] = 0.36030;
67
68
69 modetype[0]= 23;
70 modetype[1]= 12;
71 modetype[2]= 12;
72 modetype[3]= 23;
73
74/*
75 std::cout << "Initializing EvtDTopipi0pi0" << std::endl;
76 for (int i=0; i<6; i++) {
77 cout << i << " rho= " << rho[i] << " phi= " << phi[i] << endl;
78 }
79*/
80
81 width[0] = 0.1867;
82 width[1] = 0.1502;
83 width[2] = 0.4;
84 width[3] = 0.020;
85
86 mass[0] = 1.2755;
87 mass[1] = 0.7665;
88 mass[2] = 1.465;
89 mass[3] = 0.900;
90
91 mD = 1.86486;
92 mDs = 1.9683;
93 rRes = 9.0;
94 rD = 5.0;
95 metap = 0.95778;
96 mkstr = 0.89594;
97 mk0 = 0.497614;
98 mass_Kaon = 0.49368;
99 mass_Pion = 0.13957;
100 mass_Pi0 = 0.1349766;
101 math_pi = 3.1415926;
102 ma0 = 0.99;
103 Ga0 = 0.0756;
104 meta = 0.547862;
105 mass_Kaon = 0.49368;
106 mass_Eta = 0.547862;
107 mass_Etap = 0.95778;
108
109 GS1 = 0.636619783;
110 GS2 = 0.01860182466;
111 GS3 = 0.1591549458; // 1/(2*math_2pi)
112 GS4 = 0.00620060822; // mass_Pion2/math_pi
113
114 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
115 for (int i=0; i<4; i++) {
116 for (int j=0; j<4; j++) {
117 G[i][j] = GG[i][j];
118 }
119 }
120}
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 EvtDTopipi0pi0::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 122 of file EvtDTopipi0pi0.cc.

122 {
123 setProbMax(815.0);
124}
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition EvtDecayBase.cc:297
The documentation for this class was generated from the following files: