BOSS 7.1.1
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtD0Toa0enu Class Reference

#include <EvtD0Toa0enu.hh>

+ Inheritance diagram for EvtD0Toa0enu:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtD0Toa0enu()

EvtD0Toa0enu::EvtD0Toa0enu ( )
inline

Definition at line 12 of file EvtD0Toa0enu.hh.

12{}

Referenced by clone().

◆ ~EvtD0Toa0enu()

EvtD0Toa0enu::~EvtD0Toa0enu ( )
virtual

Definition at line 30 of file EvtD0Toa0enu.cc.

30{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtD0Toa0enu::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 36 of file EvtD0Toa0enu.cc.

36 {
37 return new EvtD0Toa0enu;
38}

◆ decay()

void EvtD0Toa0enu::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 86 of file EvtD0Toa0enu.cc.

86 {
87/*
88 double maxprob = 0.0;
89 for(int ir=0;ir<=60000000;ir++){
90 p->initializePhaseSpace(getNDaug(),getDaugs());
91 EvtVector4R _pi1 = p->getDaug(0)->getP4();
92 EvtVector4R _pi2 = p->getDaug(1)->getP4();
93 EvtVector4R _e = p->getDaug(2)->getP4();
94 EvtVector4R _nu = p->getDaug(3)->getP4();
95
96 int pid = EvtPDL::getStdHep(p->getDaug(2)->getId());
97 int charm;
98 if(pid == 11) charm = -1;
99 else charm = 1;
100 double m2, q2, cosV, cosL, chi;
101 KinVGen(_pi1, _pi2, _e, _nu, charm, m2, q2, cosV, cosL, chi);
102 double _prob = calPDF(m2, q2, cosV, cosL, chi);
103 if(_prob>maxprob) {
104 maxprob=_prob;
105 std::cout << "Max PDF = " << ir << " charm= " << charm << " prob= " << _prob << std::endl;
106 }
107 }
108 std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
109*/
111 EvtVector4R pi1 = p->getDaug(0)->getP4();
112 EvtVector4R pi2 = p->getDaug(1)->getP4();
113 EvtVector4R e = p->getDaug(2)->getP4();
114 EvtVector4R nu = p->getDaug(3)->getP4();
115
116 int pid = EvtPDL::getStdHep(p->getDaug(2)->getId());
117 int charm;
118 if(pid==11||pid==13) charm =-1;
119 else charm = 1;
120 double m2, q2, cosV, cosL, chi;
121 KinVGen(pi1, pi2, e, nu, charm, m2, q2, cosV, cosL, chi);
122 double prob = calPDF(m2, q2, cosV, cosL, chi);
123 setProb(prob);
124 return;
125}
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 EvtD0Toa0enu::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 32 of file EvtD0Toa0enu.cc.

32 {
33 model_name="D0Toa0enu";
34}

◆ init()

void EvtD0Toa0enu::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 40 of file EvtD0Toa0enu.cc.

40 {
41 checkNArg(0);
42 checkNDaug(4);
46
47
48 mode = 0; // D0
49 ProbMax = 15.55;
50
51 //std::cout << "Initializing EvtD0Toa0enu (D0 -> a0(980)- e+ nu_e) ProbMax = " << ProbMax << std::endl;
52
53 double m_Kaon = 0.493677;
54 m2_Kaon = m_Kaon*m_Kaon;
55 double m_K0 = 0.497614;
56 m2_K0 = m_K0*m_K0;
57 double m_eta = 0.547853;
58 m2_eta = m_eta*m_eta;
59 double m_Pion = 0.13957;
60 m2_Pion = m_Pion*m_Pion;
61 double m_Pi0 = 0.1349766;
62 m2_Pi0 = m_Pi0*m_Pi0;
63 m_D0 = 1.86486;
64 m2_D0 = m_D0*m_D0;
65 m_D = 1.86962;
66 m2_D = m_D*m_D;
67
68 ciR = EvtComplex(1.0, 0.0);
69 ciM = EvtComplex(0.0, 1.0);
70
71 ///////////// for D to a0 e nu
72 double mAD0 = 2.42;
73 m2AD0 = mAD0*mAD0;
74 double mAD = 2.42;
75 m2AD = mAD*mAD;
76 double m_a0 = 0.990;
77 m2_a0 = m_a0*m_a0;
78 flatte_g1 = 0.341;
79 flatte_g2 = 0.341*0.892;
80}
*********DOUBLE PRECISION m_pi INTEGER m_lenwt !max no of aux weights INTEGER m_phmax !maximum photon multiplicity ISR FSR *DOUBLE COMPLEX m_Pauli4 DOUBLE COMPLEX m_AmpBorn DOUBLE COMPLEX m_AmpBoxy DOUBLE COMPLEX m_AmpBorn1 DOUBLE COMPLEX m_AmpBorn2 DOUBLE COMPLEX m_AmpExpo2p DOUBLE COMPLEX m_Rmat DOUBLE COMPLEX m_BoxGZut !DOUBLE COMPLEX m_F1finPair2 !DOUBLE PRECISION m_Vcut DOUBLE PRECISION m_Alfinv DOUBLE PRECISION m_Lorin1 DOUBLE PRECISION m_Lorin2 DOUBLE PRECISION m_eta
Definition GPS.h:30
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 EvtD0Toa0enu::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 82 of file EvtD0Toa0enu.cc.

82 {
83 setProbMax(ProbMax);
84}
void setProbMax(double prbmx)

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