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

#include <EvtD0ToKpiEtap.hh>

+ Inheritance diagram for EvtD0ToKpiEtap:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtD0ToKpiEtap()

EvtD0ToKpiEtap::EvtD0ToKpiEtap ( )
inline

Definition at line 12 of file EvtD0ToKpiEtap.hh.

12{}

Referenced by clone().

◆ ~EvtD0ToKpiEtap()

EvtD0ToKpiEtap::~EvtD0ToKpiEtap ( )
virtual

Definition at line 30 of file EvtD0ToKpiEtap.cc.

30{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtD0ToKpiEtap::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 36 of file EvtD0ToKpiEtap.cc.

36 {
37 return new EvtD0ToKpiEtap;
38}

◆ decay()

void EvtD0ToKpiEtap::decay ( EvtParticle p)
virtual

Implements EvtDecayBase.

Definition at line 93 of file EvtD0ToKpiEtap.cc.

93 {
94
95 // This piece of code could in principle be used to calculate maximum
96 // probablity on fly. But as it uses high number of points and model
97 // deals with single final state, we keep hardcoded number for now rather
98 // than adapting code to work here.
99
100 /* double maxprob = 0.0;
101 for(int ir=0;ir<=60000000;ir++){
102 p->initializePhaseSpace(getNDaug(),getDaugs());
103 EvtVector4R D1 = p->getDaug(0)->getP4();
104 EvtVector4R D2 = p->getDaug(1)->getP4();
105 EvtVector4R D3 = p->getDaug(2)->getP4();
106
107 double P1[4], P2[4], P3[4];
108 P1[0] = D1.get(0); P1[1] = D1.get(1); P1[2] = D1.get(2); P1[3] = D1.get(3);
109 P2[0] = D2.get(0); P2[1] = D2.get(1); P2[2] = D2.get(2); P2[3] = D2.get(3);
110 P3[0] = D3.get(0); P3[1] = D3.get(1); P3[2] = D3.get(2); P3[3] = D3.get(3);
111
112 double value;
113 value = calDalEva(P1, P2, P3);
114 if(value>maxprob) {
115 maxprob=value;
116 cout << "ir = " << ir << " maxProb= " << value << endl;
117 }
118 }
119 cout << "maxProb = " << maxprob << endl;
120 */
122 EvtVector4R D1 = p->getDaug(0)->getP4();
123 EvtVector4R D2 = p->getDaug(1)->getP4();
124 EvtVector4R D3 = p->getDaug(2)->getP4();
125
126 double P1[4], P2[4], P3[4];
127 P1[0] = D1.get(0); P1[1] = D1.get(1); P1[2] = D1.get(2); P1[3] = D1.get(3);
128 P2[0] = D2.get(0); P2[1] = D2.get(1); P2[2] = D2.get(2); P2[3] = D2.get(3);
129 P3[0] = D3.get(0); P3[1] = D3.get(1); P3[2] = D3.get(2); P3[3] = D3.get(3);
130
131
132/* double P1[4] = {0.5891186 , 0.1178592 , -0.289217 , 0.0761997};
133 double P2[4] = {0.2860281 , -0.095573 , 0.2069696 , -0.101791};
134 double P3[4] = {1.0120198 , -0.171013 , 0.2017959 , -0.192027};
135*/
136
137 double value;
138 int g0[3]={0,1,1};
139 int spin[3]={0,1,1};
140 int nstates=3;
141 calEva(P1, P2, P3, mass, width, rho, phi, g0, spin, modetype, nstates, value);
142
143 setProb(value);
144 return ;
145
146}
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 EvtD0ToKpiEtap::getName ( std::string &  name)
virtual

Implements EvtDecayBase.

Definition at line 32 of file EvtD0ToKpiEtap.cc.

32 {
33 model_name="D0ToKpiEtap";
34}

◆ init()

void EvtD0ToKpiEtap::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 39 of file EvtD0ToKpiEtap.cc.

39 {
40 checkNArg(0);
41 checkNDaug(3);
46
47 phi[0] = 0.0; rho[0] = 1.0;
48 phi[1] = 3.5567 ; rho[1] = 2.0540e-01; //892
49 phi[2] = 0.61028; rho[2] = 1.5846e+00;//1680
50
51 //std::cout << "Initializing EvtD0ToKpiEtap" << std::endl;
52 //for (int i=0; i<3; i++) {
53 // cout << i << " rho= " << rho[i] << " phi= " << phi[i] << endl;
54 //}
55
56 modetype[0]= 132;
57 modetype[1]= 12;
58 modetype[2]= 12;
59
60 mass[2] = 1.718;
61 mass[1] = 0.89555;
62 width[2]= 0.322;
63 width[1]= 0.0473;
64 mass[0] = 1.414;
65 width[0]= 0.232;
66 const double mk0 = 0.497614;
67 const double mass_Kaon = 0.49368;
68 const double mass_Pion = 0.13957;
69 const double mass_Pi0 = 0.1349766;
70 const double meta = 0.547862;
71 const double metap = 0.95778;
72 mpi = 0.13957;
73 mD = 1.86486;
74 mDs = 1.9683;
75 rRes= 9.0;
76 rD = 5.0;
77 math_pi = 3.1415926;
78 GS1 = 0.636619783;
79 GS2 = 0.01860182466;
80 GS3 = 0.1591549458; // 1/(2*math_2pi)
81 GS4 = 0.00620060822; // mass_Pion2/math_pi
82
83 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
84 for (int i=0; i<4; i++) { for (int j=0; j<4; j++) {
85 G[i][j] = GG[i][j];
86 }
87 }
88}
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 EvtD0ToKpiEtap::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 89 of file EvtD0ToKpiEtap.cc.

89 {
90 setProbMax(3.76);
91}
void setProbMax(double prbmx)

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