BOSS 8.0.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 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 10 of file EvtD0ToKpipi0.hh.

10{}

Referenced by clone().

◆ ~EvtD0ToKpipi0()

EvtD0ToKpipi0::~EvtD0ToKpipi0 ( )
virtual

Definition at line 34 of file EvtD0ToKpipi0.cc.

34{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtD0ToKpipi0::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 38 of file EvtD0ToKpipi0.cc.

38{ return new EvtD0ToKpipi0; }

◆ decay()

void EvtD0ToKpipi0::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 133 of file EvtD0ToKpipi0.cc.

133 {
134 //-----------for max value------------------
135 /*
136 double maxprob = 0.0;
137 for(int ir=0;ir<=60000000;ir++){
138 p->initializePhaseSpace(getNDaug(),getDaugs());
139 EvtVector4R ks0 = p->getDaug(0)->getP4();
140 EvtVector4R pic = p->getDaug(1)->getP4();
141 EvtVector4R pi0 = p->getDaug(2)->getP4();
142
143 double Ks0[4],Pic[4],Pi0[4];
144 Ks0[0] = ks0.get(0); Pic[0] = pic.get(0); Pi0[0] = pi0.get(0);
145 Ks0[1] = ks0.get(1); Pic[1] = pic.get(1); Pi0[1] = pi0.get(1);
146 Ks0[2] = ks0.get(2); Pic[2] = pic.get(2); Pi0[2] = pi0.get(2);
147 Ks0[3] = ks0.get(3); Pic[3] = pic.get(3); Pi0[3] = pi0.get(3);
148 double value;
149 calPDF(Ks0, Pic, Pi0, value);
150 if(value>maxprob) {
151 maxprob=value;
152 std::cout << "Max PDF = " << ir << " prob= " << value << std::endl;
153 }
154 }
155 std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
156 return;
157 */
158 //-----------------------------------------------
160 EvtVector4R ks0 = p->getDaug( 0 )->getP4();
161 EvtVector4R pic = p->getDaug( 1 )->getP4();
162 EvtVector4R pi0 = p->getDaug( 2 )->getP4();
163
164 double Ks0[4], Pic[4], Pi0[4];
165 Ks0[0] = ks0.get( 0 );
166 Pic[0] = pic.get( 0 );
167 Pi0[0] = pi0.get( 0 );
168 Ks0[1] = ks0.get( 1 );
169 Pic[1] = pic.get( 1 );
170 Pi0[1] = pi0.get( 1 );
171 Ks0[2] = ks0.get( 2 );
172 Pic[2] = pic.get( 2 );
173 Pi0[2] = pi0.get( 2 );
174 Ks0[3] = ks0.get( 3 );
175 Pic[3] = pic.get( 3 );
176 Pi0[3] = pi0.get( 3 );
177
178 double value;
179 calPDF( Ks0, Pic, Pi0, value );
180
181 setProb( value );
182 return;
183}
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 36 of file EvtD0ToKpipi0.cc.

36{ model_name = "D0ToKpipi0"; }

◆ init()

void EvtD0ToKpipi0::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 40 of file EvtD0ToKpipi0.cc.

40 {
41 checkNArg( 0 );
42 checkNDaug( 3 );
44
45 mrho_770 = 0.76519, mrho_1450 = 1.465, mKn892 = 0.89555, mKp892 = 0.89167, mK0_1430 = 1.466,
46 mK2_1430 = 1.432, mKn_1680 = 1.717;
47 Grho_770 = 0.150, Grho_1450 = 0.310, GKn892 = 0.05, GKp892 = 0.0514, GK0_1430 = 0.174,
48 GK2_1430 = 0.109, GKn_1680 = 0.322;
49
50 rho[0] = 1; // rho770-pipi0
51 rho[1] = 1.5481; // rho1450-pipi0
52 rho[2] = 0.41186; // K*892p-Kpi
53 rho[3] = 0.38950; // K*982n-Kpi0
54 rho[4] = 0.37221; // K0*1430n-Kpi
55 rho[5] = 0.31480; // K0*1430p-Kpi0
56 rho[6] = 0.35958; // K2*1430n-Kpi
57 rho[7] = 0.36683; // K2*1430p-Kpi0
58 rho[8] = 1.0967; // K*1680n-Kpi
59 rho[9] = 0.26035; // K*1680p-Kpi
60 rho[10] = 1.0210; // res-S
61 rho[11] = 0.91605; // res-Kpi0p
62
63 phi[0] = 0; // rho770
64 phi[1] = 8.8326; // rho1450-pipi0
65 phi[2] = 2.9453; // K*892p-Kpi
66 phi[3] = 6.4385; // K*982n-Kpi0
67 phi[4] = 0.42526; // K0*1430n-Kpi
68 phi[5] = -1.8798; // K0*1430p-Kpi0
69 phi[6] = 5.7080; // K2*1430n-Kpi
70 phi[7] = 2.6428; // K2*1430p-Kpi0
71 phi[8] = 1.0485; // K*1680n-Kpi
72 phi[9] = 1.2300; // K*1680p-Kpi
73 phi[10] = 10.459; // res-S
74 phi[11] = 10.449; // res-Kpi0p
75
76 spin[0] = 1;
77 spin[1] = 1;
78 spin[2] = 1;
79 spin[3] = 1;
80 spin[4] = 0;
81 spin[5] = 0;
82 spin[6] = 2;
83 spin[7] = 2;
84 spin[8] = 1;
85 spin[9] = 1;
86 spin[10] = 0;
87 spin[11] = 1;
88
89 modetype[0] = 23;
90 modetype[1] = 23;
91 modetype[2] = 13;
92 modetype[3] = 12;
93 modetype[4] = 12;
94 modetype[5] = 13;
95 modetype[6] = 12;
96 modetype[7] = 13;
97 modetype[8] = 12;
98 modetype[9] = 13;
99 modetype[10] = 12;
100 modetype[11] = 13;
101
102 /*
103 std::cout << "EvtD0ToKpipi0 (May 01, 2023) ==> Initialization" << std::endl;
104 for (int i=0; i<13; i++) {
105 cout << i << " rho,phi = " << rho[i] << ", "<< phi[i] << endl;
106 }
107 */
108 mass_Ks = 0.4977;
109 mass_Eta = 0.547862;
110 mD = 1.864;
111 rD = 5;
112 metap = 0.95778;
113 mkstr = 0.89594;
114 mk0 = 0.497611;
115 mass_Kaon = 0.493677;
116 mass_Pion = 0.13957;
117 mass_Pi0 = 0.1349768;
118 math_pi = 3.1415926;
119 pi = 3.1415926;
120 mpi = 0.13957;
121
122 int GG[4][4] = { { 1, 0, 0, 0 }, { 0, -1, 0, 0 }, { 0, 0, -1, 0 }, { 0, 0, 0, -1 } };
123 for ( int i = 0; i < 4; i++ )
124 {
125 for ( int j = 0; j < 4; j++ ) { G[i][j] = GG[i][j]; }
126 }
127}
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 129 of file EvtD0ToKpipi0.cc.

129 {
130 setProbMax( 188.35 ); // setProbMax;
131}
void setProbMax(double prbmx)

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