BOSS 7.1.2
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 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 135 of file EvtD0ToKpipi0.cc.

135 {
136 //-----------for max value------------------
137/*
138 double maxprob = 0.0;
139 for(int ir=0;ir<=60000000;ir++){
140 p->initializePhaseSpace(getNDaug(),getDaugs());
141 EvtVector4R ks0 = p->getDaug(0)->getP4();
142 EvtVector4R pic = p->getDaug(1)->getP4();
143 EvtVector4R pi0 = p->getDaug(2)->getP4();
144
145 double Ks0[4],Pic[4],Pi0[4];
146 Ks0[0] = ks0.get(0); Pic[0] = pic.get(0); Pi0[0] = pi0.get(0);
147 Ks0[1] = ks0.get(1); Pic[1] = pic.get(1); Pi0[1] = pi0.get(1);
148 Ks0[2] = ks0.get(2); Pic[2] = pic.get(2); Pi0[2] = pi0.get(2);
149 Ks0[3] = ks0.get(3); Pic[3] = pic.get(3); Pi0[3] = pi0.get(3);
150 double value;
151 calPDF(Ks0, Pic, Pi0, value);
152 if(value>maxprob) {
153 maxprob=value;
154 std::cout << "Max PDF = " << ir << " prob= " << value << std::endl;
155 }
156 }
157 std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
158 return;
159*/
160 //-----------------------------------------------
162 EvtVector4R ks0 = p->getDaug(0)->getP4();
163 EvtVector4R pic = p->getDaug(1)->getP4();
164 EvtVector4R pi0 = p->getDaug(2)->getP4();
165
166 double Ks0[4],Pic[4],Pi0[4];
167 Ks0[0] = ks0.get(0); Pic[0] = pic.get(0); Pi0[0] = pi0.get(0);
168 Ks0[1] = ks0.get(1); Pic[1] = pic.get(1); Pi0[1] = pi0.get(1);
169 Ks0[2] = ks0.get(2); Pic[2] = pic.get(2); Pi0[2] = pi0.get(2);
170 Ks0[3] = ks0.get(3); Pic[3] = pic.get(3); Pi0[3] = pi0.get(3);
171
172 double value;
173 calPDF(Ks0, Pic, Pi0, value);
174
175 setProb(value);
176 return;
177}
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 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.05, 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.5481; //rho1450-pipi0
53 rho[2] = 0.41186; //K*892p-Kpi
54 rho[3] = 0.38950; //K*982n-Kpi0
55 rho[4] = 0.37221; //K0*1430n-Kpi
56 rho[5] = 0.31480; //K0*1430p-Kpi0
57 rho[6] = 0.35958; //K2*1430n-Kpi
58 rho[7] = 0.36683; //K2*1430p-Kpi0
59 rho[8] = 1.0967; //K*1680n-Kpi
60 rho[9] = 0.26035; //K*1680p-Kpi
61 rho[10] = 1.0210; //res-S
62 rho[11] = 0.91605; //res-Kpi0p
63
64 phi[0] = 0; //rho770
65 phi[1] = 8.8326; //rho1450-pipi0
66 phi[2] = 2.9453; //K*892p-Kpi
67 phi[3] = 6.4385; //K*982n-Kpi0
68 phi[4] = 0.42526; //K0*1430n-Kpi
69 phi[5] = -1.8798; //K0*1430p-Kpi0
70 phi[6] = 5.7080; //K2*1430n-Kpi
71 phi[7] = 2.6428; //K2*1430p-Kpi0
72 phi[8] = 1.0485; //K*1680n-Kpi
73 phi[9] = 1.2300; //K*1680p-Kpi
74 phi[10] = 10.459; //res-S
75 phi[11] = 10.449; //res-Kpi0p
76
77 spin[0] = 1;
78 spin[1] = 1;
79 spin[2] = 1;
80 spin[3] = 1;
81 spin[4] = 0;
82 spin[5] = 0;
83 spin[6] = 2;
84 spin[7] = 2;
85 spin[8] = 1;
86 spin[9] = 1;
87 spin[10] = 0;
88 spin[11] = 1;
89
90 modetype[0] = 23;
91 modetype[1] = 23;
92 modetype[2] = 13;
93 modetype[3] = 12;
94 modetype[4] = 12;
95 modetype[5] = 13;
96 modetype[6] = 12;
97 modetype[7] = 13;
98 modetype[8] = 12;
99 modetype[9] = 13;
100 modetype[10] = 12;
101 modetype[11] = 13;
102
103/*
104 std::cout << "EvtD0ToKpipi0 (May 01, 2023) ==> Initialization" << std::endl;
105 for (int i=0; i<13; i++) {
106 cout << i << " rho,phi = " << rho[i] << ", "<< phi[i] << endl;
107 }
108*/
109 mass_Ks=0.4977;
110 mass_Eta=0.547862;
111 mD = 1.864;
112 rD = 5;
113 metap = 0.95778;
114 mkstr = 0.89594;
115 mk0 = 0.497611;
116 mass_Kaon = 0.493677;
117 mass_Pion = 0.13957;
118 mass_Pi0 = 0.1349768;
119 math_pi = 3.1415926;
120 pi = 3.1415926;
121 mpi = 0.13957;
122
123 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
124 for (int i=0; i<4; i++) {
125 for (int j=0; j<4; j++) {
126 G[i][j] = GG[i][j];
127 }
128 }
129}
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 131 of file EvtD0ToKpipi0.cc.

131 {
132 setProbMax(188.35);//setProbMax;
133}
void setProbMax(double prbmx)

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