BOSS 7.1.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtPhiDalitz Class Reference

#include <EvtPhiDalitz.hh>

+ Inheritance diagram for EvtPhiDalitz:

Public Member Functions

 EvtPhiDalitz ()
 
virtual ~EvtPhiDalitz ()
 
void getName (std::string &name)
 
EvtDecayBaseclone ()
 
void init ()
 
void decay (EvtParticle *p)
 
- Public Member Functions inherited from EvtDecayAmp
void makeDecay (EvtParticle *p)
 
void setWeight (double weight)
 
void vertex (const EvtComplex &amp)
 
void vertex (int i1, const EvtComplex &amp)
 
void vertex (int i1, int i2, const EvtComplex &amp)
 
void vertex (int i1, int i2, int i3, const EvtComplex &amp)
 
void vertex (int *i1, const EvtComplex &amp)
 
virtual ~EvtDecayAmp ()
 
- Public Member Functions inherited from EvtDecayBase
virtual void getName (std::string &name)=0
 
virtual void decay (EvtParticle *p)=0
 
virtual void makeDecay (EvtParticle *p)=0
 
virtual EvtDecayBaseclone ()=0
 
virtual void init ()
 
virtual void initProbMax ()
 
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 EvtDecayAmp
EvtAmp _amp2
 
- Protected Attributes inherited from EvtDecayBase
bool _daugsDecayedByParentModel
 

Detailed Description

Definition at line 9 of file EvtPhiDalitz.hh.

Constructor & Destructor Documentation

◆ EvtPhiDalitz()

EvtPhiDalitz::EvtPhiDalitz ( )
inline

Definition at line 13 of file EvtPhiDalitz.hh.

13{}

Referenced by clone().

◆ ~EvtPhiDalitz()

EvtPhiDalitz::~EvtPhiDalitz ( )
virtual

Definition at line 13 of file EvtPhiDalitz.cc.

13{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtPhiDalitz::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 22 of file EvtPhiDalitz.cc.

22 {
23
24 return new EvtPhiDalitz;
25
26}

◆ decay()

void EvtPhiDalitz::decay ( EvtParticle p)
virtual

Implements EvtDecayBase.

Definition at line 43 of file EvtPhiDalitz.cc.

43 {
44
45 EvtId PIP=EvtPDL::getId("pi+");
46 EvtId PIM=EvtPDL::getId("pi-");
47 EvtId PIZ=EvtPDL::getId("pi0");
48 EvtId RHOP=EvtPDL::getId("rho+");
49 EvtId RHOM=EvtPDL::getId("rho-");
50 EvtId RHOZ=EvtPDL::getId("rho0");
51 EvtId OMEGA=EvtPDL::getId("omega");
52
54
55 EvtVector4R Ppip = p->getDaug(0)->getP4();
56 EvtVector4R Ppim = p->getDaug(1)->getP4();
57 EvtVector4R Ppi0 = p->getDaug(2)->getP4();
58 EvtVector4R Qp = (Ppim + Ppi0);
59 EvtVector4R Qm = (Ppip + Ppi0);
60 EvtVector4R Q0 = (Ppip + Ppim);
61 double m2_pip = pow(EvtPDL::getMeanMass(PIP),2);
62 double m2_pim = pow(EvtPDL::getMeanMass(PIM),2);
63 double m2_pi0 = pow(EvtPDL::getMeanMass(PIZ),2);
64 double M2rhop = pow(EvtPDL::getMeanMass(RHOP),2);
65 double M2rhom = pow(EvtPDL::getMeanMass(RHOM),2);
66 double M2rho0 = pow(EvtPDL::getMeanMass(RHOZ),2);
67 double M2omega = pow(EvtPDL::getMeanMass(OMEGA),2);
68
69 double Wrhop = EvtPDL::getWidth(RHOP);
70 double Wrhom = EvtPDL::getWidth(RHOM);
71 double Wrho0 = EvtPDL::getWidth(RHOZ);
72 double Womega = EvtPDL::getWidth(OMEGA);
73
74 EvtComplex Atot(0,0);
75
76 //Rho+ Risonance Amplitude
77 double Gp = Wrhop*pow(((Qp.mass2()-m2_pim-m2_pi0)/2-M2rhop/4)/(M2rhop/4-(m2_pim+m2_pi0)/2),3/2)*(M2rhop/Qp.mass2());
78 EvtComplex Drhop((Qp.mass2()-M2rhop),Qp.mass()*Gp);
79 EvtComplex A1(M2rhop/Drhop);
80
81 //Rho- Risonance Amplitude
82 double Gm = Wrhom*pow(((Qm.mass2()-m2_pip-m2_pi0)/2-M2rhom/4)/(M2rhom/4-(m2_pip+m2_pi0)/2),3/2)*(M2rhom/Qm.mass2());
83 EvtComplex Drhom((Qm.mass2()-M2rhom),Qm.mass()*Gm);
84 EvtComplex A2(M2rhom/Drhom);
85
86 //Rho0 Risonance Amplitude
87 double G0 = Wrho0*pow(((Q0.mass2()-m2_pip-m2_pim)/2-M2rho0/4)/(M2rho0/4-(m2_pip+m2_pim)/2),3/2)*(M2rho0/Q0.mass2());
88 EvtComplex Drho0((Q0.mass2()-M2rho0),Q0.mass()*G0);
89 EvtComplex A3(M2rho0/Drho0);
90
91 //Omega Risonance Amplitude
92 EvtComplex OmegaPhase(0,-0.22);
93 EvtComplex DOmega((Q0.mass2()-M2omega),Q0.mass()*Womega);
94 EvtComplex A4(0.0071*M2omega*exp(OmegaPhase)/DOmega);
95
96 //Direct Decay Amplitude
97 EvtComplex DirPhase(0,2.47);
98 EvtComplex A5(0.78*exp(DirPhase));
99
100 Atot=A1+A2+A3+A4+A5;
101
102 vertex(0,Atot);
103 vertex(1,Atot);
104 vertex(2,Atot);
105
106 return ;
107
108}
EvtComplex exp(const EvtComplex &c)
Definition: EvtComplex.hh:252
void vertex(const EvtComplex &amp)
Definition: EvtDecayAmp.hh:37
EvtId * getDaugs()
Definition: EvtDecayBase.hh:65
Definition: EvtId.hh:27
static double getWidth(EvtId i)
Definition: EvtPDL.hh:54
static double getMeanMass(EvtId i)
Definition: EvtPDL.hh:45
static EvtId getId(const std::string &name)
Definition: EvtPDL.cc:287
const EvtVector4R & getP4() const
Definition: EvtParticle.cc:120
EvtParticle * getDaug(int i)
Definition: EvtParticle.cc:84
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
double mass() const
Definition: EvtVector4R.cc:39
double mass2() const
Definition: EvtVector4R.hh:116

◆ getName()

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

Implements EvtDecayBase.

Definition at line 15 of file EvtPhiDalitz.cc.

15 {
16
17 model_name="PHI_DALITZ";
18
19}

◆ init()

void EvtPhiDalitz::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 28 of file EvtPhiDalitz.cc.

28 {
29
30 // check that there are 0 arguments
31 checkNArg(0);
32 checkNDaug(3);
33
35
39}
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)

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