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

#include <EvtD0Topipienu.hh>

Inheritance diagram for EvtD0Topipienu:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtD0Topipienu()

EvtD0Topipienu::EvtD0Topipienu ( )
inline

Definition at line 11 of file EvtD0Topipienu.hh.

11{}

Referenced by clone().

◆ ~EvtD0Topipienu()

EvtD0Topipienu::~EvtD0Topipienu ( )
virtual

Definition at line 31 of file EvtD0Topipienu.cc.

31{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtD0Topipienu::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 35 of file EvtD0Topipienu.cc.

35{ return new EvtD0Topipienu; }

◆ decay()

void EvtD0Topipienu::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 98 of file EvtD0Topipienu.cc.

98 {
99 /*
100 double maxprob = 0.0;
101 for(int ir=0;ir<=60000000;ir++){
102 p->initializePhaseSpace(getNDaug(),getDaugs());
103 EvtVector4R _pi1 = p->getDaug(0)->getP4();
104 EvtVector4R _pi2 = p->getDaug(1)->getP4();
105 EvtVector4R _e = p->getDaug(2)->getP4();
106 EvtVector4R _nu = p->getDaug(3)->getP4();
107
108 int pid = EvtPDL::getStdHep(p->getDaug(0)->getId());
109 int charm;
110 if(pid == -211) charm = 1;
111 else charm = -1;
112 double m2, q2, cosV, cosL, chi;
113 KinVGen(_pi1, _pi2, _e, _nu, charm, m2, q2, cosV, cosL, chi);
114 double _prob = calPDF(m2, q2, cosV, cosL, chi);
115 if(_prob>maxprob) {
116 maxprob=_prob;
117 std::cout << "Max PDF = " << ir << " charm= " << charm << " prob= " << _prob <<
118 std::endl;
119 }
120 }
121 std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
122 */
124 EvtVector4R pi1 = p->getDaug( 0 )->getP4();
125 EvtVector4R pi2 = p->getDaug( 1 )->getP4();
126 EvtVector4R e = p->getDaug( 2 )->getP4();
127 EvtVector4R nu = p->getDaug( 3 )->getP4();
128
129 int pid = EvtPDL::getStdHep( p->getDaug( 0 )->getId() );
130 int charm;
131 if ( pid == -211 ) charm = 1;
132 else charm = -1;
133 double m2, q2, cosV, cosL, chi;
134 KinVGen( pi1, pi2, e, nu, charm, m2, q2, cosV, cosL, chi );
135 double prob = calPDF( m2, q2, cosV, cosL, chi );
136 setProb( prob );
137 return;
138}
character *LEPTONflag integer iresonances real pi2
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 EvtD0Topipienu::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 33 of file EvtD0Topipienu.cc.

33{ model_name = "D0Topipienu"; }

◆ init()

void EvtD0Topipienu::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 37 of file EvtD0Topipienu.cc.

37 {
38 checkNArg( 0 );
39 checkNDaug( 4 );
43
44 first = 0;
45 last = 1;
46 ProbMax = 405000;
47
48 // cout << "Initializing EvtD0Topipienu: ProbMax = " << ProbMax << endl;
49
50 type[0] = 0;
51 type[1] = 1;
52 type[2] = 2;
53 type[3] = 3;
54
55 mV = 2.01;
56 mA = 2.42;
57 V_0 = 1.6948;
58 A1_0 = 1;
59 A2_0 = 0.84489;
60
61 m0 = 0.77526;
62 width0 = 0.14910;
63 rBW = 3.0;
64 rho = 1.0;
65 phi = 0.0;
66 BF = 1.0;
67
68 m0_omega = 0.78265;
69 width0_omega = 0.00849;
70 rho_omega = 0.12902;
71 phi_omega = 2.9285;
72 BF_omega = 1.0;
73
74 m0_S = 0.953;
75 rho_S = 135.27;
76 phi_S = 3.4044;
77
78 Dp_mD = 1.86962;
79 Dp_mPi1 = 0.13957;
80 Dp_mPi2 = 0.13957;
81 D0_mD = 1.86486;
82 D0_mPi1 = 0.13957;
83 D0_mPi2 = 0.1349766;
84
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
91 mKa = 0.493677;
92 mPi = 0.13957;
93 mEt = 0.547853;
94}
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 EvtD0Topipienu::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 96 of file EvtD0Topipienu.cc.

96{ setProbMax( ProbMax ); }
void setProbMax(double prbmx)

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