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

#include <EvtDTopipi0Eta.hh>

+ Inheritance diagram for EvtDTopipi0Eta:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtDTopipi0Eta()

EvtDTopipi0Eta::EvtDTopipi0Eta ( )
inline

Definition at line 11 of file EvtDTopipi0Eta.hh.

11{}

Referenced by clone().

◆ ~EvtDTopipi0Eta()

EvtDTopipi0Eta::~EvtDTopipi0Eta ( )
virtual

Definition at line 32 of file EvtDTopipi0Eta.cc.

32{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDTopipi0Eta::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 38 of file EvtDTopipi0Eta.cc.

38 {
39 return new EvtDTopipi0Eta;
40}

◆ decay()

void EvtDTopipi0Eta::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 97 of file EvtDTopipi0Eta.cc.

97 {
98
99 // This piece of code could in principle be used to calculate maximum
100 // probablity on fly. But as it uses high number of points and model
101 // deals with single final state, we keep hardcoded number for now rather
102 // than adapting code to work here.
103
104// double maxprob = 0.0;
105// for(int ir=0;ir<=60000000;ir++){
106// p->initializePhaseSpace(getNDaug(),getDaugs());
107// EvtVector4R D1 = p->getDaug(0)->getP4();
108// EvtVector4R D2 = p->getDaug(1)->getP4();
109// EvtVector4R D3 = p->getDaug(2)->getP4();
110//
111// double P1[4], P2[4], P3[4];
112// P1[0] = D1.get(0); P1[1] = D1.get(1); P1[2] = D1.get(2); P1[3] = D1.get(3);
113// P2[0] = D2.get(0); P2[1] = D2.get(1); P2[2] = D2.get(2); P2[3] = D2.get(3);
114// P3[0] = D3.get(0); P3[1] = D3.get(1); P3[2] = D3.get(2); P3[3] = D3.get(3);
115//
116// double value;
117// value = calDalEva(P1, P2, P3);
118// if(value>maxprob) {
119// maxprob=value;
120// cout << "ir = " << ir << " maxProb= " << value << endl;
121// }
122// }
123// cout << "maxProb = " << maxprob << endl;
124
126 EvtVector4R D1 = p->getDaug(0)->getP4();
127 EvtVector4R D2 = p->getDaug(1)->getP4();
128 EvtVector4R D3 = p->getDaug(2)->getP4();
129
130 double P1[4], P2[4], P3[4];
131 P1[0] = D1.get(0); P1[1] = D1.get(1); P1[2] = D1.get(2); P1[3] = D1.get(3);
132 P2[0] = D2.get(0); P2[1] = D2.get(1); P2[2] = D2.get(2); P2[3] = D2.get(3);
133 P3[0] = D3.get(0); P3[1] = D3.get(1); P3[2] = D3.get(2); P3[3] = D3.get(3);
134
135 double value;
136 value = calDalEva(P1, P2, P3);
137 setProb(value);
138 return;
139
140}
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 EvtDTopipi0Eta::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 34 of file EvtDTopipi0Eta.cc.

34 {
35 model_name="DTopipi0Eta";
36}

◆ init()

void EvtDTopipi0Eta::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 42 of file EvtDTopipi0Eta.cc.

42 {
43 // check that there are 0 arguments
44 checkNArg(0);
45 checkNDaug(3);
50
51 phi[0] = -3.3276; rho[0] = 0.31478; //rho eta
52 phi[1] = 0.0; rho[1] = 1.0; //a0+ pi0
53 phi[2] = 0.0; rho[2] = -1.0; //a00 pi+
54
55 //cout << "Initializing DTopipi0Eta" << endl;
56 //for (int i=0; i<3; i++) {
57 // cout << i << " rho= " << rho[i] << " phi= " << phi[i] << endl;
58 //}
59 mrho = 0.77511;
60 ma0 = 0.99;
61 Grho = 0.1491;
62 Ga0 = 0.0756;
63
64 const double mk0 = 0.497614;
65 const double mass_Kaon = 0.49368;
66 const double mass_Pion = 0.13957;
67 const double mass_Pi0 = 0.1349766;
68 const double meta = 0.547862;
69 mpi = 0.13957;
70 mD = 1.86966;
71 sD = mD*mD;
72 spi = mpi*mpi;
73 snk = mk0*mk0;
74 sck = mass_Kaon*mass_Kaon;
75 scpi = mass_Pion*mass_Pion;
76 snpi = mass_Pi0*mass_Pi0;
77 seta = meta*meta;
78
79 pi = 3.1415926;
80
81 ci = EvtComplex(0.0,1.0);
82 one = EvtComplex(1.0,0.0);
83
84 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
85 for (int i=0; i<4; i++) {
86 for (int j=0; j<4; j++) {
87 G[i][j] = GG[i][j];
88 }
89 }
90
91}
const double mass_Pion
double meta
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)
const double mk0
Definition inclks.cxx:46

◆ initProbMax()

void EvtDTopipi0Eta::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 93 of file EvtDTopipi0Eta.cc.

93 {
94 setProbMax(20.0);
95}
void setProbMax(double prbmx)

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