70 {
71
74
75
76
77
78
79
80
81
82
83 EvtComplex h0=EvtComplex(0.0,0.0);
84
85 EvtParticle *v1,*ph;
86
90 EvtVector4R momv1 = v1->
getP4();
91 EvtVector4R momph = ph->
getP4();
92
93 EvtTensor4C d,g;
94
96
98
100 vp=(momv1+momph)/(momv1+momph).
mass();
101
102 d=((1.0/sqrt(3.0))*(h0-(hp+hm))*(-1.0/sqrt(3.0)))*g+
103 ((1.0/sqrt(2.0))*(hp-hm)*EvtComplex(0.0,1.0)*(sqrt(1.0/2.0)))*
dual(
directProd(
v,vp))+
104 (sqrt(2.0/3.0)*(h0+0.5*(hp+hm))*sqrt(3.0/2.0))*(
directProd(
v,
v)+(1.0/3.0)*g);
105
106 EvtVector4C ep0,ep1,ep2;
107
111
112 EvtVector4C ep20,ep21,ep22;
113
116
119
122
125
126
127 return ;
128
129}
double sin(const BesAngle a)
double cos(const BesAngle a)
Evt3Rank3C directProd(const EvtVector3C &c1, const EvtVector3C &c2, const EvtVector3C &c3)
EvtTensor4C dual(const EvtTensor4C &t2)
**********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 &)
virtual EvtVector4C epsParentPhoton(int i)
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
virtual EvtVector4C eps(int i) const
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
void setdiag(double t00, double t11, double t22, double t33)
EvtVector4C cont1(const EvtVector4C &v4) const