70 {
71
73
76
82
83
85 double theta =angles.getHelAng(1);
86 double phi =angles.getHelAng(2);
87 double gamma=0;
88
91 double F10=F01;
92 double F11=F00;
93
94 vertex(0,0,0,
Djmn(2, 0, 0,phi,theta,gamma)*F00);
95 vertex(0,0,1,
Djmn(2, 0, 1,phi,theta,gamma)*F01);
96 vertex(0,1,0,
Djmn(2, 0,-1,phi,theta,gamma)*F10);
97 vertex(0,1,1,
Djmn(2, 0, 0,phi,theta,gamma)*F11);
98
99 vertex(1,0,0,
Djmn(2, 2, 0,phi,theta,gamma)*F00);
100 vertex(1,0,1,
Djmn(2, 2, 1,phi,theta,gamma)*F01);
101 vertex(1,1,0,
Djmn(2, 2,-1,phi,theta,gamma)*F10);
102 vertex(1,1,1,
Djmn(2, 2, 0,phi,theta,gamma)*F11);
103
104 vertex(2,0,0,
Djmn(2,-2, 0,phi,theta,gamma)*F00);
105 vertex(2,0,1,
Djmn(2,-2, 1,phi,theta,gamma)*F01);
106 vertex(2,1,0,
Djmn(2,-2,-1,phi,theta,gamma)*F10);
107 vertex(2,1,1,
Djmn(2,-2, 0,phi,theta,gamma)*F11);
108
109 vertex(3,0,0,
Djmn(2, 1, 0,phi,theta,gamma)*F00);
110 vertex(3,0,1,
Djmn(2, 1, 1,phi,theta,gamma)*F01);
111 vertex(3,1,0,
Djmn(2, 1,-1,phi,theta,gamma)*F10);
112 vertex(3,1,1,
Djmn(2, 1, 0,phi,theta,gamma)*F11);
113
114 vertex(4,0,0,
Djmn(2,-1, 0,phi,theta,gamma)*F00);
115 vertex(4,0,1,
Djmn(2,-1, 1,phi,theta,gamma)*F01);
116 vertex(4,1,0,
Djmn(2,-1,-1,phi,theta,gamma)*F10);
117 vertex(4,1,1,
Djmn(2,-1, 0,phi,theta,gamma)*F11);
118
119 return ;
120
121}
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 &)
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)