BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDToKSKSK Class Reference

#include <EvtDToKSKSK.hh>

Inheritance diagram for EvtDToKSKSK:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtDToKSKSK()

EvtDToKSKSK::EvtDToKSKSK ( )
inline

Definition at line 12 of file EvtDToKSKSK.hh.

12{}

Referenced by clone().

◆ ~EvtDToKSKSK()

EvtDToKSKSK::~EvtDToKSKSK ( )
virtual

Definition at line 31 of file EvtDToKSKSK.cc.

31{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDToKSKSK::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 35 of file EvtDToKSKSK.cc.

35{ return new EvtDToKSKSK; }

◆ decay()

void EvtDToKSKSK::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 93 of file EvtDToKSKSK.cc.

93 {
94
95 /* double maxprob = 0.0;
96 for(int ir=0;ir<=60000000;ir++){
97 p->initializePhaseSpace(getNDaug(),getDaugs());
98 EvtVector4R D1 = p->getDaug(0)->getP4();
99 EvtVector4R D2 = p->getDaug(1)->getP4();
100 EvtVector4R D3 = p->getDaug(2)->getP4();
101
102 double P1[4], P2[4], P3[4];
103 P1[0] = D1.get(0); P1[1] = D1.get(1); P1[2] = D1.get(2); P1[3] = D1.get(3);
104 P2[0] = D2.get(0); P2[1] = D2.get(1); P2[2] = D2.get(2); P2[3] = D2.get(3);
105 P3[0] = D3.get(0); P3[1] = D3.get(1); P3[2] = D3.get(2); P3[3] = D3.get(3);
106
107 double value;
108 int g0[3]={1,1,1};
109 int spin[3]={0,1,0};
110 int nstates=2;
111 calEva(P1, P2, P3, mass, width, rho, phi, g0, spin, modetype, nstates,
112 value); if (value<0) continue; if(value>maxprob) { maxprob=value; cout << "ir= " << ir <<
113 endl; cout << "double P1[4] = {" << P1[0] <<","<< P1[1] <<","<< P1[2] <<","<< P1[3]
114 <<"};"<< endl; cout << "double P2[4] = {" << P2[0] <<","<< P2[1] <<","<< P2[2] <<","<<
115 P2[3] <<"};"<< endl; cout << "double P3[4] = {" << P3[0] <<","<< P3[1] <<","<< P3[2]
116 <<","<< P3[3] <<"};"<< endl; cout << "MAX====> " << maxprob << endl;
117 }
118 }
119 printf("MAXprob = %.10f\n",maxprob);
120
121 */
123 EvtVector4R D1 = p->getDaug( 0 )->getP4();
124 EvtVector4R D2 = p->getDaug( 1 )->getP4();
125 EvtVector4R D3 = p->getDaug( 2 )->getP4();
126
127 double P1[4], P2[4], P3[4];
128 P1[0] = D1.get( 0 );
129 P1[1] = D1.get( 1 );
130 P1[2] = D1.get( 2 );
131 P1[3] = D1.get( 3 );
132 P2[0] = D2.get( 0 );
133 P2[1] = D2.get( 1 );
134 P2[2] = D2.get( 2 );
135 P2[3] = D2.get( 3 );
136 P3[0] = D3.get( 0 );
137 P3[1] = D3.get( 1 );
138 P3[2] = D3.get( 2 );
139 P3[3] = D3.get( 3 );
140
141 double value;
142 int g0[2] = { 1, 1 };
143 int spin[2] = { 0, 0 };
144 int nstates = 2;
145 calEva( P1, P2, P3, mass, width, rho, phi, g0, spin, modetype, nstates, value );
146 setProb( value );
147
148 return;
149}
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 EvtDToKSKSK::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 33 of file EvtDToKSKSK.cc.

33{ model_name = "DToKSKSK"; }

◆ init()

void EvtDToKSKSK::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 37 of file EvtDToKSKSK.cc.

37 {
38 checkNArg( 0 );
39 checkNDaug( 3 );
41
42 mass[0] = 0.980;
43 mass[1] = 0.980;
44 width[0] = 0.0972;
45 width[1] = 0.0972;
46
47 rho[0] = 1.0;
48 rho[1] = 1.1141e+01;
49 phi[0] = 0.00;
50 phi[1] = -5.6220e+00;
51
52 // spin[0] = 0;
53 // spin[1] = 1;
54 // spin[2] = 1;
55 modetype[0] = 12;
56 modetype[1] = 122;
57
58 // std::cout << "Initializing EvtDToKSKSK" << std::endl;
59 // for (int i=0; i<2; i++) {
60 // cout << i << " rho,phi = " << rho[i] << ", "<< phi[i] << endl;
61 // }
62
63 math_K = 0.493677;
64 GS1 = 0.636619783;
65 GS2 = 0.23459086;
66 GS3 = 0.1591549458;
67 GS4 = 0.078196955;
68
69 mD0 = 1.86966;
70 mK0 = 0.497611;
71 mKa = 0.493677;
72 mPi = 0.13957;
73 mK02 = 0.237616707;
74 mPi2 = 0.01947978;
75 mass_EtaP = 0.95778;
76 mass_Kaon = 0.49368;
77
78 math_pi = 3.1415926;
79 mass_Pion = 0.13957;
80 mass_Pion2 = 0.0194797849;
81 mass_2Pion = 0.27914;
82 math_2pi = 6.2831852;
83 rD2 = 25.0;
84 rRes2 = 9.0;
85 g2 = 0.23;
86
87 rho_omega = 0.00294;
88 phi_omega = -0.02;
89}
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 EvtDToKSKSK::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 91 of file EvtDToKSKSK.cc.

91{ setProbMax( 878.9 ); }
void setProbMax(double prbmx)

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