BOSS 7.1.2
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDsToKKpi Class Reference

#include <EvtDsToKKpi.hh>

+ Inheritance diagram for EvtDsToKKpi:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtDsToKKpi()

EvtDsToKKpi::EvtDsToKKpi ( )
inline

Definition at line 12 of file EvtDsToKKpi.hh.

12{}

Referenced by clone().

◆ ~EvtDsToKKpi()

EvtDsToKKpi::~EvtDsToKKpi ( )
virtual

Definition at line 33 of file EvtDsToKKpi.cc.

33{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDsToKKpi::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 39 of file EvtDsToKKpi.cc.

39 {
40 return new EvtDsToKKpi;
41}

◆ decay()

void EvtDsToKKpi::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 119 of file EvtDsToKKpi.cc.

119 {
120/*
121 double maxprob = 0.0;
122 cout<<" calculate maxprob for DsToKKpi: "<<endl;
123 for(int ir=0;ir<=60000000;ir++){
124 p->initializePhaseSpace(getNDaug(),getDaugs());
125 EvtVector4R D1 = p->getDaug(0)->getP4();
126 EvtVector4R D2 = p->getDaug(1)->getP4();
127 EvtVector4R D3 = p->getDaug(2)->getP4();
128 double P1[4], P2[4], P3[4];
129 P1[1] = D1.get(1); P1[2] = D1.get(2); P1[3] = D1.get(3);
130 P2[1] = D2.get(1); P2[2] = D2.get(2); P2[3] = D2.get(3);
131 P3[1] = D3.get(1); P3[2] = D3.get(2); P3[3] = D3.get(3);
132 double value;
133 int g0[6]={1,1,1,1,1,1};
134 int nstates=6;
135 calEvaMy( P1, P2, P3, mass, width, rho, phi, g0, modetype, nstates, value);
136 if(value>maxprob) {
137 maxprob=value;
138 cout << "Max PDF = " << ir << " " << maxprob << endl;
139 }
140 }
141 printf("MAXprob = %.10f\n",maxprob);
142*/
144 EvtVector4R D1 = p->getDaug(0)->getP4();
145 EvtVector4R D2 = p->getDaug(1)->getP4();
146 EvtVector4R D3 = p->getDaug(2)->getP4();
147 double P1[4], P2[4], P3[4];
148 P1[1] = D1.get(1); P1[2] = D1.get(2); P1[3] = D1.get(3);
149 P2[1] = D2.get(1); P2[2] = D2.get(2); P2[3] = D2.get(3);
150 P3[1] = D3.get(1); P3[2] = D3.get(2); P3[3] = D3.get(3);
151
152 double value;
153 int g0[6]={1,1,1,1,1,1};
154 int nstates=6;
155 calEvaMy(P1, P2, P3, mass, width, rho, phi, g0, modetype, nstates, value);
156 setProb(value);
157 return ;
158}
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 EvtDsToKKpi::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 35 of file EvtDsToKKpi.cc.

35 {
36 model_name="DsToKKpi";
37}

◆ init()

void EvtDsToKKpi::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 43 of file EvtDsToKKpi.cc.

43 {
44 // check that there are 0 arguments
45 checkNArg(0);
46 checkNDaug(3);
51
52 phi[1] = 6.1944e+00;
53 phi[2] = 4.7398e+00;
54 phi[3] = 2.9047e+00;
55 phi[4] = 1.0068e+00;
56 phi[5] = 5.8035e-01;
57
58 rho[1] = 1.0963e+00;
59 rho[2] = 2.7818e+00;
60 rho[3] = 1.2570e+00;
61 rho[4] = 7.7351e-01;
62 rho[5] = 5.6277e-01;
63
64 mass[0] = 8.9581e-01;
65 mass[1] = 1.019461e+00;
66 mass[2] = 0.919;
67 mass[3] = 1.4712e+00;
68 mass[4] = 1.7220e+00;
69 mass[5] = 1.3500e+00;
70
71 width[0] = 4.7400e-02;
72 width[1] = 4.2660e-03;
73 width[2] = 0.272;
74 width[3] = 2.7000e-01;
75 width[4] = 1.3500e-01;
76 width[5] = 2.6500e-01;
77
78 modetype[0]= 0;
79 modetype[1]= 1;
80 modetype[2]= 13;
81 modetype[3]= 3;
82 modetype[4]= 4;
83 modetype[5]= 5;
84
85
86 phi[0] = 0;
87 rho[0] = 1;
88
89 //cout << "DsToKKpi :"<< endl;
90 //for (int i=0; i<6; i++) {
91 // cout << i << " rho= " << rho[i] << " phi= " << phi[i] << endl;
92 //}
93
94 mD = 1.86486;
95 mDs = 1.9683;
96 rRes = 1.5;
97 rD = 3.0;
98 mkstr = 0.89594;
99 mk0 = 0.497614;
100 mass_Kaon = 0.49368;
101 mass_Pion = 0.13957;
102 mass_Pi0 = 0.1349766;
103 mass_EtaP = 0.95778;
104 mass_Eta = 0.547862;
105 math_pi = 3.1415926;
106 afRatio= 2.04835; //BABAR // afRatio= 1.23202; //BES2
107 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
108 for (int i=0; i<4; i++) {
109 for (int j=0; j<4; j++) {
110 G[i][j] = GG[i][j];
111 }
112 }
113}
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
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 EvtDsToKKpi::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 115 of file EvtDsToKKpi.cc.

115 {
116 setProbMax(35638.0);
117}
void setProbMax(double prbmx)

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