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

#include <EvtDToKSpieta.hh>

+ Inheritance diagram for EvtDToKSpieta:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtDToKSpieta()

EvtDToKSpieta::EvtDToKSpieta ( )
inline

Definition at line 12 of file EvtDToKSpieta.hh.

12{}

Referenced by clone().

◆ ~EvtDToKSpieta()

EvtDToKSpieta::~EvtDToKSpieta ( )
virtual

Definition at line 32 of file EvtDToKSpieta.cc.

32{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDToKSpieta::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 38 of file EvtDToKSpieta.cc.

38 {
39 return new EvtDToKSpieta;
40}

◆ decay()

void EvtDToKSpieta::decay ( EvtParticle p)
virtual

Implements EvtDecayBase.

Definition at line 121 of file EvtDToKSpieta.cc.

121 {
122 //-----------for max value------------------
123/* double maxprob = 0.0;
124 for(int ir=0;ir<=60000000;ir++){
125 p->initializePhaseSpace(getNDaug(),getDaugs());
126 EvtVector4R ks0 = p->getDaug(0)->getP4();
127 EvtVector4R pic = p->getDaug(1)->getP4();
128 EvtVector4R pi0 = p->getDaug(2)->getP4();
129
130 double Ks0[4],Pic[4],Pi0[4];
131 Ks0[0] = ks0.get(0); Pic[0] = pic.get(0); Pi0[0] = pi0.get(0);
132 Ks0[1] = ks0.get(1); Pic[1] = pic.get(1); Pi0[1] = pi0.get(1);
133 Ks0[2] = ks0.get(2); Pic[2] = pic.get(2); Pi0[2] = pi0.get(2);
134 Ks0[3] = ks0.get(3); Pic[3] = pic.get(3); Pi0[3] = pi0.get(3);
135 double value;
136 calPDF(Ks0, Pic, Pi0, 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 //-----------------------------------------------
146 EvtVector4R ks0 = p->getDaug(0)->getP4();
147 EvtVector4R pic = p->getDaug(1)->getP4();
148 EvtVector4R pi0 = p->getDaug(2)->getP4();
149
150 double Ks0[4],Pic[4],Pi0[4];
151 Ks0[0] = ks0.get(0); Pic[0] = pic.get(0); Pi0[0] = pi0.get(0);
152 Ks0[1] = ks0.get(1); Pic[1] = pic.get(1); Pi0[1] = pi0.get(1);
153 Ks0[2] = ks0.get(2); Pic[2] = pic.get(2); Pi0[2] = pi0.get(2);
154 Ks0[3] = ks0.get(3); Pic[3] = pic.get(3); Pi0[3] = pi0.get(3);
155
156 double value;
157 calPDF(Ks0, Pic, Pi0, value);
158 setProb(value);
159 return;
160}
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 EvtDToKSpieta::getName ( std::string &  name)
virtual

Implements EvtDecayBase.

Definition at line 34 of file EvtDToKSpieta.cc.

34 {
35 model_name="DToKSpieta";
36}

◆ init()

void EvtDToKSpieta::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 42 of file EvtDToKSpieta.cc.

42 {
43 checkNArg(0);
44 checkNDaug(3);
46 ma0 = 0.987; Ga0 = 0.0732;
47 mKn_1430 = 1.423; GKn_1430 = 0.1815;
48 mK1270 = 1.272; mK1400 = 1.403;
49 GK1270 = 0.09; GK1400 = 0.174;
50 phi_omega = -0.02;
51
52 rho_omega = 0.00294;
53
54 rho[0] = 1;//a0980-pieta
55 phi[0] = 0;
56
57 rho[1] = 0.26085;//K1430-kseta
58 phi[1] = 2.6200;
59
60 spin[0]=0;
61 spin[1]=0;
62
63 modetype[0]=23;
64 modetype[1]=13;
65
66 //cout << "Initializing EvtDToKSpieta" << endl;
67 //for (int i=0; i<2; i++) {
68 // cout << i << " rho,phi = " << rho[i] << ", "<< phi[i] << endl;
69 //}
70 mass_Ks=0.4977;
71 mass_Eta=0.547862;
72 mD = 1.86965;
73 rD = 5;
74 metap = 0.95778;
75 mkstr = 0.89594;
76 mk0 = 0.497611;
77 mass_Kaon = 0.493677;
78 mass_Pion = 0.13957;
79 mass_Pi0 = 0.1349768;
80 math_pi = 3.1415926;
81
82 pi = 3.1415926;
83 mpi = 0.13957;
84 g1 = 0.5468;
85 g2 = 0.23;
86
87 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
88 int EE[4][4][4][4] =
89 { { {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
90 {{0,0,0,0}, {0,0,0,0}, {0,0,0,1}, {0,0,-1,0}},
91 {{0,0,0,0}, {0,0,0,-1}, {0,0,0,0}, {0,1,0,0} },
92 {{0,0,0,0}, {0,0,1,0}, {0,-1,0,0}, {0,0,0,0} } },
93 { {{0,0,0,0}, {0,0,0,0}, {0,0,0,-1}, {0,0,1,0} },
94 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
95 {{0,0,0,1}, {0,0,0,0}, {0,0,0,0}, {-1,0,0,0}},
96 {{0,0,-1,0}, {0,0,0,0}, {1,0,0,0}, {0,0,0,0} } },
97 { {{0,0,0,0}, {0,0,0,1}, {0,0,0,0}, {0,-1,0,0}},
98 {{0,0,0,-1}, {0,0,0,0}, {0,0,0,0}, {1,0,0,0} },
99 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
100 {{0,1,0,0}, {-1,0,0,0}, {0,0,0,0}, {0,0,0,0} } },
101 { {{0,0,0,0}, {0,0,-1,0}, {0,1,0,0}, {0,0,0,0} },
102 {{0,0,1,0}, {0,0,0,0}, {-1,0,0,0}, {0,0,0,0} },
103 {{0,-1,0,0}, {1,0,0,0}, {0,0,0,0}, {0,0,0,0} },
104 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} } } };
105 for (int i=0; i<4; i++) {
106 for (int j=0; j<4; j++) {
107 G[i][j] = GG[i][j];
108 for (int k=0; k<4; k++) {
109 for (int l=0; l<4; l++) {
110 E[i][j][k][l] = EE[i][j][k][l];
111 }
112 }
113 }
114 }
115}
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 EvtDToKSpieta::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 117 of file EvtDToKSpieta.cc.

117 {
118 setProbMax(13.00);//setProbMax(1680.0);
119}
void setProbMax(double prbmx)

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