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

#include <EvtD0ToKpipi0.hh>

+ Inheritance diagram for EvtD0ToKpipi0:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtD0ToKpipi0()

EvtD0ToKpipi0::EvtD0ToKpipi0 ( )
inline

Definition at line 11 of file EvtD0ToKpipi0.hh.

11{}

Referenced by clone().

◆ ~EvtD0ToKpipi0()

EvtD0ToKpipi0::~EvtD0ToKpipi0 ( )
virtual

Definition at line 33 of file EvtD0ToKpipi0.cc.

33{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtD0ToKpipi0::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 39 of file EvtD0ToKpipi0.cc.

39 {
40 return new EvtD0ToKpipi0;
41}

◆ decay()

void EvtD0ToKpipi0::decay ( EvtParticle p)
virtual

Implements EvtDecayBase.

Definition at line 129 of file EvtD0ToKpipi0.cc.

129 {
130 //-----------for max value------------------
131 /* double maxprob = 0.0;
132 for(int ir=0;ir<=60000000;ir++){
133 p->initializePhaseSpace(getNDaug(),getDaugs());
134 EvtVector4R ks0 = p->getDaug(0)->getP4();
135 EvtVector4R pic = p->getDaug(1)->getP4();
136 EvtVector4R pi0 = p->getDaug(2)->getP4();
137
138 double Ks0[4],Pic[4],Pi0[4];
139 Ks0[0] = ks0.get(0); Pic[0] = pic.get(0); Pi0[0] = pi0.get(0);
140 Ks0[1] = ks0.get(1); Pic[1] = pic.get(1); Pi0[1] = pi0.get(1);
141 Ks0[2] = ks0.get(2); Pic[2] = pic.get(2); Pi0[2] = pi0.get(2);
142 Ks0[3] = ks0.get(3); Pic[3] = pic.get(3); Pi0[3] = pi0.get(3);
143 double value;
144 calPDF(Ks0, Pic, Pi0, value);
145 if(value>maxprob) {
146 maxprob=value;
147 std::cout << "Max PDF = " << ir << " prob= " << value << std::endl;
148 }
149 }
150 std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
151 return;*/
152 //-----------------------------------------------
154 EvtVector4R ks0 = p->getDaug(0)->getP4();
155 EvtVector4R pic = p->getDaug(1)->getP4();
156 EvtVector4R pi0 = p->getDaug(2)->getP4();
157
158 double Ks0[4],Pic[4],Pi0[4];
159 Ks0[0] = ks0.get(0); Pic[0] = pic.get(0); Pi0[0] = pi0.get(0);
160 Ks0[1] = ks0.get(1); Pic[1] = pic.get(1); Pi0[1] = pi0.get(1);
161 Ks0[2] = ks0.get(2); Pic[2] = pic.get(2); Pi0[2] = pi0.get(2);
162 Ks0[3] = ks0.get(3); Pic[3] = pic.get(3); Pi0[3] = pi0.get(3);
163
164 double value;
165 calPDF(Ks0, Pic, Pi0, value);
166 setProb(value);
167 return;
168}
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 EvtD0ToKpipi0::getName ( std::string &  name)
virtual

Implements EvtDecayBase.

Definition at line 35 of file EvtD0ToKpipi0.cc.

35 {
36 model_name="D0ToKpipi0";
37}

◆ init()

void EvtD0ToKpipi0::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 43 of file EvtD0ToKpipi0.cc.

43 {
44 checkNArg(0);
45 checkNDaug(3);
47
48 mrho_770 = 0.76519, mrho_1450 = 1.465, mKn892 = 0.89555, mKp892 = 0.89167, mK0_1430 = 1.466, mK2_1430 = 1.432, mKn_1680 = 1.717;
49 Grho_770 = 0.150, Grho_1450 = 0.310, GKn892 = 0.0473, GKp892 = 0.0514, GK0_1430 = 0.174, GK2_1430 = 0.109, GKn_1680 = 0.322;
50
51 rho[0] = 1;//rho770-pipi0
52 rho[1] = 1.32950;//rho1450-pipi0
53 rho[2] = 0.34447;//K*892n-Kpi
54 rho[3] = 0.37795;//K*982p-Kpi0
55 rho[4] = 0.15531;//K0*1430n-Kpi
56 rho[5] = 0.27334;//K0*1430p-Kpi0
57 rho[6] = 0.34413;//K2*1430n-Kpi
58 rho[7] = 0.34723;//K2*1430p-Kpi0
59 rho[8] = 0.38826;//K*1680n-Kpi
60 rho[9] = 0.86076;//KpiSW
61
62 phi[0] = 0;//rho770
63 phi[1] = -3.5053; //rho1450-pipi0
64 phi[2] = 0.12764; //K*892n-Kpi
65 phi[3] = 9.1910; //K*982p-Kpi0
66 phi[4] = -5.4878; //K0*1430n-Kpi
67 phi[5] = 4.4317; //K0*1430p-Kpi0
68 phi[6] = -0.69265;//K2*1430n-Kpi
69 phi[7] = -3.5656; //K2*1430p-Kpi0
70 phi[8] = -4.3046; //K*1680n-Kpi
71 phi[9] = 8.8127; //KpiSW
72
73 spin[0] = 1;
74 spin[1] = 1;
75 spin[2] = 1;
76 spin[3] = 1;
77 spin[4] = 0;
78 spin[5] = 0;
79 spin[6] = 2;
80 spin[7] = 2;
81 spin[8] = 1;
82 spin[9] = 0;
83
84 modetype[0] = 23;
85 modetype[1] = 23;
86 modetype[2] = 12;
87 modetype[3] = 13;
88 modetype[4] = 12;
89 modetype[5] = 13;
90 modetype[6] = 12;
91 modetype[7] = 13;
92 modetype[8] = 12;
93 modetype[9] = 132;
94
95 //std::cout << "Initializing EvtD0ToKpipi0" << std::endl;
96 //for (int i=0; i<10; i++) {
97 // cout << i << " rho,phi = " << rho[i] << ", "<< phi[i] << endl;
98 //}
99
100 mass_Ks=0.4977;
101 mass_Eta=0.547862;
102 mD = 1.86965;
103 rD = 5;
104 metap = 0.95778;
105 mkstr = 0.89594;
106 mk0 = 0.497611;
107 mass_Kaon = 0.493677;
108 mass_Pion = 0.13957;
109 mass_Pi0 = 0.1349768;
110 math_pi = 3.1415926;
111
112 pi = 3.1415926;
113 mpi = 0.13957;
114 g1 = 0.5468;
115 g2 = 0.23;
116
117 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
118 for (int i=0; i<4; i++) {
119 for (int j=0; j<4; j++) {
120 G[i][j] = GG[i][j];
121 }
122 }
123}
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 EvtD0ToKpipi0::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 125 of file EvtD0ToKpipi0.cc.

125 {
126 setProbMax(144.5);//setProbMax;
127}
void setProbMax(double prbmx)

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