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

#include <EvtDTopipienu.hh>

+ Inheritance diagram for EvtDTopipienu:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtDTopipienu()

EvtDTopipienu::EvtDTopipienu ( )
inline

Definition at line 12 of file EvtDTopipienu.hh.

12{}

Referenced by clone().

◆ ~EvtDTopipienu()

EvtDTopipienu::~EvtDTopipienu ( )
virtual

Definition at line 31 of file EvtDTopipienu.cc.

31{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDTopipienu::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 37 of file EvtDTopipienu.cc.

37 {
38 return new EvtDTopipienu;
39}

◆ decay()

void EvtDTopipienu::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 103 of file EvtDTopipienu.cc.

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

◆ getName()

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

Implements EvtDecayBase.

Definition at line 33 of file EvtDTopipienu.cc.

33 {
34 model_name="DTopipienu";
35}

◆ init()

void EvtDTopipienu::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 41 of file EvtDTopipienu.cc.

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

Reimplemented from EvtDecayBase.

Definition at line 99 of file EvtDTopipienu.cc.

99 {
100 setProbMax(ProbMax);
101}
void setProbMax(double prbmx)

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