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

#include <EvtDsToKKmunu.hh>

Inheritance diagram for EvtDsToKKmunu:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtDsToKKmunu()

EvtDsToKKmunu::EvtDsToKKmunu ( )
inline

Definition at line 11 of file EvtDsToKKmunu.hh.

11{}

Referenced by clone().

◆ ~EvtDsToKKmunu()

EvtDsToKKmunu::~EvtDsToKKmunu ( )
virtual

Definition at line 27 of file EvtDsToKKmunu.cc.

27{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDsToKKmunu::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 31 of file EvtDsToKKmunu.cc.

31{ return new EvtDsToKKmunu; }

◆ decay()

void EvtDsToKKmunu::decay ( EvtParticle * p)
virtual

test test test PDF

the end

Implements EvtDecayBase.

Definition at line 92 of file EvtDsToKKmunu.cc.

92 {
93 /*
94 double maxprob = 0.0;
95 for(int ir=0;ir<=2000000;ir++){
96 p->initializePhaseSpace(getNDaug(),getDaugs());
97 EvtVector4R _K = p->getDaug(0)->getP4();
98 EvtVector4R _pi = p->getDaug(1)->getP4();
99 EvtVector4R _e = p->getDaug(2)->getP4();
100 EvtVector4R _nu = p->getDaug(3)->getP4();
101 int pid = EvtPDL::getStdHep(p->getDaug(0)->getId());
102 int charm;
103 if(pid == -321) charm = 1;
104 else charm = -1;
105 double m2, q2, cosV, cosL, chi;
106 KinVGen(_K, _pi, _e, _nu, charm, m2, q2, cosV, cosL, chi);
107 double _prob = calPDF(m2, q2, cosV, cosL, chi);
108 if(_prob>maxprob) {
109 maxprob=_prob;
110 std::cout << "Max PDF = " << ir << " charm= " << charm << " prob= " << _prob <<
111 std::endl;
112 }
113 }
114 std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
115 */
116 /////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
117 /////////////////////////////////////////////////////////////////////////////////////////
118 /// test test test PDF
119 ///////////////////////////////////////////////////////////////////////////////////////////////////
120
121 /*
122 p->initializePhaseSpace(getNDaug(),getDaugs());
123 EvtVector4R _K = p->getDaug(0)->getP4();
124 EvtVector4R _pi = p->getDaug(1)->getP4();
125 EvtVector4R _e = p->getDaug(2)->getP4();
126 EvtVector4R _nu = p->getDaug(3)->getP4();
127 int pid = EvtPDL::getStdHep(p->getDaug(0)->getId());
128 int charm;
129 if(pid == -321) charm = 1;
130 else charm = -1;
131
132 double m2, q2, cosV, cosL, chi;
133 KinVGen(_K, _pi, _e, _nu, charm, m2, q2, cosV, cosL, chi);
134 //m2 = 0.7317376201248132;
135 //q2 = 0.7930526827833025;
136 //cosV = -0.1044604580373859;
137 //cosL = 0.0503218879658880;
138 //chi = 1.9753715152960665;
139 //double _prob = calPDF(m2, q2, cosV, cosL, chi);
140
141 double _prob = calPDF(0.7317376201248132, 0.7930526827833025, -0.1044604580373859,
142 0.0503218879658880, 1.9753715152960665); std::cout << "PDF PDF1 !!!!!!!!!! = " << _prob <<
143 std::endl; _prob = calPDF(0.8253460062251410, 0.7590531075741469, -0.5990247196671010,
144 -0.9155344197677775, -1.3296895383709690); std::cout << "PDF PDF2 !!!!!!!!!! = " << _prob
145 << std::endl; _prob = calPDF(0.8161884101984322, 0.2477364661970792, -0.2764898614626146,
146 0.2826827644931852, 0.9307086619444224); std::cout << "PDF PDF3 !!!!!!!!!! = " << _prob <<
147 std::endl; _prob = calPDF(1.0912051424949705, 0.1906228414935633, 0.8132818897807017,
148 -0.1032400648183377, 0.9686208243587908); std::cout << "PDF PDF4 !!!!!!!!!! = " << _prob
149 << std::endl; _prob = calPDF(0.8079301573422222, 0.9019596037894472, 0.7789653251154771,
150 0.5324409155087387, -2.4926364105516909); std::cout << "PDF PDF5 !!!!!!!!!! = " << _prob
151 << std::endl; _prob = calPDF(0.6698684771433824, 0.1744218225194828, 0.7286179682544673,
152 -0.1705130677785587, 1.8691658148034282); std::cout << "PDF PDF6 !!!!!!!!!! = " << _prob
153 << std::endl;
154
155 double _prob = calPDF(1.0754665082688288, 0.6798367607443577, -0.1726785245179900,
156 -0.8252479196790240, 2.8955919246019519); std::cout << "PDF PDF1 !!!!!!!!!! = " << _prob
157 << std::endl; _prob = calPDF(1.0386029533289431, 0.4389802155039476, -0.5990247196671010,
158 -0.9155344197677775, -1.3296895383709690); std::cout << "PDF PDF2 !!!!!!!!!! = " << _prob
159 << std::endl; _prob = calPDF(1.0307383656072204, 0.2894701649730593, -0.2764898614626146,
160 0.2826827644931852, 0.9307086619444224); std::cout << "PDF PDF3 !!!!!!!!!! = " << _prob <<
161 std::endl; _prob = calPDF(1.0501028363282328, 0.1987927934599780, 0.8132818897807017,
162 -0.1032400648183377, 0.9686208243587908); std::cout << "PDF PDF4 !!!!!!!!!! = " << _prob
163 << std::endl; _prob = calPDF(1.0433160107621029, 0.2206250677498859, 0.7789653251154771,
164 0.5324409155087387, -2.4926364105516909); std::cout << "PDF PDF5 !!!!!!!!!! = " << _prob
165 << std::endl; _prob = calPDF(1.0493574400231911, 0.2460442691281926, 0.7286179682544673,
166 -0.1705130677785587, 1.8691658148034282);
167
168 */
169 //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
170 ///////////////////////////////////////////////////////////////////////////////////////////////
171 /// the end
172 /////////////////////////////////////////////////////////////////////////////////////////////////////////
173
174 p->initializePhaseSpace( getNDaug(), getDaugs() );
175 EvtVector4R K = p->getDaug( 0 )->getP4();
176 EvtVector4R pi = p->getDaug( 1 )->getP4();
177 EvtVector4R e = p->getDaug( 2 )->getP4();
178 EvtVector4R nu = p->getDaug( 3 )->getP4();
179
180 int pid = EvtPDL::getStdHep( p->getDaug( 0 )->getId() );
181 int charm;
182 if ( pid == -321 ) charm = 1;
183 else charm = -1;
184 double m2, q2, cosV, cosL, chi;
185 KinVGen( K, pi, e, nu, charm, m2, q2, cosV, cosL, chi );
186 double prob = calPDF( m2, q2, cosV, cosL, chi );
187 setProb( prob );
188 return;
189}
pi
double pi
Definition EvtPhokharaDef.hh:17
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition EvtDecayBase.hh:68
EvtDecayProb::setProb
void setProb(double prob)
Definition EvtDecayProb.hh:40
EvtPDL::getStdHep
static int getStdHep(EvtId id)
Definition EvtPDL.hh:61
EvtParticle::getId
EvtId getId() const
Definition EvtParticle.cc:100
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition EvtParticle.cc:108
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition EvtParticle.cc:80
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition EvtParticle.cc:1042
m2
double double * m2
Definition qcdloop1.h:83

◆ getName()

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

Implements EvtDecayBase.

Definition at line 29 of file EvtDsToKKmunu.cc.

29{ model_name = "DsToKKmunu"; }

◆ init()

void EvtDsToKKmunu::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 33 of file EvtDsToKKmunu.cc.

33 {
34 checkNArg( 0 );
35 checkNDaug( 4 );
36 checkSpinParent( EvtSpinType::SCALAR );
37 checkSpinDaughter( 0, EvtSpinType::SCALAR );
38 checkSpinDaughter( 1, EvtSpinType::SCALAR );
39
40 // std::cout << "EvtDsToKKmunu ==> Initialization !" << std::endl;
41 nAmps = 2;
42
43 // rS = -11.27; // S-wave
44 // rS1 = 0.08;
45 // a_delta = 1.58;
46 // b_delta = -0.81;
47 // m0_1430_S = 1.425;
48 // width0_1430_S = 0.270;
49 // type[0] = 0;
50
51 mV = 2.1;
52 mA = 2.5;
53 V_0 = 1.58;
54 A1_0 = 1;
55 A2_0 = 0.71;
56
57 m0 = 1.01946; // P-wave phi
58 width0 = 0.004249;
59 rBW = 3.0;
60 rho = 1;
61 phi = 0;
62 type[1] = 1;
63
64 m0_1410 = 1.414; // P-wave K*(1410)
65 width0_1410 = 0.232;
66 rho_1410 = 0.1;
67 phi_1410 = 0.;
68 type[2] = 2;
69
70 TV_0 = 1; // D-wave K*2(1430)
71 T1_0 = 1;
72 T2_0 = 1;
73 m0_1430 = 1.4324;
74 width0_1430 = 0.109;
75 rho_1430 = 15;
76 phi_1430 = 0;
77 type[3] = 3;
78
79 mD = 1.96834; // Ds Mass
80 mPi = 0.49368; // K+ Mass
81 mK = 0.49368; // K- Mass
82
83 Pi = atan2( 0.0, -1.0 );
84 root2 = sqrt( 2. );
85 root2d3 = sqrt( 2. / 3 );
86 root1d2 = sqrt( 0.5 );
87 root3d2 = sqrt( 1.5 );
88}
EvtDecayBase::checkSpinDaughter
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:492
EvtDecayBase::checkSpinParent
void checkSpinParent(EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:479
EvtDecayBase::checkNDaug
void checkNDaug(int d1, int d2=-1)
Definition EvtDecayBase.cc:465
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition EvtDecayBase.cc:449

◆ initProbMax()

void EvtDsToKKmunu::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 90 of file EvtDsToKKmunu.cc.

90{ setProbMax( 462526.1 ); }
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition EvtDecayBase.cc:261
The documentation for this class was generated from the following files: