23 Double_t tPy = fpy1 + fpy2;
24 Double_t tPz = fpz1 + fpz2;
25 Double_t tE = fe1 + fe2;
26 Double_t tPt = tPx * tPx + tPy * tPy;
27 Double_t tMt = tE * tE - tPz * tPz;
28 Double_t tM = TMath::Sqrt(tMt - tPt);
29 tMt = TMath::Sqrt(tMt);
30 tPt = TMath::Sqrt(tPt);
37 Double_t tBeta = tPz / tE;
38 Double_t tGamma = tE / tMt;
39 fZ = tGamma * (fpz1 - tBeta * fe1);
40 Double_t tE1L = tGamma * (fe1 - tBeta * fpz1);
53 fX = (
fpx1 * tPx + fpy1 * tPy) / tPt;
54 fY = (-
fpx1 * tPy + fpy1 * tPx) / tPt;
62 fX = tMt / tM * (
fX - tPt / tMt * tE1L);
67 fT =
fX > 0. ? 1. : -1.;
73 Double_t tPy = fpy1 - fpy2;
74 Double_t tPz = fpz1 + fpz2;
75 Double_t tE = fe1 + fe2;
76 Double_t tPt = tPx * tPx + tPy * tPy;
77 Double_t tMt = tE * tE - tPz * tPz;
78 Double_t tM = TMath::Sqrt(tMt - tPt);
79 tMt = TMath::Sqrt(tMt);
80 tPt = TMath::Sqrt(tPt);
82 Double_t tBeta = tPz / tE;
83 Double_t tGamma = tE / tMt;
84 fZ = tGamma * (fpz1 - tBeta * fe1);
85 Double_t tE1L = tGamma * (fe1 - tBeta * fpz1);
88 fX = (
fpx1 * tPx + fpy1 * tPy) / tPt;
89 fY = (-
fpx1 * tPy + fpy1 * tPx) / tPt;
90 fX = tMt / tM * (
fX - tPt / tMt * tE1L);
91 fT =
fX > 0. ? 1. : -1.;
97 Double_t tPy = fpy1 - fpy2;
98 Double_t tPz = fpz1 - fpz2;
99 Double_t tE = fe1 + fe2;
100 Double_t tPt = tPx * tPx + tPy * tPy;
101 Double_t tMt = tE * tE - tPz * tPz;
102 Double_t tM = TMath::Sqrt(tMt - tPt);
103 tMt = TMath::Sqrt(tMt);
104 tPt = TMath::Sqrt(tPt);
106 Double_t tBeta = tPz / tE;
107 Double_t tGamma = tE / tMt;
108 fZ = tGamma * (fpz1 - tBeta * fe1);
109 Double_t tE1L = tGamma * (fe1 - tBeta * fpz1);
112 fX = (
fpx1 * tPx + fpy1 * tPy) / tPt;
113 fY = (-
fpx1 * tPy + fpy1 * tPx) / tPt;
114 fX = tMt / tM * (
fX - tPt / tMt * tE1L);
115 fT =
fX > 0. ? 1. : -1.;
121 Double_t tPy = fpy1 + fpy2;
122 Double_t tPz = fpz1 + fpz2;
123 Double_t tE = fe1 + fe2;
124 Double_t tPt = tPx * tPx + tPy * tPy;
125 Double_t tMt = tE * tE - tPz * tPz;
126 tMt = TMath::Sqrt(tMt);
127 tPt = TMath::Sqrt(tPt);
128 Double_t tBeta = tPz / tE;
129 Double_t tGamma = tE / tMt;
133 Double_t particle1lcms_pz = tGamma * (fpz1 - tBeta * fe1);
134 Double_t particle1lcms_e = tGamma * (fe1 - tBeta * fpz1);
135 Double_t particle2lcms_pz = tGamma * (fpz2 - tBeta * fe2);
136 Double_t particle2lcms_e = tGamma * (fe2 - tBeta * fpz2);
140 Double_t particle1lcms_px = (
fpx1 * tPx + fpy1 * tPy) / tPt;
141 Double_t particle1lcms_py = (-
fpx1 * tPy + fpy1 * tPx) / tPt;
143 Double_t particle2lcms_px = (
fpx2 * tPx + fpy2 * tPy) / tPt;
144 Double_t particle2lcms_py = (-
fpx2 * tPy + fpy2 * tPx) / tPt;
146 fX = particle1lcms_px - particle2lcms_px;
147 fY = particle1lcms_py - particle2lcms_py;
148 fZ = particle1lcms_pz - particle2lcms_pz;
149 Double_t mDE = particle1lcms_e - particle2lcms_e;
150 fT = TMath::Sqrt(TMath::Abs(
fX *
fX +
fY *
fY +
fZ *
fZ - mDE * mDE));
156 Double_t tPy = fpy1 - fpy2;
157 Double_t tPz = fpz1 + fpz2;
158 Double_t tE = fe1 + fe2;
159 Double_t tPt = tPx * tPx + tPy * tPy;
160 Double_t tMt = tE * tE - tPz * tPz;
161 tMt = TMath::Sqrt(tMt);
162 tPt = TMath::Sqrt(tPt);
163 Double_t tBeta = tPz / tE;
164 Double_t tGamma = tE / tMt;
168 Double_t particle1lcms_pz = tGamma * (fpz1 - tBeta * fe1);
169 Double_t particle1lcms_e = tGamma * (fe1 - tBeta * fpz1);
170 Double_t particle2lcms_pz = tGamma * (fpz2 - tBeta * fe2);
171 Double_t particle2lcms_e = tGamma * (fe2 - tBeta * fpz2);
175 Double_t particle1lcms_px = (
fpx1 * tPx + fpy1 * tPy) / tPt;
176 Double_t particle1lcms_py = (-
fpx1 * tPy + fpy1 * tPx) / tPt;
178 Double_t particle2lcms_px = (-
fpx2 * tPx - fpy2 * tPy) / tPt;
179 Double_t particle2lcms_py = (+
fpx2 * tPy - fpy2 * tPx) / tPt;
181 fX = particle1lcms_px - particle2lcms_px;
182 fY = particle1lcms_py - particle2lcms_py;
183 fZ = particle1lcms_pz - particle2lcms_pz;
185 Double_t mDE = particle1lcms_e - particle2lcms_e;
186 fT = TMath::Sqrt(TMath::Abs(
fX *
fX +
fY *
fY +
fZ *
fZ - mDE * mDE));
191 Double_t tPy = fpy1 - fpy2;
192 Double_t tPz = fpz1 - fpz2;
193 Double_t tE = fe1 + fe2;
194 Double_t tPt = tPx * tPx + tPy * tPy;
195 Double_t tMt = tE * tE - tPz * tPz;
196 tMt = TMath::Sqrt(tMt);
197 tPt = TMath::Sqrt(tPt);
198 Double_t tBeta = tPz / tE;
199 Double_t tGamma = tE / tMt;
203 Double_t particle1lcms_pz = tGamma * (fpz1 - tBeta * fe1);
204 Double_t particle1lcms_e = tGamma * (fe1 - tBeta * fpz1);
205 Double_t particle2lcms_pz = tGamma * (-fpz2 - tBeta * fe2);
206 Double_t particle2lcms_e = tGamma * (fe2 + tBeta * fpz2);
210 Double_t particle1lcms_px = (
fpx1 * tPx + fpy1 * tPy) / tPt;
211 Double_t particle1lcms_py = (-
fpx1 * tPy + fpy1 * tPx) / tPt;
213 Double_t particle2lcms_px = (-
fpx2 * tPx - fpy2 * tPy) / tPt;
214 Double_t particle2lcms_py = (+
fpx2 * tPy - fpy2 * tPx) / tPt;
217 fX = particle1lcms_px - particle2lcms_px;
218 fY = particle1lcms_py - particle2lcms_py;
219 fZ = particle1lcms_pz - particle2lcms_pz;
221 Double_t mDE = particle1lcms_e - particle2lcms_e;
222 fT = TMath::Sqrt(TMath::Abs(
fX *
fX +
fY *
fY +
fZ *
fZ - mDE * mDE));
227 Double_t tPy = fpy1 + fpy2;
228 Double_t tPz = fpz1 + fpz2;
229 Double_t tE = fe1 + fe2;
230 Double_t tPt = tPx * tPx + tPy * tPy;
231 Double_t tMt = tE * tE - tPz * tPz;
232 Double_t tM = TMath::Sqrt(tMt - tPt);
233 tMt = TMath::Sqrt(tMt);
234 tPt = TMath::Sqrt(tPt);
236 Double_t tBeta = tPz / tE;
237 Double_t tGamma = tE / tMt;
238 fZ = tGamma * (fpz1 - tBeta * fe1);
239 Double_t tE1L = tGamma * (fe1 - tBeta * fpz1);
241 fX = (
fpx1 * tPx + fpy1 * tPy) / tPt;
242 fY = (-
fpx1 * tPy + fpy1 * tPx) / tPt;
243 fX = tMt / tM * (
fX - tPt / tMt * tE1L);
251 Double_t tPy = fpy1 - fpy2;
252 Double_t tPz = fpz1 + fpz2;
253 Double_t tE = fe1 + fe2;
254 Double_t tPt = tPx * tPx + tPy * tPy;
255 Double_t tMt = tE * tE - tPz * tPz;
256 Double_t tM = TMath::Sqrt(tMt - tPt);
257 tMt = TMath::Sqrt(tMt);
258 tPt = TMath::Sqrt(tPt);
260 Double_t tBeta = tPz / tE;
261 Double_t tGamma = tE / tMt;
262 fZ = tGamma * (fpz1 - tBeta * fe1);
263 Double_t tE1L = tGamma * (fe1 - tBeta * fpz1);
266 fX = (
fpx1 * tPx + fpy1 * tPy) / tPt;
267 fY = (-
fpx1 * tPy + fpy1 * tPx) / tPt;
268 fX = tMt / tM * (
fX - tPt / tMt * tE1L);
269 fT =
fX > 0. ? 1. : -1.;
276 Double_t tPy = fpy1 - fpy2;
277 Double_t tPz = fpz1 - fpz2;
278 Double_t tE = fe1 + fe2;
279 Double_t tPt = tPx * tPx + tPy * tPy;
280 Double_t tMt = tE * tE - tPz * tPz;
281 Double_t tM = TMath::Sqrt(tMt - tPt);
282 tMt = TMath::Sqrt(tMt);
283 tPt = TMath::Sqrt(tPt);
285 Double_t tBeta = tPz / tE;
286 Double_t tGamma = tE / tMt;
287 fZ = tGamma * (fpz1 - tBeta * fe1);
288 Double_t tE1L = tGamma * (fe1 - tBeta * fpz1);
291 fX = (
fpx1 * tPx + fpy1 * tPy) / tPt;
292 fY = (-
fpx1 * tPy + fpy1 * tPx) / tPt;
293 fX = tMt / tM * (
fX - tPt / tMt * tE1L);
294 fT =
fX > 0. ? 1. : -1.;
301 Double_t tPy = fpy1 + fpy2;
302 Double_t tPz = fpz1 + fpz2;
303 Double_t tE = fe1 + fe2;
304 Double_t tPt = tPx * tPx + tPy * tPy;
305 Double_t tMt = tE * tE - tPz * tPz;
306 tMt = TMath::Sqrt(tMt);
307 tPt = TMath::Sqrt(tPt);
308 Double_t tBeta = tPz / tE;
309 Double_t tGamma = tE / tMt;
313 Double_t particle1lcms_pz = tGamma * (fpz1 - tBeta * fe1);
314 Double_t particle2lcms_pz = tGamma * (fpz2 - tBeta * fe2);
318 Double_t particle1lcms_px = (
fpx1 * tPx + fpy1 * tPy) / tPt;
319 Double_t particle1lcms_py = (-
fpx1 * tPy + fpy1 * tPx) / tPt;
321 Double_t particle2lcms_px = (
fpx2 * tPx + fpy2 * tPy) / tPt;
322 Double_t particle2lcms_py = (-
fpx2 * tPy + fpy2 * tPx) / tPt;
324 fX = particle1lcms_px - particle2lcms_px;
325 fY = particle1lcms_py - particle2lcms_py;
326 fZ = particle1lcms_pz - particle2lcms_pz;
334 Double_t tPy = fpy1 - fpy2;
335 Double_t tPz = fpz1 + fpz2;
336 Double_t tE = fe1 + fe2;
337 Double_t tPt = tPx * tPx + tPy * tPy;
338 Double_t tMt = tE * tE - tPz * tPz;
339 tMt = TMath::Sqrt(tMt);
340 tPt = TMath::Sqrt(tPt);
341 Double_t tBeta = tPz / tE;
342 Double_t tGamma = tE / tMt;
346 Double_t particle1lcms_pz = tGamma * (fpz1 - tBeta * fe1);
347 Double_t particle2lcms_pz = tGamma * (fpz2 - tBeta * fe2);
351 Double_t particle1lcms_px = (
fpx1 * tPx + fpy1 * tPy) / tPt;
352 Double_t particle1lcms_py = (-
fpx1 * tPy + fpy1 * tPx) / tPt;
354 Double_t particle2lcms_px = (-
fpx2 * tPx - fpy2 * tPy) / tPt;
355 Double_t particle2lcms_py = (+
fpx2 * tPy - fpy2 * tPx) / tPt;
358 fX = particle1lcms_px - particle2lcms_px;
359 fY = particle1lcms_py - particle2lcms_py;
360 fZ = particle1lcms_pz - particle2lcms_pz;
366 Double_t tPy = fpy1 - fpy2;
367 Double_t tPz = fpz1 - fpz2;
368 Double_t tE = fe1 + fe2;
369 Double_t tPt = tPx * tPx + tPy * tPy;
370 Double_t tMt = tE * tE - tPz * tPz;
371 tMt = TMath::Sqrt(tMt);
372 tPt = TMath::Sqrt(tPt);
373 Double_t tBeta = tPz / tE;
374 Double_t tGamma = tE / tMt;
378 Double_t particle1lcms_pz = tGamma * (fpz1 - tBeta * fe1);
379 Double_t particle2lcms_pz = tGamma * (-fpz2 - tBeta * fe2);
383 Double_t particle1lcms_px = (
fpx1 * tPx + fpy1 * tPy) / tPt;
384 Double_t particle1lcms_py = (-
fpx1 * tPy + fpy1 * tPx) / tPt;
386 Double_t particle2lcms_px = (-
fpx2 * tPx - fpy2 * tPy) / tPt;
387 Double_t particle2lcms_py = (+
fpx2 * tPy - fpy2 * tPx) / tPt;
390 fX = particle1lcms_px - particle2lcms_px;
391 fY = particle1lcms_py - particle2lcms_py;
392 fZ = particle1lcms_pz - particle2lcms_pz;