BOSS 7.1.2
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDsToEtapi2pi0 Class Reference

#include <EvtDsToEtapi2pi0.hh>

+ Inheritance diagram for EvtDsToEtapi2pi0:

Public Member Functions

 EvtDsToEtapi2pi0 ()
 
virtual ~EvtDsToEtapi2pi0 ()
 
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 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 8 of file EvtDsToEtapi2pi0.hh.

Constructor & Destructor Documentation

◆ EvtDsToEtapi2pi0()

EvtDsToEtapi2pi0::EvtDsToEtapi2pi0 ( )
inline

Definition at line 12 of file EvtDsToEtapi2pi0.hh.

12{}

Referenced by clone().

◆ ~EvtDsToEtapi2pi0()

EvtDsToEtapi2pi0::~EvtDsToEtapi2pi0 ( )
virtual

Definition at line 33 of file EvtDsToEtapi2pi0.cc.

33{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDsToEtapi2pi0::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 39 of file EvtDsToEtapi2pi0.cc.

39 {
40 return new EvtDsToEtapi2pi0;
41}

◆ decay()

void EvtDsToEtapi2pi0::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 110 of file EvtDsToEtapi2pi0.cc.

110 {
111/*
112 double maxprob = 0.0;
113 for(int ir=0;ir<=60000000;ir++){
114 p->initializePhaseSpace(getNDaug(),getDaugs());
115 EvtVector4R pic = p->getDaug(0)->getP4();
116 EvtVector4R pi01 = p->getDaug(1)->getP4();
117 EvtVector4R pi02 = p->getDaug(2)->getP4();
118 EvtVector4R eta = p->getDaug(3)->getP4();
119
120 double Pic[4],Pi01[4],Pi02[4],Eta[4];
121 Pic[0]=pic.get(0); Pic[1]=pic.get(1); Pic[2]=pic.get(2); Pic[3]=pic.get(3);
122 Pi01[0]=pi01.get(0); Pi01[1]=pi01.get(1); Pi01[2]=pi01.get(2); Pi01[3]=pi01.get(3);
123 Pi02[0]=pi02.get(0); Pi02[1]=pi02.get(1); Pi02[2]=pi02.get(2); Pi02[3]=pi02.get(3);
124 Eta[0] = eta.get(0); Eta[1] = eta.get(1); Eta[2] = eta.get(2); Eta[3] = eta.get(3);
125
126 double prob;
127 int g0[3]={1,1,1};
128 int g1[3]={1,1,1};
129 int g2[3]={0,1,0};
130 int nstates=3;
131
132 calEvaMy(Pic, Pi01, Pi02, Eta, mass1, mass2, width1, width2, rho, phi, g0, g1, g2, modetype, nstates, prob);
133
134 if(prob>maxprob) {
135 maxprob=prob;
136 printf("nrun = %d,maxprob = %.10f,prob = %.10f\n\n",ir,maxprob,prob);
137 }
138 }
139*/
140
142 EvtVector4R pic = p->getDaug(0)->getP4();
143 EvtVector4R pi01 = p->getDaug(1)->getP4();
144 EvtVector4R pi02 = p->getDaug(2)->getP4();
145 EvtVector4R eta = p->getDaug(3)->getP4();
146
147 double Pic[4],Pi01[4],Pi02[4],Eta[4];
148 Pic[0]=pic.get(0); Pic[1]=pic.get(1); Pic[2]=pic.get(2); Pic[3]=pic.get(3);
149 Pi01[0]=pi01.get(0); Pi01[1]=pi01.get(1); Pi01[2]=pi01.get(2); Pi01[3]=pi01.get(3);
150 Pi02[0]=pi02.get(0); Pi02[1]=pi02.get(1); Pi02[2]=pi02.get(2); Pi02[3]=pi02.get(3);
151 Eta[0]=eta.get(0); Eta[1]=eta.get(1); Eta[2]=eta.get(2); Eta[3]=eta.get(3);
152
153 double prob;
154 int g0[3]={1,1,1};
155 int g1[3]={1,1,1};
156 int g2[3]={0,1,0};
157 int nstates=3;
158
159 calEvaMy(Pic, Pi01, Pi02, Eta, mass1, mass2, width1, width2, rho, phi, g0, g1, g2, modetype, nstates, prob);
160 setProb(prob);
161
162 return ;
163}
TF1 * g1
EvtId * getDaugs()
void setProb(double prob)
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
double get(int i) const

◆ getName()

void EvtDsToEtapi2pi0::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 35 of file EvtDsToEtapi2pi0.cc.

35 {
36 model_name="DsToEtapi2pi0";
37}

◆ init()

void EvtDsToEtapi2pi0::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 43 of file EvtDsToEtapi2pi0.cc.

43 {
44
45 // check that there are 0 arguments
46 checkNArg(0);
47 checkNDaug(4);
48
54
55 int mode = 0;
56 phi[0] = 0.0000e+00;
57 phi[1] = 2.2539e+00;
58 phi[2] = -4.1144e+00;
59 rho[0] = 1.0000e+00;
60 rho[1] = 3.0726e-01;
61 rho[2] = 8.8537e-01;
62 mass1[0] = 7.7526e-01;
63 mass1[1] = 9.8000e-01;
64 mass1[2] = 9.8000e-01;
65 mass2[0] = 1.2300e+00;
66 mass2[1] = 7.7526e-01;
67 mass2[2] = 1.8000e+00;
68 width1[0] = 1.4910e-01;
69 width1[1] = 7.0000e-02;
70 width1[2] = 7.0000e-02;
71 width2[0] = 3.3049e-01;
72 width2[1] = 1.4910e-01;
73 width2[2] = 3.8900e-01;
74
75 modetype[0] = 1;
76 modetype[1] = 2;
77 modetype[2] = 3;
78
79 //cout << "DsToEtapi2pi0 :" << endl;
80 //for (int i=0; i<3; i++) {
81 // cout << i << " rho= " << rho[i] << " phi= " << phi[i] << " m1= " << mass1[i] << " m2= " << mass2[i] << " w1= " << width1[i] << " w2= " << width2[i] << endl;
82 //}
83
84 rRes = 3.0;
85 rRes2 = 9.0;
86 rD = 5.0;
87 rD2 = 25.0;
88 mass_Ks = 0.497611;
89 mass_Kaon = 0.493677;
90 mass_Pion = 0.13957;
91 mass_Pi0 = 0.1349768;
92 meta = 0.547862;
93 GS1 = 0.636619783;
94 GS2 = 0.01860182466;
95 GS3 = 0.1591549458; // 1/(2*math_2pi)
96 GS4 = 0.00620060822;
97
98 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
99 for (int a=0; a<4; a++) {
100 for (int j=0; j<4; j++) {
101 G[a][j] = GG[a][j];
102 }
103 }
104}
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
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 EvtDsToEtapi2pi0::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 106 of file EvtDsToEtapi2pi0.cc.

106 {
107 setProbMax(4400.0);
108}
void setProbMax(double prbmx)

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