BOSS 7.1.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDTopiEtaEta Class Reference

#include <EvtDTopiEtaEta.hh>

+ Inheritance diagram for EvtDTopiEtaEta:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtDTopiEtaEta()

EvtDTopiEtaEta::EvtDTopiEtaEta ( )
inline

Definition at line 11 of file EvtDTopiEtaEta.hh.

11{}

Referenced by clone().

◆ ~EvtDTopiEtaEta()

EvtDTopiEtaEta::~EvtDTopiEtaEta ( )
virtual

Definition at line 32 of file EvtDTopiEtaEta.cc.

32{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDTopiEtaEta::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 38 of file EvtDTopiEtaEta.cc.

38 {
39 return new EvtDTopiEtaEta;
40}

◆ decay()

void EvtDTopiEtaEta::decay ( EvtParticle p)
virtual

Implements EvtDecayBase.

Definition at line 98 of file EvtDTopiEtaEta.cc.

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

Implements EvtDecayBase.

Definition at line 34 of file EvtDTopiEtaEta.cc.

34 {
35 model_name="DTopiEtaEta";
36}

◆ init()

void EvtDTopiEtaEta::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 42 of file EvtDTopiEtaEta.cc.

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

Reimplemented from EvtDecayBase.

Definition at line 94 of file EvtDTopiEtaEta.cc.

94 {
95 setProbMax(62.5);
96}
void setProbMax(double prbmx)

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