73 {
74
76
79
85
86
88 double theta =angles.getHelAng(1);
89 double phi =angles.getHelAng(2);
90 double gamma =0.;
91
92 double A10=1.0,A11=2.08,A12=3.03;
93 double Af10=A10,Af1f1=A11,Af1f2=A12;
94 int j1,j2,j3;
95 for(j1=0;j1<=2;){
96 for(j2=0;j2<=1;){
97 for(j3=0;j3<=4;){
98 vertex(j1,j2,j3,0.0);j3++;}
99j2++;}
100j1++;}
101
102 vertex(0,0,1,
Djmn(1, 1,-1,phi,theta,gamma)*A12);
103 vertex(0,0,3,
Djmn(1, 1, 0,phi,theta,gamma)*A11);
104 vertex(0,0,0,
Djmn(1, 1, 1,phi,theta,gamma)*A10);
105 vertex(0,1,0,
Djmn(1, 1,-1,phi,theta,gamma)*Af10);
106 vertex(0,1,4,
Djmn(1, 1, 0,phi,theta,gamma)*Af1f1);
107 vertex(0,1,2,
Djmn(1, 1, 1,phi,theta,gamma)*Af1f2);
108
109 vertex(1,0,1,
Djmn(1,-1,-1,phi,theta,gamma)*A12);
110 vertex(1,0,3,
Djmn(1,-1, 0,phi,theta,gamma)*A11);
111 vertex(1,0,0,
Djmn(1,-1, 1,phi,theta,gamma)*A10);
112 vertex(1,1,0,
Djmn(1,-1,-1,phi,theta,gamma)*Af10);
113 vertex(1,1,4,
Djmn(1,-1, 0,phi,theta,gamma)*Af1f1);
114 vertex(1,1,2,
Djmn(1,-1, 1,phi,theta,gamma)*Af1f2);
115
116 vertex(2,0,1,
Djmn(1, 0,-1,phi,theta,gamma)*A12);
117 vertex(2,0,3,
Djmn(1, 0, 0,phi,theta,gamma)*A11);
118 vertex(2,0,0,
Djmn(1, 0, 1,phi,theta,gamma)*A10);
119 vertex(2,1,0,
Djmn(1, 0,-1,phi,theta,gamma)*Af10);
120 vertex(2,1,4,
Djmn(1, 0, 0,phi,theta,gamma)*Af1f1);
121 vertex(2,1,2,
Djmn(1, 0, 1,phi,theta,gamma)*Af1f2);
122
123 return ;
124
125}
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)