20 #include "TDirectory.h"
32 #include "CbmRichUtil.h"
36 #include "FairRootManager.h"
37 #include "KFParticle.h"
45 , fGlobalTracks(nullptr)
47 , fStsTrackMatches(nullptr)
48 , fRichProjections(nullptr)
50 , fRichRingMatches(nullptr)
52 , fArrayMvdHit(nullptr)
53 , fArrayStsHit(nullptr)
54 , fPrimVertex(nullptr)
57 , VStsTrack_minus_Outside()
58 , VMCtracks_minus_Outside()
59 , VRings_minus_Outside()
60 , VStsIndex_minus_Outside()
61 , VRichRing_minus_Outside()
62 , VMCIndex_minus_Outside()
63 , VStsTrack_plus_Outside()
64 , VMCtracks_plus_Outside()
65 , VRings_plus_Outside()
66 , VStsIndex_plus_Outside()
67 , VRichRing_plus_Outside()
68 , VMCIndex_plus_Outside()
69 , VMCtracks_minus_Target()
70 , VStsTrack_minus_Target()
71 , VMomenta_minus_Target()
72 , VRings_minus_Target()
73 , VStsIndex_minus_Target()
74 , VRichRing_minus_Target()
75 , VMCIndex_minus_Target()
76 , VMCtracks_plus_Target()
77 , VStsTrack_plus_Target()
78 , VMomenta_plus_Target()
79 , VRings_plus_Target()
80 , VStsIndex_plus_Target()
81 , VRichRing_plus_Target()
82 , VMCIndex_plus_Target()
87 , EMT_eta_gg_Event_Target()
88 , EMT_eta_gg_pair_momenta_Target()
89 , EMT_eta_gg_NofRings_Target()
90 , EMT_eta_gg_Event_Outside()
91 , EMT_eta_gg_pair_momenta_Outside()
92 , EMT_eta_gg_NofRings_Outside()
93 , EMT_eta_gg_Event_Both()
94 , EMT_eta_gg_pair_momenta_Both()
95 , EMT_eta_gg_NofRings_Both()
98 , Gammas_onetwo_Target()
99 , Gammas_stsIndex_all_Target()
100 , Gammas_stsIndex_two_Target()
101 , Gammas_stsIndex_onetwo_Target()
102 , Gammas_MCIndex_all_Target()
103 , Gammas_MCIndex_two_Target()
104 , Gammas_MCIndex_onetwo_Target()
105 , Gammas_MC_all_Target()
106 , Gammas_MC_two_Target()
107 , Gammas_MC_onetwo_Target()
108 , Gammas_all_Outside()
109 , Gammas_two_Outside()
110 , Gammas_onetwo_Outside()
111 , Gammas_stsIndex_all_Outside()
112 , Gammas_stsIndex_two_Outside()
113 , Gammas_stsIndex_onetwo_Outside()
114 , Gammas_MCIndex_all_Outside()
115 , Gammas_MCIndex_two_Outside()
116 , Gammas_MCIndex_onetwo_Outside()
117 , Gammas_MC_all_Outside()
118 , Gammas_MC_two_Outside()
119 , Gammas_MC_onetwo_Outside()
122 , Gammas_onetwo_Both()
123 , Gammas_stsIndex_all_Both()
124 , Gammas_stsIndex_two_Both()
125 , Gammas_stsIndex_onetwo_Both()
126 , Gammas_MCIndex_all_Both()
127 , Gammas_MCIndex_two_Both()
128 , Gammas_MCIndex_onetwo_Both()
129 , Gammas_MC_all_Both()
130 , Gammas_MC_two_Both()
131 , Gammas_MC_onetwo_Both()
133 , Particle_pull_px_Target(nullptr)
134 , Particle_pull_py_Target(nullptr)
135 , Particle_pull_pz_Target(nullptr)
136 , Particle_pull_px_Outside(nullptr)
137 , Particle_pull_py_Outside(nullptr)
138 , Particle_pull_pz_Outside(nullptr)
139 , Particle_pull_X_Outside(nullptr)
140 , Particle_pull_Y_Outside(nullptr)
141 , Particle_pull_Z_Outside(nullptr)
142 , fHistoList_Eta_cuts_Target()
143 , fHistoList_Eta_cuts_zeroInRich_Target()
144 , GammasInvMass_fromEta_beforeCuts_zeroInRich_Target(nullptr)
145 , GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Target(nullptr)
146 , GammasMomentum_fromEta_beforeCuts_zeroInRich_Target(nullptr)
147 , GammasRapidity_fromEta_beforeCuts_zeroInRich_Target(nullptr)
148 , PlaneAngles_last_fromEta_zeroInRich_Target(nullptr)
149 , PlaneAngles_first_fromEta_zeroInRich_Target(nullptr)
150 , GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Target(nullptr)
151 , GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Target(nullptr)
152 , GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Target(nullptr)
153 , GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Target(nullptr)
154 , PlaneAngles_last_wrongpairs_zeroInRich_Target(nullptr)
155 , PlaneAngles_first_wrongpairs_zeroInRich_Target(nullptr)
156 , fHistoList_Eta_cuts_oneInRich_Target()
157 , GammasInvMass_fromEta_beforeCuts_oneInRich_Target(nullptr)
158 , GammasOpenAngle_fromEta_beforeCuts_oneInRich_Target(nullptr)
159 , GammasMomentum_fromEta_beforeCuts_oneInRich_Target(nullptr)
160 , GammasRapidity_fromEta_beforeCuts_oneInRich_Target(nullptr)
161 , PlaneAngles_last_fromEta_oneInRich_Target(nullptr)
162 , PlaneAngles_first_fromEta_oneInRich_Target(nullptr)
163 , GammasInvMass_wrongpairs_beforeCuts_oneInRich_Target(nullptr)
164 , GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Target(nullptr)
165 , GammasMomentum_wrongpairs_beforeCuts_oneInRich_Target(nullptr)
166 , GammasRapidity_wrongpairs_beforeCuts_oneInRich_Target(nullptr)
167 , PlaneAngles_last_wrongpairs_oneInRich_Target(nullptr)
168 , PlaneAngles_first_wrongpairs_oneInRich_Target(nullptr)
169 , fHistoList_Eta_cuts_twoInRich_Target()
170 , GammasInvMass_fromEta_beforeCuts_twoInRich_Target(nullptr)
171 , GammasOpenAngle_fromEta_beforeCuts_twoInRich_Target(nullptr)
172 , GammasMomentum_fromEta_beforeCuts_twoInRich_Target(nullptr)
173 , GammasRapidity_fromEta_beforeCuts_twoInRich_Target(nullptr)
174 , PlaneAngles_last_fromEta_twoInRich_Target(nullptr)
175 , PlaneAngles_first_fromEta_twoInRich_Target(nullptr)
176 , GammasInvMass_wrongpairs_beforeCuts_twoInRich_Target(nullptr)
177 , GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Target(nullptr)
178 , GammasMomentum_wrongpairs_beforeCuts_twoInRich_Target(nullptr)
179 , GammasRapidity_wrongpairs_beforeCuts_twoInRich_Target(nullptr)
180 , PlaneAngles_last_wrongpairs_twoInRich_Target(nullptr)
181 , PlaneAngles_first_wrongpairs_twoInRich_Target(nullptr)
182 , fHistoList_Eta_cuts_Outside()
183 , fHistoList_Eta_cuts_zeroInRich_Outside()
184 , GammasInvMass_fromEta_beforeCuts_zeroInRich_Outside(nullptr)
185 , GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Outside(nullptr)
186 , GammasMomentum_fromEta_beforeCuts_zeroInRich_Outside(nullptr)
187 , GammasRapidity_fromEta_beforeCuts_zeroInRich_Outside(nullptr)
188 , PlaneAngles_last_fromEta_zeroInRich_Outside(nullptr)
189 , PlaneAngles_first_fromEta_zeroInRich_Outside(nullptr)
190 , GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Outside(nullptr)
191 , GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Outside(nullptr)
192 , GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Outside(nullptr)
193 , GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Outside(nullptr)
194 , PlaneAngles_last_wrongpairs_zeroInRich_Outside(nullptr)
195 , PlaneAngles_first_wrongpairs_zeroInRich_Outside(nullptr)
196 , fHistoList_Eta_cuts_oneInRich_Outside()
197 , GammasInvMass_fromEta_beforeCuts_oneInRich_Outside(nullptr)
198 , GammasOpenAngle_fromEta_beforeCuts_oneInRich_Outside(nullptr)
199 , GammasMomentum_fromEta_beforeCuts_oneInRich_Outside(nullptr)
200 , GammasRapidity_fromEta_beforeCuts_oneInRich_Outside(nullptr)
201 , PlaneAngles_last_fromEta_oneInRich_Outside(nullptr)
202 , PlaneAngles_first_fromEta_oneInRich_Outside(nullptr)
203 , GammasInvMass_wrongpairs_beforeCuts_oneInRich_Outside(nullptr)
204 , GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Outside(nullptr)
205 , GammasMomentum_wrongpairs_beforeCuts_oneInRich_Outside(nullptr)
206 , GammasRapidity_wrongpairs_beforeCuts_oneInRich_Outside(nullptr)
207 , PlaneAngles_last_wrongpairs_oneInRich_Outside(nullptr)
208 , PlaneAngles_first_wrongpairs_oneInRich_Outside(nullptr)
209 , fHistoList_Eta_cuts_twoInRich_Outside()
210 , GammasInvMass_fromEta_beforeCuts_twoInRich_Outside(nullptr)
211 , GammasOpenAngle_fromEta_beforeCuts_twoInRich_Outside(nullptr)
212 , GammasMomentum_fromEta_beforeCuts_twoInRich_Outside(nullptr)
213 , GammasRapidity_fromEta_beforeCuts_twoInRich_Outside(nullptr)
214 , PlaneAngles_last_fromEta_twoInRich_Outside(nullptr)
215 , PlaneAngles_first_fromEta_twoInRich_Outside(nullptr)
216 , GammasInvMass_wrongpairs_beforeCuts_twoInRich_Outside(nullptr)
217 , GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Outside(nullptr)
218 , GammasMomentum_wrongpairs_beforeCuts_twoInRich_Outside(nullptr)
219 , GammasRapidity_wrongpairs_beforeCuts_twoInRich_Outside(nullptr)
220 , PlaneAngles_last_wrongpairs_twoInRich_Outside(nullptr)
221 , PlaneAngles_first_wrongpairs_twoInRich_Outside(nullptr)
222 , fHistoList_Eta_cuts_Both()
223 , fHistoList_Eta_cuts_zeroInRich_Both()
224 , GammasInvMass_fromEta_beforeCuts_zeroInRich_Both(nullptr)
225 , GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Both(nullptr)
226 , GammasMomentum_fromEta_beforeCuts_zeroInRich_Both(nullptr)
227 , GammasRapidity_fromEta_beforeCuts_zeroInRich_Both(nullptr)
228 , PlaneAngles_last_fromEta_zeroInRich_Both(nullptr)
229 , PlaneAngles_first_fromEta_zeroInRich_Both(nullptr)
230 , GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Both(nullptr)
231 , GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Both(nullptr)
232 , GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Both(nullptr)
233 , GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Both(nullptr)
234 , PlaneAngles_last_wrongpairs_zeroInRich_Both(nullptr)
235 , PlaneAngles_first_wrongpairs_zeroInRich_Both(nullptr)
236 , fHistoList_Eta_cuts_oneInRich_Both()
237 , GammasInvMass_fromEta_beforeCuts_oneInRich_Both(nullptr)
238 , GammasOpenAngle_fromEta_beforeCuts_oneInRich_Both(nullptr)
239 , GammasMomentum_fromEta_beforeCuts_oneInRich_Both(nullptr)
240 , GammasRapidity_fromEta_beforeCuts_oneInRich_Both(nullptr)
241 , PlaneAngles_last_fromEta_oneInRich_Both(nullptr)
242 , PlaneAngles_first_fromEta_oneInRich_Both(nullptr)
243 , GammasInvMass_wrongpairs_beforeCuts_oneInRich_Both(nullptr)
244 , GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Both(nullptr)
245 , GammasMomentum_wrongpairs_beforeCuts_oneInRich_Both(nullptr)
246 , GammasRapidity_wrongpairs_beforeCuts_oneInRich_Both(nullptr)
247 , PlaneAngles_last_wrongpairs_oneInRich_Both(nullptr)
248 , PlaneAngles_first_wrongpairs_oneInRich_Both(nullptr)
249 , PlaneAngles_last_fromEta_afterCuts_oneInRich_Both(nullptr)
250 , PlaneAngles_first_fromEta_afterCuts_oneInRich_Both(nullptr)
251 , PlaneAngles_last_wrongpairs_afterCuts_oneInRich_Both(nullptr)
252 , PlaneAngles_first_wrongpairs_afterCuts_oneInRich_Both(nullptr)
253 , fHistoList_Eta_cuts_twoInRich_Both()
254 , GammasInvMass_fromEta_beforeCuts_twoInRich_Both(nullptr)
255 , GammasOpenAngle_fromEta_beforeCuts_twoInRich_Both(nullptr)
256 , GammasMomentum_fromEta_beforeCuts_twoInRich_Both(nullptr)
257 , GammasRapidity_fromEta_beforeCuts_twoInRich_Both(nullptr)
258 , PlaneAngles_last_fromEta_twoInRich_Both(nullptr)
259 , PlaneAngles_first_fromEta_twoInRich_Both(nullptr)
260 , GammasInvMass_wrongpairs_beforeCuts_twoInRich_Both(nullptr)
261 , GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Both(nullptr)
262 , GammasMomentum_wrongpairs_beforeCuts_twoInRich_Both(nullptr)
263 , GammasRapidity_wrongpairs_beforeCuts_twoInRich_Both(nullptr)
264 , PlaneAngles_last_wrongpairs_twoInRich_Both(nullptr)
265 , PlaneAngles_first_wrongpairs_twoInRich_Both(nullptr)
266 , fHistoList_Eta_all_Target()
267 , AngleBetweenGammas_trueEta_before_cuts_all_Target(nullptr)
268 , AngleBetweenGammas_wrongEta_before_cuts_all_Target(nullptr)
269 , InvMass_GammaCombinations_all_Target(nullptr)
270 , Particles_PDG_all_Target(nullptr)
271 , TrueEta_pt_vs_rap_all_Target(nullptr)
272 , WrongEta_pt_vs_rap_all_Target(nullptr)
273 , TrueEta_pt_vs_rap_est_all_Target(nullptr)
274 , WrongEta_pt_vs_rap_est_all_Target(nullptr)
275 , TrueEta_InvMass_after_cuts_all_Target(nullptr)
276 , WrongEta_InvMass_after_cuts_all_Target(nullptr)
277 , TrueEta_AngleBetweenGammas_after_cuts_all_Target(nullptr)
278 , WrongEta_AngleBetweenGammas_after_cuts_all_Target(nullptr)
279 , EMT_Eta_InvMass_all_Target(nullptr)
280 , fHistoList_Eta_onetwo_Target()
281 , AngleBetweenGammas_trueEta_before_cuts_onetwo_Target(nullptr)
282 , AngleBetweenGammas_wrongEta_before_cuts_onetwo_Target(nullptr)
283 , InvMass_GammaCombinations_onetwo_Target(nullptr)
284 , Particles_PDG_onetwo_Target(nullptr)
285 , TrueEta_pt_vs_rap_onetwo_Target(nullptr)
286 , WrongEta_pt_vs_rap_onetwo_Target(nullptr)
287 , TrueEta_pt_vs_rap_est_onetwo_Target(nullptr)
288 , WrongEta_pt_vs_rap_est_onetwo_Target(nullptr)
289 , TrueEta_InvMass_after_cuts_onetwo_Target(nullptr)
290 , WrongEta_InvMass_after_cuts_onetwo_Target(nullptr)
291 , TrueEta_AngleBetweenGammas_after_cuts_onetwo_Target(nullptr)
292 , WrongEta_AngleBetweenGammas_after_cuts_onetwo_Target(nullptr)
293 , EMT_Eta_InvMass_onetwo_Target(nullptr)
294 , fHistoList_Eta_two_Target()
295 , AngleBetweenGammas_trueEta_before_cuts_two_Target(nullptr)
296 , AngleBetweenGammas_wrongEta_before_cuts_two_Target(nullptr)
297 , InvMass_GammaCombinations_two_Target(nullptr)
298 , Particles_PDG_two_Target(nullptr)
299 , TrueEta_pt_vs_rap_two_Target(nullptr)
300 , WrongEta_pt_vs_rap_two_Target(nullptr)
301 , TrueEta_pt_vs_rap_est_two_Target(nullptr)
302 , WrongEta_pt_vs_rap_est_two_Target(nullptr)
303 , TrueEta_InvMass_after_cuts_two_Target(nullptr)
304 , WrongEta_InvMass_after_cuts_two_Target(nullptr)
305 , TrueEta_AngleBetweenGammas_after_cuts_two_Target(nullptr)
306 , WrongEta_AngleBetweenGammas_after_cuts_two_Target(nullptr)
307 , EMT_Eta_InvMass_two_Target(nullptr)
308 , fHistoList_rap_pt_Eta_all_Target()
309 , fHistoList_rap_pt_Eta_onetwo_Target()
310 , fHistoList_rap_pt_Eta_two_Target()
311 , fHistoList_Eta_all_Outside()
312 , AngleBetweenGammas_trueEta_before_cuts_all_Outside(nullptr)
313 , AngleBetweenGammas_wrongEta_before_cuts_all_Outside(nullptr)
314 , InvMass_GammaCombinations_all_Outside(nullptr)
315 , Particles_PDG_all_Outside(nullptr)
316 , TrueEta_pt_vs_rap_all_Outside(nullptr)
317 , WrongEta_pt_vs_rap_all_Outside(nullptr)
318 , TrueEta_pt_vs_rap_est_all_Outside(nullptr)
319 , WrongEta_pt_vs_rap_est_all_Outside(nullptr)
320 , TrueEta_InvMass_after_cuts_all_Outside(nullptr)
321 , WrongEta_InvMass_after_cuts_all_Outside(nullptr)
322 , TrueEta_AngleBetweenGammas_after_cuts_all_Outside(nullptr)
323 , WrongEta_AngleBetweenGammas_after_cuts_all_Outside(nullptr)
324 , EMT_Eta_InvMass_all_Outside(nullptr)
325 , fHistoList_Eta_onetwo_Outside()
326 , AngleBetweenGammas_trueEta_before_cuts_onetwo_Outside(nullptr)
327 , AngleBetweenGammas_wrongEta_before_cuts_onetwo_Outside(nullptr)
328 , InvMass_GammaCombinations_onetwo_Outside(nullptr)
329 , Particles_PDG_onetwo_Outside(nullptr)
330 , TrueEta_pt_vs_rap_onetwo_Outside(nullptr)
331 , WrongEta_pt_vs_rap_onetwo_Outside(nullptr)
332 , TrueEta_pt_vs_rap_est_onetwo_Outside(nullptr)
333 , WrongEta_pt_vs_rap_est_onetwo_Outside(nullptr)
334 , TrueEta_InvMass_after_cuts_onetwo_Outside(nullptr)
335 , WrongEta_InvMass_after_cuts_onetwo_Outside(nullptr)
336 , TrueEta_AngleBetweenGammas_after_cuts_onetwo_Outside(nullptr)
337 , WrongEta_AngleBetweenGammas_after_cuts_onetwo_Outside(nullptr)
338 , EMT_Eta_InvMass_onetwo_Outside(nullptr)
339 , fHistoList_Eta_two_Outside()
340 , AngleBetweenGammas_trueEta_before_cuts_two_Outside(nullptr)
341 , AngleBetweenGammas_wrongEta_before_cuts_two_Outside(nullptr)
342 , InvMass_GammaCombinations_two_Outside(nullptr)
343 , Particles_PDG_two_Outside(nullptr)
344 , TrueEta_pt_vs_rap_two_Outside(nullptr)
345 , WrongEta_pt_vs_rap_two_Outside(nullptr)
346 , TrueEta_pt_vs_rap_est_two_Outside(nullptr)
347 , WrongEta_pt_vs_rap_est_two_Outside(nullptr)
348 , TrueEta_InvMass_after_cuts_two_Outside(nullptr)
349 , WrongEta_InvMass_after_cuts_two_Outside(nullptr)
350 , TrueEta_AngleBetweenGammas_after_cuts_two_Outside(nullptr)
351 , WrongEta_AngleBetweenGammas_after_cuts_two_Outside(nullptr)
352 , EMT_Eta_InvMass_two_Outside(nullptr)
353 , fHistoList_rap_pt_Eta_all_Outside()
354 , fHistoList_rap_pt_Eta_onetwo_Outside()
355 , fHistoList_rap_pt_Eta_two_Outside()
356 , fHistoList_Eta_all_Both()
357 , AngleBetweenGammas_trueEta_before_cuts_all_Both(nullptr)
358 , AngleBetweenGammas_wrongEta_before_cuts_all_Both(nullptr)
359 , InvMass_GammaCombinations_all_Both(nullptr)
360 , Particles_PDG_all_Both(nullptr)
361 , TrueEta_pt_vs_rap_all_Both(nullptr)
362 , WrongEta_pt_vs_rap_all_Both(nullptr)
363 , TrueEta_pt_vs_rap_est_all_Both(nullptr)
364 , WrongEta_pt_vs_rap_est_all_Both(nullptr)
365 , TrueEta_InvMass_after_cuts_all_Both(nullptr)
366 , WrongEta_InvMass_after_cuts_all_Both(nullptr)
367 , TrueEta_AngleBetweenGammas_after_cuts_all_Both(nullptr)
368 , WrongEta_AngleBetweenGammas_after_cuts_all_Both(nullptr)
369 , EMT_Eta_InvMass_all_Both(nullptr)
370 , fHistoList_rap_pt_Eta_all_Both()
371 , multi_InvMass_Eta_all_Both_1(nullptr)
372 , multi_InvMass_Eta_all_Both_2(nullptr)
373 , multi_InvMass_Eta_all_Both_3(nullptr)
374 , multi_InvMass_Eta_all_Both_4(nullptr)
375 , multi_InvMass_Eta_all_Both_5(nullptr)
376 , multi_InvMass_Eta_all_Both_6(nullptr)
377 , multi_InvMass_Eta_all_Both_7(nullptr)
378 , multi_InvMass_Eta_all_Both_8(nullptr)
379 , multi_InvMass_Eta_all_Both_9(nullptr)
380 , multi_InvMass_Eta_all_Both_10(nullptr)
381 , multi_InvMass_Eta_all_Both_11(nullptr)
382 , multi_InvMass_Eta_all_Both_12(nullptr)
383 , multi_InvMass_Eta_all_Both_13(nullptr)
384 , multi_InvMass_Eta_all_Both_14(nullptr)
385 , multi_InvMass_Eta_all_Both_15(nullptr)
386 , multi_InvMass_Eta_all_Both_16(nullptr)
387 , multi_InvMass_Eta_all_Both_17(nullptr)
388 , multi_EMT_Eta_all_Both_1(nullptr)
389 , multi_EMT_Eta_all_Both_2(nullptr)
390 , multi_EMT_Eta_all_Both_3(nullptr)
391 , multi_EMT_Eta_all_Both_4(nullptr)
392 , multi_EMT_Eta_all_Both_5(nullptr)
393 , multi_EMT_Eta_all_Both_6(nullptr)
394 , multi_EMT_Eta_all_Both_7(nullptr)
395 , multi_EMT_Eta_all_Both_8(nullptr)
396 , multi_EMT_Eta_all_Both_9(nullptr)
397 , multi_EMT_Eta_all_Both_10(nullptr)
398 , multi_EMT_Eta_all_Both_11(nullptr)
399 , multi_EMT_Eta_all_Both_12(nullptr)
400 , multi_EMT_Eta_all_Both_13(nullptr)
401 , multi_EMT_Eta_all_Both_14(nullptr)
402 , multi_EMT_Eta_all_Both_15(nullptr)
403 , multi_EMT_Eta_all_Both_16(nullptr)
404 , multi_EMT_Eta_all_Both_17(nullptr)
405 , fHistoList_Eta_onetwo_Both()
406 , AngleBetweenGammas_trueEta_before_cuts_onetwo_Both(nullptr)
407 , AngleBetweenGammas_wrongEta_before_cuts_onetwo_Both(nullptr)
408 , InvMass_GammaCombinations_onetwo_Both(nullptr)
409 , Particles_PDG_onetwo_Both(nullptr)
410 , TrueEta_pt_vs_rap_onetwo_Both(nullptr)
411 , WrongEta_pt_vs_rap_onetwo_Both(nullptr)
412 , TrueEta_pt_vs_rap_est_onetwo_Both(nullptr)
413 , WrongEta_pt_vs_rap_est_onetwo_Both(nullptr)
414 , TrueEta_InvMass_after_cuts_onetwo_Both(nullptr)
415 , WrongEta_InvMass_after_cuts_onetwo_Both(nullptr)
416 , TrueEta_AngleBetweenGammas_after_cuts_onetwo_Both(nullptr)
417 , WrongEta_AngleBetweenGammas_after_cuts_onetwo_Both(nullptr)
418 , EMT_Eta_InvMass_onetwo_Both(nullptr)
419 , fHistoList_rap_pt_Eta_onetwo_Both()
420 , multi_InvMass_Eta_onetwo_Both_1(nullptr)
421 , multi_InvMass_Eta_onetwo_Both_2(nullptr)
422 , multi_InvMass_Eta_onetwo_Both_3(nullptr)
423 , multi_InvMass_Eta_onetwo_Both_4(nullptr)
424 , multi_InvMass_Eta_onetwo_Both_5(nullptr)
425 , multi_InvMass_Eta_onetwo_Both_6(nullptr)
426 , multi_InvMass_Eta_onetwo_Both_7(nullptr)
427 , multi_InvMass_Eta_onetwo_Both_8(nullptr)
428 , multi_InvMass_Eta_onetwo_Both_9(nullptr)
429 , multi_InvMass_Eta_onetwo_Both_10(nullptr)
430 , multi_InvMass_Eta_onetwo_Both_11(nullptr)
431 , multi_InvMass_Eta_onetwo_Both_12(nullptr)
432 , multi_InvMass_Eta_onetwo_Both_13(nullptr)
433 , multi_InvMass_Eta_onetwo_Both_14(nullptr)
434 , multi_InvMass_Eta_onetwo_Both_15(nullptr)
435 , multi_InvMass_Eta_onetwo_Both_16(nullptr)
436 , multi_InvMass_Eta_onetwo_Both_17(nullptr)
437 , multi_EMT_Eta_onetwo_Both_1(nullptr)
438 , multi_EMT_Eta_onetwo_Both_2(nullptr)
439 , multi_EMT_Eta_onetwo_Both_3(nullptr)
440 , multi_EMT_Eta_onetwo_Both_4(nullptr)
441 , multi_EMT_Eta_onetwo_Both_5(nullptr)
442 , multi_EMT_Eta_onetwo_Both_6(nullptr)
443 , multi_EMT_Eta_onetwo_Both_7(nullptr)
444 , multi_EMT_Eta_onetwo_Both_8(nullptr)
445 , multi_EMT_Eta_onetwo_Both_9(nullptr)
446 , multi_EMT_Eta_onetwo_Both_10(nullptr)
447 , multi_EMT_Eta_onetwo_Both_11(nullptr)
448 , multi_EMT_Eta_onetwo_Both_12(nullptr)
449 , multi_EMT_Eta_onetwo_Both_13(nullptr)
450 , multi_EMT_Eta_onetwo_Both_14(nullptr)
451 , multi_EMT_Eta_onetwo_Both_15(nullptr)
452 , multi_EMT_Eta_onetwo_Both_16(nullptr)
453 , multi_EMT_Eta_onetwo_Both_17(nullptr)
454 , fHistoList_Eta_two_Both()
455 , AngleBetweenGammas_trueEta_before_cuts_two_Both(nullptr)
456 , AngleBetweenGammas_wrongEta_before_cuts_two_Both(nullptr)
457 , InvMass_GammaCombinations_two_Both(nullptr)
458 , Particles_PDG_two_Both(nullptr)
459 , TrueEta_pt_vs_rap_two_Both(nullptr)
460 , WrongEta_pt_vs_rap_two_Both(nullptr)
461 , TrueEta_pt_vs_rap_est_two_Both(nullptr)
462 , WrongEta_pt_vs_rap_est_two_Both(nullptr)
463 , TrueEta_InvMass_after_cuts_two_Both(nullptr)
464 , WrongEta_InvMass_after_cuts_two_Both(nullptr)
465 , TrueEta_AngleBetweenGammas_after_cuts_two_Both(nullptr)
466 , WrongEta_AngleBetweenGammas_after_cuts_two_Both(nullptr)
467 , EMT_Eta_InvMass_two_Both(nullptr)
468 , fHistoList_rap_pt_Eta_two_Both()
469 , multi_InvMass_Eta_two_Both_1(nullptr)
470 , multi_InvMass_Eta_two_Both_2(nullptr)
471 , multi_InvMass_Eta_two_Both_3(nullptr)
472 , multi_InvMass_Eta_two_Both_4(nullptr)
473 , multi_InvMass_Eta_two_Both_5(nullptr)
474 , multi_InvMass_Eta_two_Both_6(nullptr)
475 , multi_InvMass_Eta_two_Both_7(nullptr)
476 , multi_InvMass_Eta_two_Both_8(nullptr)
477 , multi_InvMass_Eta_two_Both_9(nullptr)
478 , multi_InvMass_Eta_two_Both_10(nullptr)
479 , multi_InvMass_Eta_two_Both_11(nullptr)
480 , multi_InvMass_Eta_two_Both_12(nullptr)
481 , multi_InvMass_Eta_two_Both_13(nullptr)
482 , multi_InvMass_Eta_two_Both_14(nullptr)
483 , multi_InvMass_Eta_two_Both_15(nullptr)
484 , multi_InvMass_Eta_two_Both_16(nullptr)
485 , multi_InvMass_Eta_two_Both_17(nullptr)
486 , multi_EMT_Eta_two_Both_1(nullptr)
487 , multi_EMT_Eta_two_Both_2(nullptr)
488 , multi_EMT_Eta_two_Both_3(nullptr)
489 , multi_EMT_Eta_two_Both_4(nullptr)
490 , multi_EMT_Eta_two_Both_5(nullptr)
491 , multi_EMT_Eta_two_Both_6(nullptr)
492 , multi_EMT_Eta_two_Both_7(nullptr)
493 , multi_EMT_Eta_two_Both_8(nullptr)
494 , multi_EMT_Eta_two_Both_9(nullptr)
495 , multi_EMT_Eta_two_Both_10(nullptr)
496 , multi_EMT_Eta_two_Both_11(nullptr)
497 , multi_EMT_Eta_two_Both_12(nullptr)
498 , multi_EMT_Eta_two_Both_13(nullptr)
499 , multi_EMT_Eta_two_Both_14(nullptr)
500 , multi_EMT_Eta_two_Both_15(nullptr)
501 , multi_EMT_Eta_two_Both_16(nullptr)
502 , multi_EMT_Eta_two_Both_17(nullptr) {}
507 FairRootManager* ioman = FairRootManager::Instance();
508 if (
nullptr == ioman) {
509 Fatal(
"CbmKresEta::Init",
"RootManager not instantised!");
512 fMcTracks = (TClonesArray*) ioman->GetObject(
"MCTrack");
513 if (
nullptr ==
fMcTracks) { Fatal(
"CbmKresEta::Init",
"No MCTrack array!"); }
520 LOG(fatal) <<
"CbmKresEta::Init No PrimaryVertex array!";
523 fGlobalTracks = (TClonesArray*) ioman->GetObject(
"GlobalTrack");
525 Fatal(
"CbmKresEta::Init",
"No GlobalTrack array!");
528 fStsTracks = (TClonesArray*) ioman->GetObject(
"StsTrack");
530 Fatal(
"CbmKresEta::Init",
"No StsTrack array!");
535 Fatal(
"CbmKresEta::Init",
"No StsTrackMatch array!");
540 Fatal(
"CbmKresEta::Init",
"No RichProjection array!");
543 fRichRings = (TClonesArray*) ioman->GetObject(
"RichRing");
545 Fatal(
"CbmKresEta::Init",
"No RichRing array!");
550 Fatal(
"CbmKresEta::Init",
"No RichRingMatch array!");
553 fRichHits = (TClonesArray*) ioman->GetObject(
"RichHit");
554 if (
nullptr ==
fRichHits) { Fatal(
"CbmKresEta::Init",
"No RichHit array!"); }
556 fArrayMvdHit = (TClonesArray*) ioman->GetObject(
"MvdHit");
558 Fatal(
"CbmKresEta::Init",
"No MvdHit array!");
561 fArrayStsHit = (TClonesArray*) ioman->GetObject(
"StsHit");
563 Fatal(
"CbmKresEta::Init",
"No StsHit array!");
572 double OpeningAngleCut,
573 double GammaInvMassCut,
580 Fatal(
"CbmKresConversionManual::Exec",
"No PrimaryVertex array!");
649 for (Int_t
i = 0;
i < ngTracks;
i++) {
651 if (
nullptr == gTrack)
continue;
655 if (stsInd < 0)
continue;
657 if (stsTrack ==
nullptr)
continue;
660 if (stsMatch ==
nullptr)
continue;
663 if (stsMcTrackId < 0)
continue;
665 if (mcTrack ==
nullptr)
continue;
668 if (richInd < 0 && proj->GetX() > -115 && proj->GetX() < 115
669 && ((proj->GetY() < -120 && proj->GetY() > -200)
670 || (proj->GetY() > 120 && proj->GetY() < 200)))
686 double charge = track_par->GetQp();
688 if (chi2 != chi2)
continue;
689 if (chi2 == 0)
continue;
693 mcTrack, stsTrack, charge, stsInd, richInd, stsMcTrackId, Ring);
695 if (chi2 > 3)
continue;
698 mcTrack, stsTrack, Momentum, charge, stsInd, richInd, stsMcTrackId, Ring);
815 if (fEventNumEta % numformix == 0) {
842 if (fEventNumEta % numformix == 0) {
849 if (fEventNumEta % numformix == 0) {
865 int InRich =
FindInRich(richInd, stsMcTrackId);
893 int InRich =
FindInRich(richInd, stsMcTrackId);
919 if (richMatch !=
nullptr && richMatch->
GetNofLinks() > 0) {
921 if (richMcTrackId > 0) {
926 if (mcTrack2 !=
nullptr) {
928 if (TMath::Abs(pdgRICH) == 11) RingsInRich++;
940 if (
nullptr != ring) {
943 for (
int i = 0;
i < nofHits;
i++) {
946 if (
nullptr == hit)
continue;
962 double FinalAngle = 400;
977 Xpart1 = stsHit1->
GetX();
978 Ypart1 = stsHit1->
GetY();
979 Zpart1 = stsHit1->
GetZ();
983 Xpart1 = mvdHit1->
GetX();
984 Ypart1 = mvdHit1->
GetY();
985 Zpart1 = mvdHit1->
GetZ();
991 Xpart2 = stsHit2->
GetX();
992 Ypart2 = stsHit2->
GetY();
993 Zpart2 = stsHit2->
GetZ();
997 Xpart2 = mvdHit2->
GetX();
998 Ypart2 = mvdHit2->
GetY();
999 Zpart2 = mvdHit2->
GetZ();
1004 TMath::Abs(Zpart1 - Zpart2) > 2
1007 for (
int i = hits1sts - 2;
i > -1;
i--) {
1008 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
1011 Xpart1 = stsHit->
GetX();
1012 Ypart1 = stsHit->
GetY();
1013 Zpart1 = stsHit->
GetZ();
1015 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 > Zpart2) {
1016 for (
int i = hits1mvd - 2;
i > -1;
i--) {
1017 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
1020 Xpart1 = mvdHit->
GetX();
1021 Ypart1 = mvdHit->
GetY();
1022 Zpart1 = mvdHit->
GetZ();
1028 TMath::Abs(Zpart1 - Zpart2) > 2
1031 for (
int i = hits2sts - 2;
i > -1;
i--) {
1032 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
1035 Xpart2 = stsHit->
GetX();
1036 Ypart2 = stsHit->
GetY();
1037 Zpart2 = stsHit->
GetZ();
1039 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 < Zpart2) {
1040 for (
int i = hits2mvd - 2;
i > -1;
i--) {
1041 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
1044 Xpart2 = mvdHit->
GetX();
1045 Ypart2 = mvdHit->
GetY();
1046 Zpart2 = mvdHit->
GetZ();
1052 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
1054 TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 /
TMath::Pi());
1057 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
1062 double FinalAngle = 400;
1077 Xpart1 = mvdHit1->
GetX();
1078 Ypart1 = mvdHit1->
GetY();
1079 Zpart1 = mvdHit1->
GetZ();
1083 Xpart1 = stsHit1->
GetX();
1084 Ypart1 = stsHit1->
GetY();
1085 Zpart1 = stsHit1->
GetZ();
1091 Xpart2 = mvdHit2->
GetX();
1092 Ypart2 = mvdHit2->
GetY();
1093 Zpart2 = mvdHit2->
GetZ();
1097 Xpart2 = stsHit2->
GetX();
1098 Ypart2 = stsHit2->
GetY();
1099 Zpart2 = stsHit2->
GetZ();
1104 TMath::Abs(Zpart1 - Zpart2) > 2
1107 for (
int i = 1;
i < hits1mvd;
i++) {
1108 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
1111 Xpart1 = mvdHit->
GetX();
1112 Ypart1 = mvdHit->
GetY();
1113 Zpart1 = mvdHit->
GetZ();
1115 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 < Zpart2) {
1116 for (
int i = 0;
i < hits1sts;
i++) {
1117 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
1120 Xpart1 = stsHit->
GetX();
1121 Ypart1 = stsHit->
GetY();
1122 Zpart1 = stsHit->
GetZ();
1128 TMath::Abs(Zpart1 - Zpart2) > 2
1131 for (
int i = 1;
i < hits2mvd;
i++) {
1132 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
1135 Xpart2 = mvdHit->
GetX();
1136 Ypart2 = mvdHit->
GetY();
1137 Zpart2 = mvdHit->
GetZ();
1139 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 > Zpart2) {
1140 for (
int i = 0;
i < hits2sts;
i++) {
1141 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
1144 Xpart2 = stsHit->
GetX();
1145 Ypart2 = stsHit->
GetY();
1146 Zpart2 = stsHit->
GetZ();
1152 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
1154 TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 /
TMath::Pi());
1157 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
1165 vector<CbmMCTrack*> MCtracks_minus,
1166 vector<CbmMCTrack*> MCtracks_plus,
1167 vector<CbmStsTrack*> StsTrack_minus,
1168 vector<CbmStsTrack*> StsTrack_plus,
1169 vector<TVector3> Momenta_minus,
1170 vector<TVector3> Momenta_plus,
1171 std::vector<int> Rings_minus,
1172 std::vector<int> Rings_plus,
1173 std::vector<int> stsIndex_minus,
1174 std::vector<int> stsIndex_plus,
1175 vector<CbmRichRing*> richRing_minus,
1176 vector<CbmRichRing*> richRing_plus,
1177 vector<Int_t> MCIndex_minus,
1178 vector<Int_t> MCIndex_plus) {
1179 for (
size_t i = 0;
i < Momenta_minus.size();
i++) {
1180 for (
size_t j = 0; j < Momenta_plus.size(); j++) {
1185 TVector3 part1 = Momenta_minus[
i];
1186 TVector3 part2 = Momenta_plus[j];
1190 int richcheck_0 = 0;
1191 int richcheck_1 = 0;
1196 richcheck_1 =
CheckIfElectron(richRing_plus[j], Momenta_plus[j].Mag());
1199 richcheck_0 = Rings_minus[
i];
1200 richcheck_1 = Rings_plus[j];
1202 int richcheck = richcheck_0 + richcheck_1;
1204 double InvmassReco =
1206 double OpeningAngle =
1223 if (mcTrackgrmama->
GetPdgCode() == 221) { TruePair = 1; }
1227 if (TruePair == 1) {
1228 if (richcheck == 0) {
1249 if (richcheck == 1) {
1269 if (richcheck == 2) {
1296 if (richcheck == 0) {
1320 if (richcheck == 1) {
1342 if (richcheck == 2) {
1368 if (TMath::Abs(OpeningAngle) > AngleCut)
continue;
1369 if (TMath::Abs(InvmassReco) > InvMassCut)
continue;
1371 if (richcheck == 1 && TruePair == 1) {
1377 if (richcheck == 1 && TruePair == 0) {
1388 frefId.push_back(stsIndex_minus[
i]);
1389 frefId.push_back(stsIndex_plus[j]);
1391 fMCId.push_back(MCIndex_minus[
i]);
1392 fMCId.push_back(MCIndex_plus[j]);
1404 if (richcheck == 0 || richcheck == 1 || richcheck == 2) {
1412 if (richcheck == 2) {
1420 if (richcheck == 1 || richcheck == 2) {
1435 vector<CbmMCTrack*> MCtracks_minus_Outside,
1436 vector<CbmMCTrack*> MCtracks_plus_Outside,
1437 vector<CbmStsTrack*> StsTrack_minus_Outside,
1438 vector<CbmStsTrack*> StsTrack_plus_Outside,
1439 std::vector<int> Rings_minus_Outside,
1440 std::vector<int> Rings_plus_Outside,
1441 std::vector<int> stsIndex_minus_Outside,
1442 std::vector<int> stsIndex_plus_Outside,
1443 vector<CbmRichRing*> richRing_minus_Outside,
1444 vector<CbmRichRing*> richRing_plus_Outside,
1445 vector<Int_t> MCIndex_minus_Outside,
1446 vector<Int_t> MCIndex_plus_Outside) {
1447 for (
size_t i = 0;
i < StsTrack_minus_Outside.size();
i++) {
1448 for (
size_t j = 0; j < StsTrack_plus_Outside.size(); j++) {
1454 CbmMCTrack* part2MC = MCtracks_plus_Outside[j];
1456 KFParticle electron;
1458 KFParticle positron;
1460 const KFParticle* daughters[2] = {&electron, &positron};
1461 KFParticle intersection;
1462 intersection.Construct(daughters, 2);
1464 if (intersection.GetZ() > 75 || intersection.GetZ() < -5)
1471 double chi_electron = 0;
1472 double chi_positron = 0;
1475 intersection.GetX(),
1476 intersection.GetY(),
1477 intersection.GetZ(),
1481 intersection.GetX(),
1482 intersection.GetY(),
1483 intersection.GetZ(),
1485 if (chi_electron != chi_electron)
continue;
1486 if (chi_electron == 0)
continue;
1487 if (chi_electron > 3)
continue;
1488 if (chi_positron != chi_positron)
continue;
1489 if (chi_positron == 0)
continue;
1490 if (chi_positron > 3)
continue;
1493 int richcheck_0 = 0;
1494 int richcheck_1 = 0;
1499 richcheck = richcheck_0 + richcheck_1;
1502 richcheck = Rings_minus_Outside[
i] + Rings_plus_Outside[j];
1505 double InvmassReco =
1507 double OpeningAngle =
1524 if (mcTrackgrmama->
GetPdgCode() == 221) { TruePair = 1; }
1528 if (TruePair == 1) {
1529 if (richcheck == 0) {
1551 if (richcheck == 1) {
1571 if (richcheck == 2) {
1598 - intersection.GetX());
1600 - intersection.GetY());
1602 - intersection.GetZ());
1604 - intersection.GetX());
1606 - intersection.GetY());
1608 - intersection.GetZ());
1610 if (richcheck == 0) {
1634 if (richcheck == 1) {
1657 if (richcheck == 2) {
1684 if (TMath::Abs(OpeningAngle) > AngleCut)
continue;
1685 if (TMath::Abs(InvmassReco) > InvMassCut)
continue;
1688 if (richcheck == 1 && TruePair == 1) {
1694 if (richcheck == 1 && TruePair == 0) {
1706 frefId.push_back(stsIndex_minus_Outside[
i]);
1707 frefId.push_back(stsIndex_plus_Outside[j]);
1709 fMCId.push_back(MCIndex_minus_Outside[
i]);
1710 fMCId.push_back(MCIndex_plus_Outside[j]);
1722 if (richcheck == 0 || richcheck == 1 || richcheck == 2) {
1730 if (richcheck == 2) {
1738 if (richcheck == 1 || richcheck == 2) {
1827 vector<vector<TVector3>> Gammas,
1828 vector<vector<int>> StsIndex,
1829 vector<vector<int>> MCIndex,
1830 vector<vector<CbmMCTrack*>> GammasMC,
1832 vector<TH1*> rap_pt_separation) {
1833 if (Gammas.size() < 2)
return;
1834 for (
size_t gamma1 = 0; gamma1 < Gammas.size() - 1; gamma1++) {
1835 for (
size_t gamma2 = gamma1 + 1; gamma2 < Gammas.size(); gamma2++) {
1837 TVector3 e11 = Gammas[gamma1][0];
1838 TVector3 e12 = Gammas[gamma1][1];
1839 TVector3 e21 = Gammas[gamma2][0];
1840 TVector3 e22 = Gammas[gamma2][1];
1848 if (StsIndex[gamma1][0] == StsIndex[gamma2][0]
1849 || StsIndex[gamma1][0] == StsIndex[gamma2][1]
1850 || StsIndex[gamma1][1] == StsIndex[gamma2][0]
1851 || StsIndex[gamma1][1] == StsIndex[gamma2][1])
1856 e11, e12, e21, e22);
1857 double openingAngleBetweenGammasReco =
1859 e11, e12, e21, e22);
1862 int STSmcId1 = MCIndex[gamma1][0];
1863 int STSmcId2 = MCIndex[gamma1][1];
1864 int STSmcId3 = MCIndex[gamma2][0];
1865 int STSmcId4 = MCIndex[gamma2][1];
1866 if (STSmcId1 != STSmcId2 && STSmcId1 != STSmcId3 && STSmcId1 != STSmcId4
1867 && STSmcId2 != STSmcId3 && STSmcId2 != STSmcId4
1868 && STSmcId3 != STSmcId4) {
1872 && TMath::Abs(mcTrack4->
GetPdgCode()) == 11) {
1914 if (TrueEta == 1) gg[0]->Fill(openingAngleBetweenGammasReco);
1915 if (TrueEta == 0) gg[1]->Fill(openingAngleBetweenGammasReco);
1917 if (openingAngleBetweenGammasReco < 10
1918 || openingAngleBetweenGammasReco > 40)
1921 gg[2]->Fill(params.
fMinv);
1923 gg[3]->Fill(TMath::Abs(mcTrack1->
GetPdgCode()));
1924 gg[3]->Fill(TMath::Abs(mcTrack2->
GetPdgCode()));
1925 gg[3]->Fill(TMath::Abs(mcTrack3->
GetPdgCode()));
1926 gg[3]->Fill(TMath::Abs(mcTrack4->
GetPdgCode()));
1929 if (position ==
"Both") {
1931 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
1932 rap_pt_separation[1]->Fill(params.
fMinv);
1933 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
1934 rap_pt_separation[2]->Fill(params.
fMinv);
1935 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
1936 rap_pt_separation[3]->Fill(params.
fMinv);
1937 if (params.
fPt > 1.2 && params.
fPt <= 1.6)
1938 rap_pt_separation[4]->Fill(params.
fMinv);
1941 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
1942 rap_pt_separation[5]->Fill(params.
fMinv);
1943 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
1944 rap_pt_separation[6]->Fill(params.
fMinv);
1945 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
1946 rap_pt_separation[7]->Fill(params.
fMinv);
1947 if (params.
fPt > 1.2 && params.
fPt <= 1.6)
1948 rap_pt_separation[8]->Fill(params.
fMinv);
1951 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
1952 rap_pt_separation[9]->Fill(params.
fMinv);
1953 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
1954 rap_pt_separation[10]->Fill(params.
fMinv);
1955 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
1956 rap_pt_separation[11]->Fill(params.
fMinv);
1957 if (params.
fPt > 1.2 && params.
fPt <= 1.6)
1958 rap_pt_separation[12]->Fill(params.
fMinv);
1961 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
1962 rap_pt_separation[13]->Fill(params.
fMinv);
1963 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
1964 rap_pt_separation[14]->Fill(params.
fMinv);
1965 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
1966 rap_pt_separation[15]->Fill(params.
fMinv);
1969 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
1970 rap_pt_separation[16]->Fill(params.
fMinv);
1971 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
1972 rap_pt_separation[17]->Fill(params.
fMinv);
1979 gg[8]->Fill(params.
fMinv);
1980 gg[10]->Fill(openingAngleBetweenGammasReco);
1986 gg[9]->Fill(params.
fMinv);
1987 gg[11]->Fill(openingAngleBetweenGammasReco);
1997 cout <<
"Mixing for Eta(Target) - nof entries " << nof_Target << endl;
1998 for (Int_t a = 0; a < nof_Target - 1; a++) {
1999 for (Int_t b = a + 1; b < nof_Target; b++) {
2008 e11, e12, e21, e22);
2009 double openingAngleBetweenGammasReco =
2011 e11, e12, e21, e22);
2012 if (openingAngleBetweenGammasReco < 10
2013 || openingAngleBetweenGammasReco > 40)
2034 cout <<
"Mixing for Eta(Outside) - nof entries " << nof_Outside << endl;
2035 for (Int_t a = 0; a < nof_Outside - 1; a++) {
2036 for (Int_t b = a + 1; b < nof_Outside; b++) {
2045 e11, e12, e21, e22);
2046 double openingAngleBetweenGammasReco =
2048 e11, e12, e21, e22);
2049 if (openingAngleBetweenGammasReco < 10
2050 || openingAngleBetweenGammasReco > 40)
2067 vector<TH1*> rap_pt_separation_onetwo,
2068 vector<TH1*> rap_pt_separation_two)
2072 cout <<
"Mixing for Eta(Both) - nof entries " << nof_Both << endl;
2073 for (Int_t a = 0; a < nof_Both - 1; a++) {
2074 for (Int_t b = a + 1; b < nof_Both; b++) {
2083 e11, e12, e21, e22);
2084 double openingAngleBetweenGammasReco =
2086 e11, e12, e21, e22);
2087 if (openingAngleBetweenGammasReco < 10
2088 || openingAngleBetweenGammasReco > 40)
2101 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2102 rap_pt_separation_all[18]->Fill(params.
fMinv);
2103 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2104 rap_pt_separation_all[19]->Fill(params.
fMinv);
2105 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
2106 rap_pt_separation_all[20]->Fill(params.
fMinv);
2107 if (params.
fPt > 1.2 && params.
fPt <= 1.6)
2108 rap_pt_separation_all[21]->Fill(params.
fMinv);
2111 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2112 rap_pt_separation_all[22]->Fill(params.
fMinv);
2113 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2114 rap_pt_separation_all[23]->Fill(params.
fMinv);
2115 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
2116 rap_pt_separation_all[24]->Fill(params.
fMinv);
2117 if (params.
fPt > 1.2 && params.
fPt <= 1.6)
2118 rap_pt_separation_all[25]->Fill(params.
fMinv);
2121 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2122 rap_pt_separation_all[26]->Fill(params.
fMinv);
2123 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2124 rap_pt_separation_all[27]->Fill(params.
fMinv);
2125 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
2126 rap_pt_separation_all[28]->Fill(params.
fMinv);
2127 if (params.
fPt > 1.2 && params.
fPt <= 1.6)
2128 rap_pt_separation_all[29]->Fill(params.
fMinv);
2131 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2132 rap_pt_separation_all[30]->Fill(params.
fMinv);
2133 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2134 rap_pt_separation_all[31]->Fill(params.
fMinv);
2135 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
2136 rap_pt_separation_all[32]->Fill(params.
fMinv);
2139 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2140 rap_pt_separation_all[33]->Fill(params.
fMinv);
2141 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2142 rap_pt_separation_all[34]->Fill(params.
fMinv);
2150 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2151 rap_pt_separation_onetwo[18]->Fill(params.
fMinv);
2152 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2153 rap_pt_separation_onetwo[19]->Fill(params.
fMinv);
2154 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
2155 rap_pt_separation_onetwo[20]->Fill(params.
fMinv);
2156 if (params.
fPt > 1.2 && params.
fPt <= 1.6)
2157 rap_pt_separation_onetwo[21]->Fill(params.
fMinv);
2160 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2161 rap_pt_separation_onetwo[22]->Fill(params.
fMinv);
2162 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2163 rap_pt_separation_onetwo[23]->Fill(params.
fMinv);
2164 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
2165 rap_pt_separation_onetwo[24]->Fill(params.
fMinv);
2166 if (params.
fPt > 1.2 && params.
fPt <= 1.6)
2167 rap_pt_separation_onetwo[25]->Fill(params.
fMinv);
2170 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2171 rap_pt_separation_onetwo[26]->Fill(params.
fMinv);
2172 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2173 rap_pt_separation_onetwo[27]->Fill(params.
fMinv);
2174 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
2175 rap_pt_separation_onetwo[28]->Fill(params.
fMinv);
2176 if (params.
fPt > 1.2 && params.
fPt <= 1.6)
2177 rap_pt_separation_onetwo[29]->Fill(params.
fMinv);
2180 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2181 rap_pt_separation_onetwo[30]->Fill(params.
fMinv);
2182 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2183 rap_pt_separation_onetwo[31]->Fill(params.
fMinv);
2184 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
2185 rap_pt_separation_onetwo[32]->Fill(params.
fMinv);
2188 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2189 rap_pt_separation_onetwo[33]->Fill(params.
fMinv);
2190 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2191 rap_pt_separation_onetwo[34]->Fill(params.
fMinv);
2199 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2200 rap_pt_separation_two[18]->Fill(params.
fMinv);
2201 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2202 rap_pt_separation_two[19]->Fill(params.
fMinv);
2203 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
2204 rap_pt_separation_two[20]->Fill(params.
fMinv);
2205 if (params.
fPt > 1.2 && params.
fPt <= 1.6)
2206 rap_pt_separation_two[21]->Fill(params.
fMinv);
2209 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2210 rap_pt_separation_two[22]->Fill(params.
fMinv);
2211 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2212 rap_pt_separation_two[23]->Fill(params.
fMinv);
2213 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
2214 rap_pt_separation_two[24]->Fill(params.
fMinv);
2215 if (params.
fPt > 1.2 && params.
fPt <= 1.6)
2216 rap_pt_separation_two[25]->Fill(params.
fMinv);
2219 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2220 rap_pt_separation_two[26]->Fill(params.
fMinv);
2221 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2222 rap_pt_separation_two[27]->Fill(params.
fMinv);
2223 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
2224 rap_pt_separation_two[28]->Fill(params.
fMinv);
2225 if (params.
fPt > 1.2 && params.
fPt <= 1.6)
2226 rap_pt_separation_two[29]->Fill(params.
fMinv);
2229 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2230 rap_pt_separation_two[30]->Fill(params.
fMinv);
2231 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2232 rap_pt_separation_two[31]->Fill(params.
fMinv);
2233 if (params.
fPt > 0.8 && params.
fPt <= 1.2)
2234 rap_pt_separation_two[32]->Fill(params.
fMinv);
2237 if (params.
fPt > 0.0 && params.
fPt <= 0.4)
2238 rap_pt_separation_two[33]->Fill(params.
fMinv);
2239 if (params.
fPt > 0.4 && params.
fPt <= 0.8)
2240 rap_pt_separation_two[34]->Fill(params.
fMinv);
2249 gDirectory->mkdir(
"Eta");
2250 gDirectory->cd(
"Eta");
2252 gDirectory->mkdir(
"Cuts_for_Eta_gg");
2253 gDirectory->cd(
"Cuts_for_Eta_gg");
2255 gDirectory->mkdir(
"Cuts_for_Eta_gg_Target");
2256 gDirectory->cd(
"Cuts_for_Eta_gg_Target");
2257 gDirectory->mkdir(
"Cuts_for_Eta_gg_zeroInRich_Target");
2258 gDirectory->cd(
"Cuts_for_Eta_gg_zeroInRich_Target");
2262 gDirectory->cd(
"..");
2263 gDirectory->mkdir(
"Cuts_for_Eta_gg_oneInRich_Target");
2264 gDirectory->cd(
"Cuts_for_Eta_gg_oneInRich_Target");
2268 gDirectory->cd(
"..");
2269 gDirectory->mkdir(
"Cuts_for_Eta_gg_twoInRich_Target");
2270 gDirectory->cd(
"Cuts_for_Eta_gg_twoInRich_Target");
2274 gDirectory->cd(
"..");
2278 gDirectory->cd(
"..");
2280 gDirectory->mkdir(
"Cuts_for_Eta_gg_Outside");
2281 gDirectory->cd(
"Cuts_for_Eta_gg_Outside");
2282 gDirectory->mkdir(
"Cuts_for_Eta_gg_zeroInRich_Outside");
2283 gDirectory->cd(
"Cuts_for_Eta_gg_zeroInRich_Outside");
2287 gDirectory->cd(
"..");
2288 gDirectory->mkdir(
"Cuts_for_Eta_gg_oneInRich_Outside");
2289 gDirectory->cd(
"Cuts_for_Eta_gg_oneInRich_Outside");
2293 gDirectory->cd(
"..");
2294 gDirectory->mkdir(
"Cuts_for_Eta_gg_twoInRich_Outside");
2295 gDirectory->cd(
"Cuts_for_Eta_gg_twoInRich_Outside");
2299 gDirectory->cd(
"..");
2303 gDirectory->cd(
"..");
2305 gDirectory->mkdir(
"Cuts_for_Eta_gg_Both");
2306 gDirectory->cd(
"Cuts_for_Eta_gg_Both");
2307 gDirectory->mkdir(
"Cuts_for_Eta_gg_zeroInRich_Both");
2308 gDirectory->cd(
"Cuts_for_Eta_gg_zeroInRich_Both");
2312 gDirectory->cd(
"..");
2313 gDirectory->mkdir(
"Cuts_for_Eta_gg_oneInRich_Both");
2314 gDirectory->cd(
"Cuts_for_Eta_gg_oneInRich_Both");
2318 gDirectory->cd(
"..");
2319 gDirectory->mkdir(
"Cuts_for_Eta_gg_twoInRich_Both");
2320 gDirectory->cd(
"Cuts_for_Eta_gg_twoInRich_Both");
2324 gDirectory->cd(
"..");
2328 gDirectory->cd(
"..");
2330 gDirectory->cd(
"..");
2333 gDirectory->mkdir(
"Eta_gg_Target");
2334 gDirectory->cd(
"Eta_gg_Target");
2336 gDirectory->mkdir(
"all");
2337 gDirectory->cd(
"all");
2341 gDirectory->cd(
"..");
2343 gDirectory->mkdir(
"onetwo");
2344 gDirectory->cd(
"onetwo");
2348 gDirectory->cd(
"..");
2350 gDirectory->mkdir(
"two");
2351 gDirectory->cd(
"two");
2355 gDirectory->cd(
"..");
2357 gDirectory->cd(
"..");
2360 gDirectory->mkdir(
"Eta_gg_Outside");
2361 gDirectory->cd(
"Eta_gg_Outside");
2363 gDirectory->mkdir(
"all");
2364 gDirectory->cd(
"all");
2368 gDirectory->cd(
"..");
2370 gDirectory->mkdir(
"onetwo");
2371 gDirectory->cd(
"onetwo");
2375 gDirectory->cd(
"..");
2377 gDirectory->mkdir(
"two");
2378 gDirectory->cd(
"two");
2382 gDirectory->cd(
"..");
2384 gDirectory->cd(
"..");
2387 gDirectory->mkdir(
"Eta_gg_Both");
2388 gDirectory->cd(
"Eta_gg_Both");
2390 gDirectory->mkdir(
"all");
2391 gDirectory->cd(
"all");
2392 gDirectory->mkdir(
"multi_all_Both");
2393 gDirectory->cd(
"multi_all_Both");
2397 gDirectory->cd(
"..");
2401 gDirectory->cd(
"..");
2403 gDirectory->mkdir(
"onetwo");
2404 gDirectory->cd(
"onetwo");
2405 gDirectory->mkdir(
"multi_onetwo_Both");
2406 gDirectory->cd(
"multi_onetwo_Both");
2410 gDirectory->cd(
"..");
2414 gDirectory->cd(
"..");
2416 gDirectory->mkdir(
"two");
2417 gDirectory->cd(
"two");
2418 gDirectory->mkdir(
"multi_two_Both");
2419 gDirectory->cd(
"multi_two_Both");
2423 gDirectory->cd(
"..");
2427 gDirectory->cd(
"..");
2429 gDirectory->cd(
"..");
2437 gDirectory->cd(
"..");
2442 new TH1D(
"Particle_pull_px_Target",
2443 "Target track, pull Px: mc - reco; P_{x}(mc) - P_{x}(reco)",
2449 new TH1D(
"Particle_pull_py_Target",
2450 "Target track, pull Px: mc - reco; P_{y}(mc) - P_{y}(reco)",
2456 new TH1D(
"Particle_pull_pz_Target",
2457 "Target track, pull Px: mc - reco; P_{z}(mc) - P_{z}(reco)",
2463 new TH1D(
"Particle_pull_px_Outside",
2464 "Outside track, pull Px: mc - reco; P_{x}(mc) - P_{x}(reco)",
2470 new TH1D(
"Particle_pull_py_Outside",
2471 "Outside track, pull Px: mc - reco; P_{y}(mc) - P_{y}(reco)",
2477 new TH1D(
"Particle_pull_pz_Outside",
2478 "Outside track, pull Px: mc - reco; P_{z}(mc) - P_{z}(reco)",
2484 new TH1D(
"Particle_pull_X_Outside",
2485 "Outside track, pull Px: mc - reco; X(mc) - X(reco)",
2491 new TH1D(
"Particle_pull_Y_Outside",
2492 "Outside track, pull Px: mc - reco; Y(mc) - Y(reco)",
2498 new TH1D(
"Particle_pull_Z_Outside",
2499 "Outside track, pull Px: mc - reco; Z(mc) - Z(reco)",
2507 new TH1D(
"GammasInvMass_fromEta_beforeCuts_zeroInRich_Target",
2508 "GammasInvMass_fromEta_beforeCuts_zeroInRich_Target; invariant "
2509 "mass in GeV/c^{2}",
2516 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Target",
2517 "GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Target; opening "
2525 new TH1D(
"GammasMomentum_fromEta_beforeCuts_zeroInRich_Target",
2526 "GammasMomentum_fromEta_beforeCuts_zeroInRich_Target; p in GeV/c",
2533 new TH1D(
"GammasRapidity_fromEta_beforeCuts_zeroInRich_Target",
2534 "GammasRapidity_fromEta_beforeCuts_zeroInRich_Target; rapidity",
2541 "PlaneAngles_last_fromEta_zeroInRich_Target",
2542 "PlaneAngles_last_fromEta_zeroInRich_Target; #theta angle in degree;#",
2549 "PlaneAngles_first_fromEta_zeroInRich_Target",
2550 "PlaneAngles_first_fromEta_zeroInRich_Target; #theta angle in degree;#",
2557 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Target",
2558 "GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Target; invariant "
2559 "mass in GeV/c^{2}",
2566 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Target",
2567 "GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Target; opening "
2575 "GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Target",
2576 "GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Target; p in GeV/c",
2583 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Target",
2584 "GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Target; rapidity",
2591 "PlaneAngles_last_wrongpairs_zeroInRich_Target",
2592 "PlaneAngles_last_wrongpairs_zeroInRich_Target; #theta angle in degree;#",
2599 "PlaneAngles_first_wrongpairs_zeroInRich_Target",
2600 "PlaneAngles_first_wrongpairs_zeroInRich_Target; #theta angle in degree;#",
2608 new TH1D(
"GammasInvMass_fromEta_beforeCuts_oneInRich_Target",
2609 "GammasInvMass_fromEta_beforeCuts_oneInRich_Target; invariant "
2610 "mass in GeV/c^{2}",
2617 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_oneInRich_Target",
2618 "GammasOpenAngle_fromEta_beforeCuts_oneInRich_Target; opening "
2626 new TH1D(
"GammasMomentum_fromEta_beforeCuts_oneInRich_Target",
2627 "GammasMomentum_fromEta_beforeCuts_oneInRich_Target; p in GeV/c",
2634 new TH1D(
"GammasRapidity_fromEta_beforeCuts_oneInRich_Target",
2635 "GammasRapidity_fromEta_beforeCuts_oneInRich_Target; rapidity",
2642 "PlaneAngles_last_fromEta_oneInRich_Target",
2643 "PlaneAngles_last_fromEta_oneInRich_Target; #theta angle in degree;#",
2650 "PlaneAngles_first_fromEta_oneInRich_Target",
2651 "PlaneAngles_first_fromEta_oneInRich_Target; #theta angle in degree;#",
2658 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_oneInRich_Target",
2659 "GammasInvMass_wrongpairs_beforeCuts_oneInRich_Target; invariant "
2660 "mass in GeV/c^{2}",
2667 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Target",
2668 "GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Target; opening "
2676 "GammasMomentum_wrongpairs_beforeCuts_oneInRich_Target",
2677 "GammasMomentum_wrongpairs_beforeCuts_oneInRich_Target; p in GeV/c",
2684 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_oneInRich_Target",
2685 "GammasRapidity_wrongpairs_beforeCuts_oneInRich_Target; rapidity",
2692 "PlaneAngles_last_wrongpairs_oneInRich_Target",
2693 "PlaneAngles_last_wrongpairs_oneInRich_Target; #theta angle in degree;#",
2700 "PlaneAngles_first_wrongpairs_oneInRich_Target",
2701 "PlaneAngles_first_wrongpairs_oneInRich_Target; #theta angle in degree;#",
2709 new TH1D(
"GammasInvMass_fromEta_beforeCuts_twoInRich_Target",
2710 "GammasInvMass_fromEta_beforeCuts_twoInRich_Target; invariant "
2711 "mass in GeV/c^{2}",
2718 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_twoInRich_Target",
2719 "GammasOpenAngle_fromEta_beforeCuts_twoInRich_Target; opening "
2727 new TH1D(
"GammasMomentum_fromEta_beforeCuts_twoInRich_Target",
2728 "GammasMomentum_fromEta_beforeCuts_twoInRich_Target; p in GeV/c",
2735 new TH1D(
"GammasRapidity_fromEta_beforeCuts_twoInRich_Target",
2736 "GammasRapidity_fromEta_beforeCuts_twoInRich_Target; rapidity",
2743 "PlaneAngles_last_fromEta_twoInRich_Target",
2744 "PlaneAngles_last_fromEta_twoInRich_Target; #theta angle in degree;#",
2751 "PlaneAngles_first_fromEta_twoInRich_Target",
2752 "PlaneAngles_first_fromEta_twoInRich_Target; #theta angle in degree;#",
2759 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_twoInRich_Target",
2760 "GammasInvMass_wrongpairs_beforeCuts_twoInRich_Target; invariant "
2761 "mass in GeV/c^{2}",
2768 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Target",
2769 "GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Target; opening "
2777 "GammasMomentum_wrongpairs_beforeCuts_twoInRich_Target",
2778 "GammasMomentum_wrongpairs_beforeCuts_twoInRich_Target; p in GeV/c",
2785 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_twoInRich_Target",
2786 "GammasRapidity_wrongpairs_beforeCuts_twoInRich_Target; rapidity",
2793 "PlaneAngles_last_wrongpairs_twoInRich_Target",
2794 "PlaneAngles_last_wrongpairs_twoInRich_Target; #theta angle in degree;#",
2801 "PlaneAngles_first_wrongpairs_twoInRich_Target",
2802 "PlaneAngles_first_wrongpairs_twoInRich_Target; #theta angle in degree;#",
2812 new TH1D(
"GammasInvMass_fromEta_beforeCuts_zeroInRich_Outside",
2813 "GammasInvMass_fromEta_beforeCuts_zeroInRich_Outside; invariant "
2814 "mass in GeV/c^{2}",
2821 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Outside",
2822 "GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Outside; opening "
2830 new TH1D(
"GammasMomentum_fromEta_beforeCuts_zeroInRich_Outside",
2831 "GammasMomentum_fromEta_beforeCuts_zeroInRich_Outside; p in GeV/c",
2838 new TH1D(
"GammasRapidity_fromEta_beforeCuts_zeroInRich_Outside",
2839 "GammasRapidity_fromEta_beforeCuts_zeroInRich_Outside; rapidity",
2846 "PlaneAngles_last_fromEta_zeroInRich_Outside",
2847 "PlaneAngles_last_fromEta_zeroInRich_Outside; #theta angle in degree;#",
2854 "PlaneAngles_first_fromEta_zeroInRich_Outside",
2855 "PlaneAngles_first_fromEta_zeroInRich_Outside; #theta angle in degree;#",
2862 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Outside",
2863 "GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Outside; "
2864 "invariant mass in GeV/c^{2}",
2871 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Outside",
2872 "GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Outside; "
2873 "opening angle in degree",
2880 "GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Outside",
2881 "GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Outside; p in GeV/c",
2888 "GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Outside",
2889 "GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Outside; rapidity",
2896 "PlaneAngles_last_wrongpairs_zeroInRich_Outside",
2897 "PlaneAngles_last_wrongpairs_zeroInRich_Outside; #theta angle in degree;#",
2904 "PlaneAngles_first_wrongpairs_zeroInRich_Outside",
2905 "PlaneAngles_first_wrongpairs_zeroInRich_Outside; #theta angle in degree;#",
2913 new TH1D(
"GammasInvMass_fromEta_beforeCuts_oneInRich_Outside",
2914 "GammasInvMass_fromEta_beforeCuts_oneInRich_Outside; invariant "
2915 "mass in GeV/c^{2}",
2922 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_oneInRich_Outside",
2923 "GammasOpenAngle_fromEta_beforeCuts_oneInRich_Outside; opening "
2931 new TH1D(
"GammasMomentum_fromEta_beforeCuts_oneInRich_Outside",
2932 "GammasMomentum_fromEta_beforeCuts_oneInRich_Outside; p in GeV/c",
2939 new TH1D(
"GammasRapidity_fromEta_beforeCuts_oneInRich_Outside",
2940 "GammasRapidity_fromEta_beforeCuts_oneInRich_Outside; rapidity",
2947 "PlaneAngles_last_fromEta_oneInRich_Outside",
2948 "PlaneAngles_last_fromEta_oneInRich_Outside; #theta angle in degree;#",
2955 "PlaneAngles_first_fromEta_oneInRich_Outside",
2956 "PlaneAngles_first_fromEta_oneInRich_Outside; #theta angle in degree;#",
2963 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_oneInRich_Outside",
2964 "GammasInvMass_wrongpairs_beforeCuts_oneInRich_Outside; invariant "
2965 "mass in GeV/c^{2}",
2972 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Outside",
2973 "GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Outside; opening "
2981 "GammasMomentum_wrongpairs_beforeCuts_oneInRich_Outside",
2982 "GammasMomentum_wrongpairs_beforeCuts_oneInRich_Outside; p in GeV/c",
2989 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_oneInRich_Outside",
2990 "GammasRapidity_wrongpairs_beforeCuts_oneInRich_Outside; rapidity",
2997 "PlaneAngles_last_wrongpairs_oneInRich_Outside",
2998 "PlaneAngles_last_wrongpairs_oneInRich_Outside; #theta angle in degree;#",
3005 "PlaneAngles_first_wrongpairs_oneInRich_Outside",
3006 "PlaneAngles_first_wrongpairs_oneInRich_Outside; #theta angle in degree;#",
3014 new TH1D(
"GammasInvMass_fromEta_beforeCuts_twoInRich_Outside",
3015 "GammasInvMass_fromEta_beforeCuts_twoInRich_Outside; invariant "
3016 "mass in GeV/c^{2}",
3023 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_twoInRich_Outside",
3024 "GammasOpenAngle_fromEta_beforeCuts_twoInRich_Outside; opening "
3032 new TH1D(
"GammasMomentum_fromEta_beforeCuts_twoInRich_Outside",
3033 "GammasMomentum_fromEta_beforeCuts_twoInRich_Outside; p in GeV/c",
3040 new TH1D(
"GammasRapidity_fromEta_beforeCuts_twoInRich_Outside",
3041 "GammasRapidity_fromEta_beforeCuts_twoInRich_Outside; rapidity",
3048 "PlaneAngles_last_fromEta_twoInRich_Outside",
3049 "PlaneAngles_last_fromEta_twoInRich_Outside; #theta angle in degree;#",
3056 "PlaneAngles_first_fromEta_twoInRich_Outside",
3057 "PlaneAngles_first_fromEta_twoInRich_Outside; #theta angle in degree;#",
3064 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_twoInRich_Outside",
3065 "GammasInvMass_wrongpairs_beforeCuts_twoInRich_Outside; invariant "
3066 "mass in GeV/c^{2}",
3073 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Outside",
3074 "GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Outside; opening "
3082 "GammasMomentum_wrongpairs_beforeCuts_twoInRich_Outside",
3083 "GammasMomentum_wrongpairs_beforeCuts_twoInRich_Outside; p in GeV/c",
3090 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_twoInRich_Outside",
3091 "GammasRapidity_wrongpairs_beforeCuts_twoInRich_Outside; rapidity",
3098 "PlaneAngles_last_wrongpairs_twoInRich_Outside",
3099 "PlaneAngles_last_wrongpairs_twoInRich_Outside; #theta angle in degree;#",
3106 "PlaneAngles_first_wrongpairs_twoInRich_Outside",
3107 "PlaneAngles_first_wrongpairs_twoInRich_Outside; #theta angle in degree;#",
3117 new TH1D(
"GammasInvMass_fromEta_beforeCuts_zeroInRich_Both",
3118 "GammasInvMass_fromEta_beforeCuts_zeroInRich_Both; invariant mass "
3126 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Both",
3127 "GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Both; opening "
3135 new TH1D(
"GammasMomentum_fromEta_beforeCuts_zeroInRich_Both",
3136 "GammasMomentum_fromEta_beforeCuts_zeroInRich_Both; p in GeV/c",
3143 new TH1D(
"GammasRapidity_fromEta_beforeCuts_zeroInRich_Both",
3144 "GammasRapidity_fromEta_beforeCuts_zeroInRich_Both; rapidity",
3151 "PlaneAngles_last_fromEta_zeroInRich_Both",
3152 "PlaneAngles_last_fromEta_zeroInRich_Both; #theta angle in degree;#",
3159 "PlaneAngles_first_fromEta_zeroInRich_Both",
3160 "PlaneAngles_first_fromEta_zeroInRich_Both; #theta angle in degree;#",
3167 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Both",
3168 "GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Both; invariant "
3169 "mass in GeV/c^{2}",
3176 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Both",
3177 "GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Both; opening "
3185 new TH1D(
"GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Both",
3186 "GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Both; p in GeV/c",
3193 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Both",
3194 "GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Both; rapidity",
3201 "PlaneAngles_last_wrongpairs_zeroInRich_Both",
3202 "PlaneAngles_last_wrongpairs_zeroInRich_Both; #theta angle in degree;#",
3209 "PlaneAngles_first_wrongpairs_zeroInRich_Both",
3210 "PlaneAngles_first_wrongpairs_zeroInRich_Both; #theta angle in degree;#",
3218 new TH1D(
"GammasInvMass_fromEta_beforeCuts_oneInRich_Both",
3219 "GammasInvMass_fromEta_beforeCuts_oneInRich_Both; invariant mass "
3227 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_oneInRich_Both",
3228 "GammasOpenAngle_fromEta_beforeCuts_oneInRich_Both; opening angle "
3236 new TH1D(
"GammasMomentum_fromEta_beforeCuts_oneInRich_Both",
3237 "GammasMomentum_fromEta_beforeCuts_oneInRich_Both; p in GeV/c",
3244 new TH1D(
"GammasRapidity_fromEta_beforeCuts_oneInRich_Both",
3245 "GammasRapidity_fromEta_beforeCuts_oneInRich_Both; rapidity",
3252 "PlaneAngles_last_fromEta_oneInRich_Both",
3253 "PlaneAngles_last_fromEta_oneInRich_Both; #theta angle in degree;#",
3260 "PlaneAngles_first_fromEta_oneInRich_Both",
3261 "PlaneAngles_first_fromEta_oneInRich_Both; #theta angle in degree;#",
3268 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_oneInRich_Both",
3269 "GammasInvMass_wrongpairs_beforeCuts_oneInRich_Both; invariant "
3270 "mass in GeV/c^{2}",
3277 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Both",
3278 "GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Both; opening "
3286 new TH1D(
"GammasMomentum_wrongpairs_beforeCuts_oneInRich_Both",
3287 "GammasMomentum_wrongpairs_beforeCuts_oneInRich_Both; p in GeV/c",
3294 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_oneInRich_Both",
3295 "GammasRapidity_wrongpairs_beforeCuts_oneInRich_Both; rapidity",
3302 "PlaneAngles_last_wrongpairs_oneInRich_Both",
3303 "PlaneAngles_last_wrongpairs_oneInRich_Both; #theta angle in degree;#",
3310 "PlaneAngles_first_wrongpairs_oneInRich_Both",
3311 "PlaneAngles_first_wrongpairs_oneInRich_Both; #theta angle in degree;#",
3318 new TH1D(
"PlaneAngles_last_fromEta_afterCuts_oneInRich_Both",
3319 "PlaneAngles_last_fromEta_afterCuts_oneInRich_Both; #theta angle "
3327 new TH1D(
"PlaneAngles_first_fromEta_afterCuts_oneInRich_Both",
3328 "PlaneAngles_first_fromEta_afterCuts_oneInRich_Both; #theta angle "
3336 new TH1D(
"PlaneAngles_last_wrongpairs_afterCuts_oneInRich_Both",
3337 "PlaneAngles_last_wrongpairs_afterCuts_oneInRich_Both; #theta "
3338 "angle in degree;#",
3345 new TH1D(
"PlaneAngles_first_wrongpairs_afterCuts_oneInRich_Both",
3346 "PlaneAngles_first_wrongpairs_afterCuts_oneInRich_Both; #theta "
3347 "angle in degree;#",
3355 new TH1D(
"GammasInvMass_fromEta_beforeCuts_twoInRich_Both",
3356 "GammasInvMass_fromEta_beforeCuts_twoInRich_Both; invariant mass "
3364 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_twoInRich_Both",
3365 "GammasOpenAngle_fromEta_beforeCuts_twoInRich_Both; opening angle "
3373 new TH1D(
"GammasMomentum_fromEta_beforeCuts_twoInRich_Both",
3374 "GammasMomentum_fromEta_beforeCuts_twoInRich_Both; p in GeV/c",
3381 new TH1D(
"GammasRapidity_fromEta_beforeCuts_twoInRich_Both",
3382 "GammasRapidity_fromEta_beforeCuts_twoInRich_Both; rapidity",
3389 "PlaneAngles_last_fromEta_twoInRich_Both",
3390 "PlaneAngles_last_fromEta_twoInRich_Both; #theta angle in degree;#",
3397 "PlaneAngles_first_fromEta_twoInRich_Both",
3398 "PlaneAngles_first_fromEta_twoInRich_Both; #theta angle in degree;#",
3405 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_twoInRich_Both",
3406 "GammasInvMass_wrongpairs_beforeCuts_twoInRich_Both; invariant "
3407 "mass in GeV/c^{2}",
3414 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Both",
3415 "GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Both; opening "
3423 new TH1D(
"GammasMomentum_wrongpairs_beforeCuts_twoInRich_Both",
3424 "GammasMomentum_wrongpairs_beforeCuts_twoInRich_Both; p in GeV/c",
3431 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_twoInRich_Both",
3432 "GammasRapidity_wrongpairs_beforeCuts_twoInRich_Both; rapidity",
3439 "PlaneAngles_last_wrongpairs_twoInRich_Both",
3440 "PlaneAngles_last_wrongpairs_twoInRich_Both; #theta angle in degree;#",
3447 "PlaneAngles_first_wrongpairs_twoInRich_Both",
3448 "PlaneAngles_first_wrongpairs_twoInRich_Both; #theta angle in degree;#",
3459 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_all_Target",
3460 "AngleBetweenGammas_trueEta_before_cuts_all_Target (between "
3461 "#gamma#gamma from #eta); angle [deg]",
3468 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_all_Target",
3469 "AngleBetweenGammas_wrongEta_before_cuts_all_Target (between "
3470 "#gamma#gamma); angle [deg]",
3477 "InvMass_GammaCombinations_all_Target",
3478 "InvMass_GammaCombinations_all_Target; invariant mass in GeV/c^{2};#",
3484 "Particles_PDG_all_Target ; PDG",
3490 new TH2D(
"TrueEta_pt_vs_rap_all_Target",
3491 "TrueEta_pt_vs_rap_all_Target; rapidity y; p_{t} in GeV/c ",
3500 new TH2D(
"WrongEta_pt_vs_rap_all_Target",
3501 "WrongEta_pt_vs_rap_all_Target; rapidity y; p_{t} in GeV/c ",
3510 new TH2D(
"TrueEta_pt_vs_rap_est_all_Target",
3511 "TrueEta_pt_vs_rap_est_all_Target; rapidity y; p_{t} in GeV/c ",
3520 new TH2D(
"WrongEta_pt_vs_rap_est_all_Target",
3521 "WrongEta_pt_vs_rap_est_all_Target; rapidity y; p_{t} in GeV/c ",
3530 "TrueEta_InvMass_after_cuts_all_Target",
3531 "TrueEta_InvMass_after_cuts_all_Target; invariant mass in GeV/c^{2};#",
3537 "WrongEta_InvMass_after_cuts_all_Target",
3538 "WrongEta_InvMass_after_cuts_all_Target; invariant mass in GeV/c^{2};#",
3544 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_all_Target",
3545 "TrueEta_AngleBetweenGammas_after_cuts_all_Target (between "
3546 "#gamma#gamma from #eta); angle [deg]",
3553 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_all_Target",
3554 "WrongEta_AngleBetweenGammas_after_cuts_all_Target (between "
3555 "#gamma#gamma); angle [deg]",
3562 new TH1D(
"EMT_Eta_InvMass_all_Target",
3563 "EMT_Eta_InvMass_all_Target; invariant mass in GeV/c^{2};#",
3570 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_onetwo_Target",
3571 "AngleBetweenGammas_trueEta_before_cuts_onetwo_Target (between "
3572 "#gamma#gamma from #eta); angle [deg]",
3579 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_onetwo_Target",
3580 "AngleBetweenGammas_wrongEta_before_cuts_onetwo_Target (between "
3581 "#gamma#gamma); angle [deg]",
3588 "InvMass_GammaCombinations_onetwo_Target",
3589 "InvMass_GammaCombinations_onetwo_Target; invariant mass in GeV/c^{2};#",
3596 "Particles_PDG_onetwo_Target ; PDG",
3602 new TH2D(
"TrueEta_pt_vs_rap_onetwo_Target",
3603 "TrueEta_pt_vs_rap_onetwo_Target; rapidity y; p_{t} in GeV/c ",
3612 new TH2D(
"WrongEta_pt_vs_rap_onetwo_Target",
3613 "WrongEta_pt_vs_rap_onetwo_Target; rapidity y; p_{t} in GeV/c ",
3622 new TH2D(
"TrueEta_pt_vs_rap_est_onetwo_Target",
3623 "TrueEta_pt_vs_rap_est_onetwo_Target; rapidity y; p_{t} in GeV/c ",
3632 "WrongEta_pt_vs_rap_est_onetwo_Target",
3633 "WrongEta_pt_vs_rap_est_onetwo_Target; rapidity y; p_{t} in GeV/c ",
3642 "TrueEta_InvMass_after_cuts_onetwo_Target",
3643 "TrueEta_InvMass_after_cuts_onetwo_Target; invariant mass in GeV/c^{2};#",
3650 "WrongEta_InvMass_after_cuts_onetwo_Target",
3651 "WrongEta_InvMass_after_cuts_onetwo_Target; invariant mass in GeV/c^{2};#",
3658 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_onetwo_Target",
3659 "TrueEta_AngleBetweenGammas_after_cuts_onetwo_Target (between "
3660 "#gamma#gamma from #eta); angle [deg]",
3667 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_onetwo_Target",
3668 "WrongEta_AngleBetweenGammas_after_cuts_onetwo_Target (between "
3669 "#gamma#gamma); angle [deg]",
3676 new TH1D(
"EMT_Eta_InvMass_onetwo_Target",
3677 "EMT_Eta_InvMass_onetwo_Target; invariant mass in GeV/c^{2};#",
3684 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_two_Target",
3685 "AngleBetweenGammas_trueEta_before_cuts_two_Target (between "
3686 "#gamma#gamma from #eta); angle [deg]",
3693 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_two_Target",
3694 "AngleBetweenGammas_wrongEta_before_cuts_two_Target (between "
3695 "#gamma#gamma); angle [deg]",
3702 "InvMass_GammaCombinations_two_Target",
3703 "InvMass_GammaCombinations_two_Target; invariant mass in GeV/c^{2};#",
3709 "Particles_PDG_two_Target ; PDG",
3715 new TH2D(
"TrueEta_pt_vs_rap_two_Target",
3716 "TrueEta_pt_vs_rap_two_Target; rapidity y; p_{t} in GeV/c ",
3725 new TH2D(
"WrongEta_pt_vs_rap_two_Target",
3726 "WrongEta_pt_vs_rap_two_Target; rapidity y; p_{t} in GeV/c ",
3735 new TH2D(
"TrueEta_pt_vs_rap_est_two_Target",
3736 "TrueEta_pt_vs_rap_est_two_Target; rapidity y; p_{t} in GeV/c ",
3745 new TH2D(
"WrongEta_pt_vs_rap_est_two_Target",
3746 "WrongEta_pt_vs_rap_est_two_Target; rapidity y; p_{t} in GeV/c ",
3755 "TrueEta_InvMass_after_cuts_two_Target",
3756 "TrueEta_InvMass_after_cuts_two_Target; invariant mass in GeV/c^{2};#",
3762 "WrongEta_InvMass_after_cuts_two_Target",
3763 "WrongEta_InvMass_after_cuts_two_Target; invariant mass in GeV/c^{2};#",
3769 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_two_Target",
3770 "TrueEta_AngleBetweenGammas_after_cuts_two_Target (between "
3771 "#gamma#gamma from #eta); angle [deg]",
3778 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_two_Target",
3779 "WrongEta_AngleBetweenGammas_after_cuts_two_Target (between "
3780 "#gamma#gamma); angle [deg]",
3787 new TH1D(
"EMT_Eta_InvMass_two_Target",
3788 "EMT_Eta_InvMass_two_Target; invariant mass in GeV/c^{2};#",
3799 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_all_Outside",
3800 "AngleBetweenGammas_trueEta_before_cuts_all_Outside (between "
3801 "#gamma#gamma from #eta); angle [deg]",
3808 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_all_Outside",
3809 "AngleBetweenGammas_wrongEta_before_cuts_all_Outside (between "
3810 "#gamma#gamma); angle [deg]",
3817 "InvMass_GammaCombinations_all_Outside",
3818 "InvMass_GammaCombinations_all_Outside; invariant mass in GeV/c^{2};#",
3824 "Particles_PDG_all_Outside ; PDG",
3830 new TH2D(
"TrueEta_pt_vs_rap_all_Outside",
3831 "TrueEta_pt_vs_rap_all_Outside; rapidity y; p_{t} in GeV/c ",
3840 new TH2D(
"WrongEta_pt_vs_rap_all_Outside",
3841 "WrongEta_pt_vs_rap_all_Outside; rapidity y; p_{t} in GeV/c ",
3850 new TH2D(
"TrueEta_pt_vs_rap_est_all_Outside",
3851 "TrueEta_pt_vs_rap_est_all_Outside; rapidity y; p_{t} in GeV/c ",
3860 new TH2D(
"WrongEta_pt_vs_rap_est_all_Outside",
3861 "WrongEta_pt_vs_rap_est_all_Outside; rapidity y; p_{t} in GeV/c ",
3870 "TrueEta_InvMass_after_cuts_all_Outside",
3871 "TrueEta_InvMass_after_cuts_all_Outside; invariant mass in GeV/c^{2};#",
3877 "WrongEta_InvMass_after_cuts_all_Outside",
3878 "WrongEta_InvMass_after_cuts_all_Outside; invariant mass in GeV/c^{2};#",
3884 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_all_Outside",
3885 "TrueEta_AngleBetweenGammas_after_cuts_all_Outside (between "
3886 "#gamma#gamma from #eta); angle [deg]",
3893 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_all_Outside",
3894 "WrongEta_AngleBetweenGammas_after_cuts_all_Outside (between "
3895 "#gamma#gamma); angle [deg]",
3902 new TH1D(
"EMT_Eta_InvMass_all_Outside",
3903 "EMT_Eta_InvMass_all_Outside; invariant mass in GeV/c^{2};#",
3910 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_onetwo_Outside",
3911 "AngleBetweenGammas_trueEta_before_cuts_onetwo_Outside (between "
3912 "#gamma#gamma from #eta); angle [deg]",
3919 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_onetwo_Outside",
3920 "AngleBetweenGammas_wrongEta_before_cuts_onetwo_Outside (between "
3921 "#gamma#gamma); angle [deg]",
3928 "InvMass_GammaCombinations_onetwo_Outside",
3929 "InvMass_GammaCombinations_onetwo_Outside; invariant mass in GeV/c^{2};#",
3936 "Particles_PDG_onetwo_Outside ; PDG",
3942 new TH2D(
"TrueEta_pt_vs_rap_onetwo_Outside",
3943 "TrueEta_pt_vs_rap_onetwo_Outside; rapidity y; p_{t} in GeV/c ",
3952 new TH2D(
"WrongEta_pt_vs_rap_onetwo_Outside",
3953 "WrongEta_pt_vs_rap_onetwo_Outside; rapidity y; p_{t} in GeV/c ",
3962 "TrueEta_pt_vs_rap_est_onetwo_Outside",
3963 "TrueEta_pt_vs_rap_est_onetwo_Outside; rapidity y; p_{t} in GeV/c ",
3972 "WrongEta_pt_vs_rap_est_onetwo_Outside",
3973 "WrongEta_pt_vs_rap_est_onetwo_Outside; rapidity y; p_{t} in GeV/c ",
3983 "TrueEta_InvMass_after_cuts_onetwo_Outside",
3984 "TrueEta_InvMass_after_cuts_onetwo_Outside; invariant mass in GeV/c^{2};#",
3991 "WrongEta_InvMass_after_cuts_onetwo_Outside",
3992 "WrongEta_InvMass_after_cuts_onetwo_Outside; invariant mass in GeV/c^{2};#",
3999 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_onetwo_Outside",
4000 "TrueEta_AngleBetweenGammas_after_cuts_onetwo_Outside (between "
4001 "#gamma#gamma from #eta); angle [deg]",
4008 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_onetwo_Outside",
4009 "WrongEta_AngleBetweenGammas_after_cuts_onetwo_Outside (between "
4010 "#gamma#gamma); angle [deg]",
4017 new TH1D(
"EMT_Eta_InvMass_onetwo_Outside",
4018 "EMT_Eta_InvMass_onetwo_Outside; invariant mass in GeV/c^{2};#",
4025 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_two_Outside",
4026 "AngleBetweenGammas_trueEta_before_cuts_two_Outside (between "
4027 "#gamma#gamma from #eta); angle [deg]",
4034 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_two_Outside",
4035 "AngleBetweenGammas_wrongEta_before_cuts_two_Outside (between "
4036 "#gamma#gamma); angle [deg]",
4043 "InvMass_GammaCombinations_two_Outside",
4044 "InvMass_GammaCombinations_two_Outside; invariant mass in GeV/c^{2};#",
4050 "Particles_PDG_two_Outside ; PDG",
4056 new TH2D(
"TrueEta_pt_vs_rap_two_Outside",
4057 "TrueEta_pt_vs_rap_two_Outside; rapidity y; p_{t} in GeV/c ",
4066 new TH2D(
"WrongEta_pt_vs_rap_two_Outside",
4067 "WrongEta_pt_vs_rap_two_Outside; rapidity y; p_{t} in GeV/c ",
4076 new TH2D(
"TrueEta_pt_vs_rap_est_two_Outside",
4077 "TrueEta_pt_vs_rap_est_two_Outside; rapidity y; p_{t} in GeV/c ",
4086 new TH2D(
"WrongEta_pt_vs_rap_est_two_Outside",
4087 "WrongEta_pt_vs_rap_est_two_Outside; rapidity y; p_{t} in GeV/c ",
4096 "TrueEta_InvMass_after_cuts_two_Outside",
4097 "TrueEta_InvMass_after_cuts_two_Outside; invariant mass in GeV/c^{2};#",
4103 "WrongEta_InvMass_after_cuts_two_Outside",
4104 "WrongEta_InvMass_after_cuts_two_Outside; invariant mass in GeV/c^{2};#",
4110 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_two_Outside",
4111 "TrueEta_AngleBetweenGammas_after_cuts_two_Outside (between "
4112 "#gamma#gamma from #eta); angle [deg]",
4119 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_two_Outside",
4120 "WrongEta_AngleBetweenGammas_after_cuts_two_Outside (between "
4121 "#gamma#gamma); angle [deg]",
4128 new TH1D(
"EMT_Eta_InvMass_two_Outside",
4129 "EMT_Eta_InvMass_two_Outside; invariant mass in GeV/c^{2};#",
4140 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_all_Both",
4141 "AngleBetweenGammas_trueEta_before_cuts_all_Both (between "
4142 "#gamma#gamma from #eta); angle [deg]",
4149 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_all_Both",
4150 "AngleBetweenGammas_wrongEta_before_cuts_all_Both (between "
4151 "#gamma#gamma); angle [deg]",
4158 "InvMass_GammaCombinations_all_Both",
4159 "InvMass_GammaCombinations_all_Both; invariant mass in GeV/c^{2};#",
4165 "Particles_PDG_all_Both",
"Particles_PDG_all_Both ; PDG", 2300, 0, 2300);
4168 new TH2D(
"TrueEta_pt_vs_rap_all_Both",
4169 "TrueEta_pt_vs_rap_all_Both; rapidity y; p_{t} in GeV/c ",
4178 new TH2D(
"WrongEta_pt_vs_rap_all_Both",
4179 "WrongEta_pt_vs_rap_all_Both; rapidity y; p_{t} in GeV/c ",
4188 new TH2D(
"TrueEta_pt_vs_rap_est_all_Both",
4189 "TrueEta_pt_vs_rap_est_all_Both; rapidity y; p_{t} in GeV/c ",
4198 new TH2D(
"WrongEta_pt_vs_rap_est_all_Both",
4199 "WrongEta_pt_vs_rap_est_all_Both; rapidity y; p_{t} in GeV/c ",
4208 "TrueEta_InvMass_after_cuts_all_Both",
4209 "TrueEta_InvMass_after_cuts_all_Both; invariant mass in GeV/c^{2};#",
4215 "WrongEta_InvMass_after_cuts_all_Both",
4216 "WrongEta_InvMass_after_cuts_all_Both; invariant mass in GeV/c^{2};#",
4222 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_all_Both",
4223 "TrueEta_AngleBetweenGammas_after_cuts_all_Both (between "
4224 "#gamma#gamma from #eta); angle [deg]",
4231 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_all_Both",
4232 "WrongEta_AngleBetweenGammas_after_cuts_all_Both (between "
4233 "#gamma#gamma); angle [deg]",
4240 new TH1D(
"EMT_Eta_InvMass_all_Both",
4241 "EMT_Eta_InvMass_all_Both; invariant mass in GeV/c^{2};#",
4249 new TH1D(
"multi_InvMass_Eta_all_Both_1",
4250 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4251 "invariant mass in GeV/c^{2};#",
4257 new TH1D(
"multi_InvMass_Eta_all_Both_2",
4258 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4259 "invariant mass in GeV/c^{2};#",
4265 new TH1D(
"multi_InvMass_Eta_all_Both_3",
4266 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4267 "invariant mass in GeV/c^{2};#",
4273 new TH1D(
"multi_InvMass_Eta_all_Both_4",
4274 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
4275 "invariant mass in GeV/c^{2};#",
4281 new TH1D(
"multi_InvMass_Eta_all_Both_5",
4282 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4283 "invariant mass in GeV/c^{2};#",
4289 new TH1D(
"multi_InvMass_Eta_all_Both_6",
4290 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4291 "invariant mass in GeV/c^{2};#",
4297 new TH1D(
"multi_InvMass_Eta_all_Both_7",
4298 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4299 "invariant mass in GeV/c^{2};#",
4305 new TH1D(
"multi_InvMass_Eta_all_Both_8",
4306 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
4307 "invariant mass in GeV/c^{2};#",
4313 new TH1D(
"multi_InvMass_Eta_all_Both_9",
4314 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4315 "invariant mass in GeV/c^{2};#",
4321 new TH1D(
"multi_InvMass_Eta_all_Both_10",
4322 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4323 "invariant mass in GeV/c^{2};#",
4329 new TH1D(
"multi_InvMass_Eta_all_Both_11",
4330 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4331 "invariant mass in GeV/c^{2};#",
4337 new TH1D(
"multi_InvMass_Eta_all_Both_12",
4338 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
4339 "invariant mass in GeV/c^{2};#",
4345 new TH1D(
"multi_InvMass_Eta_all_Both_13",
4346 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4347 "invariant mass in GeV/c^{2};#",
4353 new TH1D(
"multi_InvMass_Eta_all_Both_14",
4354 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4355 "invariant mass in GeV/c^{2};#",
4361 new TH1D(
"multi_InvMass_Eta_all_Both_15",
4362 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4363 "invariant mass in GeV/c^{2};#",
4369 new TH1D(
"multi_InvMass_Eta_all_Both_16",
4370 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4371 "invariant mass in GeV/c^{2};#",
4377 new TH1D(
"multi_InvMass_Eta_all_Both_17",
4378 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4379 "invariant mass in GeV/c^{2};#",
4385 new TH1D(
"multi_EMT_Eta_all_Both_1",
4386 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4387 "invariant mass in GeV/c^{2};#",
4393 new TH1D(
"multi_EMT_Eta_all_Both_2",
4394 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4395 "invariant mass in GeV/c^{2};#",
4401 new TH1D(
"multi_EMT_Eta_all_Both_3",
4402 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4403 "invariant mass in GeV/c^{2};#",
4409 new TH1D(
"multi_EMT_Eta_all_Both_4",
4410 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
4411 "invariant mass in GeV/c^{2};#",
4417 new TH1D(
"multi_EMT_Eta_all_Both_5",
4418 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4419 "invariant mass in GeV/c^{2};#",
4425 new TH1D(
"multi_EMT_Eta_all_Both_6",
4426 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4427 "invariant mass in GeV/c^{2};#",
4433 new TH1D(
"multi_EMT_Eta_all_Both_7",
4434 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4435 "invariant mass in GeV/c^{2};#",
4441 new TH1D(
"multi_EMT_Eta_all_Both_8",
4442 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
4443 "invariant mass in GeV/c^{2};#",
4449 new TH1D(
"multi_EMT_Eta_all_Both_9",
4450 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4451 "invariant mass in GeV/c^{2};#",
4457 new TH1D(
"multi_EMT_Eta_all_Both_10",
4458 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4459 "invariant mass in GeV/c^{2};#",
4465 new TH1D(
"multi_EMT_Eta_all_Both_11",
4466 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4467 "invariant mass in GeV/c^{2};#",
4473 new TH1D(
"multi_EMT_Eta_all_Both_12",
4474 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
4475 "invariant mass in GeV/c^{2};#",
4481 new TH1D(
"multi_EMT_Eta_all_Both_13",
4482 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4483 "invariant mass in GeV/c^{2};#",
4489 new TH1D(
"multi_EMT_Eta_all_Both_14",
4490 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4491 "invariant mass in GeV/c^{2};#",
4497 new TH1D(
"multi_EMT_Eta_all_Both_15",
4498 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4499 "invariant mass in GeV/c^{2};#",
4505 new TH1D(
"multi_EMT_Eta_all_Both_16",
4506 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4507 "invariant mass in GeV/c^{2};#",
4513 new TH1D(
"multi_EMT_Eta_all_Both_17",
4514 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4515 "invariant mass in GeV/c^{2};#",
4524 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_onetwo_Both",
4525 "AngleBetweenGammas_trueEta_before_cuts_onetwo_Both (between "
4526 "#gamma#gamma from #eta); angle [deg]",
4533 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_onetwo_Both",
4534 "AngleBetweenGammas_wrongEta_before_cuts_onetwo_Both (between "
4535 "#gamma#gamma); angle [deg]",
4542 "InvMass_GammaCombinations_onetwo_Both",
4543 "InvMass_GammaCombinations_onetwo_Both; invariant mass in GeV/c^{2};#",
4549 "Particles_PDG_onetwo_Both ; PDG",
4555 new TH2D(
"TrueEta_pt_vs_rap_onetwo_Both",
4556 "TrueEta_pt_vs_rap_onetwo_Both; rapidity y; p_{t} in GeV/c ",
4565 new TH2D(
"WrongEta_pt_vs_rap_onetwo_Both",
4566 "WrongEta_pt_vs_rap_onetwo_Both; rapidity y; p_{t} in GeV/c ",
4575 new TH2D(
"TrueEta_pt_vs_rap_est_onetwo_Both",
4576 "TrueEta_pt_vs_rap_est_onetwo_Both; rapidity y; p_{t} in GeV/c ",
4585 new TH2D(
"WrongEta_pt_vs_rap_est_onetwo_Both",
4586 "WrongEta_pt_vs_rap_est_onetwo_Both; rapidity y; p_{t} in GeV/c ",
4595 "TrueEta_InvMass_after_cuts_onetwo_Both",
4596 "TrueEta_InvMass_after_cuts_onetwo_Both; invariant mass in GeV/c^{2};#",
4602 "WrongEta_InvMass_after_cuts_onetwo_Both",
4603 "WrongEta_InvMass_after_cuts_onetwo_Both; invariant mass in GeV/c^{2};#",
4609 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_onetwo_Both",
4610 "TrueEta_AngleBetweenGammas_after_cuts_onetwo_Both (between "
4611 "#gamma#gamma from #eta); angle [deg]",
4618 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_onetwo_Both",
4619 "WrongEta_AngleBetweenGammas_after_cuts_onetwo_Both (between "
4620 "#gamma#gamma); angle [deg]",
4627 new TH1D(
"EMT_Eta_InvMass_onetwo_Both",
4628 "EMT_Eta_InvMass_onetwo_Both; invariant mass in GeV/c^{2};#",
4637 new TH1D(
"multi_InvMass_Eta_onetwo_Both_1",
4638 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4639 "invariant mass in GeV/c^{2};#",
4645 new TH1D(
"multi_InvMass_Eta_onetwo_Both_2",
4646 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4647 "invariant mass in GeV/c^{2};#",
4653 new TH1D(
"multi_InvMass_Eta_onetwo_Both_3",
4654 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4655 "invariant mass in GeV/c^{2};#",
4661 new TH1D(
"multi_InvMass_Eta_onetwo_Both_4",
4662 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
4663 "invariant mass in GeV/c^{2};#",
4669 new TH1D(
"multi_InvMass_Eta_onetwo_Both_5",
4670 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4671 "invariant mass in GeV/c^{2};#",
4677 new TH1D(
"multi_InvMass_Eta_onetwo_Both_6",
4678 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4679 "invariant mass in GeV/c^{2};#",
4685 new TH1D(
"multi_InvMass_Eta_onetwo_Both_7",
4686 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4687 "invariant mass in GeV/c^{2};#",
4693 new TH1D(
"multi_InvMass_Eta_onetwo_Both_8",
4694 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
4695 "invariant mass in GeV/c^{2};#",
4701 new TH1D(
"multi_InvMass_Eta_onetwo_Both_9",
4702 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4703 "invariant mass in GeV/c^{2};#",
4709 new TH1D(
"multi_InvMass_Eta_onetwo_Both_10",
4710 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4711 "invariant mass in GeV/c^{2};#",
4717 new TH1D(
"multi_InvMass_Eta_onetwo_Both_11",
4718 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4719 "invariant mass in GeV/c^{2};#",
4725 new TH1D(
"multi_InvMass_Eta_onetwo_Both_12",
4726 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
4727 "invariant mass in GeV/c^{2};#",
4733 new TH1D(
"multi_InvMass_Eta_onetwo_Both_13",
4734 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4735 "invariant mass in GeV/c^{2};#",
4741 new TH1D(
"multi_InvMass_Eta_onetwo_Both_14",
4742 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4743 "invariant mass in GeV/c^{2};#",
4749 new TH1D(
"multi_InvMass_Eta_onetwo_Both_15",
4750 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4751 "invariant mass in GeV/c^{2};#",
4757 new TH1D(
"multi_InvMass_Eta_onetwo_Both_16",
4758 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4759 "invariant mass in GeV/c^{2};#",
4765 new TH1D(
"multi_InvMass_Eta_onetwo_Both_17",
4766 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4767 "invariant mass in GeV/c^{2};#",
4773 new TH1D(
"multi_EMT_Eta_onetwo_Both_1",
4774 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4775 "invariant mass in GeV/c^{2};#",
4781 new TH1D(
"multi_EMT_Eta_onetwo_Both_2",
4782 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4783 "invariant mass in GeV/c^{2};#",
4789 new TH1D(
"multi_EMT_Eta_onetwo_Both_3",
4790 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4791 "invariant mass in GeV/c^{2};#",
4797 new TH1D(
"multi_EMT_Eta_onetwo_Both_4",
4798 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
4799 "invariant mass in GeV/c^{2};#",
4805 new TH1D(
"multi_EMT_Eta_onetwo_Both_5",
4806 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4807 "invariant mass in GeV/c^{2};#",
4813 new TH1D(
"multi_EMT_Eta_onetwo_Both_6",
4814 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4815 "invariant mass in GeV/c^{2};#",
4821 new TH1D(
"multi_EMT_Eta_onetwo_Both_7",
4822 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4823 "invariant mass in GeV/c^{2};#",
4829 new TH1D(
"multi_EMT_Eta_onetwo_Both_8",
4830 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
4831 "invariant mass in GeV/c^{2};#",
4837 new TH1D(
"multi_EMT_Eta_onetwo_Both_9",
4838 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4839 "invariant mass in GeV/c^{2};#",
4845 new TH1D(
"multi_EMT_Eta_onetwo_Both_10",
4846 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4847 "invariant mass in GeV/c^{2};#",
4853 new TH1D(
"multi_EMT_Eta_onetwo_Both_11",
4854 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4855 "invariant mass in GeV/c^{2};#",
4861 new TH1D(
"multi_EMT_Eta_onetwo_Both_12",
4862 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
4863 "invariant mass in GeV/c^{2};#",
4869 new TH1D(
"multi_EMT_Eta_onetwo_Both_13",
4870 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4871 "invariant mass in GeV/c^{2};#",
4877 new TH1D(
"multi_EMT_Eta_onetwo_Both_14",
4878 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4879 "invariant mass in GeV/c^{2};#",
4885 new TH1D(
"multi_EMT_Eta_onetwo_Both_15",
4886 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
4887 "invariant mass in GeV/c^{2};#",
4893 new TH1D(
"multi_EMT_Eta_onetwo_Both_16",
4894 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
4895 "invariant mass in GeV/c^{2};#",
4901 new TH1D(
"multi_EMT_Eta_onetwo_Both_17",
4902 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
4903 "invariant mass in GeV/c^{2};#",
4912 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_two_Both",
4913 "AngleBetweenGammas_trueEta_before_cuts_two_Both (between "
4914 "#gamma#gamma from #eta); angle [deg]",
4921 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_two_Both",
4922 "AngleBetweenGammas_wrongEta_before_cuts_two_Both (between "
4923 "#gamma#gamma); angle [deg]",
4930 "InvMass_GammaCombinations_two_Both",
4931 "InvMass_GammaCombinations_two_Both; invariant mass in GeV/c^{2};#",
4937 "Particles_PDG_two_Both",
"Particles_PDG_two_Both ; PDG", 2300, 0, 2300);
4940 new TH2D(
"TrueEta_pt_vs_rap_two_Both",
4941 "TrueEta_pt_vs_rap_two_Both; rapidity y; p_{t} in GeV/c ",
4950 new TH2D(
"WrongEta_pt_vs_rap_two_Both",
4951 "WrongEta_pt_vs_rap_two_Both; rapidity y; p_{t} in GeV/c ",
4960 new TH2D(
"TrueEta_pt_vs_rap_est_two_Both",
4961 "TrueEta_pt_vs_rap_est_two_Both; rapidity y; p_{t} in GeV/c ",
4970 new TH2D(
"WrongEta_pt_vs_rap_est_two_Both",
4971 "WrongEta_pt_vs_rap_est_two_Both; rapidity y; p_{t} in GeV/c ",
4980 "TrueEta_InvMass_after_cuts_two_Both",
4981 "TrueEta_InvMass_after_cuts_two_Both; invariant mass in GeV/c^{2};#",
4987 "WrongEta_InvMass_after_cuts_two_Both",
4988 "WrongEta_InvMass_after_cuts_two_Both; invariant mass in GeV/c^{2};#",
4994 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_two_Both",
4995 "TrueEta_AngleBetweenGammas_after_cuts_two_Both (between "
4996 "#gamma#gamma from #eta); angle [deg]",
5003 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_two_Both",
5004 "WrongEta_AngleBetweenGammas_after_cuts_two_Both (between "
5005 "#gamma#gamma); angle [deg]",
5012 new TH1D(
"EMT_Eta_InvMass_two_Both",
5013 "EMT_Eta_InvMass_two_Both; invariant mass in GeV/c^{2};#",
5021 new TH1D(
"multi_InvMass_Eta_two_Both_1",
5022 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
5023 "invariant mass in GeV/c^{2};#",
5029 new TH1D(
"multi_InvMass_Eta_two_Both_2",
5030 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
5031 "invariant mass in GeV/c^{2};#",
5037 new TH1D(
"multi_InvMass_Eta_two_Both_3",
5038 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
5039 "invariant mass in GeV/c^{2};#",
5045 new TH1D(
"multi_InvMass_Eta_two_Both_4",
5046 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
5047 "invariant mass in GeV/c^{2};#",
5053 new TH1D(
"multi_InvMass_Eta_two_Both_5",
5054 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
5055 "invariant mass in GeV/c^{2};#",
5061 new TH1D(
"multi_InvMass_Eta_two_Both_6",
5062 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
5063 "invariant mass in GeV/c^{2};#",
5069 new TH1D(
"multi_InvMass_Eta_two_Both_7",
5070 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
5071 "invariant mass in GeV/c^{2};#",
5077 new TH1D(
"multi_InvMass_Eta_two_Both_8",
5078 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
5079 "invariant mass in GeV/c^{2};#",
5085 new TH1D(
"multi_InvMass_Eta_two_Both_9",
5086 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
5087 "invariant mass in GeV/c^{2};#",
5093 new TH1D(
"multi_InvMass_Eta_two_Both_10",
5094 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
5095 "invariant mass in GeV/c^{2};#",
5101 new TH1D(
"multi_InvMass_Eta_two_Both_11",
5102 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
5103 "invariant mass in GeV/c^{2};#",
5109 new TH1D(
"multi_InvMass_Eta_two_Both_12",
5110 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
5111 "invariant mass in GeV/c^{2};#",
5117 new TH1D(
"multi_InvMass_Eta_two_Both_13",
5118 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
5119 "invariant mass in GeV/c^{2};#",
5125 new TH1D(
"multi_InvMass_Eta_two_Both_14",
5126 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
5127 "invariant mass in GeV/c^{2};#",
5133 new TH1D(
"multi_InvMass_Eta_two_Both_15",
5134 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
5135 "invariant mass in GeV/c^{2};#",
5141 new TH1D(
"multi_InvMass_Eta_two_Both_16",
5142 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
5143 "invariant mass in GeV/c^{2};#",
5149 new TH1D(
"multi_InvMass_Eta_two_Both_17",
5150 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
5151 "invariant mass in GeV/c^{2};#",
5157 new TH1D(
"multi_EMT_Eta_two_Both_1",
5158 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
5159 "invariant mass in GeV/c^{2};#",
5165 new TH1D(
"multi_EMT_Eta_two_Both_2",
5166 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
5167 "invariant mass in GeV/c^{2};#",
5173 new TH1D(
"multi_EMT_Eta_two_Both_3",
5174 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
5175 "invariant mass in GeV/c^{2};#",
5181 new TH1D(
"multi_EMT_Eta_two_Both_4",
5182 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
5183 "invariant mass in GeV/c^{2};#",
5189 new TH1D(
"multi_EMT_Eta_two_Both_5",
5190 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
5191 "invariant mass in GeV/c^{2};#",
5197 new TH1D(
"multi_EMT_Eta_two_Both_6",
5198 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
5199 "invariant mass in GeV/c^{2};#",
5205 new TH1D(
"multi_EMT_Eta_two_Both_7",
5206 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
5207 "invariant mass in GeV/c^{2};#",
5213 new TH1D(
"multi_EMT_Eta_two_Both_8",
5214 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
5215 "invariant mass in GeV/c^{2};#",
5221 new TH1D(
"multi_EMT_Eta_two_Both_9",
5222 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
5223 "invariant mass in GeV/c^{2};#",
5229 new TH1D(
"multi_EMT_Eta_two_Both_10",
5230 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
5231 "invariant mass in GeV/c^{2};#",
5237 new TH1D(
"multi_EMT_Eta_two_Both_11",
5238 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
5239 "invariant mass in GeV/c^{2};#",
5245 new TH1D(
"multi_EMT_Eta_two_Both_12",
5246 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
5247 "invariant mass in GeV/c^{2};#",
5253 new TH1D(
"multi_EMT_Eta_two_Both_13",
5254 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
5255 "invariant mass in GeV/c^{2};#",
5261 new TH1D(
"multi_EMT_Eta_two_Both_14",
5262 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
5263 "invariant mass in GeV/c^{2};#",
5269 new TH1D(
"multi_EMT_Eta_two_Both_15",
5270 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
5271 "invariant mass in GeV/c^{2};#",
5277 new TH1D(
"multi_EMT_Eta_two_Both_16",
5278 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
5279 "invariant mass in GeV/c^{2};#",
5285 new TH1D(
"multi_EMT_Eta_two_Both_17",
5286 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
5287 "invariant mass in GeV/c^{2};#",