66 {
67
69
71 EvtVector4R pv, ps, ppr;
72
78 EvtVector4R pvcm =
v->getP4();
79
80
81 EvtHelSys angles( ppr, pv );
82 double theta = angles.getHelAng( 1 );
83 double phi = angles.getHelAng( 2 );
84 double gamma = 0;
85
87 double Mpsip = fabs( mass_M - 3.686 );
88 double Mpsipp = fabs( mass_M - 3.771 );
89 double FP = pvcm.
d3mag();
90
91 if ( Mpsip < 0.1 ) { FP = FP * FP; }
92 else { FP = 1; }
93
94 double FM = -FP;
96 Djmn( 1, 1, 1, phi, theta, gamma ) * FP );
97 vertex( 0, 1,
Djmn( 1, 1, -1, phi, theta, gamma ) * FM );
98 vertex( 1, 0,
Djmn( 1, -1, 1, phi, theta, gamma ) * FP );
99
100 vertex( 1, 1,
Djmn( 1, -1, -1, phi, theta, gamma ) * FM );
102 Djmn( 1, 0, 1, phi, theta, gamma ) * FP );
103 vertex( 2, 1,
Djmn( 1, 0, -1, phi, theta, gamma ) * FM );
104
105 return;
106}
EvtComplex Djmn(int j, int m, int n, double phi, double theta, double gamma)
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyWtm INTEGER m_idyfs DOUBLE PRECISION m_zini DOUBLE PRECISION m_q2 DOUBLE PRECISION m_Wt_KF DOUBLE PRECISION m_WtCut INTEGER m_KFfin *COMMON c_KarLud $ !Input CMS energy[GeV] $ !CMS energy after beam spread beam strahlung[GeV] $ !Beam energy spread[GeV] $ !z boost due to beam spread $ !electron beam mass *ff pair spectrum $ !minimum v
void vertex(const EvtComplex &)
static double getMass(EvtId i)
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)