32 {
33
41
43
47 EvtVector4R Qp = ( Ppim + Ppi0 );
48 EvtVector4R Qm = ( Ppip + Ppi0 );
49 EvtVector4R Q0 = ( Ppip + Ppim );
57
62
63 EvtComplex Atot( 0, 0 );
64
65
66 double Gp = Wrhop *
67 pow( ( ( Qp.
mass2() - m2_pim - m2_pi0 ) / 2 - M2rhop / 4 ) /
68 ( M2rhop / 4 - ( m2_pim + m2_pi0 ) / 2 ),
69 3 / 2 ) *
70 ( M2rhop / Qp.
mass2() );
71 EvtComplex Drhop( ( Qp.
mass2() - M2rhop ), Qp.
mass() * Gp );
72 EvtComplex A1( M2rhop / Drhop );
73
74
75 double Gm = Wrhom *
76 pow( ( ( Qm.
mass2() - m2_pip - m2_pi0 ) / 2 - M2rhom / 4 ) /
77 ( M2rhom / 4 - ( m2_pip + m2_pi0 ) / 2 ),
78 3 / 2 ) *
79 ( M2rhom / Qm.
mass2() );
80 EvtComplex Drhom( ( Qm.
mass2() - M2rhom ), Qm.
mass() * Gm );
81 EvtComplex A2( M2rhom / Drhom );
82
83
84 double G0 = Wrho0 *
85 pow( ( ( Q0.
mass2() - m2_pip - m2_pim ) / 2 - M2rho0 / 4 ) /
86 ( M2rho0 / 4 - ( m2_pip + m2_pim ) / 2 ),
87 3 / 2 ) *
88 ( M2rho0 / Q0.
mass2() );
89 EvtComplex Drho0( ( Q0.
mass2() - M2rho0 ), Q0.
mass() * G0 );
90 EvtComplex A3( M2rho0 / Drho0 );
91
92
93 EvtComplex OmegaPhase( 0, -0.22 );
94 EvtComplex DOmega( ( Q0.
mass2() - M2omega ), Q0.
mass() * Womega );
95 EvtComplex A4( 0.0071 * M2omega *
exp( OmegaPhase ) / DOmega );
96
97
98 EvtComplex DirPhase( 0, 2.47 );
99 EvtComplex A5( 0.78 *
exp( DirPhase ) );
100
101 Atot = A1 + A2 + A3 + A4 + A5;
102
106
107 return;
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)