75 {
76
78
80 EvtVector4R pv,ps,ppr;
81
87 EvtVector4R pvcm=
v->getP4();
88
89
90 EvtHelSys angles(ppr,pv);
91 double theta =angles.getHelAng(1);
92 double phi =angles.getHelAng(2);
93 double gamma=0;
94
96 double Mpsip=fabs(mass_M-3.686);
97 double Mpsipp=fabs(mass_M-3.771);
98 double FP=pvcm.
d3mag();
99 if(Mpsip<0.1) {FP=FP*FP;} else{FP=1;}
100
101
102 double FM=-FP;
109
110
111 return ;
112
113}
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)