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

#include <EvtDToKpienu.hh>

Inheritance diagram for EvtDToKpienu:

Public Member Functions

 EvtDToKpienu ()
virtual ~EvtDToKpienu ()
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 EvtDToKpienu.hh.

Constructor & Destructor Documentation

◆ EvtDToKpienu()

EvtDToKpienu::EvtDToKpienu ( )
inline

Definition at line 11 of file EvtDToKpienu.hh.

11{}

Referenced by clone().

◆ ~EvtDToKpienu()

EvtDToKpienu::~EvtDToKpienu ( )
virtual

Definition at line 30 of file EvtDToKpienu.cc.

30{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDToKpienu::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 34 of file EvtDToKpienu.cc.

34{ return new EvtDToKpienu; }

◆ decay()

void EvtDToKpienu::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 94 of file EvtDToKpienu.cc.

94 {
95 /*
96 double maxprob = 0.0;
97 for(int ir=0;ir<=60000000;ir++){
98 p->initializePhaseSpace(getNDaug(),getDaugs());
99 EvtVector4R _K = p->getDaug(0)->getP4();
100 EvtVector4R _pi = p->getDaug(1)->getP4();
101 EvtVector4R _e = p->getDaug(2)->getP4();
102 EvtVector4R _nu = p->getDaug(3)->getP4();
103 int pid = EvtPDL::getStdHep(p->getDaug(0)->getId());
104 int charm;
105 if(pid == -321) charm = 1;
106 else charm = -1;
107 double m2, q2, cosV, cosL, chi;
108 KinVGen(_K, _pi, _e, _nu, charm, m2, q2, cosV, cosL, chi);
109 double _prob = calPDF(m2, q2, cosV, cosL, chi);
110 if(_prob>maxprob) {
111 maxprob=_prob;
112 std::cout << "Max PDF = " << ir << " charm= " << charm << " prob= " << _prob <<
113 std::endl;
114 }
115 }
116 std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
117 */
119 EvtVector4R K = p->getDaug( 0 )->getP4();
120 EvtVector4R pi = p->getDaug( 1 )->getP4();
121 EvtVector4R e = p->getDaug( 2 )->getP4();
122 EvtVector4R nu = p->getDaug( 3 )->getP4();
123
124 int pid = EvtPDL::getStdHep( p->getDaug( 0 )->getId() );
125 int charm;
126 if ( pid == -321 ) charm = 1;
127 else charm = -1;
128 double m2, q2, cosV, cosL, chi;
129 KinVGen( K, pi, e, nu, charm, m2, q2, cosV, cosL, chi );
130 double prob = calPDF( m2, q2, cosV, cosL, chi );
131 setProb( prob );
132 return;
133}
double pi
EvtId * getDaugs()
void setProb(double prob)
static int getStdHep(EvtId id)
Definition EvtPDL.hh:61
EvtId getId() const
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
double double * m2
Definition qcdloop1.h:83

◆ getName()

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

Implements EvtDecayBase.

Definition at line 32 of file EvtDToKpienu.cc.

32{ model_name = "DToKpienu"; }

◆ init()

void EvtDToKpienu::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 36 of file EvtDToKpienu.cc.

36 {
37 checkNArg( 0 );
38 checkNDaug( 4 );
42
43 // std::cout << "Initializing EvtDToKpienu" << std::endl;
44 nAmps = 2;
45
46 rS = -11.57; // S-wave
47 rS1 = 0.08;
48 a_delta = 1.94;
49 b_delta = -0.81;
50 m0_1430_S = 1.425;
51 width0_1430_S = 0.270;
52 type[0] = 0;
53
54 mV = 1.81;
55 mA = 2.61;
56 V_0 = 1.411;
57 A1_0 = 1;
58 A2_0 = 0.788;
59
60 m0 = 0.8946; // P-wave K*
61 width0 = 0.04642;
62 rBW = 3.07;
63 rho = 1;
64 phi = 0;
65 type[1] = 1;
66
67 m0_1410 = 1.414; // P-wave K*(1410)
68 width0_1410 = 0.232;
69 rho_1410 = 0.1;
70 phi_1410 = 0.;
71 type[2] = 2;
72
73 TV_0 = 1; // D-wave K*2(1430)
74 T1_0 = 1;
75 T2_0 = 1;
76 m0_1430 = 1.4324;
77 width0_1430 = 0.109;
78 rho_1430 = 15;
79 phi_1430 = 0;
80 type[3] = 3;
81
82 mD = 1.86962;
83 mPi = 0.13957;
84 mK = 0.49368;
85 Pi = atan2( 0.0, -1.0 );
86 root2 = sqrt( 2. );
87 root2d3 = sqrt( 2. / 3 );
88 root1d2 = sqrt( 0.5 );
89 root3d2 = sqrt( 1.5 );
90}
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 EvtDToKpienu::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 92 of file EvtDToKpienu.cc.

92{ setProbMax( 22750.0 ); }
void setProbMax(double prbmx)

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