43 {
44
52
54
58 EvtVector4R Qp = (Ppim + Ppi0);
59 EvtVector4R Qm = (Ppip + Ppi0);
60 EvtVector4R Q0 = (Ppip + Ppim);
68
73
74 EvtComplex Atot(0,0);
75
76
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
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
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
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
97 EvtComplex DirPhase(0,2.47);
98 EvtComplex A5(0.78*
exp(DirPhase));
99
100 Atot=A1+A2+A3+A4+A5;
101
105
106 return ;
107
108}
EvtComplex exp(const EvtComplex &c)
void vertex(const EvtComplex &)
static double getWidth(EvtId i)
static double getMeanMass(EvtId i)
static EvtId getId(const std::string &name)
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)