69 {
70
72
75
81
82
84 double theta =angles.getHelAng(1);
85 double phi =angles.getHelAng(2);
86 double gamma=0;
87 double F00,F01,F02,F03,F10,F11,F13,F20,F22,F23,F30,F31,F32,F33;
93
94 F33=F00; F32=F01; F31=F02; F30=F03; F22=F11;
95 F10=F01; F20=F02; F23=F10; F13=F20;
96
97
98
99
100 vertex(0,0,0,
Djmn(2, 0, 0,phi,theta,gamma)*F00);
101 vertex(0,0,1,
Djmn(2, 0,-1,phi,theta,gamma)*F01);
102 vertex(0,0,2,
Djmn(2, 0, 2,phi,theta,gamma)*F02);
103 vertex(0,0,3,
Djmn(2, 0, 1,phi,theta,gamma)*F03);
104 vertex(0,1,0,
Djmn(2, 0, 1,phi,theta,gamma)*F10);
105 vertex(0,1,1,
Djmn(2, 0, 0,phi,theta,gamma)*F11);
107 vertex(0,1,3,
Djmn(2, 0, 2,phi,theta,gamma)*F13);
108 vertex(0,2,0,
Djmn(2, 0,-2,phi,theta,gamma)*F20);
110 vertex(0,2,2,
Djmn(2, 0, 0,phi,theta,gamma)*F22);
111 vertex(0,2,3,
Djmn(2, 0,-1,phi,theta,gamma)*F23);
112 vertex(0,3,0,
Djmn(2, 0,-1,phi,theta,gamma)*F30);
113 vertex(0,3,1,
Djmn(2, 0,-2,phi,theta,gamma)*F31);
114 vertex(0,3,2,
Djmn(2, 0, 1,phi,theta,gamma)*F32);
115 vertex(0,3,3,
Djmn(2, 0, 0,phi,theta,gamma)*F33);
116
117 vertex(1,0,0,
Djmn(2, 2, 0,phi,theta,gamma)*F00);
118 vertex(1,0,1,
Djmn(2, 2,-1,phi,theta,gamma)*F01);
119 vertex(1,0,2,
Djmn(2, 2, 2,phi,theta,gamma)*F02);
120 vertex(1,0,3,
Djmn(2, 2, 1,phi,theta,gamma)*F03);
121 vertex(1,1,0,
Djmn(2, 2, 1,phi,theta,gamma)*F10);
122 vertex(1,1,1,
Djmn(2, 2, 0,phi,theta,gamma)*F11);
124 vertex(1,1,3,
Djmn(2, 2, 2,phi,theta,gamma)*F13);
125 vertex(1,2,0,
Djmn(2, 2,-2,phi,theta,gamma)*F20);
127 vertex(1,2,2,
Djmn(2, 2, 0,phi,theta,gamma)*F22);
128 vertex(1,2,3,
Djmn(2, 2,-1,phi,theta,gamma)*F23);
129 vertex(1,3,0,
Djmn(2, 2,-1,phi,theta,gamma)*F30);
130 vertex(1,3,1,
Djmn(2, 2,-2,phi,theta,gamma)*F31);
131 vertex(1,3,2,
Djmn(2, 2, 1,phi,theta,gamma)*F32);
132 vertex(1,3,3,
Djmn(2, 2, 0,phi,theta,gamma)*F33);
133
134 vertex(2,0,0,
Djmn(2,-2, 0,phi,theta,gamma)*F00);
135 vertex(2,0,1,
Djmn(2,-2,-1,phi,theta,gamma)*F01);
136 vertex(2,0,2,
Djmn(2,-2, 2,phi,theta,gamma)*F02);
137 vertex(2,0,3,
Djmn(2,-2, 1,phi,theta,gamma)*F03);
138 vertex(2,1,0,
Djmn(2,-2, 1,phi,theta,gamma)*F10);
139 vertex(2,1,1,
Djmn(2,-2, 0,phi,theta,gamma)*F11);
141 vertex(2,1,3,
Djmn(2,-2, 2,phi,theta,gamma)*F13);
142 vertex(2,2,0,
Djmn(2,-2,-2,phi,theta,gamma)*F20);
144 vertex(2,2,2,
Djmn(2,-2, 0,phi,theta,gamma)*F22);
145 vertex(2,2,3,
Djmn(2,-2,-1,phi,theta,gamma)*F23);
146 vertex(2,3,0,
Djmn(2,-2,-1,phi,theta,gamma)*F30);
147 vertex(2,3,1,
Djmn(2,-2,-2,phi,theta,gamma)*F31);
148 vertex(2,3,2,
Djmn(2,-2, 1,phi,theta,gamma)*F32);
149 vertex(2,3,3,
Djmn(2,-2, 0,phi,theta,gamma)*F33);
150
151 vertex(3,0,0,
Djmn(2, 1, 0,phi,theta,gamma)*F00);
152 vertex(3,0,1,
Djmn(2, 1,-1,phi,theta,gamma)*F01);
153 vertex(3,0,2,
Djmn(2, 1, 2,phi,theta,gamma)*F02);
154 vertex(3,0,3,
Djmn(2, 1, 1,phi,theta,gamma)*F03);
155 vertex(3,1,0,
Djmn(2, 1, 1,phi,theta,gamma)*F10);
156 vertex(3,1,1,
Djmn(2, 1, 0,phi,theta,gamma)*F11);
158 vertex(3,1,3,
Djmn(2, 1, 2,phi,theta,gamma)*F13);
159 vertex(3,2,0,
Djmn(2, 1,-2,phi,theta,gamma)*F20);
161 vertex(3,2,2,
Djmn(2, 1, 0,phi,theta,gamma)*F22);
162 vertex(3,2,3,
Djmn(2, 1,-1,phi,theta,gamma)*F23);
163 vertex(3,3,0,
Djmn(2, 1,-1,phi,theta,gamma)*F30);
164 vertex(3,3,1,
Djmn(2, 1,-2,phi,theta,gamma)*F31);
165 vertex(3,3,2,
Djmn(2, 1, 1,phi,theta,gamma)*F32);
166 vertex(3,3,3,
Djmn(2, 1, 0,phi,theta,gamma)*F33);
167
168 vertex(4,0,0,
Djmn(2,-1, 0,phi,theta,gamma)*F00);
169 vertex(4,0,1,
Djmn(2,-1,-1,phi,theta,gamma)*F01);
170 vertex(4,0,2,
Djmn(2,-1, 2,phi,theta,gamma)*F02);
171 vertex(4,0,3,
Djmn(2,-1, 1,phi,theta,gamma)*F03);
172 vertex(4,1,0,
Djmn(2,-1, 1,phi,theta,gamma)*F10);
173 vertex(4,1,1,
Djmn(2,-1, 0,phi,theta,gamma)*F11);
175 vertex(4,1,3,
Djmn(2,-1, 2,phi,theta,gamma)*F13);
176 vertex(4,2,0,
Djmn(2,-1,-2,phi,theta,gamma)*F20);
178 vertex(4,2,2,
Djmn(2,-1, 0,phi,theta,gamma)*F22);
179 vertex(4,2,3,
Djmn(2,-1,-1,phi,theta,gamma)*F23);
180 vertex(4,3,0,
Djmn(2,-1,-1,phi,theta,gamma)*F30);
181 vertex(4,3,1,
Djmn(2,-1,-2,phi,theta,gamma)*F31);
182 vertex(4,3,2,
Djmn(2,-1, 1,phi,theta,gamma)*F32);
183 vertex(4,3,3,
Djmn(2,-1, 0,phi,theta,gamma)*F33);
184
185 return ;
186
187}
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)