BOSS 7.1.1
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | List of all members
EvtD0ToKSpi0pi0pi0 Class Reference

#include <EvtD0ToKSpi0pi0pi0.hh>

+ Inheritance diagram for EvtD0ToKSpi0pi0pi0:

Public Member Functions

 EvtD0ToKSpi0pi0pi0 ()
 
virtual ~EvtD0ToKSpi0pi0pi0 ()
 
void getName (std::string &name)
 
EvtDecayBase * clone ()
 
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 ()
 
EvtId * getDaugs ()
 
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 EvtD0ToKSpi0pi0pi0.hh.

Constructor & Destructor Documentation

◆ EvtD0ToKSpi0pi0pi0()

EvtD0ToKSpi0pi0pi0::EvtD0ToKSpi0pi0pi0 ( )
inline

Definition at line 11 of file EvtD0ToKSpi0pi0pi0.hh.

11{}

Referenced by clone().

◆ ~EvtD0ToKSpi0pi0pi0()

EvtD0ToKSpi0pi0pi0::~EvtD0ToKSpi0pi0pi0 ( )
virtual

Definition at line 32 of file EvtD0ToKSpi0pi0pi0.cc.

32{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtD0ToKSpi0pi0pi0::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 38 of file EvtD0ToKSpi0pi0pi0.cc.

38 {
39 return new EvtD0ToKSpi0pi0pi0;
40}

◆ decay()

void EvtD0ToKSpi0pi0pi0::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 114 of file EvtD0ToKSpi0pi0pi0.cc.

114 {
115 //-----------for max value------------------
116/* double maxprob = 0.0;
117 for(int ir=0;ir<=60000000;ir++){
118 p->initializePhaseSpace(getNDaug(),getDaugs());
119 EvtVector4R ks = p->getDaug(0)->getP4();
120 EvtVector4R kp = p->getDaug(1)->getP4();
121 EvtVector4R pim = p->getDaug(2)->getP4();
122 EvtVector4R pi0 = p->getDaug(3)->getP4();
123 mother_c=EvtPDL::getStdHep(p->getId());
124 int cc;
125 if(mother_c==421){
126 cc=1;
127 } else if(mother_c==-421){
128 cc=-1;
129 }
130 double value;
131 double KS[4],Pi01[4],Pi02[4],Pi03[4];
132 KS[0] = ks.get(0); Pi01[0] = kp.get(0); Pi02[0] = pim.get(0); Pi03[0] = pi0.get(0);
133 KS[1] = cc*ks.get(1); Pi01[1] = cc*kp.get(1); Pi02[1] = cc*pim.get(1); Pi03[1] = cc*pi0.get(1);
134 KS[2] = cc*ks.get(2); Pi01[2] = cc*kp.get(2); Pi02[2] = cc*pim.get(2); Pi03[2] = cc*pi0.get(2);
135 KS[3] = cc*ks.get(3); Pi01[3] = cc*kp.get(3); Pi02[3] = cc*pim.get(3); Pi03[3] = cc*pi0.get(3);
136 calPDF(KS, Pi01, Pi02, Pi03, value);
137 if(value>maxprob) {
138 maxprob=value;
139 std::cout << "Max PDF = " << ir << " prob= " << value << std::endl;
140 }
141 }
142 std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
143 return;*/
144 //-----------------------------------------------
145 p->initializePhaseSpace(getNDaug(),getDaugs());
146 EvtVector4R ks = p->getDaug(0)->getP4();
147 EvtVector4R pi01 = p->getDaug(1)->getP4();
148 EvtVector4R pi02 = p->getDaug(2)->getP4();
149 EvtVector4R pi03 = p->getDaug(3)->getP4();
150
151 mother_c=EvtPDL::getStdHep(p->getId());
152 int cc;
153 if(mother_c==421){
154 cc=1;
155 } else if(mother_c==-421){
156 cc=-1;
157 }
158
159 double KS[4],Pi01[4],Pi02[4],Pi03[4];
160 KS[0] = ks.get(0); Pi01[0] = pi01.get(0); Pi02[0] = pi02.get(0); Pi03[0] = pi03.get(0);
161 KS[1] = cc*ks.get(1); Pi01[1] = cc*pi01.get(1); Pi02[1] = cc*pi02.get(1); Pi03[1] = cc*pi03.get(1);
162 KS[2] = cc*ks.get(2); Pi01[2] = cc*pi01.get(2); Pi02[2] = cc*pi02.get(2); Pi03[2] = cc*pi03.get(2);
163 KS[3] = cc*ks.get(3); Pi01[3] = cc*pi01.get(3); Pi02[3] = cc*pi02.get(3); Pi03[3] = cc*pi03.get(3);
164 //KS[0] = 0.916636911124628; Pi01[0] = 0.450839902755463; Pi02[0] = 0.255815397731456; Pi03[0] = 0.269747834431058;
165 //KS[1] = -0.039197614059022; Pi01[1] = 0.002321391987197; Pi02[1] = 0.134785309537075; Pi03[1] = -0.227415473128566;
166 //KS[2] = 0.742055044382617; Pi01[2] = -0.397791163697056; Pi02[2] = -0.058650300320192; Pi03[2] = -0.039227963851607;
167 //KS[3] = -0.201058144869571; Pi01[3] = -0.163690811604007; Pi02[3] = 0.160049489145264; Pi03[3] = 0.035895968170634;
168
169 double value;
170 calPDF(KS, Pi01, Pi02, Pi03, value);
171 setProb(value);
172 return;
173
174}
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition EvtDecayBase.hh:65
EvtDecayProb::setProb
void setProb(double prob)
Definition EvtDecayProb.hh:34
EvtPDL::getStdHep
static int getStdHep(EvtId id)
Definition EvtPDL.hh:56
EvtParticle::getId
EvtId getId() const
Definition EvtParticle.cc:112
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition EvtParticle.cc:120
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition EvtParticle.cc:84
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition EvtParticle.cc:1070
EvtVector4R::get
double get(int i) const
Definition EvtVector4R.hh:179

◆ getName()

void EvtD0ToKSpi0pi0pi0::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 34 of file EvtD0ToKSpi0pi0pi0.cc.

34 {
35 model_name="D0ToKSpi0pi0pi0";
36}

◆ init()

void EvtD0ToKSpi0pi0pi0::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 42 of file EvtD0ToKSpi0pi0pi0.cc.

42 {
43 checkNArg(0);
44 checkNDaug(4);
45 checkSpinParent(EvtSpinType::SCALAR);
46 //cout << "Initializing EvtD0ToKSpi0pi0pi0" << endl;
47
48 mD = 1.86484; mPi = 0.13957; mKa = 0.493677;
49 math_pi = 3.1415926;
50 mKstr0 = 0.89555; mrhop = 0.77511;
51 GKstr0 = 0.0473; Grhop = 0.1491;
52 mKstrp = 0.89167; mrho0 = 0.77526;
53 GKstrp = 0.0514; Grho0 = 0.1474;
54 ma980 = 0.980;
55 Ga980 = 0.0972;
56 ma1_1260 = 1.230; mK1_1270 = 1.272;
57 Ga1_1260 = 0.420; GK1_1270 = 0.087;
58 meta1475 = 1.475; mK1_1400 = 1.463;
59 Geta1475 = 0.09; GK1_1400 = 0.3;
60
61 phi[1] = 2.3000e+00; //09 1270
62 phi[4] = 2.9764e+00; //126
63 phi[3] = 2.7663e+00; //21 kpiS
64 phi[0] = -2.8780e+00; //22 f0
65 rho[1] = 9.7790e-02;
66 rho[4] = 2.0000e+01;
67 rho[3] = 2.5101e+00;
68 rho[0] = 4.2006e-01;
69
70 phi[2] = 0.0; //1400
71 rho[2] = 1.0;
72
73
74 modetype[0]=1001;
75 modetype[1]=2006;
76 modetype[2]=2006;
77 modetype[3]=2007;
78 modetype[4]=126;
79
80 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
81 int EE[4][4][4][4] =
82 { { {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
83 {{0,0,0,0}, {0,0,0,0}, {0,0,0,1}, {0,0,-1,0}},
84 {{0,0,0,0}, {0,0,0,-1}, {0,0,0,0}, {0,1,0,0} },
85 {{0,0,0,0}, {0,0,1,0}, {0,-1,0,0}, {0,0,0,0} } },
86 { {{0,0,0,0}, {0,0,0,0}, {0,0,0,-1}, {0,0,1,0} },
87 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
88 {{0,0,0,1}, {0,0,0,0}, {0,0,0,0}, {-1,0,0,0}},
89 {{0,0,-1,0}, {0,0,0,0}, {1,0,0,0}, {0,0,0,0} } },
90 { {{0,0,0,0}, {0,0,0,1}, {0,0,0,0}, {0,-1,0,0}},
91 {{0,0,0,-1}, {0,0,0,0}, {0,0,0,0}, {1,0,0,0} },
92 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
93 {{0,1,0,0}, {-1,0,0,0}, {0,0,0,0}, {0,0,0,0} } },
94 { {{0,0,0,0}, {0,0,-1,0}, {0,1,0,0}, {0,0,0,0} },
95 {{0,0,1,0}, {0,0,0,0}, {-1,0,0,0}, {0,0,0,0} },
96 {{0,-1,0,0}, {1,0,0,0}, {0,0,0,0}, {0,0,0,0} },
97 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} } } };
98 for (int i=0; i<4; i++) {
99 for (int j=0; j<4; j++) {
100 G[i][j] = GG[i][j];
101 for (int k=0; k<4; k++) {
102 for (int l=0; l<4; l++) {
103 E[i][j][k][l] = EE[i][j][k][l];
104 }
105 }
106 }
107 }
108}
EvtDecayBase::checkSpinParent
void checkSpinParent(EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:520
EvtDecayBase::checkNDaug
void checkNDaug(int d1, int d2=-1)
Definition EvtDecayBase.cc:504
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition EvtDecayBase.cc:482

◆ initProbMax()

void EvtD0ToKSpi0pi0pi0::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 110 of file EvtD0ToKSpi0pi0pi0.cc.

110 {
111 setProbMax(8272);
112}
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition EvtDecayBase.cc:297
The documentation for this class was generated from the following files: