CbmRoot
Loading...
Searching...
No Matches
CbmKresEta.cxx
Go to the documentation of this file.
1/* Copyright (C) 2018-2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
2 SPDX-License-Identifier: GPL-3.0-only
3 Authors: Ievgenii Kres, Florian Uhlig [committer] */
4
21#include "CbmKresEta.h"
22
23#include "CbmGlobalTrack.h"
25#include "CbmKresFunctions.h"
26#include "CbmMCTrack.h"
27#include "CbmMvdHit.h"
28#include "CbmRichHit.h"
29#include "CbmRichRing.h"
31#include "CbmRichRingLight.h"
32#include "CbmRichUtil.h"
33#include "CbmStsHit.h"
34#include "CbmStsTrack.h"
35#include "CbmTrackMatchNew.h"
36
37#include "FairRootManager.h"
38
39#include "TDirectory.h"
40
41#include <iostream>
42
43#include "KFParticle.h"
44
45
46using namespace std;
47
49 : fMcTracks(nullptr)
50 , fGlobalTracks(nullptr)
51 , fStsTracks(nullptr)
52 , fStsTrackMatches(nullptr)
53 , fRichProjections(nullptr)
54 , fRichRings(nullptr)
55 , fRichRingMatches(nullptr)
56 , fRichHits(nullptr)
57 , fArrayMvdHit(nullptr)
58 , fArrayStsHit(nullptr)
59 , fPrimVertex(nullptr)
60 , fKFVertex()
61 , fTauFit(nullptr)
62 , VStsTrack_minus_Outside()
63 , VMCtracks_minus_Outside()
64 , VRings_minus_Outside()
65 , VStsIndex_minus_Outside()
66 , VRichRing_minus_Outside()
67 , VMCIndex_minus_Outside()
68 , VStsTrack_plus_Outside()
69 , VMCtracks_plus_Outside()
70 , VRings_plus_Outside()
71 , VStsIndex_plus_Outside()
72 , VRichRing_plus_Outside()
73 , VMCIndex_plus_Outside()
74 , VMCtracks_minus_Target()
75 , VStsTrack_minus_Target()
76 , VMomenta_minus_Target()
77 , VRings_minus_Target()
78 , VStsIndex_minus_Target()
79 , VRichRing_minus_Target()
80 , VMCIndex_minus_Target()
81 , VMCtracks_plus_Target()
82 , VStsTrack_plus_Target()
83 , VMomenta_plus_Target()
84 , VRings_plus_Target()
85 , VStsIndex_plus_Target()
86 , VRichRing_plus_Target()
87 , VMCIndex_plus_Target()
88 , frefmomenta()
89 , frefId()
90 , fMCId()
91 , fMCtracks()
92 , EMT_eta_gg_Event_Target()
93 , EMT_eta_gg_pair_momenta_Target()
94 , EMT_eta_gg_NofRings_Target()
95 , EMT_eta_gg_Event_Outside()
96 , EMT_eta_gg_pair_momenta_Outside()
97 , EMT_eta_gg_NofRings_Outside()
98 , EMT_eta_gg_Event_Both()
99 , EMT_eta_gg_pair_momenta_Both()
100 , EMT_eta_gg_NofRings_Both()
101 , Gammas_all_Target()
102 , Gammas_two_Target()
103 , Gammas_onetwo_Target()
104 , Gammas_stsIndex_all_Target()
105 , Gammas_stsIndex_two_Target()
106 , Gammas_stsIndex_onetwo_Target()
107 , Gammas_MCIndex_all_Target()
108 , Gammas_MCIndex_two_Target()
109 , Gammas_MCIndex_onetwo_Target()
110 , Gammas_MC_all_Target()
111 , Gammas_MC_two_Target()
112 , Gammas_MC_onetwo_Target()
113 , Gammas_all_Outside()
114 , Gammas_two_Outside()
115 , Gammas_onetwo_Outside()
116 , Gammas_stsIndex_all_Outside()
117 , Gammas_stsIndex_two_Outside()
118 , Gammas_stsIndex_onetwo_Outside()
119 , Gammas_MCIndex_all_Outside()
120 , Gammas_MCIndex_two_Outside()
121 , Gammas_MCIndex_onetwo_Outside()
122 , Gammas_MC_all_Outside()
123 , Gammas_MC_two_Outside()
124 , Gammas_MC_onetwo_Outside()
125 , Gammas_all_Both()
126 , Gammas_two_Both()
127 , Gammas_onetwo_Both()
128 , Gammas_stsIndex_all_Both()
129 , Gammas_stsIndex_two_Both()
130 , Gammas_stsIndex_onetwo_Both()
131 , Gammas_MCIndex_all_Both()
132 , Gammas_MCIndex_two_Both()
133 , Gammas_MCIndex_onetwo_Both()
134 , Gammas_MC_all_Both()
135 , Gammas_MC_two_Both()
136 , Gammas_MC_onetwo_Both()
137 , fHistoList_Eta()
138 , Particle_pull_px_Target(nullptr)
139 , Particle_pull_py_Target(nullptr)
140 , Particle_pull_pz_Target(nullptr)
141 , Particle_pull_px_Outside(nullptr)
142 , Particle_pull_py_Outside(nullptr)
143 , Particle_pull_pz_Outside(nullptr)
144 , Particle_pull_X_Outside(nullptr)
145 , Particle_pull_Y_Outside(nullptr)
146 , Particle_pull_Z_Outside(nullptr)
147 , fHistoList_Eta_cuts_Target()
148 , fHistoList_Eta_cuts_zeroInRich_Target()
149 , GammasInvMass_fromEta_beforeCuts_zeroInRich_Target(nullptr)
150 , GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Target(nullptr)
151 , GammasMomentum_fromEta_beforeCuts_zeroInRich_Target(nullptr)
152 , GammasRapidity_fromEta_beforeCuts_zeroInRich_Target(nullptr)
153 , PlaneAngles_last_fromEta_zeroInRich_Target(nullptr)
154 , PlaneAngles_first_fromEta_zeroInRich_Target(nullptr)
155 , GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Target(nullptr)
156 , GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Target(nullptr)
157 , GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Target(nullptr)
158 , GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Target(nullptr)
159 , PlaneAngles_last_wrongpairs_zeroInRich_Target(nullptr)
160 , PlaneAngles_first_wrongpairs_zeroInRich_Target(nullptr)
161 , fHistoList_Eta_cuts_oneInRich_Target()
162 , GammasInvMass_fromEta_beforeCuts_oneInRich_Target(nullptr)
163 , GammasOpenAngle_fromEta_beforeCuts_oneInRich_Target(nullptr)
164 , GammasMomentum_fromEta_beforeCuts_oneInRich_Target(nullptr)
165 , GammasRapidity_fromEta_beforeCuts_oneInRich_Target(nullptr)
166 , PlaneAngles_last_fromEta_oneInRich_Target(nullptr)
167 , PlaneAngles_first_fromEta_oneInRich_Target(nullptr)
168 , GammasInvMass_wrongpairs_beforeCuts_oneInRich_Target(nullptr)
169 , GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Target(nullptr)
170 , GammasMomentum_wrongpairs_beforeCuts_oneInRich_Target(nullptr)
171 , GammasRapidity_wrongpairs_beforeCuts_oneInRich_Target(nullptr)
172 , PlaneAngles_last_wrongpairs_oneInRich_Target(nullptr)
173 , PlaneAngles_first_wrongpairs_oneInRich_Target(nullptr)
174 , fHistoList_Eta_cuts_twoInRich_Target()
175 , GammasInvMass_fromEta_beforeCuts_twoInRich_Target(nullptr)
176 , GammasOpenAngle_fromEta_beforeCuts_twoInRich_Target(nullptr)
177 , GammasMomentum_fromEta_beforeCuts_twoInRich_Target(nullptr)
178 , GammasRapidity_fromEta_beforeCuts_twoInRich_Target(nullptr)
179 , PlaneAngles_last_fromEta_twoInRich_Target(nullptr)
180 , PlaneAngles_first_fromEta_twoInRich_Target(nullptr)
181 , GammasInvMass_wrongpairs_beforeCuts_twoInRich_Target(nullptr)
182 , GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Target(nullptr)
183 , GammasMomentum_wrongpairs_beforeCuts_twoInRich_Target(nullptr)
184 , GammasRapidity_wrongpairs_beforeCuts_twoInRich_Target(nullptr)
185 , PlaneAngles_last_wrongpairs_twoInRich_Target(nullptr)
186 , PlaneAngles_first_wrongpairs_twoInRich_Target(nullptr)
187 , fHistoList_Eta_cuts_Outside()
188 , fHistoList_Eta_cuts_zeroInRich_Outside()
189 , GammasInvMass_fromEta_beforeCuts_zeroInRich_Outside(nullptr)
190 , GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Outside(nullptr)
191 , GammasMomentum_fromEta_beforeCuts_zeroInRich_Outside(nullptr)
192 , GammasRapidity_fromEta_beforeCuts_zeroInRich_Outside(nullptr)
193 , PlaneAngles_last_fromEta_zeroInRich_Outside(nullptr)
194 , PlaneAngles_first_fromEta_zeroInRich_Outside(nullptr)
195 , GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Outside(nullptr)
196 , GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Outside(nullptr)
197 , GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Outside(nullptr)
198 , GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Outside(nullptr)
199 , PlaneAngles_last_wrongpairs_zeroInRich_Outside(nullptr)
200 , PlaneAngles_first_wrongpairs_zeroInRich_Outside(nullptr)
201 , fHistoList_Eta_cuts_oneInRich_Outside()
202 , GammasInvMass_fromEta_beforeCuts_oneInRich_Outside(nullptr)
203 , GammasOpenAngle_fromEta_beforeCuts_oneInRich_Outside(nullptr)
204 , GammasMomentum_fromEta_beforeCuts_oneInRich_Outside(nullptr)
205 , GammasRapidity_fromEta_beforeCuts_oneInRich_Outside(nullptr)
206 , PlaneAngles_last_fromEta_oneInRich_Outside(nullptr)
207 , PlaneAngles_first_fromEta_oneInRich_Outside(nullptr)
208 , GammasInvMass_wrongpairs_beforeCuts_oneInRich_Outside(nullptr)
209 , GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Outside(nullptr)
210 , GammasMomentum_wrongpairs_beforeCuts_oneInRich_Outside(nullptr)
211 , GammasRapidity_wrongpairs_beforeCuts_oneInRich_Outside(nullptr)
212 , PlaneAngles_last_wrongpairs_oneInRich_Outside(nullptr)
213 , PlaneAngles_first_wrongpairs_oneInRich_Outside(nullptr)
214 , fHistoList_Eta_cuts_twoInRich_Outside()
215 , GammasInvMass_fromEta_beforeCuts_twoInRich_Outside(nullptr)
216 , GammasOpenAngle_fromEta_beforeCuts_twoInRich_Outside(nullptr)
217 , GammasMomentum_fromEta_beforeCuts_twoInRich_Outside(nullptr)
218 , GammasRapidity_fromEta_beforeCuts_twoInRich_Outside(nullptr)
219 , PlaneAngles_last_fromEta_twoInRich_Outside(nullptr)
220 , PlaneAngles_first_fromEta_twoInRich_Outside(nullptr)
221 , GammasInvMass_wrongpairs_beforeCuts_twoInRich_Outside(nullptr)
222 , GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Outside(nullptr)
223 , GammasMomentum_wrongpairs_beforeCuts_twoInRich_Outside(nullptr)
224 , GammasRapidity_wrongpairs_beforeCuts_twoInRich_Outside(nullptr)
225 , PlaneAngles_last_wrongpairs_twoInRich_Outside(nullptr)
226 , PlaneAngles_first_wrongpairs_twoInRich_Outside(nullptr)
227 , fHistoList_Eta_cuts_Both()
228 , fHistoList_Eta_cuts_zeroInRich_Both()
229 , GammasInvMass_fromEta_beforeCuts_zeroInRich_Both(nullptr)
230 , GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Both(nullptr)
231 , GammasMomentum_fromEta_beforeCuts_zeroInRich_Both(nullptr)
232 , GammasRapidity_fromEta_beforeCuts_zeroInRich_Both(nullptr)
233 , PlaneAngles_last_fromEta_zeroInRich_Both(nullptr)
234 , PlaneAngles_first_fromEta_zeroInRich_Both(nullptr)
235 , GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Both(nullptr)
236 , GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Both(nullptr)
237 , GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Both(nullptr)
238 , GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Both(nullptr)
239 , PlaneAngles_last_wrongpairs_zeroInRich_Both(nullptr)
240 , PlaneAngles_first_wrongpairs_zeroInRich_Both(nullptr)
241 , fHistoList_Eta_cuts_oneInRich_Both()
242 , GammasInvMass_fromEta_beforeCuts_oneInRich_Both(nullptr)
243 , GammasOpenAngle_fromEta_beforeCuts_oneInRich_Both(nullptr)
244 , GammasMomentum_fromEta_beforeCuts_oneInRich_Both(nullptr)
245 , GammasRapidity_fromEta_beforeCuts_oneInRich_Both(nullptr)
246 , PlaneAngles_last_fromEta_oneInRich_Both(nullptr)
247 , PlaneAngles_first_fromEta_oneInRich_Both(nullptr)
248 , GammasInvMass_wrongpairs_beforeCuts_oneInRich_Both(nullptr)
249 , GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Both(nullptr)
250 , GammasMomentum_wrongpairs_beforeCuts_oneInRich_Both(nullptr)
251 , GammasRapidity_wrongpairs_beforeCuts_oneInRich_Both(nullptr)
252 , PlaneAngles_last_wrongpairs_oneInRich_Both(nullptr)
253 , PlaneAngles_first_wrongpairs_oneInRich_Both(nullptr)
254 , PlaneAngles_last_fromEta_afterCuts_oneInRich_Both(nullptr)
255 , PlaneAngles_first_fromEta_afterCuts_oneInRich_Both(nullptr)
256 , PlaneAngles_last_wrongpairs_afterCuts_oneInRich_Both(nullptr)
257 , PlaneAngles_first_wrongpairs_afterCuts_oneInRich_Both(nullptr)
258 , fHistoList_Eta_cuts_twoInRich_Both()
259 , GammasInvMass_fromEta_beforeCuts_twoInRich_Both(nullptr)
260 , GammasOpenAngle_fromEta_beforeCuts_twoInRich_Both(nullptr)
261 , GammasMomentum_fromEta_beforeCuts_twoInRich_Both(nullptr)
262 , GammasRapidity_fromEta_beforeCuts_twoInRich_Both(nullptr)
263 , PlaneAngles_last_fromEta_twoInRich_Both(nullptr)
264 , PlaneAngles_first_fromEta_twoInRich_Both(nullptr)
265 , GammasInvMass_wrongpairs_beforeCuts_twoInRich_Both(nullptr)
266 , GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Both(nullptr)
267 , GammasMomentum_wrongpairs_beforeCuts_twoInRich_Both(nullptr)
268 , GammasRapidity_wrongpairs_beforeCuts_twoInRich_Both(nullptr)
269 , PlaneAngles_last_wrongpairs_twoInRich_Both(nullptr)
270 , PlaneAngles_first_wrongpairs_twoInRich_Both(nullptr)
271 , fHistoList_Eta_all_Target()
272 , AngleBetweenGammas_trueEta_before_cuts_all_Target(nullptr)
273 , AngleBetweenGammas_wrongEta_before_cuts_all_Target(nullptr)
274 , InvMass_GammaCombinations_all_Target(nullptr)
275 , Particles_PDG_all_Target(nullptr)
276 , TrueEta_pt_vs_rap_all_Target(nullptr)
277 , WrongEta_pt_vs_rap_all_Target(nullptr)
278 , TrueEta_pt_vs_rap_est_all_Target(nullptr)
279 , WrongEta_pt_vs_rap_est_all_Target(nullptr)
280 , TrueEta_InvMass_after_cuts_all_Target(nullptr)
281 , WrongEta_InvMass_after_cuts_all_Target(nullptr)
282 , TrueEta_AngleBetweenGammas_after_cuts_all_Target(nullptr)
283 , WrongEta_AngleBetweenGammas_after_cuts_all_Target(nullptr)
284 , EMT_Eta_InvMass_all_Target(nullptr)
285 , fHistoList_Eta_onetwo_Target()
286 , AngleBetweenGammas_trueEta_before_cuts_onetwo_Target(nullptr)
287 , AngleBetweenGammas_wrongEta_before_cuts_onetwo_Target(nullptr)
288 , InvMass_GammaCombinations_onetwo_Target(nullptr)
289 , Particles_PDG_onetwo_Target(nullptr)
290 , TrueEta_pt_vs_rap_onetwo_Target(nullptr)
291 , WrongEta_pt_vs_rap_onetwo_Target(nullptr)
292 , TrueEta_pt_vs_rap_est_onetwo_Target(nullptr)
293 , WrongEta_pt_vs_rap_est_onetwo_Target(nullptr)
294 , TrueEta_InvMass_after_cuts_onetwo_Target(nullptr)
295 , WrongEta_InvMass_after_cuts_onetwo_Target(nullptr)
296 , TrueEta_AngleBetweenGammas_after_cuts_onetwo_Target(nullptr)
297 , WrongEta_AngleBetweenGammas_after_cuts_onetwo_Target(nullptr)
298 , EMT_Eta_InvMass_onetwo_Target(nullptr)
299 , fHistoList_Eta_two_Target()
300 , AngleBetweenGammas_trueEta_before_cuts_two_Target(nullptr)
301 , AngleBetweenGammas_wrongEta_before_cuts_two_Target(nullptr)
302 , InvMass_GammaCombinations_two_Target(nullptr)
303 , Particles_PDG_two_Target(nullptr)
304 , TrueEta_pt_vs_rap_two_Target(nullptr)
305 , WrongEta_pt_vs_rap_two_Target(nullptr)
306 , TrueEta_pt_vs_rap_est_two_Target(nullptr)
307 , WrongEta_pt_vs_rap_est_two_Target(nullptr)
308 , TrueEta_InvMass_after_cuts_two_Target(nullptr)
309 , WrongEta_InvMass_after_cuts_two_Target(nullptr)
310 , TrueEta_AngleBetweenGammas_after_cuts_two_Target(nullptr)
311 , WrongEta_AngleBetweenGammas_after_cuts_two_Target(nullptr)
312 , EMT_Eta_InvMass_two_Target(nullptr)
313 , fHistoList_rap_pt_Eta_all_Target()
314 , fHistoList_rap_pt_Eta_onetwo_Target()
315 , fHistoList_rap_pt_Eta_two_Target()
316 , fHistoList_Eta_all_Outside()
317 , AngleBetweenGammas_trueEta_before_cuts_all_Outside(nullptr)
318 , AngleBetweenGammas_wrongEta_before_cuts_all_Outside(nullptr)
319 , InvMass_GammaCombinations_all_Outside(nullptr)
320 , Particles_PDG_all_Outside(nullptr)
321 , TrueEta_pt_vs_rap_all_Outside(nullptr)
322 , WrongEta_pt_vs_rap_all_Outside(nullptr)
323 , TrueEta_pt_vs_rap_est_all_Outside(nullptr)
324 , WrongEta_pt_vs_rap_est_all_Outside(nullptr)
325 , TrueEta_InvMass_after_cuts_all_Outside(nullptr)
326 , WrongEta_InvMass_after_cuts_all_Outside(nullptr)
327 , TrueEta_AngleBetweenGammas_after_cuts_all_Outside(nullptr)
328 , WrongEta_AngleBetweenGammas_after_cuts_all_Outside(nullptr)
329 , EMT_Eta_InvMass_all_Outside(nullptr)
330 , fHistoList_Eta_onetwo_Outside()
331 , AngleBetweenGammas_trueEta_before_cuts_onetwo_Outside(nullptr)
332 , AngleBetweenGammas_wrongEta_before_cuts_onetwo_Outside(nullptr)
333 , InvMass_GammaCombinations_onetwo_Outside(nullptr)
334 , Particles_PDG_onetwo_Outside(nullptr)
335 , TrueEta_pt_vs_rap_onetwo_Outside(nullptr)
336 , WrongEta_pt_vs_rap_onetwo_Outside(nullptr)
337 , TrueEta_pt_vs_rap_est_onetwo_Outside(nullptr)
338 , WrongEta_pt_vs_rap_est_onetwo_Outside(nullptr)
339 , TrueEta_InvMass_after_cuts_onetwo_Outside(nullptr)
340 , WrongEta_InvMass_after_cuts_onetwo_Outside(nullptr)
341 , TrueEta_AngleBetweenGammas_after_cuts_onetwo_Outside(nullptr)
342 , WrongEta_AngleBetweenGammas_after_cuts_onetwo_Outside(nullptr)
343 , EMT_Eta_InvMass_onetwo_Outside(nullptr)
344 , fHistoList_Eta_two_Outside()
345 , AngleBetweenGammas_trueEta_before_cuts_two_Outside(nullptr)
346 , AngleBetweenGammas_wrongEta_before_cuts_two_Outside(nullptr)
347 , InvMass_GammaCombinations_two_Outside(nullptr)
348 , Particles_PDG_two_Outside(nullptr)
349 , TrueEta_pt_vs_rap_two_Outside(nullptr)
350 , WrongEta_pt_vs_rap_two_Outside(nullptr)
351 , TrueEta_pt_vs_rap_est_two_Outside(nullptr)
352 , WrongEta_pt_vs_rap_est_two_Outside(nullptr)
353 , TrueEta_InvMass_after_cuts_two_Outside(nullptr)
354 , WrongEta_InvMass_after_cuts_two_Outside(nullptr)
355 , TrueEta_AngleBetweenGammas_after_cuts_two_Outside(nullptr)
356 , WrongEta_AngleBetweenGammas_after_cuts_two_Outside(nullptr)
357 , EMT_Eta_InvMass_two_Outside(nullptr)
358 , fHistoList_rap_pt_Eta_all_Outside()
359 , fHistoList_rap_pt_Eta_onetwo_Outside()
360 , fHistoList_rap_pt_Eta_two_Outside()
361 , fHistoList_Eta_all_Both()
362 , AngleBetweenGammas_trueEta_before_cuts_all_Both(nullptr)
363 , AngleBetweenGammas_wrongEta_before_cuts_all_Both(nullptr)
364 , InvMass_GammaCombinations_all_Both(nullptr)
365 , Particles_PDG_all_Both(nullptr)
366 , TrueEta_pt_vs_rap_all_Both(nullptr)
367 , WrongEta_pt_vs_rap_all_Both(nullptr)
368 , TrueEta_pt_vs_rap_est_all_Both(nullptr)
369 , WrongEta_pt_vs_rap_est_all_Both(nullptr)
370 , TrueEta_InvMass_after_cuts_all_Both(nullptr)
371 , WrongEta_InvMass_after_cuts_all_Both(nullptr)
372 , TrueEta_AngleBetweenGammas_after_cuts_all_Both(nullptr)
373 , WrongEta_AngleBetweenGammas_after_cuts_all_Both(nullptr)
374 , EMT_Eta_InvMass_all_Both(nullptr)
375 , fHistoList_rap_pt_Eta_all_Both()
376 , multi_InvMass_Eta_all_Both_1(nullptr)
377 , multi_InvMass_Eta_all_Both_2(nullptr)
378 , multi_InvMass_Eta_all_Both_3(nullptr)
379 , multi_InvMass_Eta_all_Both_4(nullptr)
380 , multi_InvMass_Eta_all_Both_5(nullptr)
381 , multi_InvMass_Eta_all_Both_6(nullptr)
382 , multi_InvMass_Eta_all_Both_7(nullptr)
383 , multi_InvMass_Eta_all_Both_8(nullptr)
384 , multi_InvMass_Eta_all_Both_9(nullptr)
385 , multi_InvMass_Eta_all_Both_10(nullptr)
386 , multi_InvMass_Eta_all_Both_11(nullptr)
387 , multi_InvMass_Eta_all_Both_12(nullptr)
388 , multi_InvMass_Eta_all_Both_13(nullptr)
389 , multi_InvMass_Eta_all_Both_14(nullptr)
390 , multi_InvMass_Eta_all_Both_15(nullptr)
391 , multi_InvMass_Eta_all_Both_16(nullptr)
392 , multi_InvMass_Eta_all_Both_17(nullptr)
393 , multi_EMT_Eta_all_Both_1(nullptr)
394 , multi_EMT_Eta_all_Both_2(nullptr)
395 , multi_EMT_Eta_all_Both_3(nullptr)
396 , multi_EMT_Eta_all_Both_4(nullptr)
397 , multi_EMT_Eta_all_Both_5(nullptr)
398 , multi_EMT_Eta_all_Both_6(nullptr)
399 , multi_EMT_Eta_all_Both_7(nullptr)
400 , multi_EMT_Eta_all_Both_8(nullptr)
401 , multi_EMT_Eta_all_Both_9(nullptr)
402 , multi_EMT_Eta_all_Both_10(nullptr)
403 , multi_EMT_Eta_all_Both_11(nullptr)
404 , multi_EMT_Eta_all_Both_12(nullptr)
405 , multi_EMT_Eta_all_Both_13(nullptr)
406 , multi_EMT_Eta_all_Both_14(nullptr)
407 , multi_EMT_Eta_all_Both_15(nullptr)
408 , multi_EMT_Eta_all_Both_16(nullptr)
409 , multi_EMT_Eta_all_Both_17(nullptr)
410 , fHistoList_Eta_onetwo_Both()
411 , AngleBetweenGammas_trueEta_before_cuts_onetwo_Both(nullptr)
412 , AngleBetweenGammas_wrongEta_before_cuts_onetwo_Both(nullptr)
413 , InvMass_GammaCombinations_onetwo_Both(nullptr)
414 , Particles_PDG_onetwo_Both(nullptr)
415 , TrueEta_pt_vs_rap_onetwo_Both(nullptr)
416 , WrongEta_pt_vs_rap_onetwo_Both(nullptr)
417 , TrueEta_pt_vs_rap_est_onetwo_Both(nullptr)
418 , WrongEta_pt_vs_rap_est_onetwo_Both(nullptr)
419 , TrueEta_InvMass_after_cuts_onetwo_Both(nullptr)
420 , WrongEta_InvMass_after_cuts_onetwo_Both(nullptr)
421 , TrueEta_AngleBetweenGammas_after_cuts_onetwo_Both(nullptr)
422 , WrongEta_AngleBetweenGammas_after_cuts_onetwo_Both(nullptr)
423 , EMT_Eta_InvMass_onetwo_Both(nullptr)
424 , fHistoList_rap_pt_Eta_onetwo_Both()
425 , multi_InvMass_Eta_onetwo_Both_1(nullptr)
426 , multi_InvMass_Eta_onetwo_Both_2(nullptr)
427 , multi_InvMass_Eta_onetwo_Both_3(nullptr)
428 , multi_InvMass_Eta_onetwo_Both_4(nullptr)
429 , multi_InvMass_Eta_onetwo_Both_5(nullptr)
430 , multi_InvMass_Eta_onetwo_Both_6(nullptr)
431 , multi_InvMass_Eta_onetwo_Both_7(nullptr)
432 , multi_InvMass_Eta_onetwo_Both_8(nullptr)
433 , multi_InvMass_Eta_onetwo_Both_9(nullptr)
434 , multi_InvMass_Eta_onetwo_Both_10(nullptr)
435 , multi_InvMass_Eta_onetwo_Both_11(nullptr)
436 , multi_InvMass_Eta_onetwo_Both_12(nullptr)
437 , multi_InvMass_Eta_onetwo_Both_13(nullptr)
438 , multi_InvMass_Eta_onetwo_Both_14(nullptr)
439 , multi_InvMass_Eta_onetwo_Both_15(nullptr)
440 , multi_InvMass_Eta_onetwo_Both_16(nullptr)
441 , multi_InvMass_Eta_onetwo_Both_17(nullptr)
442 , multi_EMT_Eta_onetwo_Both_1(nullptr)
443 , multi_EMT_Eta_onetwo_Both_2(nullptr)
444 , multi_EMT_Eta_onetwo_Both_3(nullptr)
445 , multi_EMT_Eta_onetwo_Both_4(nullptr)
446 , multi_EMT_Eta_onetwo_Both_5(nullptr)
447 , multi_EMT_Eta_onetwo_Both_6(nullptr)
448 , multi_EMT_Eta_onetwo_Both_7(nullptr)
449 , multi_EMT_Eta_onetwo_Both_8(nullptr)
450 , multi_EMT_Eta_onetwo_Both_9(nullptr)
451 , multi_EMT_Eta_onetwo_Both_10(nullptr)
452 , multi_EMT_Eta_onetwo_Both_11(nullptr)
453 , multi_EMT_Eta_onetwo_Both_12(nullptr)
454 , multi_EMT_Eta_onetwo_Both_13(nullptr)
455 , multi_EMT_Eta_onetwo_Both_14(nullptr)
456 , multi_EMT_Eta_onetwo_Both_15(nullptr)
457 , multi_EMT_Eta_onetwo_Both_16(nullptr)
458 , multi_EMT_Eta_onetwo_Both_17(nullptr)
459 , fHistoList_Eta_two_Both()
460 , AngleBetweenGammas_trueEta_before_cuts_two_Both(nullptr)
461 , AngleBetweenGammas_wrongEta_before_cuts_two_Both(nullptr)
462 , InvMass_GammaCombinations_two_Both(nullptr)
463 , Particles_PDG_two_Both(nullptr)
464 , TrueEta_pt_vs_rap_two_Both(nullptr)
465 , WrongEta_pt_vs_rap_two_Both(nullptr)
466 , TrueEta_pt_vs_rap_est_two_Both(nullptr)
467 , WrongEta_pt_vs_rap_est_two_Both(nullptr)
468 , TrueEta_InvMass_after_cuts_two_Both(nullptr)
469 , WrongEta_InvMass_after_cuts_two_Both(nullptr)
470 , TrueEta_AngleBetweenGammas_after_cuts_two_Both(nullptr)
471 , WrongEta_AngleBetweenGammas_after_cuts_two_Both(nullptr)
472 , EMT_Eta_InvMass_two_Both(nullptr)
473 , fHistoList_rap_pt_Eta_two_Both()
474 , multi_InvMass_Eta_two_Both_1(nullptr)
475 , multi_InvMass_Eta_two_Both_2(nullptr)
476 , multi_InvMass_Eta_two_Both_3(nullptr)
477 , multi_InvMass_Eta_two_Both_4(nullptr)
478 , multi_InvMass_Eta_two_Both_5(nullptr)
479 , multi_InvMass_Eta_two_Both_6(nullptr)
480 , multi_InvMass_Eta_two_Both_7(nullptr)
481 , multi_InvMass_Eta_two_Both_8(nullptr)
482 , multi_InvMass_Eta_two_Both_9(nullptr)
483 , multi_InvMass_Eta_two_Both_10(nullptr)
484 , multi_InvMass_Eta_two_Both_11(nullptr)
485 , multi_InvMass_Eta_two_Both_12(nullptr)
486 , multi_InvMass_Eta_two_Both_13(nullptr)
487 , multi_InvMass_Eta_two_Both_14(nullptr)
488 , multi_InvMass_Eta_two_Both_15(nullptr)
489 , multi_InvMass_Eta_two_Both_16(nullptr)
490 , multi_InvMass_Eta_two_Both_17(nullptr)
491 , multi_EMT_Eta_two_Both_1(nullptr)
492 , multi_EMT_Eta_two_Both_2(nullptr)
493 , multi_EMT_Eta_two_Both_3(nullptr)
494 , multi_EMT_Eta_two_Both_4(nullptr)
495 , multi_EMT_Eta_two_Both_5(nullptr)
496 , multi_EMT_Eta_two_Both_6(nullptr)
497 , multi_EMT_Eta_two_Both_7(nullptr)
498 , multi_EMT_Eta_two_Both_8(nullptr)
499 , multi_EMT_Eta_two_Both_9(nullptr)
500 , multi_EMT_Eta_two_Both_10(nullptr)
501 , multi_EMT_Eta_two_Both_11(nullptr)
502 , multi_EMT_Eta_two_Both_12(nullptr)
503 , multi_EMT_Eta_two_Both_13(nullptr)
504 , multi_EMT_Eta_two_Both_14(nullptr)
505 , multi_EMT_Eta_two_Both_15(nullptr)
506 , multi_EMT_Eta_two_Both_16(nullptr)
507 , multi_EMT_Eta_two_Both_17(nullptr)
508{
509}
510
512
514{
515 FairRootManager* ioman = FairRootManager::Instance();
516 if (nullptr == ioman) { Fatal("CbmKresEta::Init", "RootManager not instantised!"); }
517
518 fMcTracks = (TClonesArray*) ioman->GetObject("MCTrack");
519 if (nullptr == fMcTracks) { Fatal("CbmKresEta::Init", "No MCTrack array!"); }
520
521 fPrimVertex = dynamic_cast<CbmVertex*>(ioman->GetObject("PrimaryVertex."));
522 if (nullptr == fPrimVertex) { fPrimVertex = dynamic_cast<CbmVertex*>(ioman->GetObject("PrimaryVertex")); }
523 if (nullptr == fPrimVertex) { LOG(fatal) << "CbmKresEta::Init No PrimaryVertex array!"; }
524
525 fGlobalTracks = (TClonesArray*) ioman->GetObject("GlobalTrack");
526 if (nullptr == fGlobalTracks) { Fatal("CbmKresEta::Init", "No GlobalTrack array!"); }
527
528 fStsTracks = (TClonesArray*) ioman->GetObject("StsTrack");
529 if (nullptr == fStsTracks) { Fatal("CbmKresEta::Init", "No StsTrack array!"); }
530
531 fStsTrackMatches = (TClonesArray*) ioman->GetObject("StsTrackMatch");
532 if (nullptr == fStsTrackMatches) { Fatal("CbmKresEta::Init", "No StsTrackMatch array!"); }
533
534 fRichProjections = (TClonesArray*) ioman->GetObject("RichProjection");
535 if (nullptr == fRichProjections) { Fatal("CbmKresEta::Init", "No RichProjection array!"); }
536
537 fRichRings = (TClonesArray*) ioman->GetObject("RichRing");
538 if (nullptr == fRichRings) { Fatal("CbmKresEta::Init", "No RichRing array!"); }
539
540 fRichRingMatches = (TClonesArray*) ioman->GetObject("RichRingMatch");
541 if (nullptr == fRichRingMatches) { Fatal("CbmKresEta::Init", "No RichRingMatch array!"); }
542
543 fRichHits = (TClonesArray*) ioman->GetObject("RichHit");
544 if (nullptr == fRichHits) { Fatal("CbmKresEta::Init", "No RichHit array!"); }
545
546 fArrayMvdHit = (TClonesArray*) ioman->GetObject("MvdHit");
547 if (nullptr == fArrayMvdHit) { Fatal("CbmKresEta::Init", "No MvdHit array!"); }
548
549 fArrayStsHit = (TClonesArray*) ioman->GetObject("StsHit");
550 if (nullptr == fArrayStsHit) { Fatal("CbmKresEta::Init", "No StsHit array!"); }
551
553
555}
556
557void CbmKresEta::Exec(int fEventNumEta, double OpeningAngleCut, double GammaInvMassCut, int RealPID)
558{
559 // cout << "CbmKresEta, event No. " << fEventNumEta << endl;
560
561 if (fPrimVertex != nullptr) { fKFVertex = CbmKFVertex(*fPrimVertex); }
562 else {
563 Fatal("CbmKresConversionManual::Exec", "No PrimaryVertex array!");
564 }
565
566 Gammas_all_Target.clear();
567 Gammas_two_Target.clear();
568 Gammas_onetwo_Target.clear();
575 Gammas_MC_all_Target.clear();
576 Gammas_MC_two_Target.clear();
578 Gammas_all_Outside.clear();
579 Gammas_two_Outside.clear();
580 Gammas_onetwo_Outside.clear();
587 Gammas_MC_all_Outside.clear();
588 Gammas_MC_two_Outside.clear();
590 Gammas_all_Both.clear();
591 Gammas_two_Both.clear();
592 Gammas_onetwo_Both.clear();
599 Gammas_MC_all_Both.clear();
600 Gammas_MC_two_Both.clear();
601 Gammas_MC_onetwo_Both.clear();
602
605 VRings_minus_Outside.clear();
611 VRings_plus_Outside.clear();
614 VMCIndex_plus_Outside.clear();
615
617 VMomenta_minus_Target.clear();
619 VRings_minus_Target.clear();
622 VMCIndex_minus_Target.clear();
623 VMCtracks_plus_Target.clear();
624 VStsTrack_plus_Target.clear();
625 VMomenta_plus_Target.clear();
626 VRings_plus_Target.clear();
627 VStsIndex_plus_Target.clear();
628 VRichRing_plus_Target.clear();
629 VMCIndex_plus_Target.clear();
630
631 Int_t ngTracks = fGlobalTracks->GetEntriesFast();
632 for (Int_t i = 0; i < ngTracks; i++) {
633 CbmGlobalTrack* gTrack = (CbmGlobalTrack*) fGlobalTracks->At(i);
634 if (nullptr == gTrack) continue;
635 Int_t stsInd = gTrack->GetStsTrackIndex();
636 Int_t richInd = gTrack->GetRichRingIndex();
637
638 if (stsInd < 0) continue;
639 CbmStsTrack* stsTrack = (CbmStsTrack*) fStsTracks->At(stsInd);
640 if (stsTrack == nullptr) continue;
641 CbmTrackMatchNew* stsMatch = (CbmTrackMatchNew*) fStsTrackMatches->At(stsInd);
642 if (stsMatch == nullptr) continue;
643 if (stsMatch->GetNofLinks() <= 0) continue;
644 Int_t stsMcTrackId = stsMatch->GetMatchedLink().GetIndex();
645 if (stsMcTrackId < 0) continue;
646 CbmMCTrack* mcTrack = (CbmMCTrack*) fMcTracks->At(stsMcTrackId);
647 if (mcTrack == nullptr) continue;
648
649 FairTrackParam* proj = (FairTrackParam*) fRichProjections->At(i);
650 if (richInd < 0 && proj->GetX() > -115 && proj->GetX() < 115
651 && ((proj->GetY() < -120 && proj->GetY() > -200) || (proj->GetY() > 120 && proj->GetY() < 200)))
652 continue;
653 CbmRichRing* Ring = nullptr;
654 if (richInd > -1) { Ring = static_cast<CbmRichRing*>(fRichRings->At(richInd)); }
655
656 // Doing fit with Fit To primary Vertex and calculate chi2 to primary vertex
657 double chi2 = 0;
659 fKFVertex.GetRefZ(), chi2);
660 const FairTrackParam* track_par = stsTrack->GetParamFirst();
661 double charge = track_par->GetQp();
662
663 if (chi2 != chi2) continue;
664 if (chi2 == 0) continue;
665
666 if (chi2 > 3) { SaveOutsideTracks(mcTrack, stsTrack, charge, stsInd, richInd, stsMcTrackId, Ring); }
667 if (chi2 > 3) continue;
668
669 SaveTargetTracks(mcTrack, stsTrack, Momentum, charge, stsInd, richInd, stsMcTrackId, Ring);
670 }
671
672 FindGammasTarget(fEventNumEta, OpeningAngleCut, GammaInvMassCut, RealPID, VMCtracks_minus_Target,
677
678 FindGammasOutside(fEventNumEta, OpeningAngleCut, GammaInvMassCut, RealPID, VMCtracks_minus_Outside,
682
684
691
698
705
706 int numformix = 500;
707 if (fEventNumEta % numformix == 0) {
721 EMT_eta_gg_Event_Both.clear();
724 }
725
726 if (fEventNumEta % numformix == 0) {
731 }
732
733 if (fEventNumEta % numformix == 0) {
738 }
739}
740
741
742void CbmKresEta::SaveOutsideTracks(CbmMCTrack* mcTrack1, CbmStsTrack* stsTrack, double charge, int stsInd, int richInd,
743 int stsMcTrackId, CbmRichRing* RING)
744{
745 int InRich = FindInRich(richInd, stsMcTrackId);
746 if (charge < 0) {
747 VMCtracks_minus_Outside.push_back(mcTrack1);
748 VStsTrack_minus_Outside.push_back(stsTrack);
749 VRings_minus_Outside.push_back(InRich);
750 VStsIndex_minus_Outside.push_back(stsInd);
751 VRichRing_minus_Outside.push_back(RING);
752 VMCIndex_minus_Outside.push_back(stsMcTrackId);
753 }
754 if (charge > 0) {
755 VMCtracks_plus_Outside.push_back(mcTrack1);
756 VStsTrack_plus_Outside.push_back(stsTrack);
757 VRings_plus_Outside.push_back(InRich);
758 VStsIndex_plus_Outside.push_back(stsInd);
759 VRichRing_plus_Outside.push_back(RING);
760 VMCIndex_plus_Outside.push_back(stsMcTrackId);
761 }
762}
763
764
765void CbmKresEta::SaveTargetTracks(CbmMCTrack* mcTrack1, CbmStsTrack* stsTrack, TVector3 refmom, double charge,
766 int stsInd, int richInd, int stsMcTrackId, CbmRichRing* RING)
767{
768 int InRich = FindInRich(richInd, stsMcTrackId);
769 if (charge < 0) {
770 VMCtracks_minus_Target.push_back(mcTrack1);
771 VStsTrack_minus_Target.push_back(stsTrack);
772 VMomenta_minus_Target.push_back(refmom);
773 VRings_minus_Target.push_back(InRich);
774 VStsIndex_minus_Target.push_back(stsInd);
775 VRichRing_minus_Target.push_back(RING);
776 VMCIndex_minus_Target.push_back(stsMcTrackId);
777 }
778 if (charge > 0) {
779 VMCtracks_plus_Target.push_back(mcTrack1);
780 VStsTrack_plus_Target.push_back(stsTrack);
781 VMomenta_plus_Target.push_back(refmom);
782 VRings_plus_Target.push_back(InRich);
783 VStsIndex_plus_Target.push_back(stsInd);
784 VRichRing_plus_Target.push_back(RING);
785 VMCIndex_plus_Target.push_back(stsMcTrackId);
786 }
787}
788
789int CbmKresEta::FindInRich(int richInd, int stsMcTrackId)
790{
791 int RingsInRich = 0;
792 if (richInd > -1) {
793 CbmTrackMatchNew* richMatch = (CbmTrackMatchNew*) fRichRingMatches->At(richInd);
794 if (richMatch != nullptr && richMatch->GetNofLinks() > 0) {
795 int richMcTrackId = richMatch->GetMatchedLink().GetIndex();
796 if (richMcTrackId > 0) {
797 if (stsMcTrackId == richMcTrackId) { // check that global track was matched correctly for STS and RICH together
798 CbmMCTrack* mcTrack2 = (CbmMCTrack*) fMcTracks->At(richMcTrackId);
799 if (mcTrack2 != nullptr) {
800 int pdgRICH = mcTrack2->GetPdgCode();
801 if (TMath::Abs(pdgRICH) == 11) RingsInRich++;
802 }
803 }
804 }
805 }
806 }
807 return RingsInRich;
808}
809
810int CbmKresEta::CheckIfElectron(CbmRichRing* ring, double momentum)
811{
812 int identified = 0;
813
814 if (nullptr != ring) {
815 CbmRichRingLight ringHit;
816 int nofHits = ring->GetNofHits();
817 for (int i = 0; i < nofHits; i++) {
818 int hitInd = ring->GetHit(i);
819 CbmRichHit* hit = (CbmRichHit*) fRichHits->At(hitInd);
820 if (nullptr == hit) continue;
821 CbmRichHitLight hl(hit->GetX(), hit->GetY());
822 ringHit.AddHit(hl);
823 }
824 fTauFit->DoFit(&ringHit);
825 if (ringHit.GetAaxis() > 4 && ringHit.GetAaxis() < 6 && ringHit.GetBaxis() > 4 && ringHit.GetBaxis() < 6
826 && momentum > 0.2 && momentum < 4.)
827 identified++;
828 }
829
830 return identified;
831}
832
834{
835 double FinalAngle = 400;
836 Int_t hits1sts = Sts_1->GetNofStsHits();
837 Int_t hits2sts = Sts_2->GetNofStsHits();
838 Int_t hits1mvd = Sts_1->GetNofMvdHits();
839 Int_t hits2mvd = Sts_2->GetNofMvdHits();
840
841 double Xpart1 = 0;
842 double Ypart1 = 0;
843 double Zpart1 = 0;
844 double Xpart2 = 0;
845 double Ypart2 = 0;
846 double Zpart2 = 0;
847 if (hits1sts > 0) {
848 Int_t stsHitIndex1 = Sts_1->GetStsHitIndex(hits1sts - 1);
849 CbmStsHit* stsHit1 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex1);
850 Xpart1 = stsHit1->GetX();
851 Ypart1 = stsHit1->GetY();
852 Zpart1 = stsHit1->GetZ();
853 }
854 else {
855 Int_t mvdHitIndex1 = Sts_1->GetMvdHitIndex(hits1mvd - 1);
856 CbmMvdHit* mvdHit1 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex1);
857 Xpart1 = mvdHit1->GetX();
858 Ypart1 = mvdHit1->GetY();
859 Zpart1 = mvdHit1->GetZ();
860 }
861
862 if (hits2sts > 0) {
863 Int_t stsHitIndex2 = Sts_2->GetStsHitIndex(hits2sts - 1);
864 CbmStsHit* stsHit2 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex2);
865 Xpart2 = stsHit2->GetX();
866 Ypart2 = stsHit2->GetY();
867 Zpart2 = stsHit2->GetZ();
868 }
869 else {
870 Int_t mvdHitIndex2 = Sts_2->GetMvdHitIndex(hits2mvd - 1);
871 CbmMvdHit* mvdHit2 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex2);
872 Xpart2 = mvdHit2->GetX();
873 Ypart2 = mvdHit2->GetY();
874 Zpart2 = mvdHit2->GetZ();
875 }
876
877 // check difference in 2 cm, because of two slices of every STS and Mvd stations
878 if (TMath::Abs(Zpart1 - Zpart2) > 2
879 && Zpart1 > Zpart2) { // if last hits are in different stations --> try to find the latest common station
880 for (int i = hits1sts - 2; i > -1; i--) { // start from second last station
881 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
882 Int_t stsHitIndex = Sts_1->GetStsHitIndex(i);
883 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
884 Xpart1 = stsHit->GetX();
885 Ypart1 = stsHit->GetY();
886 Zpart1 = stsHit->GetZ();
887 }
888 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 > Zpart2) {
889 for (int i = hits1mvd - 2; i > -1; i--) {
890 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
891 Int_t mvdHitIndex = Sts_1->GetMvdHitIndex(i);
892 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
893 Xpart1 = mvdHit->GetX();
894 Ypart1 = mvdHit->GetY();
895 Zpart1 = mvdHit->GetZ();
896 }
897 }
898 }
899
900 if (TMath::Abs(Zpart1 - Zpart2) > 2
901 && Zpart1 < Zpart2) { // if last hits are in different stations --> try to find the latest common station
902 for (int i = hits2sts - 2; i > -1; i--) { // start from second last station
903 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
904 Int_t stsHitIndex = Sts_2->GetStsHitIndex(i);
905 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
906 Xpart2 = stsHit->GetX();
907 Ypart2 = stsHit->GetY();
908 Zpart2 = stsHit->GetZ();
909 }
910 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 < Zpart2) {
911 for (int i = hits2mvd - 2; i > -1; i--) {
912 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
913 Int_t mvdHitIndex = Sts_2->GetMvdHitIndex(i);
914 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
915 Xpart2 = mvdHit->GetX();
916 Ypart2 = mvdHit->GetY();
917 Zpart2 = mvdHit->GetZ();
918 }
919 }
920 }
921
922 // calculate angle if we found common station
923 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
924 FinalAngle = TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 / TMath::Pi());
925 }
926
927 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
928}
929
931{
932 double FinalAngle = 400;
933 Int_t hits1sts = Sts_1->GetNofStsHits();
934 Int_t hits2sts = Sts_2->GetNofStsHits();
935 Int_t hits1mvd = Sts_1->GetNofMvdHits();
936 Int_t hits2mvd = Sts_2->GetNofMvdHits();
937
938 double Xpart1 = 0;
939 double Ypart1 = 0;
940 double Zpart1 = 0;
941 double Xpart2 = 0;
942 double Ypart2 = 0;
943 double Zpart2 = 0;
944 if (hits1mvd > 0) {
945 Int_t mvdHitIndex1 = Sts_1->GetMvdHitIndex(0);
946 CbmMvdHit* mvdHit1 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex1);
947 Xpart1 = mvdHit1->GetX();
948 Ypart1 = mvdHit1->GetY();
949 Zpart1 = mvdHit1->GetZ();
950 }
951 else {
952 Int_t stsHitIndex1 = Sts_1->GetStsHitIndex(0);
953 CbmStsHit* stsHit1 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex1);
954 Xpart1 = stsHit1->GetX();
955 Ypart1 = stsHit1->GetY();
956 Zpart1 = stsHit1->GetZ();
957 }
958
959 if (hits2mvd > 0) {
960 Int_t mvdHitIndex2 = Sts_2->GetMvdHitIndex(0);
961 CbmMvdHit* mvdHit2 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex2);
962 Xpart2 = mvdHit2->GetX();
963 Ypart2 = mvdHit2->GetY();
964 Zpart2 = mvdHit2->GetZ();
965 }
966 else {
967 Int_t stsHitIndex2 = Sts_2->GetStsHitIndex(0);
968 CbmStsHit* stsHit2 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex2);
969 Xpart2 = stsHit2->GetX();
970 Ypart2 = stsHit2->GetY();
971 Zpart2 = stsHit2->GetZ();
972 }
973
974 // check difference in 2 cm, because of two slices of every STS and Mvd stations
975 if (TMath::Abs(Zpart1 - Zpart2) > 2
976 && Zpart1 < Zpart2) { // if first hits are in different stations --> try to find the earliest common station
977 for (int i = 1; i < hits1mvd; i++) { // start from second hit
978 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
979 Int_t mvdHitIndex = Sts_1->GetMvdHitIndex(i);
980 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
981 Xpart1 = mvdHit->GetX();
982 Ypart1 = mvdHit->GetY();
983 Zpart1 = mvdHit->GetZ();
984 }
985 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 < Zpart2) {
986 for (int i = 0; i < hits1sts; i++) {
987 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
988 Int_t stsHitIndex = Sts_1->GetStsHitIndex(i);
989 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
990 Xpart1 = stsHit->GetX();
991 Ypart1 = stsHit->GetY();
992 Zpart1 = stsHit->GetZ();
993 }
994 }
995 }
996
997 if (TMath::Abs(Zpart1 - Zpart2) > 2
998 && Zpart1 > Zpart2) { // if first hits are in different stations --> try to find the earliest common station
999 for (int i = 1; i < hits2mvd; i++) { // start from second hit
1000 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1001 Int_t mvdHitIndex = Sts_2->GetMvdHitIndex(i);
1002 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
1003 Xpart2 = mvdHit->GetX();
1004 Ypart2 = mvdHit->GetY();
1005 Zpart2 = mvdHit->GetZ();
1006 }
1007 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 > Zpart2) {
1008 for (int i = 0; i < hits2sts; i++) {
1009 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1010 Int_t stsHitIndex = Sts_2->GetStsHitIndex(i);
1011 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
1012 Xpart2 = stsHit->GetX();
1013 Ypart2 = stsHit->GetY();
1014 Zpart2 = stsHit->GetZ();
1015 }
1016 }
1017 }
1018
1019 // calculate angle if we found common station
1020 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
1021 FinalAngle = TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 / TMath::Pi());
1022 }
1023
1024 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
1025}
1026
1027
1028void CbmKresEta::FindGammasTarget(int EventNumEta, double AngleCut, double InvMassCut, int RealPID,
1029 vector<CbmMCTrack*> MCtracks_minus, vector<CbmMCTrack*> MCtracks_plus,
1030 vector<CbmStsTrack*> StsTrack_minus, vector<CbmStsTrack*> StsTrack_plus,
1031 vector<TVector3> Momenta_minus, vector<TVector3> Momenta_plus,
1032 std::vector<int> Rings_minus, std::vector<int> Rings_plus,
1033 std::vector<int> stsIndex_minus, std::vector<int> stsIndex_plus,
1034 vector<CbmRichRing*> richRing_minus, vector<CbmRichRing*> richRing_plus,
1035 vector<Int_t> MCIndex_minus, vector<Int_t> MCIndex_plus)
1036{
1037 for (size_t i = 0; i < Momenta_minus.size(); i++) {
1038 for (size_t j = 0; j < Momenta_plus.size(); j++) {
1039
1040 CbmStsTrack* sts1 = StsTrack_minus.at(i);
1041 CbmStsTrack* sts2 = StsTrack_plus.at(j);
1042
1043 TVector3 part1 = Momenta_minus[i];
1044 TVector3 part2 = Momenta_plus[j];
1045 CbmMCTrack* part1MC = MCtracks_minus[i];
1046 CbmMCTrack* part2MC = MCtracks_plus[j];
1047
1048 int richcheck_0 = 0;
1049 int richcheck_1 = 0;
1050 if (RealPID == 1) {
1051 // Real PID
1052 richcheck_0 = CheckIfElectron(richRing_minus[i], Momenta_minus[i].Mag());
1053 richcheck_1 = CheckIfElectron(richRing_plus[j], Momenta_plus[j].Mag());
1054 }
1055 else {
1056 // MC PID
1057 richcheck_0 = Rings_minus[i];
1058 richcheck_1 = Rings_plus[j];
1059 }
1060 int richcheck = richcheck_0 + richcheck_1;
1061
1062 double InvmassReco = CbmKresFunctions::Invmass_2particles_RECO(part1, part2);
1063 double OpeningAngle = CbmKresFunctions::CalculateOpeningAngle_Reco(part1, part2);
1064 double PlaneAngle_last = CalculatePlaneAngle_last(sts1, sts2);
1065 double PlaneAngle_first = CalculatePlaneAngle_first(sts1, sts2);
1067
1068 int TruePair =
1069 0; //correctly reconstructed photon from eta meson: 0 - means wrong pair combination; 1 - means correct pair
1070 if (part1MC->GetMotherId() == part2MC->GetMotherId() && part1MC->GetMotherId() != -1) {
1071 CbmMCTrack* mcTrackmama = (CbmMCTrack*) fMcTracks->At(part1MC->GetMotherId());
1072 if (mcTrackmama->GetMotherId() != -1 && mcTrackmama->GetPdgCode() == 22) {
1073 CbmMCTrack* mcTrackgrmama = (CbmMCTrack*) fMcTracks->At(mcTrackmama->GetMotherId());
1074 if (mcTrackgrmama->GetPdgCode() == 221) { TruePair = 1; }
1075 }
1076 }
1077
1078 if (TruePair == 1) {
1079 if (richcheck == 0) {
1084 PlaneAngles_last_fromEta_zeroInRich_Target->Fill(PlaneAngle_last);
1085 PlaneAngles_first_fromEta_zeroInRich_Target->Fill(PlaneAngle_first);
1086
1091 PlaneAngles_first_fromEta_zeroInRich_Both->Fill(PlaneAngle_first);
1092 PlaneAngles_last_fromEta_zeroInRich_Both->Fill(PlaneAngle_last);
1093 }
1094 if (richcheck == 1) {
1099 PlaneAngles_last_fromEta_oneInRich_Target->Fill(PlaneAngle_last);
1100 PlaneAngles_first_fromEta_oneInRich_Target->Fill(PlaneAngle_first);
1101
1106 PlaneAngles_first_fromEta_oneInRich_Both->Fill(PlaneAngle_first);
1107 PlaneAngles_last_fromEta_oneInRich_Both->Fill(PlaneAngle_last);
1108 }
1109 if (richcheck == 2) {
1114 PlaneAngles_last_fromEta_twoInRich_Target->Fill(PlaneAngle_last);
1115 PlaneAngles_first_fromEta_twoInRich_Target->Fill(PlaneAngle_first);
1116
1121 PlaneAngles_first_fromEta_twoInRich_Both->Fill(PlaneAngle_first);
1122 PlaneAngles_last_fromEta_twoInRich_Both->Fill(PlaneAngle_last);
1123 }
1124 Particle_pull_px_Target->Fill(part1MC->GetPx() - part1.X());
1125 Particle_pull_py_Target->Fill(part1MC->GetPy() - part1.Y());
1126 Particle_pull_pz_Target->Fill(part1MC->GetPz() - part1.Z());
1127 Particle_pull_px_Target->Fill(part2MC->GetPx() - part2.X());
1128 Particle_pull_py_Target->Fill(part2MC->GetPy() - part2.Y());
1129 Particle_pull_pz_Target->Fill(part2MC->GetPz() - part2.Z());
1130 }
1131 else {
1132 if (richcheck == 0) {
1137 PlaneAngles_last_wrongpairs_zeroInRich_Target->Fill(PlaneAngle_last);
1138 PlaneAngles_first_wrongpairs_zeroInRich_Target->Fill(PlaneAngle_first);
1139
1144 PlaneAngles_first_wrongpairs_zeroInRich_Both->Fill(PlaneAngle_first);
1145 PlaneAngles_last_wrongpairs_zeroInRich_Both->Fill(PlaneAngle_last);
1146 }
1147 if (richcheck == 1) {
1152 PlaneAngles_last_wrongpairs_oneInRich_Target->Fill(PlaneAngle_last);
1153 PlaneAngles_first_wrongpairs_oneInRich_Target->Fill(PlaneAngle_first);
1154
1159 PlaneAngles_first_wrongpairs_oneInRich_Both->Fill(PlaneAngle_first);
1160 PlaneAngles_last_wrongpairs_oneInRich_Both->Fill(PlaneAngle_last);
1161 }
1162 if (richcheck == 2) {
1167 PlaneAngles_last_wrongpairs_twoInRich_Target->Fill(PlaneAngle_last);
1168 PlaneAngles_first_wrongpairs_twoInRich_Target->Fill(PlaneAngle_first);
1169
1174 PlaneAngles_first_wrongpairs_twoInRich_Both->Fill(PlaneAngle_first);
1175 PlaneAngles_last_wrongpairs_twoInRich_Both->Fill(PlaneAngle_last);
1176 }
1177 }
1178
1179
1180 // cuts for gamma reconstruction
1181 if (TMath::Abs(OpeningAngle) > AngleCut) continue;
1182 if (TMath::Abs(InvmassReco) > InvMassCut) continue;
1183
1184 if (richcheck == 1 && TruePair == 1) {
1187 }
1188 if (richcheck == 1 && TruePair == 0) {
1191 }
1192
1193 frefmomenta.clear();
1194 frefmomenta.push_back(part1);
1195 frefmomenta.push_back(part2);
1196 frefId.clear();
1197 frefId.push_back(stsIndex_minus[i]);
1198 frefId.push_back(stsIndex_plus[j]);
1199 fMCId.clear();
1200 fMCId.push_back(MCIndex_minus[i]);
1201 fMCId.push_back(MCIndex_plus[j]);
1202 fMCtracks.clear();
1203 fMCtracks.push_back(part1MC);
1204 fMCtracks.push_back(part2MC);
1205
1206 // for event mixing
1207 EMT_eta_gg_Event_Target.push_back(EventNumEta);
1209 EMT_eta_gg_NofRings_Target.push_back(richcheck);
1210
1211
1212 // everything (RICH == 0, RICH == 1, RICH == 2) together
1213 if (richcheck == 0 || richcheck == 1 || richcheck == 2) {
1214 Gammas_all_Target.push_back(frefmomenta);
1218 }
1219
1220 // only cases, when RICH == 2
1221 if (richcheck == 2) {
1222 Gammas_two_Target.push_back(frefmomenta);
1226 }
1227
1228 // cases, when RICH == 1 or RICH == 2 together
1229 if (richcheck == 1 || richcheck == 2) {
1234 }
1235 }
1236 }
1237}
1238
1239
1240void CbmKresEta::FindGammasOutside(int EventNumEta, double AngleCut, double InvMassCut, int RealPID,
1241 vector<CbmMCTrack*> MCtracks_minus_Outside,
1242 vector<CbmMCTrack*> MCtracks_plus_Outside,
1243 vector<CbmStsTrack*> StsTrack_minus_Outside,
1244 vector<CbmStsTrack*> StsTrack_plus_Outside, std::vector<int> Rings_minus_Outside,
1245 std::vector<int> Rings_plus_Outside, std::vector<int> stsIndex_minus_Outside,
1246 std::vector<int> stsIndex_plus_Outside, vector<CbmRichRing*> richRing_minus_Outside,
1247 vector<CbmRichRing*> richRing_plus_Outside, vector<Int_t> MCIndex_minus_Outside,
1248 vector<Int_t> MCIndex_plus_Outside)
1249{
1250 for (size_t i = 0; i < StsTrack_minus_Outside.size(); i++) {
1251 for (size_t j = 0; j < StsTrack_plus_Outside.size(); j++) {
1252
1253 CbmStsTrack* sts1 = StsTrack_minus_Outside[i];
1254 CbmStsTrack* sts2 = StsTrack_plus_Outside[j];
1255
1256 CbmMCTrack* part1MC = MCtracks_minus_Outside[i];
1257 CbmMCTrack* part2MC = MCtracks_plus_Outside[j];
1258
1259 KFParticle electron;
1261 KFParticle positron;
1263 const KFParticle* daughters[2] = {&electron, &positron};
1264 KFParticle intersection;
1265 intersection.Construct(daughters, 2);
1266
1267 if (intersection.GetZ() > 75 || intersection.GetZ() < -5) continue; // kick weird intersections
1268
1269 // fit to the vertex fitter
1270 // TVector3 part1 = CbmKresFunctions::FitToVertex(sts1, intersection.GetX(), intersection.GetY(), intersection.GetZ());
1271 // TVector3 part2 = CbmKresFunctions::FitToVertex(sts2, intersection.GetX(), intersection.GetY(), intersection.GetZ());
1272
1273 double chi_electron = 0;
1274 double chi_positron = 0;
1275 TVector3 part1 = CbmKresFunctions::FitToVertexAndGetChi(sts1, intersection.GetX(), intersection.GetY(),
1276 intersection.GetZ(), chi_electron);
1277 TVector3 part2 = CbmKresFunctions::FitToVertexAndGetChi(sts2, intersection.GetX(), intersection.GetY(),
1278 intersection.GetZ(), chi_positron);
1279 if (chi_electron != chi_electron) continue;
1280 if (chi_electron == 0) continue;
1281 if (chi_electron > 3) continue;
1282 if (chi_positron != chi_positron) continue;
1283 if (chi_positron == 0) continue;
1284 if (chi_positron > 3) continue;
1285
1286 int richcheck = 0;
1287 int richcheck_0 = 0;
1288 int richcheck_1 = 0;
1289 if (RealPID == 1) {
1290 // Real PID
1291 richcheck_0 = CheckIfElectron(richRing_minus_Outside[i], part1.Mag());
1292 richcheck_1 = CheckIfElectron(richRing_plus_Outside[j], part2.Mag());
1293 richcheck = richcheck_0 + richcheck_1;
1294 }
1295 else {
1296 // MC PID
1297 richcheck = Rings_minus_Outside[i] + Rings_plus_Outside[j];
1298 }
1299
1300 double InvmassReco = CbmKresFunctions::Invmass_2particles_RECO(part1, part2);
1301 double OpeningAngle = CbmKresFunctions::CalculateOpeningAngle_Reco(part1, part2);
1302 double PlaneAngle_last = CalculatePlaneAngle_last(sts1, sts2);
1303 double PlaneAngle_first = CalculatePlaneAngle_first(sts1, sts2);
1305
1306 int TruePair =
1307 0; //correctly reconstructed photon from eta meson: 0 - means wrong pair combination; 1 - means correct pair
1308 if (part1MC->GetMotherId() == part2MC->GetMotherId() && part1MC->GetMotherId() != -1) {
1309 CbmMCTrack* mcTrackmama = (CbmMCTrack*) fMcTracks->At(part1MC->GetMotherId());
1310 if (mcTrackmama->GetMotherId() != -1 && mcTrackmama->GetPdgCode() == 22) {
1311 CbmMCTrack* mcTrackgrmama = (CbmMCTrack*) fMcTracks->At(mcTrackmama->GetMotherId());
1312 if (mcTrackgrmama->GetPdgCode() == 221) { TruePair = 1; }
1313 }
1314 }
1315
1316 if (TruePair == 1) {
1317 if (richcheck == 0) {
1322 PlaneAngles_last_fromEta_zeroInRich_Outside->Fill(PlaneAngle_last);
1323 PlaneAngles_first_fromEta_zeroInRich_Outside->Fill(PlaneAngle_first);
1324
1329 PlaneAngles_first_fromEta_zeroInRich_Both->Fill(PlaneAngle_first);
1330 PlaneAngles_last_fromEta_zeroInRich_Both->Fill(PlaneAngle_last);
1331 }
1332 if (richcheck == 1) {
1337 PlaneAngles_last_fromEta_oneInRich_Outside->Fill(PlaneAngle_last);
1338 PlaneAngles_first_fromEta_oneInRich_Outside->Fill(PlaneAngle_first);
1339
1344 PlaneAngles_first_fromEta_oneInRich_Both->Fill(PlaneAngle_first);
1345 PlaneAngles_last_fromEta_oneInRich_Both->Fill(PlaneAngle_last);
1346 }
1347 if (richcheck == 2) {
1352 PlaneAngles_last_fromEta_twoInRich_Outside->Fill(PlaneAngle_last);
1353 PlaneAngles_first_fromEta_twoInRich_Outside->Fill(PlaneAngle_first);
1354
1359 PlaneAngles_first_fromEta_twoInRich_Both->Fill(PlaneAngle_first);
1360 PlaneAngles_last_fromEta_twoInRich_Both->Fill(PlaneAngle_last);
1361 }
1362 Particle_pull_px_Outside->Fill(part1MC->GetPx() - part1.X());
1363 Particle_pull_py_Outside->Fill(part1MC->GetPy() - part1.Y());
1364 Particle_pull_pz_Outside->Fill(part1MC->GetPz() - part1.Z());
1365 Particle_pull_px_Outside->Fill(part2MC->GetPx() - part2.X());
1366 Particle_pull_py_Outside->Fill(part2MC->GetPy() - part2.Y());
1367 Particle_pull_pz_Outside->Fill(part2MC->GetPz() - part2.Z());
1368 Particle_pull_X_Outside->Fill(part1MC->GetStartX() - intersection.GetX());
1369 Particle_pull_Y_Outside->Fill(part1MC->GetStartY() - intersection.GetY());
1370 Particle_pull_Z_Outside->Fill(part1MC->GetStartZ() - intersection.GetZ());
1371 Particle_pull_X_Outside->Fill(part2MC->GetStartX() - intersection.GetX());
1372 Particle_pull_Y_Outside->Fill(part2MC->GetStartY() - intersection.GetY());
1373 Particle_pull_Z_Outside->Fill(part2MC->GetStartZ() - intersection.GetZ());
1374 }
1375 else {
1376 if (richcheck == 0) {
1383
1388 PlaneAngles_first_wrongpairs_zeroInRich_Both->Fill(PlaneAngle_first);
1389 PlaneAngles_last_wrongpairs_zeroInRich_Both->Fill(PlaneAngle_last);
1390 }
1391 if (richcheck == 1) {
1396 PlaneAngles_last_wrongpairs_oneInRich_Outside->Fill(PlaneAngle_last);
1397 PlaneAngles_first_wrongpairs_oneInRich_Outside->Fill(PlaneAngle_first);
1398
1403 PlaneAngles_first_wrongpairs_oneInRich_Both->Fill(PlaneAngle_first);
1404 PlaneAngles_last_wrongpairs_oneInRich_Both->Fill(PlaneAngle_last);
1405 }
1406 if (richcheck == 2) {
1411 PlaneAngles_last_wrongpairs_twoInRich_Outside->Fill(PlaneAngle_last);
1412 PlaneAngles_first_wrongpairs_twoInRich_Outside->Fill(PlaneAngle_first);
1413
1418 PlaneAngles_first_wrongpairs_twoInRich_Both->Fill(PlaneAngle_first);
1419 PlaneAngles_last_wrongpairs_twoInRich_Both->Fill(PlaneAngle_last);
1420 }
1421 }
1422
1423
1424 // cuts for gamma reconstruction
1425 if (TMath::Abs(OpeningAngle) > AngleCut) continue;
1426 if (TMath::Abs(InvmassReco) > InvMassCut) continue;
1427
1428
1429 if (richcheck == 1 && TruePair == 1) {
1432 }
1433 if (richcheck == 1 && TruePair == 0) {
1436 }
1437
1438
1439 frefmomenta.clear();
1440 frefmomenta.push_back(part1);
1441 frefmomenta.push_back(part2);
1442 frefId.clear();
1443 frefId.push_back(stsIndex_minus_Outside[i]);
1444 frefId.push_back(stsIndex_plus_Outside[j]);
1445 fMCId.clear();
1446 fMCId.push_back(MCIndex_minus_Outside[i]);
1447 fMCId.push_back(MCIndex_plus_Outside[j]);
1448 fMCtracks.clear();
1449 fMCtracks.push_back(part1MC);
1450 fMCtracks.push_back(part2MC);
1451
1452 // for event mixing
1453 EMT_eta_gg_Event_Outside.push_back(EventNumEta);
1455 EMT_eta_gg_NofRings_Outside.push_back(richcheck);
1456
1457
1458 // everything (RICH == 0, RICH == 1, RICH == 2) together
1459 if (richcheck == 0 || richcheck == 1 || richcheck == 2) {
1464 }
1465
1466 // only cases, when RICH == 2
1467 if (richcheck == 2) {
1472 }
1473
1474 // cases, when RICH == 1 or RICH == 2 together
1475 if (richcheck == 1 || richcheck == 2) {
1480 }
1481 }
1482 }
1483}
1484
1486{
1493
1506
1519
1528}
1529
1530
1531void CbmKresEta::FindEta(TString /*mod*/, TString position, vector<vector<TVector3>> Gammas,
1532 vector<vector<int>> StsIndex, vector<vector<int>> MCIndex,
1533 vector<vector<CbmMCTrack*>> GammasMC, vector<TH1*> gg, vector<TH1*> rap_pt_separation)
1534{
1535 if (Gammas.size() < 2) return; // min 2 gammas to form pi0 are required
1536 for (size_t gamma1 = 0; gamma1 < Gammas.size() - 1; gamma1++) {
1537 for (size_t gamma2 = gamma1 + 1; gamma2 < Gammas.size(); gamma2++) {
1538 // 4 reconstructed particles from gammas
1539 TVector3 e11 = Gammas[gamma1][0];
1540 TVector3 e12 = Gammas[gamma1][1];
1541 TVector3 e21 = Gammas[gamma2][0];
1542 TVector3 e22 = Gammas[gamma2][1];
1543
1544 // MC true data for this particles
1545 CbmMCTrack* mcTrack1 = GammasMC[gamma1][0];
1546 CbmMCTrack* mcTrack2 = GammasMC[gamma1][1];
1547 CbmMCTrack* mcTrack3 = GammasMC[gamma2][0];
1548 CbmMCTrack* mcTrack4 = GammasMC[gamma2][1];
1549
1550 if (StsIndex[gamma1][0] == StsIndex[gamma2][0] || StsIndex[gamma1][0] == StsIndex[gamma2][1]
1551 || StsIndex[gamma1][1] == StsIndex[gamma2][0] || StsIndex[gamma1][1] == StsIndex[gamma2][1])
1552 continue; // particles are not used twice --> different
1553
1555 double openingAngleBetweenGammasReco =
1557
1558 int TrueEta = 0;
1559 int STSmcId1 = MCIndex[gamma1][0];
1560 int STSmcId2 = MCIndex[gamma1][1];
1561 int STSmcId3 = MCIndex[gamma2][0];
1562 int STSmcId4 = MCIndex[gamma2][1];
1563 if (STSmcId1 != STSmcId2 && STSmcId1 != STSmcId3 && STSmcId1 != STSmcId4 && STSmcId2 != STSmcId3
1564 && STSmcId2 != STSmcId4 && STSmcId3 != STSmcId4) {
1565 if (TMath::Abs(mcTrack1->GetPdgCode()) == 11 && TMath::Abs(mcTrack2->GetPdgCode()) == 11
1566 && TMath::Abs(mcTrack3->GetPdgCode()) == 11 && TMath::Abs(mcTrack4->GetPdgCode()) == 11) {
1567 if ((mcTrack1->GetMotherId() == mcTrack2->GetMotherId() && mcTrack3->GetMotherId() == mcTrack4->GetMotherId())
1568 || (mcTrack1->GetMotherId() == mcTrack3->GetMotherId()
1569 && mcTrack2->GetMotherId() == mcTrack4->GetMotherId())
1570 || (mcTrack1->GetMotherId() == mcTrack4->GetMotherId()
1571 && mcTrack2->GetMotherId() == mcTrack3->GetMotherId())) {
1572 CbmMCTrack* MothTrack1 = (CbmMCTrack*) fMcTracks->At(mcTrack1->GetMotherId());
1573 CbmMCTrack* MothTrack2 = (CbmMCTrack*) fMcTracks->At(mcTrack2->GetMotherId());
1574 CbmMCTrack* MothTrack3 = (CbmMCTrack*) fMcTracks->At(mcTrack3->GetMotherId());
1575 CbmMCTrack* MothTrack4 = (CbmMCTrack*) fMcTracks->At(mcTrack4->GetMotherId());
1576 if (MothTrack1->GetPdgCode() == 22 && MothTrack2->GetPdgCode() == 22 && MothTrack3->GetPdgCode() == 22
1577 && MothTrack4->GetPdgCode() == 22 && MothTrack1->GetMotherId() != -1
1578 && MothTrack1->GetMotherId() == MothTrack2->GetMotherId()
1579 && MothTrack1->GetMotherId() == MothTrack3->GetMotherId()
1580 && MothTrack1->GetMotherId() == MothTrack4->GetMotherId()) {
1581 CbmMCTrack* GrandMothTrack1 = (CbmMCTrack*) fMcTracks->At(MothTrack1->GetMotherId());
1582 if (GrandMothTrack1->GetPdgCode() == 221) {
1583 TrueEta = 1;
1584 // cout << " " << mod << "; " << position << " ***** CbmKresEta ***** " << endl;
1585 // cout << "Decay eta -> gamma gamma -> e+e- e+e- detected!\t\t reco mass: " << params.fMinv << endl;
1586 // cout << "motherids: " << mcTrack1->GetMotherId() << "/" << mcTrack2->GetMotherId() << "/" << mcTrack3->GetMotherId() << "/" << mcTrack4->GetMotherId() << endl;
1587 // cout << "grandmotherid: " << MothTrack1->GetMotherId() << "/" << MothTrack2->GetMotherId() << "/" << MothTrack3->GetMotherId() << "/" << MothTrack4->GetMotherId() << endl;
1588 // cout << "pdgs: " << GrandMothTrack1->GetPdgCode() << "-->" << MothTrack1->GetPdgCode() << "/" << MothTrack2->GetPdgCode() << "/" << MothTrack3->GetPdgCode() << "/" << MothTrack4->GetPdgCode() << endl;
1589 // cout << "particle 1: \t" << mcTrack1->GetPdgCode() << ";\t pt = " << mcTrack1->GetPt() << ";\t X = " << mcTrack1->GetStartX() << ";\t Y = " << mcTrack1->GetStartY() << ";\t Z = " << mcTrack1->GetStartZ() << ";\t E = " << mcTrack1->GetEnergy() << endl;
1590 // cout << "particle 2: \t" << mcTrack2->GetPdgCode() << ";\t pt = " << mcTrack2->GetPt() << ";\t X = " << mcTrack2->GetStartX() << ";\t Y = " << mcTrack2->GetStartY() << ";\t Z = " << mcTrack2->GetStartZ() << ";\t E = " << mcTrack2->GetEnergy() << endl;
1591 // cout << "particle 3: \t" << mcTrack3->GetPdgCode() << ";\t pt = " << mcTrack3->GetPt() << ";\t X = " << mcTrack3->GetStartX() << ";\t Y = " << mcTrack3->GetStartY() << ";\t Z = " << mcTrack3->GetStartZ() << ";\t E = " << mcTrack3->GetEnergy() << endl;
1592 // cout << "particle 4: \t" << mcTrack4->GetPdgCode() << ";\t pt = " << mcTrack4->GetPt() << ";\t X = " << mcTrack4->GetStartX() << ";\t Y = " << mcTrack4->GetStartY() << ";\t Z = " << mcTrack4->GetStartZ() << ";\t E = " << mcTrack4->GetEnergy() << endl;
1593 // cout << " ***** CbmKresEta (End) ***** " << endl;
1594 }
1595 }
1596 }
1597 }
1598 }
1599
1600 if (TrueEta == 1) gg[0]->Fill(openingAngleBetweenGammasReco);
1601 if (TrueEta == 0) gg[1]->Fill(openingAngleBetweenGammasReco);
1602
1603 if (openingAngleBetweenGammasReco < 10 || openingAngleBetweenGammasReco > 40) continue;
1604
1605 gg[2]->Fill(params.fMinv);
1606
1607 gg[3]->Fill(TMath::Abs(mcTrack1->GetPdgCode()));
1608 gg[3]->Fill(TMath::Abs(mcTrack2->GetPdgCode()));
1609 gg[3]->Fill(TMath::Abs(mcTrack3->GetPdgCode()));
1610 gg[3]->Fill(TMath::Abs(mcTrack4->GetPdgCode()));
1611
1612 // separate by rap and Pt only for case "Both"
1613 if (position == "Both") {
1614 if (params.fRapidity > 1.2 && params.fRapidity <= 1.6) {
1615 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation[1]->Fill(params.fMinv);
1616 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation[2]->Fill(params.fMinv);
1617 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation[3]->Fill(params.fMinv);
1618 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_pt_separation[4]->Fill(params.fMinv);
1619 }
1620 if (params.fRapidity > 1.6 && params.fRapidity <= 2.0) {
1621 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation[5]->Fill(params.fMinv);
1622 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation[6]->Fill(params.fMinv);
1623 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation[7]->Fill(params.fMinv);
1624 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_pt_separation[8]->Fill(params.fMinv);
1625 }
1626 if (params.fRapidity > 2.0 && params.fRapidity <= 2.4) {
1627 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation[9]->Fill(params.fMinv);
1628 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation[10]->Fill(params.fMinv);
1629 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation[11]->Fill(params.fMinv);
1630 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_pt_separation[12]->Fill(params.fMinv);
1631 }
1632 if (params.fRapidity > 2.4 && params.fRapidity <= 2.8) {
1633 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation[13]->Fill(params.fMinv);
1634 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation[14]->Fill(params.fMinv);
1635 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation[15]->Fill(params.fMinv);
1636 }
1637 if (params.fRapidity > 2.8 && params.fRapidity <= 3.2) {
1638 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation[16]->Fill(params.fMinv);
1639 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation[17]->Fill(params.fMinv);
1640 }
1641 }
1642
1643 if (TrueEta == 1) {
1644 gg[4]->Fill(params.fRapidity, params.fPt);
1645 gg[6]->Fill(params.fRapidity, params.fPt);
1646 gg[8]->Fill(params.fMinv);
1647 gg[10]->Fill(openingAngleBetweenGammasReco);
1648 }
1649
1650 if (TrueEta == 0) {
1651 gg[5]->Fill(params.fRapidity, params.fPt);
1652 gg[7]->Fill(params.fRapidity, params.fPt);
1653 gg[9]->Fill(params.fMinv);
1654 gg[11]->Fill(openingAngleBetweenGammasReco);
1655 }
1656 }
1657 }
1658}
1659
1661// TARGET
1662{
1663 int nof_Target = EMT_eta_gg_Event_Target.size();
1664 cout << "Mixing for Eta(Target) - nof entries " << nof_Target << endl;
1665 for (Int_t a = 0; a < nof_Target - 1; a++) {
1666 for (Int_t b = a + 1; b < nof_Target; b++) {
1668 continue; // to make sure that the photons are from two different events
1669 TVector3 e11 = EMT_eta_gg_pair_momenta_Target[a][0];
1670 TVector3 e12 = EMT_eta_gg_pair_momenta_Target[a][1];
1671 TVector3 e21 = EMT_eta_gg_pair_momenta_Target[b][0];
1672 TVector3 e22 = EMT_eta_gg_pair_momenta_Target[b][1];
1674 double openingAngleBetweenGammasReco =
1676 if (openingAngleBetweenGammasReco < 10 || openingAngleBetweenGammasReco > 40) continue;
1677
1678 EMT_Eta_InvMass_all_Target->Fill(params.fMinv);
1680 EMT_Eta_InvMass_two_Target->Fill(params.fMinv);
1684 }
1685 }
1686}
1687
1688
1690// OUTSIDE
1691{
1692 Int_t nof_Outside = EMT_eta_gg_Event_Outside.size();
1693 cout << "Mixing for Eta(Outside) - nof entries " << nof_Outside << endl;
1694 for (Int_t a = 0; a < nof_Outside - 1; a++) {
1695 for (Int_t b = a + 1; b < nof_Outside; b++) {
1697 continue; // to make sure that the photons are from two different events
1698 TVector3 e11 = EMT_eta_gg_pair_momenta_Outside[a][0];
1699 TVector3 e12 = EMT_eta_gg_pair_momenta_Outside[a][1];
1700 TVector3 e21 = EMT_eta_gg_pair_momenta_Outside[b][0];
1701 TVector3 e22 = EMT_eta_gg_pair_momenta_Outside[b][1];
1703 double openingAngleBetweenGammasReco =
1705 if (openingAngleBetweenGammasReco < 10 || openingAngleBetweenGammasReco > 40) continue;
1706
1707 EMT_Eta_InvMass_all_Outside->Fill(params.fMinv);
1709 EMT_Eta_InvMass_two_Outside->Fill(params.fMinv);
1713 }
1714 }
1715}
1716
1717void CbmKresEta::Mixing_Both(vector<TH1*> rap_pt_separation_all, vector<TH1*> rap_pt_separation_onetwo,
1718 vector<TH1*> rap_pt_separation_two)
1719// BOTH
1720{
1721 Int_t nof_Both = EMT_eta_gg_Event_Both.size();
1722 cout << "Mixing for Eta(Both) - nof entries " << nof_Both << endl;
1723 for (Int_t a = 0; a < nof_Both - 1; a++) {
1724 for (Int_t b = a + 1; b < nof_Both; b++) {
1726 continue; // to make sure that the photons are from two different events
1727 TVector3 e11 = EMT_eta_gg_pair_momenta_Both[a][0];
1728 TVector3 e12 = EMT_eta_gg_pair_momenta_Both[a][1];
1729 TVector3 e21 = EMT_eta_gg_pair_momenta_Both[b][0];
1730 TVector3 e22 = EMT_eta_gg_pair_momenta_Both[b][1];
1732 double openingAngleBetweenGammasReco =
1734 if (openingAngleBetweenGammasReco < 10 || openingAngleBetweenGammasReco > 40) continue;
1735
1736 EMT_Eta_InvMass_all_Both->Fill(params.fMinv);
1738 EMT_Eta_InvMass_two_Both->Fill(params.fMinv);
1739 if ((EMT_eta_gg_NofRings_Both[a] == 1 || EMT_eta_gg_NofRings_Both[a] == 2)
1740 && (EMT_eta_gg_NofRings_Both[b] == 1 || EMT_eta_gg_NofRings_Both[b] == 2))
1741 EMT_Eta_InvMass_onetwo_Both->Fill(params.fMinv);
1742
1743 // separate by rap and Pt only for all
1744 if (params.fRapidity > 1.2 && params.fRapidity <= 1.6) {
1745 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_all[18]->Fill(params.fMinv);
1746 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_all[19]->Fill(params.fMinv);
1747 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation_all[20]->Fill(params.fMinv);
1748 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_pt_separation_all[21]->Fill(params.fMinv);
1749 }
1750 if (params.fRapidity > 1.6 && params.fRapidity <= 2.0) {
1751 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_all[22]->Fill(params.fMinv);
1752 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_all[23]->Fill(params.fMinv);
1753 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation_all[24]->Fill(params.fMinv);
1754 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_pt_separation_all[25]->Fill(params.fMinv);
1755 }
1756 if (params.fRapidity > 2.0 && params.fRapidity <= 2.4) {
1757 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_all[26]->Fill(params.fMinv);
1758 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_all[27]->Fill(params.fMinv);
1759 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation_all[28]->Fill(params.fMinv);
1760 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_pt_separation_all[29]->Fill(params.fMinv);
1761 }
1762 if (params.fRapidity > 2.4 && params.fRapidity <= 2.8) {
1763 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_all[30]->Fill(params.fMinv);
1764 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_all[31]->Fill(params.fMinv);
1765 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation_all[32]->Fill(params.fMinv);
1766 }
1767 if (params.fRapidity > 2.8 && params.fRapidity <= 3.2) {
1768 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_all[33]->Fill(params.fMinv);
1769 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_all[34]->Fill(params.fMinv);
1770 }
1771
1772 // separate by rap and Pt only for onetwo
1773 if ((EMT_eta_gg_NofRings_Both[a] == 1 || EMT_eta_gg_NofRings_Both[a] == 2)
1774 && (EMT_eta_gg_NofRings_Both[b] == 1 || EMT_eta_gg_NofRings_Both[b] == 2)) {
1775 if (params.fRapidity > 1.2 && params.fRapidity <= 1.6) {
1776 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_onetwo[18]->Fill(params.fMinv);
1777 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_onetwo[19]->Fill(params.fMinv);
1778 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation_onetwo[20]->Fill(params.fMinv);
1779 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_pt_separation_onetwo[21]->Fill(params.fMinv);
1780 }
1781 if (params.fRapidity > 1.6 && params.fRapidity <= 2.0) {
1782 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_onetwo[22]->Fill(params.fMinv);
1783 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_onetwo[23]->Fill(params.fMinv);
1784 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation_onetwo[24]->Fill(params.fMinv);
1785 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_pt_separation_onetwo[25]->Fill(params.fMinv);
1786 }
1787 if (params.fRapidity > 2.0 && params.fRapidity <= 2.4) {
1788 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_onetwo[26]->Fill(params.fMinv);
1789 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_onetwo[27]->Fill(params.fMinv);
1790 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation_onetwo[28]->Fill(params.fMinv);
1791 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_pt_separation_onetwo[29]->Fill(params.fMinv);
1792 }
1793 if (params.fRapidity > 2.4 && params.fRapidity <= 2.8) {
1794 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_onetwo[30]->Fill(params.fMinv);
1795 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_onetwo[31]->Fill(params.fMinv);
1796 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation_onetwo[32]->Fill(params.fMinv);
1797 }
1798 if (params.fRapidity > 2.8 && params.fRapidity <= 3.2) {
1799 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_onetwo[33]->Fill(params.fMinv);
1800 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_onetwo[34]->Fill(params.fMinv);
1801 }
1802 }
1803
1804 // separate by rap and Pt only for two
1805 if (EMT_eta_gg_NofRings_Both[a] == 2 && EMT_eta_gg_NofRings_Both[b] == 2) {
1806 if (params.fRapidity > 1.2 && params.fRapidity <= 1.6) {
1807 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_two[18]->Fill(params.fMinv);
1808 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_two[19]->Fill(params.fMinv);
1809 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation_two[20]->Fill(params.fMinv);
1810 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_pt_separation_two[21]->Fill(params.fMinv);
1811 }
1812 if (params.fRapidity > 1.6 && params.fRapidity <= 2.0) {
1813 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_two[22]->Fill(params.fMinv);
1814 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_two[23]->Fill(params.fMinv);
1815 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation_two[24]->Fill(params.fMinv);
1816 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_pt_separation_two[25]->Fill(params.fMinv);
1817 }
1818 if (params.fRapidity > 2.0 && params.fRapidity <= 2.4) {
1819 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_two[26]->Fill(params.fMinv);
1820 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_two[27]->Fill(params.fMinv);
1821 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation_two[28]->Fill(params.fMinv);
1822 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_pt_separation_two[29]->Fill(params.fMinv);
1823 }
1824 if (params.fRapidity > 2.4 && params.fRapidity <= 2.8) {
1825 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_two[30]->Fill(params.fMinv);
1826 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_two[31]->Fill(params.fMinv);
1827 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_pt_separation_two[32]->Fill(params.fMinv);
1828 }
1829 if (params.fRapidity > 2.8 && params.fRapidity <= 3.2) {
1830 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_pt_separation_two[33]->Fill(params.fMinv);
1831 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_pt_separation_two[34]->Fill(params.fMinv);
1832 }
1833 }
1834 }
1835 }
1836}
1837
1838
1840{
1841 gDirectory->mkdir("Eta");
1842 gDirectory->cd("Eta");
1843
1844 gDirectory->mkdir("Cuts_for_Eta_gg");
1845 gDirectory->cd("Cuts_for_Eta_gg");
1846
1847 gDirectory->mkdir("Cuts_for_Eta_gg_Target");
1848 gDirectory->cd("Cuts_for_Eta_gg_Target");
1849 gDirectory->mkdir("Cuts_for_Eta_gg_zeroInRich_Target");
1850 gDirectory->cd("Cuts_for_Eta_gg_zeroInRich_Target");
1851 for (UInt_t i = 0; i < fHistoList_Eta_cuts_zeroInRich_Target.size(); i++) {
1853 }
1854 gDirectory->cd("..");
1855 gDirectory->mkdir("Cuts_for_Eta_gg_oneInRich_Target");
1856 gDirectory->cd("Cuts_for_Eta_gg_oneInRich_Target");
1857 for (UInt_t i = 0; i < fHistoList_Eta_cuts_oneInRich_Target.size(); i++) {
1859 }
1860 gDirectory->cd("..");
1861 gDirectory->mkdir("Cuts_for_Eta_gg_twoInRich_Target");
1862 gDirectory->cd("Cuts_for_Eta_gg_twoInRich_Target");
1863 for (UInt_t i = 0; i < fHistoList_Eta_cuts_twoInRich_Target.size(); i++) {
1865 }
1866 gDirectory->cd("..");
1867 for (UInt_t i = 0; i < fHistoList_Eta_cuts_Target.size(); i++) {
1868 fHistoList_Eta_cuts_Target[i]->Write();
1869 }
1870 gDirectory->cd("..");
1871
1872 gDirectory->mkdir("Cuts_for_Eta_gg_Outside");
1873 gDirectory->cd("Cuts_for_Eta_gg_Outside");
1874 gDirectory->mkdir("Cuts_for_Eta_gg_zeroInRich_Outside");
1875 gDirectory->cd("Cuts_for_Eta_gg_zeroInRich_Outside");
1876 for (UInt_t i = 0; i < fHistoList_Eta_cuts_zeroInRich_Outside.size(); i++) {
1878 }
1879 gDirectory->cd("..");
1880 gDirectory->mkdir("Cuts_for_Eta_gg_oneInRich_Outside");
1881 gDirectory->cd("Cuts_for_Eta_gg_oneInRich_Outside");
1882 for (UInt_t i = 0; i < fHistoList_Eta_cuts_oneInRich_Outside.size(); i++) {
1884 }
1885 gDirectory->cd("..");
1886 gDirectory->mkdir("Cuts_for_Eta_gg_twoInRich_Outside");
1887 gDirectory->cd("Cuts_for_Eta_gg_twoInRich_Outside");
1888 for (UInt_t i = 0; i < fHistoList_Eta_cuts_twoInRich_Outside.size(); i++) {
1890 }
1891 gDirectory->cd("..");
1892 for (UInt_t i = 0; i < fHistoList_Eta_cuts_Outside.size(); i++) {
1893 fHistoList_Eta_cuts_Outside[i]->Write();
1894 }
1895 gDirectory->cd("..");
1896
1897 gDirectory->mkdir("Cuts_for_Eta_gg_Both");
1898 gDirectory->cd("Cuts_for_Eta_gg_Both");
1899 gDirectory->mkdir("Cuts_for_Eta_gg_zeroInRich_Both");
1900 gDirectory->cd("Cuts_for_Eta_gg_zeroInRich_Both");
1901 for (UInt_t i = 0; i < fHistoList_Eta_cuts_zeroInRich_Both.size(); i++) {
1903 }
1904 gDirectory->cd("..");
1905 gDirectory->mkdir("Cuts_for_Eta_gg_oneInRich_Both");
1906 gDirectory->cd("Cuts_for_Eta_gg_oneInRich_Both");
1907 for (UInt_t i = 0; i < fHistoList_Eta_cuts_oneInRich_Both.size(); i++) {
1909 }
1910 gDirectory->cd("..");
1911 gDirectory->mkdir("Cuts_for_Eta_gg_twoInRich_Both");
1912 gDirectory->cd("Cuts_for_Eta_gg_twoInRich_Both");
1913 for (UInt_t i = 0; i < fHistoList_Eta_cuts_twoInRich_Both.size(); i++) {
1915 }
1916 gDirectory->cd("..");
1917 for (UInt_t i = 0; i < fHistoList_Eta_cuts_Both.size(); i++) {
1918 fHistoList_Eta_cuts_Both[i]->Write();
1919 }
1920 gDirectory->cd("..");
1921
1922 gDirectory->cd("..");
1923
1924
1925 gDirectory->mkdir("Eta_gg_Target");
1926 gDirectory->cd("Eta_gg_Target");
1927
1928 gDirectory->mkdir("all");
1929 gDirectory->cd("all");
1930 for (UInt_t i = 0; i < fHistoList_Eta_all_Target.size(); i++) {
1931 fHistoList_Eta_all_Target[i]->Write();
1932 }
1933 gDirectory->cd("..");
1934
1935 gDirectory->mkdir("onetwo");
1936 gDirectory->cd("onetwo");
1937 for (UInt_t i = 0; i < fHistoList_Eta_onetwo_Target.size(); i++) {
1938 fHistoList_Eta_onetwo_Target[i]->Write();
1939 }
1940 gDirectory->cd("..");
1941
1942 gDirectory->mkdir("two");
1943 gDirectory->cd("two");
1944 for (UInt_t i = 0; i < fHistoList_Eta_two_Target.size(); i++) {
1945 fHistoList_Eta_two_Target[i]->Write();
1946 }
1947 gDirectory->cd("..");
1948
1949 gDirectory->cd("..");
1950
1951
1952 gDirectory->mkdir("Eta_gg_Outside");
1953 gDirectory->cd("Eta_gg_Outside");
1954
1955 gDirectory->mkdir("all");
1956 gDirectory->cd("all");
1957 for (UInt_t i = 0; i < fHistoList_Eta_all_Outside.size(); i++) {
1958 fHistoList_Eta_all_Outside[i]->Write();
1959 }
1960 gDirectory->cd("..");
1961
1962 gDirectory->mkdir("onetwo");
1963 gDirectory->cd("onetwo");
1964 for (UInt_t i = 0; i < fHistoList_Eta_onetwo_Outside.size(); i++) {
1966 }
1967 gDirectory->cd("..");
1968
1969 gDirectory->mkdir("two");
1970 gDirectory->cd("two");
1971 for (UInt_t i = 0; i < fHistoList_Eta_two_Outside.size(); i++) {
1972 fHistoList_Eta_two_Outside[i]->Write();
1973 }
1974 gDirectory->cd("..");
1975
1976 gDirectory->cd("..");
1977
1978
1979 gDirectory->mkdir("Eta_gg_Both");
1980 gDirectory->cd("Eta_gg_Both");
1981
1982 gDirectory->mkdir("all");
1983 gDirectory->cd("all");
1984 gDirectory->mkdir("multi_all_Both");
1985 gDirectory->cd("multi_all_Both");
1986 for (UInt_t i = 0; i < fHistoList_rap_pt_Eta_all_Both.size(); i++) {
1988 }
1989 gDirectory->cd("..");
1990 for (UInt_t i = 0; i < fHistoList_Eta_all_Both.size(); i++) {
1991 fHistoList_Eta_all_Both[i]->Write();
1992 }
1993 gDirectory->cd("..");
1994
1995 gDirectory->mkdir("onetwo");
1996 gDirectory->cd("onetwo");
1997 gDirectory->mkdir("multi_onetwo_Both");
1998 gDirectory->cd("multi_onetwo_Both");
1999 for (UInt_t i = 0; i < fHistoList_rap_pt_Eta_onetwo_Both.size(); i++) {
2001 }
2002 gDirectory->cd("..");
2003 for (UInt_t i = 0; i < fHistoList_Eta_onetwo_Both.size(); i++) {
2004 fHistoList_Eta_onetwo_Both[i]->Write();
2005 }
2006 gDirectory->cd("..");
2007
2008 gDirectory->mkdir("two");
2009 gDirectory->cd("two");
2010 gDirectory->mkdir("multi_two_Both");
2011 gDirectory->cd("multi_two_Both");
2012 for (UInt_t i = 0; i < fHistoList_rap_pt_Eta_two_Both.size(); i++) {
2014 }
2015 gDirectory->cd("..");
2016 for (UInt_t i = 0; i < fHistoList_Eta_two_Both.size(); i++) {
2017 fHistoList_Eta_two_Both[i]->Write();
2018 }
2019 gDirectory->cd("..");
2020
2021 gDirectory->cd("..");
2022
2023
2024 for (UInt_t i = 0; i < fHistoList_Eta.size(); i++) {
2025 fHistoList_Eta[i]->Write();
2026 }
2027
2028
2029 gDirectory->cd("..");
2030}
2031
2033{
2035 new TH1D("Particle_pull_px_Target", "Target track, pull Px: mc - reco; P_{x}(mc) - P_{x}(reco)", 100, -0.5, 0.5);
2038 new TH1D("Particle_pull_py_Target", "Target track, pull Px: mc - reco; P_{y}(mc) - P_{y}(reco)", 100, -0.5, 0.5);
2041 new TH1D("Particle_pull_pz_Target", "Target track, pull Px: mc - reco; P_{z}(mc) - P_{z}(reco)", 100, -0.5, 0.5);
2044 new TH1D("Particle_pull_px_Outside", "Outside track, pull Px: mc - reco; P_{x}(mc) - P_{x}(reco)", 100, -0.5, 0.5);
2047 new TH1D("Particle_pull_py_Outside", "Outside track, pull Px: mc - reco; P_{y}(mc) - P_{y}(reco)", 100, -0.5, 0.5);
2050 new TH1D("Particle_pull_pz_Outside", "Outside track, pull Px: mc - reco; P_{z}(mc) - P_{z}(reco)", 100, -0.5, 0.5);
2053 new TH1D("Particle_pull_X_Outside", "Outside track, pull Px: mc - reco; X(mc) - X(reco)", 500, -5, 5);
2056 new TH1D("Particle_pull_Y_Outside", "Outside track, pull Px: mc - reco; Y(mc) - Y(reco)", 500, -5, 5);
2059 new TH1D("Particle_pull_Z_Outside", "Outside track, pull Px: mc - reco; Z(mc) - Z(reco)", 500, -5, 5);
2061
2064 new TH1D("GammasInvMass_fromEta_beforeCuts_zeroInRich_Target",
2065 "GammasInvMass_fromEta_beforeCuts_zeroInRich_Target; invariant "
2066 "mass in GeV/c^{2}",
2067 510, -0.01, 0.5);
2070 new TH1D("GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Target",
2071 "GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Target; opening "
2072 "angle in degree",
2073 300, -0.1, 29.9);
2076 new TH1D("GammasMomentum_fromEta_beforeCuts_zeroInRich_Target",
2077 "GammasMomentum_fromEta_beforeCuts_zeroInRich_Target; p in GeV/c", 400, 0, 4);
2080 new TH1D("GammasRapidity_fromEta_beforeCuts_zeroInRich_Target",
2081 "GammasRapidity_fromEta_beforeCuts_zeroInRich_Target; rapidity", 400, 0., 4.);
2084 new TH1D("PlaneAngles_last_fromEta_zeroInRich_Target",
2085 "PlaneAngles_last_fromEta_zeroInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2088 new TH1D("PlaneAngles_first_fromEta_zeroInRich_Target",
2089 "PlaneAngles_first_fromEta_zeroInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2092 new TH1D("GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Target",
2093 "GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Target; invariant "
2094 "mass in GeV/c^{2}",
2095 510, -0.01, 0.5);
2098 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Target",
2099 "GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Target; opening "
2100 "angle in degree",
2101 300, -0.1, 29.9);
2104 new TH1D("GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Target",
2105 "GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Target; p in GeV/c", 400, 0, 4);
2108 new TH1D("GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Target",
2109 "GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Target; rapidity", 400, 0., 4.);
2112 new TH1D("PlaneAngles_last_wrongpairs_zeroInRich_Target",
2113 "PlaneAngles_last_wrongpairs_zeroInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2116 new TH1D("PlaneAngles_first_wrongpairs_zeroInRich_Target",
2117 "PlaneAngles_first_wrongpairs_zeroInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2121 new TH1D("GammasInvMass_fromEta_beforeCuts_oneInRich_Target",
2122 "GammasInvMass_fromEta_beforeCuts_oneInRich_Target; invariant "
2123 "mass in GeV/c^{2}",
2124 510, -0.01, 0.5);
2127 new TH1D("GammasOpenAngle_fromEta_beforeCuts_oneInRich_Target",
2128 "GammasOpenAngle_fromEta_beforeCuts_oneInRich_Target; opening "
2129 "angle in degree",
2130 300, -0.1, 29.9);
2133 new TH1D("GammasMomentum_fromEta_beforeCuts_oneInRich_Target",
2134 "GammasMomentum_fromEta_beforeCuts_oneInRich_Target; p in GeV/c", 400, 0, 4);
2137 new TH1D("GammasRapidity_fromEta_beforeCuts_oneInRich_Target",
2138 "GammasRapidity_fromEta_beforeCuts_oneInRich_Target; rapidity", 400, 0., 4.);
2141 new TH1D("PlaneAngles_last_fromEta_oneInRich_Target",
2142 "PlaneAngles_last_fromEta_oneInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2145 new TH1D("PlaneAngles_first_fromEta_oneInRich_Target",
2146 "PlaneAngles_first_fromEta_oneInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2149 new TH1D("GammasInvMass_wrongpairs_beforeCuts_oneInRich_Target",
2150 "GammasInvMass_wrongpairs_beforeCuts_oneInRich_Target; invariant "
2151 "mass in GeV/c^{2}",
2152 510, -0.01, 0.5);
2155 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Target",
2156 "GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Target; opening "
2157 "angle in degree",
2158 300, -0.1, 29.9);
2161 new TH1D("GammasMomentum_wrongpairs_beforeCuts_oneInRich_Target",
2162 "GammasMomentum_wrongpairs_beforeCuts_oneInRich_Target; p in GeV/c", 400, 0, 4);
2165 new TH1D("GammasRapidity_wrongpairs_beforeCuts_oneInRich_Target",
2166 "GammasRapidity_wrongpairs_beforeCuts_oneInRich_Target; rapidity", 400, 0., 4.);
2169 new TH1D("PlaneAngles_last_wrongpairs_oneInRich_Target",
2170 "PlaneAngles_last_wrongpairs_oneInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2173 new TH1D("PlaneAngles_first_wrongpairs_oneInRich_Target",
2174 "PlaneAngles_first_wrongpairs_oneInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2178 new TH1D("GammasInvMass_fromEta_beforeCuts_twoInRich_Target",
2179 "GammasInvMass_fromEta_beforeCuts_twoInRich_Target; invariant "
2180 "mass in GeV/c^{2}",
2181 510, -0.01, 0.5);
2184 new TH1D("GammasOpenAngle_fromEta_beforeCuts_twoInRich_Target",
2185 "GammasOpenAngle_fromEta_beforeCuts_twoInRich_Target; opening "
2186 "angle in degree",
2187 300, -0.1, 29.9);
2190 new TH1D("GammasMomentum_fromEta_beforeCuts_twoInRich_Target",
2191 "GammasMomentum_fromEta_beforeCuts_twoInRich_Target; p in GeV/c", 400, 0, 4);
2194 new TH1D("GammasRapidity_fromEta_beforeCuts_twoInRich_Target",
2195 "GammasRapidity_fromEta_beforeCuts_twoInRich_Target; rapidity", 400, 0., 4.);
2198 new TH1D("PlaneAngles_last_fromEta_twoInRich_Target",
2199 "PlaneAngles_last_fromEta_twoInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2202 new TH1D("PlaneAngles_first_fromEta_twoInRich_Target",
2203 "PlaneAngles_first_fromEta_twoInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2206 new TH1D("GammasInvMass_wrongpairs_beforeCuts_twoInRich_Target",
2207 "GammasInvMass_wrongpairs_beforeCuts_twoInRich_Target; invariant "
2208 "mass in GeV/c^{2}",
2209 510, -0.01, 0.5);
2212 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Target",
2213 "GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Target; opening "
2214 "angle in degree",
2215 300, -0.1, 29.9);
2218 new TH1D("GammasMomentum_wrongpairs_beforeCuts_twoInRich_Target",
2219 "GammasMomentum_wrongpairs_beforeCuts_twoInRich_Target; p in GeV/c", 400, 0, 4);
2222 new TH1D("GammasRapidity_wrongpairs_beforeCuts_twoInRich_Target",
2223 "GammasRapidity_wrongpairs_beforeCuts_twoInRich_Target; rapidity", 400, 0., 4.);
2226 new TH1D("PlaneAngles_last_wrongpairs_twoInRich_Target",
2227 "PlaneAngles_last_wrongpairs_twoInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2230 new TH1D("PlaneAngles_first_wrongpairs_twoInRich_Target",
2231 "PlaneAngles_first_wrongpairs_twoInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2233
2234
2237 new TH1D("GammasInvMass_fromEta_beforeCuts_zeroInRich_Outside",
2238 "GammasInvMass_fromEta_beforeCuts_zeroInRich_Outside; invariant "
2239 "mass in GeV/c^{2}",
2240 510, -0.01, 0.5);
2243 new TH1D("GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Outside",
2244 "GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Outside; opening "
2245 "angle in degree",
2246 300, -0.1, 29.9);
2249 new TH1D("GammasMomentum_fromEta_beforeCuts_zeroInRich_Outside",
2250 "GammasMomentum_fromEta_beforeCuts_zeroInRich_Outside; p in GeV/c", 400, 0, 4);
2253 new TH1D("GammasRapidity_fromEta_beforeCuts_zeroInRich_Outside",
2254 "GammasRapidity_fromEta_beforeCuts_zeroInRich_Outside; rapidity", 400, 0., 4.);
2257 new TH1D("PlaneAngles_last_fromEta_zeroInRich_Outside",
2258 "PlaneAngles_last_fromEta_zeroInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2261 new TH1D("PlaneAngles_first_fromEta_zeroInRich_Outside",
2262 "PlaneAngles_first_fromEta_zeroInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2265 new TH1D("GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Outside",
2266 "GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Outside; "
2267 "invariant mass in GeV/c^{2}",
2268 510, -0.01, 0.5);
2271 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Outside",
2272 "GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Outside; "
2273 "opening angle in degree",
2274 300, -0.1, 29.9);
2277 new TH1D("GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Outside",
2278 "GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Outside; p in GeV/c", 400, 0, 4);
2281 new TH1D("GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Outside",
2282 "GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Outside; rapidity", 400, 0., 4.);
2285 new TH1D("PlaneAngles_last_wrongpairs_zeroInRich_Outside",
2286 "PlaneAngles_last_wrongpairs_zeroInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2289 new TH1D("PlaneAngles_first_wrongpairs_zeroInRich_Outside",
2290 "PlaneAngles_first_wrongpairs_zeroInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2294 new TH1D("GammasInvMass_fromEta_beforeCuts_oneInRich_Outside",
2295 "GammasInvMass_fromEta_beforeCuts_oneInRich_Outside; invariant "
2296 "mass in GeV/c^{2}",
2297 510, -0.01, 0.5);
2300 new TH1D("GammasOpenAngle_fromEta_beforeCuts_oneInRich_Outside",
2301 "GammasOpenAngle_fromEta_beforeCuts_oneInRich_Outside; opening "
2302 "angle in degree",
2303 300, -0.1, 29.9);
2306 new TH1D("GammasMomentum_fromEta_beforeCuts_oneInRich_Outside",
2307 "GammasMomentum_fromEta_beforeCuts_oneInRich_Outside; p in GeV/c", 400, 0, 4);
2310 new TH1D("GammasRapidity_fromEta_beforeCuts_oneInRich_Outside",
2311 "GammasRapidity_fromEta_beforeCuts_oneInRich_Outside; rapidity", 400, 0., 4.);
2314 new TH1D("PlaneAngles_last_fromEta_oneInRich_Outside",
2315 "PlaneAngles_last_fromEta_oneInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2318 new TH1D("PlaneAngles_first_fromEta_oneInRich_Outside",
2319 "PlaneAngles_first_fromEta_oneInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2322 new TH1D("GammasInvMass_wrongpairs_beforeCuts_oneInRich_Outside",
2323 "GammasInvMass_wrongpairs_beforeCuts_oneInRich_Outside; invariant "
2324 "mass in GeV/c^{2}",
2325 510, -0.01, 0.5);
2328 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Outside",
2329 "GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Outside; opening "
2330 "angle in degree",
2331 300, -0.1, 29.9);
2334 new TH1D("GammasMomentum_wrongpairs_beforeCuts_oneInRich_Outside",
2335 "GammasMomentum_wrongpairs_beforeCuts_oneInRich_Outside; p in GeV/c", 400, 0, 4);
2338 new TH1D("GammasRapidity_wrongpairs_beforeCuts_oneInRich_Outside",
2339 "GammasRapidity_wrongpairs_beforeCuts_oneInRich_Outside; rapidity", 400, 0., 4.);
2342 new TH1D("PlaneAngles_last_wrongpairs_oneInRich_Outside",
2343 "PlaneAngles_last_wrongpairs_oneInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2346 new TH1D("PlaneAngles_first_wrongpairs_oneInRich_Outside",
2347 "PlaneAngles_first_wrongpairs_oneInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2351 new TH1D("GammasInvMass_fromEta_beforeCuts_twoInRich_Outside",
2352 "GammasInvMass_fromEta_beforeCuts_twoInRich_Outside; invariant "
2353 "mass in GeV/c^{2}",
2354 510, -0.01, 0.5);
2357 new TH1D("GammasOpenAngle_fromEta_beforeCuts_twoInRich_Outside",
2358 "GammasOpenAngle_fromEta_beforeCuts_twoInRich_Outside; opening "
2359 "angle in degree",
2360 300, -0.1, 29.9);
2363 new TH1D("GammasMomentum_fromEta_beforeCuts_twoInRich_Outside",
2364 "GammasMomentum_fromEta_beforeCuts_twoInRich_Outside; p in GeV/c", 400, 0, 4);
2367 new TH1D("GammasRapidity_fromEta_beforeCuts_twoInRich_Outside",
2368 "GammasRapidity_fromEta_beforeCuts_twoInRich_Outside; rapidity", 400, 0., 4.);
2371 new TH1D("PlaneAngles_last_fromEta_twoInRich_Outside",
2372 "PlaneAngles_last_fromEta_twoInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2375 new TH1D("PlaneAngles_first_fromEta_twoInRich_Outside",
2376 "PlaneAngles_first_fromEta_twoInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2379 new TH1D("GammasInvMass_wrongpairs_beforeCuts_twoInRich_Outside",
2380 "GammasInvMass_wrongpairs_beforeCuts_twoInRich_Outside; invariant "
2381 "mass in GeV/c^{2}",
2382 510, -0.01, 0.5);
2385 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Outside",
2386 "GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Outside; opening "
2387 "angle in degree",
2388 300, -0.1, 29.9);
2391 new TH1D("GammasMomentum_wrongpairs_beforeCuts_twoInRich_Outside",
2392 "GammasMomentum_wrongpairs_beforeCuts_twoInRich_Outside; p in GeV/c", 400, 0, 4);
2395 new TH1D("GammasRapidity_wrongpairs_beforeCuts_twoInRich_Outside",
2396 "GammasRapidity_wrongpairs_beforeCuts_twoInRich_Outside; rapidity", 400, 0., 4.);
2399 new TH1D("PlaneAngles_last_wrongpairs_twoInRich_Outside",
2400 "PlaneAngles_last_wrongpairs_twoInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2403 new TH1D("PlaneAngles_first_wrongpairs_twoInRich_Outside",
2404 "PlaneAngles_first_wrongpairs_twoInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2406
2407
2410 new TH1D("GammasInvMass_fromEta_beforeCuts_zeroInRich_Both",
2411 "GammasInvMass_fromEta_beforeCuts_zeroInRich_Both; invariant mass "
2412 "in GeV/c^{2}",
2413 510, -0.01, 0.5);
2416 new TH1D("GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Both",
2417 "GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Both; opening "
2418 "angle in degree",
2419 300, -0.1, 29.9);
2422 new TH1D("GammasMomentum_fromEta_beforeCuts_zeroInRich_Both",
2423 "GammasMomentum_fromEta_beforeCuts_zeroInRich_Both; p in GeV/c", 400, 0, 4);
2426 new TH1D("GammasRapidity_fromEta_beforeCuts_zeroInRich_Both",
2427 "GammasRapidity_fromEta_beforeCuts_zeroInRich_Both; rapidity", 400, 0., 4.);
2430 new TH1D("PlaneAngles_last_fromEta_zeroInRich_Both",
2431 "PlaneAngles_last_fromEta_zeroInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2434 new TH1D("PlaneAngles_first_fromEta_zeroInRich_Both",
2435 "PlaneAngles_first_fromEta_zeroInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2438 new TH1D("GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Both",
2439 "GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Both; invariant "
2440 "mass in GeV/c^{2}",
2441 510, -0.01, 0.5);
2444 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Both",
2445 "GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Both; opening "
2446 "angle in degree",
2447 300, -0.1, 29.9);
2450 new TH1D("GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Both",
2451 "GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Both; p in GeV/c", 400, 0, 4);
2454 new TH1D("GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Both",
2455 "GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Both; rapidity", 400, 0., 4.);
2458 new TH1D("PlaneAngles_last_wrongpairs_zeroInRich_Both",
2459 "PlaneAngles_last_wrongpairs_zeroInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2462 new TH1D("PlaneAngles_first_wrongpairs_zeroInRich_Both",
2463 "PlaneAngles_first_wrongpairs_zeroInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2467 new TH1D("GammasInvMass_fromEta_beforeCuts_oneInRich_Both",
2468 "GammasInvMass_fromEta_beforeCuts_oneInRich_Both; invariant mass "
2469 "in GeV/c^{2}",
2470 510, -0.01, 0.5);
2473 new TH1D("GammasOpenAngle_fromEta_beforeCuts_oneInRich_Both",
2474 "GammasOpenAngle_fromEta_beforeCuts_oneInRich_Both; opening angle "
2475 "in degree",
2476 300, -0.1, 29.9);
2479 new TH1D("GammasMomentum_fromEta_beforeCuts_oneInRich_Both",
2480 "GammasMomentum_fromEta_beforeCuts_oneInRich_Both; p in GeV/c", 400, 0, 4);
2483 new TH1D("GammasRapidity_fromEta_beforeCuts_oneInRich_Both",
2484 "GammasRapidity_fromEta_beforeCuts_oneInRich_Both; rapidity", 400, 0., 4.);
2487 new TH1D("PlaneAngles_last_fromEta_oneInRich_Both",
2488 "PlaneAngles_last_fromEta_oneInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2491 new TH1D("PlaneAngles_first_fromEta_oneInRich_Both",
2492 "PlaneAngles_first_fromEta_oneInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2495 new TH1D("GammasInvMass_wrongpairs_beforeCuts_oneInRich_Both",
2496 "GammasInvMass_wrongpairs_beforeCuts_oneInRich_Both; invariant "
2497 "mass in GeV/c^{2}",
2498 510, -0.01, 0.5);
2501 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Both",
2502 "GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Both; opening "
2503 "angle in degree",
2504 300, -0.1, 29.9);
2507 new TH1D("GammasMomentum_wrongpairs_beforeCuts_oneInRich_Both",
2508 "GammasMomentum_wrongpairs_beforeCuts_oneInRich_Both; p in GeV/c", 400, 0, 4);
2511 new TH1D("GammasRapidity_wrongpairs_beforeCuts_oneInRich_Both",
2512 "GammasRapidity_wrongpairs_beforeCuts_oneInRich_Both; rapidity", 400, 0., 4.);
2515 new TH1D("PlaneAngles_last_wrongpairs_oneInRich_Both",
2516 "PlaneAngles_last_wrongpairs_oneInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2519 new TH1D("PlaneAngles_first_wrongpairs_oneInRich_Both",
2520 "PlaneAngles_first_wrongpairs_oneInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2523 new TH1D("PlaneAngles_last_fromEta_afterCuts_oneInRich_Both",
2524 "PlaneAngles_last_fromEta_afterCuts_oneInRich_Both; #theta angle "
2525 "in degree;#",
2526 720, -1., 179.);
2529 new TH1D("PlaneAngles_first_fromEta_afterCuts_oneInRich_Both",
2530 "PlaneAngles_first_fromEta_afterCuts_oneInRich_Both; #theta angle "
2531 "in degree;#",
2532 720, -1., 179.);
2535 new TH1D("PlaneAngles_last_wrongpairs_afterCuts_oneInRich_Both",
2536 "PlaneAngles_last_wrongpairs_afterCuts_oneInRich_Both; #theta "
2537 "angle in degree;#",
2538 720, -1., 179.);
2541 new TH1D("PlaneAngles_first_wrongpairs_afterCuts_oneInRich_Both",
2542 "PlaneAngles_first_wrongpairs_afterCuts_oneInRich_Both; #theta "
2543 "angle in degree;#",
2544 720, -1., 179.);
2548 new TH1D("GammasInvMass_fromEta_beforeCuts_twoInRich_Both",
2549 "GammasInvMass_fromEta_beforeCuts_twoInRich_Both; invariant mass "
2550 "in GeV/c^{2}",
2551 510, -0.01, 0.5);
2554 new TH1D("GammasOpenAngle_fromEta_beforeCuts_twoInRich_Both",
2555 "GammasOpenAngle_fromEta_beforeCuts_twoInRich_Both; opening angle "
2556 "in degree",
2557 300, -0.1, 29.9);
2560 new TH1D("GammasMomentum_fromEta_beforeCuts_twoInRich_Both",
2561 "GammasMomentum_fromEta_beforeCuts_twoInRich_Both; p in GeV/c", 400, 0, 4);
2564 new TH1D("GammasRapidity_fromEta_beforeCuts_twoInRich_Both",
2565 "GammasRapidity_fromEta_beforeCuts_twoInRich_Both; rapidity", 400, 0., 4.);
2568 new TH1D("PlaneAngles_last_fromEta_twoInRich_Both",
2569 "PlaneAngles_last_fromEta_twoInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2572 new TH1D("PlaneAngles_first_fromEta_twoInRich_Both",
2573 "PlaneAngles_first_fromEta_twoInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2576 new TH1D("GammasInvMass_wrongpairs_beforeCuts_twoInRich_Both",
2577 "GammasInvMass_wrongpairs_beforeCuts_twoInRich_Both; invariant "
2578 "mass in GeV/c^{2}",
2579 510, -0.01, 0.5);
2582 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Both",
2583 "GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Both; opening "
2584 "angle in degree",
2585 300, -0.1, 29.9);
2588 new TH1D("GammasMomentum_wrongpairs_beforeCuts_twoInRich_Both",
2589 "GammasMomentum_wrongpairs_beforeCuts_twoInRich_Both; p in GeV/c", 400, 0, 4);
2592 new TH1D("GammasRapidity_wrongpairs_beforeCuts_twoInRich_Both",
2593 "GammasRapidity_wrongpairs_beforeCuts_twoInRich_Both; rapidity", 400, 0., 4.);
2596 new TH1D("PlaneAngles_last_wrongpairs_twoInRich_Both",
2597 "PlaneAngles_last_wrongpairs_twoInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2600 new TH1D("PlaneAngles_first_wrongpairs_twoInRich_Both",
2601 "PlaneAngles_first_wrongpairs_twoInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2603
2604
2608 new TH1D("AngleBetweenGammas_trueEta_before_cuts_all_Target",
2609 "AngleBetweenGammas_trueEta_before_cuts_all_Target (between "
2610 "#gamma#gamma from #eta); angle [deg]",
2611 900, -0.1, 89.9);
2614 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_all_Target",
2615 "AngleBetweenGammas_wrongEta_before_cuts_all_Target (between "
2616 "#gamma#gamma); angle [deg]",
2617 900, -0.1, 89.9);
2620 new TH1D("InvMass_GammaCombinations_all_Target",
2621 "InvMass_GammaCombinations_all_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2623 Particles_PDG_all_Target = new TH1D("Particles_PDG_all_Target", "Particles_PDG_all_Target ; PDG", 2300, 0, 2300);
2626 new TH2D("TrueEta_pt_vs_rap_all_Target", "TrueEta_pt_vs_rap_all_Target; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2627 60, -1., 5.);
2630 new TH2D("WrongEta_pt_vs_rap_all_Target", "WrongEta_pt_vs_rap_all_Target; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2631 60, -1., 5.);
2634 new TH2D("TrueEta_pt_vs_rap_est_all_Target", "TrueEta_pt_vs_rap_est_all_Target; rapidity y; p_{t} in GeV/c ", 10,
2635 0., 4., 5, 0., 2.);
2638 new TH2D("WrongEta_pt_vs_rap_est_all_Target", "WrongEta_pt_vs_rap_est_all_Target; rapidity y; p_{t} in GeV/c ", 10,
2639 0., 4., 5, 0., 2.);
2642 new TH1D("TrueEta_InvMass_after_cuts_all_Target",
2643 "TrueEta_InvMass_after_cuts_all_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2646 new TH1D("WrongEta_InvMass_after_cuts_all_Target",
2647 "WrongEta_InvMass_after_cuts_all_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2650 new TH1D("TrueEta_AngleBetweenGammas_after_cuts_all_Target",
2651 "TrueEta_AngleBetweenGammas_after_cuts_all_Target (between "
2652 "#gamma#gamma from #eta); angle [deg]",
2653 900, -0.1, 89.9);
2656 new TH1D("WrongEta_AngleBetweenGammas_after_cuts_all_Target",
2657 "WrongEta_AngleBetweenGammas_after_cuts_all_Target (between "
2658 "#gamma#gamma); angle [deg]",
2659 900, -0.1, 89.9);
2662 new TH1D("EMT_Eta_InvMass_all_Target", "EMT_Eta_InvMass_all_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2666 new TH1D("AngleBetweenGammas_trueEta_before_cuts_onetwo_Target",
2667 "AngleBetweenGammas_trueEta_before_cuts_onetwo_Target (between "
2668 "#gamma#gamma from #eta); angle [deg]",
2669 900, -0.1, 89.9);
2672 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_onetwo_Target",
2673 "AngleBetweenGammas_wrongEta_before_cuts_onetwo_Target (between "
2674 "#gamma#gamma); angle [deg]",
2675 900, -0.1, 89.9);
2678 new TH1D("InvMass_GammaCombinations_onetwo_Target",
2679 "InvMass_GammaCombinations_onetwo_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2682 new TH1D("Particles_PDG_onetwo_Target", "Particles_PDG_onetwo_Target ; PDG", 2300, 0, 2300);
2685 new TH2D("TrueEta_pt_vs_rap_onetwo_Target", "TrueEta_pt_vs_rap_onetwo_Target; rapidity y; p_{t} in GeV/c ", 90, -2.,
2686 7., 60, -1., 5.);
2689 new TH2D("WrongEta_pt_vs_rap_onetwo_Target", "WrongEta_pt_vs_rap_onetwo_Target; rapidity y; p_{t} in GeV/c ", 90,
2690 -2., 7., 60, -1., 5.);
2693 new TH2D("TrueEta_pt_vs_rap_est_onetwo_Target", "TrueEta_pt_vs_rap_est_onetwo_Target; rapidity y; p_{t} in GeV/c ",
2694 10, 0., 4., 5, 0., 2.);
2697 new TH2D("WrongEta_pt_vs_rap_est_onetwo_Target",
2698 "WrongEta_pt_vs_rap_est_onetwo_Target; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 5, 0., 2.);
2701 new TH1D("TrueEta_InvMass_after_cuts_onetwo_Target",
2702 "TrueEta_InvMass_after_cuts_onetwo_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2705 new TH1D("WrongEta_InvMass_after_cuts_onetwo_Target",
2706 "WrongEta_InvMass_after_cuts_onetwo_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2709 new TH1D("TrueEta_AngleBetweenGammas_after_cuts_onetwo_Target",
2710 "TrueEta_AngleBetweenGammas_after_cuts_onetwo_Target (between "
2711 "#gamma#gamma from #eta); angle [deg]",
2712 900, -0.1, 89.9);
2715 new TH1D("WrongEta_AngleBetweenGammas_after_cuts_onetwo_Target",
2716 "WrongEta_AngleBetweenGammas_after_cuts_onetwo_Target (between "
2717 "#gamma#gamma); angle [deg]",
2718 900, -0.1, 89.9);
2721 "EMT_Eta_InvMass_onetwo_Target", "EMT_Eta_InvMass_onetwo_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2725 new TH1D("AngleBetweenGammas_trueEta_before_cuts_two_Target",
2726 "AngleBetweenGammas_trueEta_before_cuts_two_Target (between "
2727 "#gamma#gamma from #eta); angle [deg]",
2728 900, -0.1, 89.9);
2731 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_two_Target",
2732 "AngleBetweenGammas_wrongEta_before_cuts_two_Target (between "
2733 "#gamma#gamma); angle [deg]",
2734 900, -0.1, 89.9);
2737 new TH1D("InvMass_GammaCombinations_two_Target",
2738 "InvMass_GammaCombinations_two_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2740 Particles_PDG_two_Target = new TH1D("Particles_PDG_two_Target", "Particles_PDG_two_Target ; PDG", 2300, 0, 2300);
2743 new TH2D("TrueEta_pt_vs_rap_two_Target", "TrueEta_pt_vs_rap_two_Target; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2744 60, -1., 5.);
2747 new TH2D("WrongEta_pt_vs_rap_two_Target", "WrongEta_pt_vs_rap_two_Target; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2748 60, -1., 5.);
2751 new TH2D("TrueEta_pt_vs_rap_est_two_Target", "TrueEta_pt_vs_rap_est_two_Target; rapidity y; p_{t} in GeV/c ", 10,
2752 0., 4., 5, 0., 2.);
2755 new TH2D("WrongEta_pt_vs_rap_est_two_Target", "WrongEta_pt_vs_rap_est_two_Target; rapidity y; p_{t} in GeV/c ", 10,
2756 0., 4., 5, 0., 2.);
2759 new TH1D("TrueEta_InvMass_after_cuts_two_Target",
2760 "TrueEta_InvMass_after_cuts_two_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2763 new TH1D("WrongEta_InvMass_after_cuts_two_Target",
2764 "WrongEta_InvMass_after_cuts_two_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2767 new TH1D("TrueEta_AngleBetweenGammas_after_cuts_two_Target",
2768 "TrueEta_AngleBetweenGammas_after_cuts_two_Target (between "
2769 "#gamma#gamma from #eta); angle [deg]",
2770 900, -0.1, 89.9);
2773 new TH1D("WrongEta_AngleBetweenGammas_after_cuts_two_Target",
2774 "WrongEta_AngleBetweenGammas_after_cuts_two_Target (between "
2775 "#gamma#gamma); angle [deg]",
2776 900, -0.1, 89.9);
2779 new TH1D("EMT_Eta_InvMass_two_Target", "EMT_Eta_InvMass_two_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2782
2783
2787 new TH1D("AngleBetweenGammas_trueEta_before_cuts_all_Outside",
2788 "AngleBetweenGammas_trueEta_before_cuts_all_Outside (between "
2789 "#gamma#gamma from #eta); angle [deg]",
2790 900, -0.1, 89.9);
2793 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_all_Outside",
2794 "AngleBetweenGammas_wrongEta_before_cuts_all_Outside (between "
2795 "#gamma#gamma); angle [deg]",
2796 900, -0.1, 89.9);
2799 new TH1D("InvMass_GammaCombinations_all_Outside",
2800 "InvMass_GammaCombinations_all_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2802 Particles_PDG_all_Outside = new TH1D("Particles_PDG_all_Outside", "Particles_PDG_all_Outside ; PDG", 2300, 0, 2300);
2805 new TH2D("TrueEta_pt_vs_rap_all_Outside", "TrueEta_pt_vs_rap_all_Outside; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2806 60, -1., 5.);
2809 new TH2D("WrongEta_pt_vs_rap_all_Outside", "WrongEta_pt_vs_rap_all_Outside; rapidity y; p_{t} in GeV/c ", 90, -2.,
2810 7., 60, -1., 5.);
2813 new TH2D("TrueEta_pt_vs_rap_est_all_Outside", "TrueEta_pt_vs_rap_est_all_Outside; rapidity y; p_{t} in GeV/c ", 10,
2814 0., 4., 5, 0., 2.);
2817 new TH2D("WrongEta_pt_vs_rap_est_all_Outside", "WrongEta_pt_vs_rap_est_all_Outside; rapidity y; p_{t} in GeV/c ",
2818 10, 0., 4., 5, 0., 2.);
2821 new TH1D("TrueEta_InvMass_after_cuts_all_Outside",
2822 "TrueEta_InvMass_after_cuts_all_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2825 new TH1D("WrongEta_InvMass_after_cuts_all_Outside",
2826 "WrongEta_InvMass_after_cuts_all_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2829 new TH1D("TrueEta_AngleBetweenGammas_after_cuts_all_Outside",
2830 "TrueEta_AngleBetweenGammas_after_cuts_all_Outside (between "
2831 "#gamma#gamma from #eta); angle [deg]",
2832 900, -0.1, 89.9);
2835 new TH1D("WrongEta_AngleBetweenGammas_after_cuts_all_Outside",
2836 "WrongEta_AngleBetweenGammas_after_cuts_all_Outside (between "
2837 "#gamma#gamma); angle [deg]",
2838 900, -0.1, 89.9);
2840 EMT_Eta_InvMass_all_Outside = new TH1D("EMT_Eta_InvMass_all_Outside",
2841 "EMT_Eta_InvMass_all_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2845 new TH1D("AngleBetweenGammas_trueEta_before_cuts_onetwo_Outside",
2846 "AngleBetweenGammas_trueEta_before_cuts_onetwo_Outside (between "
2847 "#gamma#gamma from #eta); angle [deg]",
2848 900, -0.1, 89.9);
2851 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_onetwo_Outside",
2852 "AngleBetweenGammas_wrongEta_before_cuts_onetwo_Outside (between "
2853 "#gamma#gamma); angle [deg]",
2854 900, -0.1, 89.9);
2857 new TH1D("InvMass_GammaCombinations_onetwo_Outside",
2858 "InvMass_GammaCombinations_onetwo_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2861 new TH1D("Particles_PDG_onetwo_Outside", "Particles_PDG_onetwo_Outside ; PDG", 2300, 0, 2300);
2864 new TH2D("TrueEta_pt_vs_rap_onetwo_Outside", "TrueEta_pt_vs_rap_onetwo_Outside; rapidity y; p_{t} in GeV/c ", 90,
2865 -2., 7., 60, -1., 5.);
2868 new TH2D("WrongEta_pt_vs_rap_onetwo_Outside", "WrongEta_pt_vs_rap_onetwo_Outside; rapidity y; p_{t} in GeV/c ", 90,
2869 -2., 7., 60, -1., 5.);
2872 new TH2D("TrueEta_pt_vs_rap_est_onetwo_Outside",
2873 "TrueEta_pt_vs_rap_est_onetwo_Outside; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 5, 0., 2.);
2876 new TH2D("WrongEta_pt_vs_rap_est_onetwo_Outside",
2877 "WrongEta_pt_vs_rap_est_onetwo_Outside; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 5, 0., 2.);
2880 new TH1D("TrueEta_InvMass_after_cuts_onetwo_Outside",
2881 "TrueEta_InvMass_after_cuts_onetwo_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2884 new TH1D("WrongEta_InvMass_after_cuts_onetwo_Outside",
2885 "WrongEta_InvMass_after_cuts_onetwo_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2888 new TH1D("TrueEta_AngleBetweenGammas_after_cuts_onetwo_Outside",
2889 "TrueEta_AngleBetweenGammas_after_cuts_onetwo_Outside (between "
2890 "#gamma#gamma from #eta); angle [deg]",
2891 900, -0.1, 89.9);
2894 new TH1D("WrongEta_AngleBetweenGammas_after_cuts_onetwo_Outside",
2895 "WrongEta_AngleBetweenGammas_after_cuts_onetwo_Outside (between "
2896 "#gamma#gamma); angle [deg]",
2897 900, -0.1, 89.9);
2900 "EMT_Eta_InvMass_onetwo_Outside", "EMT_Eta_InvMass_onetwo_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2904 new TH1D("AngleBetweenGammas_trueEta_before_cuts_two_Outside",
2905 "AngleBetweenGammas_trueEta_before_cuts_two_Outside (between "
2906 "#gamma#gamma from #eta); angle [deg]",
2907 900, -0.1, 89.9);
2910 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_two_Outside",
2911 "AngleBetweenGammas_wrongEta_before_cuts_two_Outside (between "
2912 "#gamma#gamma); angle [deg]",
2913 900, -0.1, 89.9);
2916 new TH1D("InvMass_GammaCombinations_two_Outside",
2917 "InvMass_GammaCombinations_two_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2919 Particles_PDG_two_Outside = new TH1D("Particles_PDG_two_Outside", "Particles_PDG_two_Outside ; PDG", 2300, 0, 2300);
2922 new TH2D("TrueEta_pt_vs_rap_two_Outside", "TrueEta_pt_vs_rap_two_Outside; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2923 60, -1., 5.);
2926 new TH2D("WrongEta_pt_vs_rap_two_Outside", "WrongEta_pt_vs_rap_two_Outside; rapidity y; p_{t} in GeV/c ", 90, -2.,
2927 7., 60, -1., 5.);
2930 new TH2D("TrueEta_pt_vs_rap_est_two_Outside", "TrueEta_pt_vs_rap_est_two_Outside; rapidity y; p_{t} in GeV/c ", 10,
2931 0., 4., 5, 0., 2.);
2934 new TH2D("WrongEta_pt_vs_rap_est_two_Outside", "WrongEta_pt_vs_rap_est_two_Outside; rapidity y; p_{t} in GeV/c ",
2935 10, 0., 4., 5, 0., 2.);
2938 new TH1D("TrueEta_InvMass_after_cuts_two_Outside",
2939 "TrueEta_InvMass_after_cuts_two_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2942 new TH1D("WrongEta_InvMass_after_cuts_two_Outside",
2943 "WrongEta_InvMass_after_cuts_two_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2946 new TH1D("TrueEta_AngleBetweenGammas_after_cuts_two_Outside",
2947 "TrueEta_AngleBetweenGammas_after_cuts_two_Outside (between "
2948 "#gamma#gamma from #eta); angle [deg]",
2949 900, -0.1, 89.9);
2952 new TH1D("WrongEta_AngleBetweenGammas_after_cuts_two_Outside",
2953 "WrongEta_AngleBetweenGammas_after_cuts_two_Outside (between "
2954 "#gamma#gamma); angle [deg]",
2955 900, -0.1, 89.9);
2957 EMT_Eta_InvMass_two_Outside = new TH1D("EMT_Eta_InvMass_two_Outside",
2958 "EMT_Eta_InvMass_two_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2961
2962
2965 AngleBetweenGammas_trueEta_before_cuts_all_Both = new TH1D("AngleBetweenGammas_trueEta_before_cuts_all_Both",
2966 "AngleBetweenGammas_trueEta_before_cuts_all_Both (between "
2967 "#gamma#gamma from #eta); angle [deg]",
2968 900, -0.1, 89.9);
2971 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_all_Both",
2972 "AngleBetweenGammas_wrongEta_before_cuts_all_Both (between "
2973 "#gamma#gamma); angle [deg]",
2974 900, -0.1, 89.9);
2977 new TH1D("InvMass_GammaCombinations_all_Both", "InvMass_GammaCombinations_all_Both; invariant mass in GeV/c^{2};#",
2978 200, 0.0, 2.0);
2980 Particles_PDG_all_Both = new TH1D("Particles_PDG_all_Both", "Particles_PDG_all_Both ; PDG", 2300, 0, 2300);
2982 TrueEta_pt_vs_rap_all_Both = new TH2D(
2983 "TrueEta_pt_vs_rap_all_Both", "TrueEta_pt_vs_rap_all_Both; rapidity y; p_{t} in GeV/c ", 90, -2., 7., 60, -1., 5.);
2986 new TH2D("WrongEta_pt_vs_rap_all_Both", "WrongEta_pt_vs_rap_all_Both; rapidity y; p_{t} in GeV/c ", 90, -2., 7., 60,
2987 -1., 5.);
2990 new TH2D("TrueEta_pt_vs_rap_est_all_Both", "TrueEta_pt_vs_rap_est_all_Both; rapidity y; p_{t} in GeV/c ", 10, 0.,
2991 4., 5, 0., 2.);
2994 new TH2D("WrongEta_pt_vs_rap_est_all_Both", "WrongEta_pt_vs_rap_est_all_Both; rapidity y; p_{t} in GeV/c ", 10, 0.,
2995 4., 5, 0., 2.);
2998 new TH1D("TrueEta_InvMass_after_cuts_all_Both",
2999 "TrueEta_InvMass_after_cuts_all_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3002 new TH1D("WrongEta_InvMass_after_cuts_all_Both",
3003 "WrongEta_InvMass_after_cuts_all_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3005 TrueEta_AngleBetweenGammas_after_cuts_all_Both = new TH1D("TrueEta_AngleBetweenGammas_after_cuts_all_Both",
3006 "TrueEta_AngleBetweenGammas_after_cuts_all_Both (between "
3007 "#gamma#gamma from #eta); angle [deg]",
3008 900, -0.1, 89.9);
3010 WrongEta_AngleBetweenGammas_after_cuts_all_Both = new TH1D("WrongEta_AngleBetweenGammas_after_cuts_all_Both",
3011 "WrongEta_AngleBetweenGammas_after_cuts_all_Both (between "
3012 "#gamma#gamma); angle [deg]",
3013 900, -0.1, 89.9);
3016 new TH1D("EMT_Eta_InvMass_all_Both", "EMT_Eta_InvMass_all_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3020 multi_InvMass_Eta_all_Both_1 = new TH1D("multi_InvMass_Eta_all_Both_1",
3021 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3022 "invariant mass in GeV/c^{2};#",
3023 200, 0.0, 2.0);
3025 multi_InvMass_Eta_all_Both_2 = new TH1D("multi_InvMass_Eta_all_Both_2",
3026 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3027 "invariant mass in GeV/c^{2};#",
3028 200, 0.0, 2.0);
3030 multi_InvMass_Eta_all_Both_3 = new TH1D("multi_InvMass_Eta_all_Both_3",
3031 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3032 "invariant mass in GeV/c^{2};#",
3033 200, 0.0, 2.0);
3035 multi_InvMass_Eta_all_Both_4 = new TH1D("multi_InvMass_Eta_all_Both_4",
3036 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3037 "invariant mass in GeV/c^{2};#",
3038 200, 0.0, 2.0);
3040 multi_InvMass_Eta_all_Both_5 = new TH1D("multi_InvMass_Eta_all_Both_5",
3041 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3042 "invariant mass in GeV/c^{2};#",
3043 200, 0.0, 2.0);
3045 multi_InvMass_Eta_all_Both_6 = new TH1D("multi_InvMass_Eta_all_Both_6",
3046 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3047 "invariant mass in GeV/c^{2};#",
3048 200, 0.0, 2.0);
3050 multi_InvMass_Eta_all_Both_7 = new TH1D("multi_InvMass_Eta_all_Both_7",
3051 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3052 "invariant mass in GeV/c^{2};#",
3053 200, 0.0, 2.0);
3055 multi_InvMass_Eta_all_Both_8 = new TH1D("multi_InvMass_Eta_all_Both_8",
3056 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3057 "invariant mass in GeV/c^{2};#",
3058 200, 0.0, 2.0);
3060 multi_InvMass_Eta_all_Both_9 = new TH1D("multi_InvMass_Eta_all_Both_9",
3061 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3062 "invariant mass in GeV/c^{2};#",
3063 200, 0.0, 2.0);
3065 multi_InvMass_Eta_all_Both_10 = new TH1D("multi_InvMass_Eta_all_Both_10",
3066 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3067 "invariant mass in GeV/c^{2};#",
3068 200, 0.0, 2.0);
3070 multi_InvMass_Eta_all_Both_11 = new TH1D("multi_InvMass_Eta_all_Both_11",
3071 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3072 "invariant mass in GeV/c^{2};#",
3073 200, 0.0, 2.0);
3075 multi_InvMass_Eta_all_Both_12 = new TH1D("multi_InvMass_Eta_all_Both_12",
3076 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3077 "invariant mass in GeV/c^{2};#",
3078 200, 0.0, 2.0);
3080 multi_InvMass_Eta_all_Both_13 = new TH1D("multi_InvMass_Eta_all_Both_13",
3081 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3082 "invariant mass in GeV/c^{2};#",
3083 200, 0.0, 2.0);
3085 multi_InvMass_Eta_all_Both_14 = new TH1D("multi_InvMass_Eta_all_Both_14",
3086 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3087 "invariant mass in GeV/c^{2};#",
3088 200, 0.0, 2.0);
3090 multi_InvMass_Eta_all_Both_15 = new TH1D("multi_InvMass_Eta_all_Both_15",
3091 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3092 "invariant mass in GeV/c^{2};#",
3093 200, 0.0, 2.0);
3095 multi_InvMass_Eta_all_Both_16 = new TH1D("multi_InvMass_Eta_all_Both_16",
3096 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3097 "invariant mass in GeV/c^{2};#",
3098 200, 0.0, 2.0);
3100 multi_InvMass_Eta_all_Both_17 = new TH1D("multi_InvMass_Eta_all_Both_17",
3101 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3102 "invariant mass in GeV/c^{2};#",
3103 200, 0.0, 2.0);
3105 multi_EMT_Eta_all_Both_1 = new TH1D("multi_EMT_Eta_all_Both_1",
3106 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3107 "invariant mass in GeV/c^{2};#",
3108 200, 0.0, 2.0);
3110 multi_EMT_Eta_all_Both_2 = new TH1D("multi_EMT_Eta_all_Both_2",
3111 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3112 "invariant mass in GeV/c^{2};#",
3113 200, 0.0, 2.0);
3115 multi_EMT_Eta_all_Both_3 = new TH1D("multi_EMT_Eta_all_Both_3",
3116 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3117 "invariant mass in GeV/c^{2};#",
3118 200, 0.0, 2.0);
3120 multi_EMT_Eta_all_Both_4 = new TH1D("multi_EMT_Eta_all_Both_4",
3121 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3122 "invariant mass in GeV/c^{2};#",
3123 200, 0.0, 2.0);
3125 multi_EMT_Eta_all_Both_5 = new TH1D("multi_EMT_Eta_all_Both_5",
3126 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3127 "invariant mass in GeV/c^{2};#",
3128 200, 0.0, 2.0);
3130 multi_EMT_Eta_all_Both_6 = new TH1D("multi_EMT_Eta_all_Both_6",
3131 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3132 "invariant mass in GeV/c^{2};#",
3133 200, 0.0, 2.0);
3135 multi_EMT_Eta_all_Both_7 = new TH1D("multi_EMT_Eta_all_Both_7",
3136 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3137 "invariant mass in GeV/c^{2};#",
3138 200, 0.0, 2.0);
3140 multi_EMT_Eta_all_Both_8 = new TH1D("multi_EMT_Eta_all_Both_8",
3141 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3142 "invariant mass in GeV/c^{2};#",
3143 200, 0.0, 2.0);
3145 multi_EMT_Eta_all_Both_9 = new TH1D("multi_EMT_Eta_all_Both_9",
3146 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3147 "invariant mass in GeV/c^{2};#",
3148 200, 0.0, 2.0);
3150 multi_EMT_Eta_all_Both_10 = new TH1D("multi_EMT_Eta_all_Both_10",
3151 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3152 "invariant mass in GeV/c^{2};#",
3153 200, 0.0, 2.0);
3155 multi_EMT_Eta_all_Both_11 = new TH1D("multi_EMT_Eta_all_Both_11",
3156 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3157 "invariant mass in GeV/c^{2};#",
3158 200, 0.0, 2.0);
3160 multi_EMT_Eta_all_Both_12 = new TH1D("multi_EMT_Eta_all_Both_12",
3161 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3162 "invariant mass in GeV/c^{2};#",
3163 200, 0.0, 2.0);
3165 multi_EMT_Eta_all_Both_13 = new TH1D("multi_EMT_Eta_all_Both_13",
3166 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3167 "invariant mass in GeV/c^{2};#",
3168 200, 0.0, 2.0);
3170 multi_EMT_Eta_all_Both_14 = new TH1D("multi_EMT_Eta_all_Both_14",
3171 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3172 "invariant mass in GeV/c^{2};#",
3173 200, 0.0, 2.0);
3175 multi_EMT_Eta_all_Both_15 = new TH1D("multi_EMT_Eta_all_Both_15",
3176 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3177 "invariant mass in GeV/c^{2};#",
3178 200, 0.0, 2.0);
3180 multi_EMT_Eta_all_Both_16 = new TH1D("multi_EMT_Eta_all_Both_16",
3181 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3182 "invariant mass in GeV/c^{2};#",
3183 200, 0.0, 2.0);
3185 multi_EMT_Eta_all_Both_17 = new TH1D("multi_EMT_Eta_all_Both_17",
3186 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3187 "invariant mass in GeV/c^{2};#",
3188 200, 0.0, 2.0);
3190
3191
3194 new TH1D("AngleBetweenGammas_trueEta_before_cuts_onetwo_Both",
3195 "AngleBetweenGammas_trueEta_before_cuts_onetwo_Both (between "
3196 "#gamma#gamma from #eta); angle [deg]",
3197 900, -0.1, 89.9);
3200 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_onetwo_Both",
3201 "AngleBetweenGammas_wrongEta_before_cuts_onetwo_Both (between "
3202 "#gamma#gamma); angle [deg]",
3203 900, -0.1, 89.9);
3206 new TH1D("InvMass_GammaCombinations_onetwo_Both",
3207 "InvMass_GammaCombinations_onetwo_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3209 Particles_PDG_onetwo_Both = new TH1D("Particles_PDG_onetwo_Both", "Particles_PDG_onetwo_Both ; PDG", 2300, 0, 2300);
3212 new TH2D("TrueEta_pt_vs_rap_onetwo_Both", "TrueEta_pt_vs_rap_onetwo_Both; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
3213 60, -1., 5.);
3216 new TH2D("WrongEta_pt_vs_rap_onetwo_Both", "WrongEta_pt_vs_rap_onetwo_Both; rapidity y; p_{t} in GeV/c ", 90, -2.,
3217 7., 60, -1., 5.);
3220 new TH2D("TrueEta_pt_vs_rap_est_onetwo_Both", "TrueEta_pt_vs_rap_est_onetwo_Both; rapidity y; p_{t} in GeV/c ", 10,
3221 0., 4., 5, 0., 2.);
3224 new TH2D("WrongEta_pt_vs_rap_est_onetwo_Both", "WrongEta_pt_vs_rap_est_onetwo_Both; rapidity y; p_{t} in GeV/c ",
3225 10, 0., 4., 5, 0., 2.);
3228 new TH1D("TrueEta_InvMass_after_cuts_onetwo_Both",
3229 "TrueEta_InvMass_after_cuts_onetwo_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3232 new TH1D("WrongEta_InvMass_after_cuts_onetwo_Both",
3233 "WrongEta_InvMass_after_cuts_onetwo_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3236 new TH1D("TrueEta_AngleBetweenGammas_after_cuts_onetwo_Both",
3237 "TrueEta_AngleBetweenGammas_after_cuts_onetwo_Both (between "
3238 "#gamma#gamma from #eta); angle [deg]",
3239 900, -0.1, 89.9);
3242 new TH1D("WrongEta_AngleBetweenGammas_after_cuts_onetwo_Both",
3243 "WrongEta_AngleBetweenGammas_after_cuts_onetwo_Both (between "
3244 "#gamma#gamma); angle [deg]",
3245 900, -0.1, 89.9);
3247 EMT_Eta_InvMass_onetwo_Both = new TH1D("EMT_Eta_InvMass_onetwo_Both",
3248 "EMT_Eta_InvMass_onetwo_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3252 multi_InvMass_Eta_onetwo_Both_1 = new TH1D("multi_InvMass_Eta_onetwo_Both_1",
3253 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3254 "invariant mass in GeV/c^{2};#",
3255 200, 0.0, 2.0);
3257 multi_InvMass_Eta_onetwo_Both_2 = new TH1D("multi_InvMass_Eta_onetwo_Both_2",
3258 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3259 "invariant mass in GeV/c^{2};#",
3260 200, 0.0, 2.0);
3262 multi_InvMass_Eta_onetwo_Both_3 = new TH1D("multi_InvMass_Eta_onetwo_Both_3",
3263 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3264 "invariant mass in GeV/c^{2};#",
3265 200, 0.0, 2.0);
3267 multi_InvMass_Eta_onetwo_Both_4 = new TH1D("multi_InvMass_Eta_onetwo_Both_4",
3268 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3269 "invariant mass in GeV/c^{2};#",
3270 200, 0.0, 2.0);
3272 multi_InvMass_Eta_onetwo_Both_5 = new TH1D("multi_InvMass_Eta_onetwo_Both_5",
3273 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3274 "invariant mass in GeV/c^{2};#",
3275 200, 0.0, 2.0);
3277 multi_InvMass_Eta_onetwo_Both_6 = new TH1D("multi_InvMass_Eta_onetwo_Both_6",
3278 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3279 "invariant mass in GeV/c^{2};#",
3280 200, 0.0, 2.0);
3282 multi_InvMass_Eta_onetwo_Both_7 = new TH1D("multi_InvMass_Eta_onetwo_Both_7",
3283 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3284 "invariant mass in GeV/c^{2};#",
3285 200, 0.0, 2.0);
3287 multi_InvMass_Eta_onetwo_Both_8 = new TH1D("multi_InvMass_Eta_onetwo_Both_8",
3288 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3289 "invariant mass in GeV/c^{2};#",
3290 200, 0.0, 2.0);
3292 multi_InvMass_Eta_onetwo_Both_9 = new TH1D("multi_InvMass_Eta_onetwo_Both_9",
3293 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3294 "invariant mass in GeV/c^{2};#",
3295 200, 0.0, 2.0);
3297 multi_InvMass_Eta_onetwo_Both_10 = new TH1D("multi_InvMass_Eta_onetwo_Both_10",
3298 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3299 "invariant mass in GeV/c^{2};#",
3300 200, 0.0, 2.0);
3302 multi_InvMass_Eta_onetwo_Both_11 = new TH1D("multi_InvMass_Eta_onetwo_Both_11",
3303 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3304 "invariant mass in GeV/c^{2};#",
3305 200, 0.0, 2.0);
3307 multi_InvMass_Eta_onetwo_Both_12 = new TH1D("multi_InvMass_Eta_onetwo_Both_12",
3308 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3309 "invariant mass in GeV/c^{2};#",
3310 200, 0.0, 2.0);
3312 multi_InvMass_Eta_onetwo_Both_13 = new TH1D("multi_InvMass_Eta_onetwo_Both_13",
3313 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3314 "invariant mass in GeV/c^{2};#",
3315 200, 0.0, 2.0);
3317 multi_InvMass_Eta_onetwo_Both_14 = new TH1D("multi_InvMass_Eta_onetwo_Both_14",
3318 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3319 "invariant mass in GeV/c^{2};#",
3320 200, 0.0, 2.0);
3322 multi_InvMass_Eta_onetwo_Both_15 = new TH1D("multi_InvMass_Eta_onetwo_Both_15",
3323 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3324 "invariant mass in GeV/c^{2};#",
3325 200, 0.0, 2.0);
3327 multi_InvMass_Eta_onetwo_Both_16 = new TH1D("multi_InvMass_Eta_onetwo_Both_16",
3328 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3329 "invariant mass in GeV/c^{2};#",
3330 200, 0.0, 2.0);
3332 multi_InvMass_Eta_onetwo_Both_17 = new TH1D("multi_InvMass_Eta_onetwo_Both_17",
3333 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3334 "invariant mass in GeV/c^{2};#",
3335 200, 0.0, 2.0);
3337 multi_EMT_Eta_onetwo_Both_1 = new TH1D("multi_EMT_Eta_onetwo_Both_1",
3338 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3339 "invariant mass in GeV/c^{2};#",
3340 200, 0.0, 2.0);
3342 multi_EMT_Eta_onetwo_Both_2 = new TH1D("multi_EMT_Eta_onetwo_Both_2",
3343 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3344 "invariant mass in GeV/c^{2};#",
3345 200, 0.0, 2.0);
3347 multi_EMT_Eta_onetwo_Both_3 = new TH1D("multi_EMT_Eta_onetwo_Both_3",
3348 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3349 "invariant mass in GeV/c^{2};#",
3350 200, 0.0, 2.0);
3352 multi_EMT_Eta_onetwo_Both_4 = new TH1D("multi_EMT_Eta_onetwo_Both_4",
3353 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3354 "invariant mass in GeV/c^{2};#",
3355 200, 0.0, 2.0);
3357 multi_EMT_Eta_onetwo_Both_5 = new TH1D("multi_EMT_Eta_onetwo_Both_5",
3358 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3359 "invariant mass in GeV/c^{2};#",
3360 200, 0.0, 2.0);
3362 multi_EMT_Eta_onetwo_Both_6 = new TH1D("multi_EMT_Eta_onetwo_Both_6",
3363 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3364 "invariant mass in GeV/c^{2};#",
3365 200, 0.0, 2.0);
3367 multi_EMT_Eta_onetwo_Both_7 = new TH1D("multi_EMT_Eta_onetwo_Both_7",
3368 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3369 "invariant mass in GeV/c^{2};#",
3370 200, 0.0, 2.0);
3372 multi_EMT_Eta_onetwo_Both_8 = new TH1D("multi_EMT_Eta_onetwo_Both_8",
3373 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3374 "invariant mass in GeV/c^{2};#",
3375 200, 0.0, 2.0);
3377 multi_EMT_Eta_onetwo_Both_9 = new TH1D("multi_EMT_Eta_onetwo_Both_9",
3378 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3379 "invariant mass in GeV/c^{2};#",
3380 200, 0.0, 2.0);
3382 multi_EMT_Eta_onetwo_Both_10 = new TH1D("multi_EMT_Eta_onetwo_Both_10",
3383 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3384 "invariant mass in GeV/c^{2};#",
3385 200, 0.0, 2.0);
3387 multi_EMT_Eta_onetwo_Both_11 = new TH1D("multi_EMT_Eta_onetwo_Both_11",
3388 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3389 "invariant mass in GeV/c^{2};#",
3390 200, 0.0, 2.0);
3392 multi_EMT_Eta_onetwo_Both_12 = new TH1D("multi_EMT_Eta_onetwo_Both_12",
3393 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3394 "invariant mass in GeV/c^{2};#",
3395 200, 0.0, 2.0);
3397 multi_EMT_Eta_onetwo_Both_13 = new TH1D("multi_EMT_Eta_onetwo_Both_13",
3398 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3399 "invariant mass in GeV/c^{2};#",
3400 200, 0.0, 2.0);
3402 multi_EMT_Eta_onetwo_Both_14 = new TH1D("multi_EMT_Eta_onetwo_Both_14",
3403 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3404 "invariant mass in GeV/c^{2};#",
3405 200, 0.0, 2.0);
3407 multi_EMT_Eta_onetwo_Both_15 = new TH1D("multi_EMT_Eta_onetwo_Both_15",
3408 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3409 "invariant mass in GeV/c^{2};#",
3410 200, 0.0, 2.0);
3412 multi_EMT_Eta_onetwo_Both_16 = new TH1D("multi_EMT_Eta_onetwo_Both_16",
3413 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3414 "invariant mass in GeV/c^{2};#",
3415 200, 0.0, 2.0);
3417 multi_EMT_Eta_onetwo_Both_17 = new TH1D("multi_EMT_Eta_onetwo_Both_17",
3418 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3419 "invariant mass in GeV/c^{2};#",
3420 200, 0.0, 2.0);
3422
3423
3425 AngleBetweenGammas_trueEta_before_cuts_two_Both = new TH1D("AngleBetweenGammas_trueEta_before_cuts_two_Both",
3426 "AngleBetweenGammas_trueEta_before_cuts_two_Both (between "
3427 "#gamma#gamma from #eta); angle [deg]",
3428 900, -0.1, 89.9);
3431 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_two_Both",
3432 "AngleBetweenGammas_wrongEta_before_cuts_two_Both (between "
3433 "#gamma#gamma); angle [deg]",
3434 900, -0.1, 89.9);
3437 new TH1D("InvMass_GammaCombinations_two_Both", "InvMass_GammaCombinations_two_Both; invariant mass in GeV/c^{2};#",
3438 200, 0.0, 2.0);
3440 Particles_PDG_two_Both = new TH1D("Particles_PDG_two_Both", "Particles_PDG_two_Both ; PDG", 2300, 0, 2300);
3442 TrueEta_pt_vs_rap_two_Both = new TH2D(
3443 "TrueEta_pt_vs_rap_two_Both", "TrueEta_pt_vs_rap_two_Both; rapidity y; p_{t} in GeV/c ", 90, -2., 7., 60, -1., 5.);
3446 new TH2D("WrongEta_pt_vs_rap_two_Both", "WrongEta_pt_vs_rap_two_Both; rapidity y; p_{t} in GeV/c ", 90, -2., 7., 60,
3447 -1., 5.);
3450 new TH2D("TrueEta_pt_vs_rap_est_two_Both", "TrueEta_pt_vs_rap_est_two_Both; rapidity y; p_{t} in GeV/c ", 10, 0.,
3451 4., 5, 0., 2.);
3454 new TH2D("WrongEta_pt_vs_rap_est_two_Both", "WrongEta_pt_vs_rap_est_two_Both; rapidity y; p_{t} in GeV/c ", 10, 0.,
3455 4., 5, 0., 2.);
3458 new TH1D("TrueEta_InvMass_after_cuts_two_Both",
3459 "TrueEta_InvMass_after_cuts_two_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3462 new TH1D("WrongEta_InvMass_after_cuts_two_Both",
3463 "WrongEta_InvMass_after_cuts_two_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3465 TrueEta_AngleBetweenGammas_after_cuts_two_Both = new TH1D("TrueEta_AngleBetweenGammas_after_cuts_two_Both",
3466 "TrueEta_AngleBetweenGammas_after_cuts_two_Both (between "
3467 "#gamma#gamma from #eta); angle [deg]",
3468 900, -0.1, 89.9);
3470 WrongEta_AngleBetweenGammas_after_cuts_two_Both = new TH1D("WrongEta_AngleBetweenGammas_after_cuts_two_Both",
3471 "WrongEta_AngleBetweenGammas_after_cuts_two_Both (between "
3472 "#gamma#gamma); angle [deg]",
3473 900, -0.1, 89.9);
3476 new TH1D("EMT_Eta_InvMass_two_Both", "EMT_Eta_InvMass_two_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3480 multi_InvMass_Eta_two_Both_1 = new TH1D("multi_InvMass_Eta_two_Both_1",
3481 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3482 "invariant mass in GeV/c^{2};#",
3483 200, 0.0, 2.0);
3485 multi_InvMass_Eta_two_Both_2 = new TH1D("multi_InvMass_Eta_two_Both_2",
3486 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3487 "invariant mass in GeV/c^{2};#",
3488 200, 0.0, 2.0);
3490 multi_InvMass_Eta_two_Both_3 = new TH1D("multi_InvMass_Eta_two_Both_3",
3491 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3492 "invariant mass in GeV/c^{2};#",
3493 200, 0.0, 2.0);
3495 multi_InvMass_Eta_two_Both_4 = new TH1D("multi_InvMass_Eta_two_Both_4",
3496 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3497 "invariant mass in GeV/c^{2};#",
3498 200, 0.0, 2.0);
3500 multi_InvMass_Eta_two_Both_5 = new TH1D("multi_InvMass_Eta_two_Both_5",
3501 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3502 "invariant mass in GeV/c^{2};#",
3503 200, 0.0, 2.0);
3505 multi_InvMass_Eta_two_Both_6 = new TH1D("multi_InvMass_Eta_two_Both_6",
3506 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3507 "invariant mass in GeV/c^{2};#",
3508 200, 0.0, 2.0);
3510 multi_InvMass_Eta_two_Both_7 = new TH1D("multi_InvMass_Eta_two_Both_7",
3511 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3512 "invariant mass in GeV/c^{2};#",
3513 200, 0.0, 2.0);
3515 multi_InvMass_Eta_two_Both_8 = new TH1D("multi_InvMass_Eta_two_Both_8",
3516 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3517 "invariant mass in GeV/c^{2};#",
3518 200, 0.0, 2.0);
3520 multi_InvMass_Eta_two_Both_9 = new TH1D("multi_InvMass_Eta_two_Both_9",
3521 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3522 "invariant mass in GeV/c^{2};#",
3523 200, 0.0, 2.0);
3525 multi_InvMass_Eta_two_Both_10 = new TH1D("multi_InvMass_Eta_two_Both_10",
3526 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3527 "invariant mass in GeV/c^{2};#",
3528 200, 0.0, 2.0);
3530 multi_InvMass_Eta_two_Both_11 = new TH1D("multi_InvMass_Eta_two_Both_11",
3531 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3532 "invariant mass in GeV/c^{2};#",
3533 200, 0.0, 2.0);
3535 multi_InvMass_Eta_two_Both_12 = new TH1D("multi_InvMass_Eta_two_Both_12",
3536 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3537 "invariant mass in GeV/c^{2};#",
3538 200, 0.0, 2.0);
3540 multi_InvMass_Eta_two_Both_13 = new TH1D("multi_InvMass_Eta_two_Both_13",
3541 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3542 "invariant mass in GeV/c^{2};#",
3543 200, 0.0, 2.0);
3545 multi_InvMass_Eta_two_Both_14 = new TH1D("multi_InvMass_Eta_two_Both_14",
3546 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3547 "invariant mass in GeV/c^{2};#",
3548 200, 0.0, 2.0);
3550 multi_InvMass_Eta_two_Both_15 = new TH1D("multi_InvMass_Eta_two_Both_15",
3551 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3552 "invariant mass in GeV/c^{2};#",
3553 200, 0.0, 2.0);
3555 multi_InvMass_Eta_two_Both_16 = new TH1D("multi_InvMass_Eta_two_Both_16",
3556 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3557 "invariant mass in GeV/c^{2};#",
3558 200, 0.0, 2.0);
3560 multi_InvMass_Eta_two_Both_17 = new TH1D("multi_InvMass_Eta_two_Both_17",
3561 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3562 "invariant mass in GeV/c^{2};#",
3563 200, 0.0, 2.0);
3565 multi_EMT_Eta_two_Both_1 = new TH1D("multi_EMT_Eta_two_Both_1",
3566 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3567 "invariant mass in GeV/c^{2};#",
3568 200, 0.0, 2.0);
3570 multi_EMT_Eta_two_Both_2 = new TH1D("multi_EMT_Eta_two_Both_2",
3571 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3572 "invariant mass in GeV/c^{2};#",
3573 200, 0.0, 2.0);
3575 multi_EMT_Eta_two_Both_3 = new TH1D("multi_EMT_Eta_two_Both_3",
3576 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3577 "invariant mass in GeV/c^{2};#",
3578 200, 0.0, 2.0);
3580 multi_EMT_Eta_two_Both_4 = new TH1D("multi_EMT_Eta_two_Both_4",
3581 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3582 "invariant mass in GeV/c^{2};#",
3583 200, 0.0, 2.0);
3585 multi_EMT_Eta_two_Both_5 = new TH1D("multi_EMT_Eta_two_Both_5",
3586 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3587 "invariant mass in GeV/c^{2};#",
3588 200, 0.0, 2.0);
3590 multi_EMT_Eta_two_Both_6 = new TH1D("multi_EMT_Eta_two_Both_6",
3591 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3592 "invariant mass in GeV/c^{2};#",
3593 200, 0.0, 2.0);
3595 multi_EMT_Eta_two_Both_7 = new TH1D("multi_EMT_Eta_two_Both_7",
3596 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3597 "invariant mass in GeV/c^{2};#",
3598 200, 0.0, 2.0);
3600 multi_EMT_Eta_two_Both_8 = new TH1D("multi_EMT_Eta_two_Both_8",
3601 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3602 "invariant mass in GeV/c^{2};#",
3603 200, 0.0, 2.0);
3605 multi_EMT_Eta_two_Both_9 = new TH1D("multi_EMT_Eta_two_Both_9",
3606 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3607 "invariant mass in GeV/c^{2};#",
3608 200, 0.0, 2.0);
3610 multi_EMT_Eta_two_Both_10 = new TH1D("multi_EMT_Eta_two_Both_10",
3611 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3612 "invariant mass in GeV/c^{2};#",
3613 200, 0.0, 2.0);
3615 multi_EMT_Eta_two_Both_11 = new TH1D("multi_EMT_Eta_two_Both_11",
3616 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3617 "invariant mass in GeV/c^{2};#",
3618 200, 0.0, 2.0);
3620 multi_EMT_Eta_two_Both_12 = new TH1D("multi_EMT_Eta_two_Both_12",
3621 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3622 "invariant mass in GeV/c^{2};#",
3623 200, 0.0, 2.0);
3625 multi_EMT_Eta_two_Both_13 = new TH1D("multi_EMT_Eta_two_Both_13",
3626 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3627 "invariant mass in GeV/c^{2};#",
3628 200, 0.0, 2.0);
3630 multi_EMT_Eta_two_Both_14 = new TH1D("multi_EMT_Eta_two_Both_14",
3631 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3632 "invariant mass in GeV/c^{2};#",
3633 200, 0.0, 2.0);
3635 multi_EMT_Eta_two_Both_15 = new TH1D("multi_EMT_Eta_two_Both_15",
3636 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3637 "invariant mass in GeV/c^{2};#",
3638 200, 0.0, 2.0);
3640 multi_EMT_Eta_two_Both_16 = new TH1D("multi_EMT_Eta_two_Both_16",
3641 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3642 "invariant mass in GeV/c^{2};#",
3643 200, 0.0, 2.0);
3645 multi_EMT_Eta_two_Both_17 = new TH1D("multi_EMT_Eta_two_Both_17",
3646 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3647 "invariant mass in GeV/c^{2};#",
3648 200, 0.0, 2.0);
3651}
Here the ring is fitted with Taubin algorithm from A. Ayriyan, G. Ososkov, N. Chernov.
Data class for a reconstructed hit in the STS.
Data class for STS tracks.
int32_t GetStsTrackIndex() const
int32_t GetRichRingIndex() const
double GetZ() const
Definition CbmHit.h:71
static void SetKFParticleFromStsTrack(CbmStsTrack *track, KFParticle *particle, Int_t pdg=211, Bool_t firstPoint=kTRUE)
Double_t & GetRefZ()
Definition CbmKFVertex.h:27
Double_t & GetRefX()
Definition CbmKFVertex.h:25
Double_t & GetRefY()
Definition CbmKFVertex.h:26
TClonesArray * fRichHits
Definition CbmKresEta.h:85
TH1D * multi_EMT_Eta_onetwo_Both_5
Definition CbmKresEta.h:553
std::vector< std::vector< int > > Gammas_MCIndex_all_Outside
Definition CbmKresEta.h:176
TH1D * multi_InvMass_Eta_onetwo_Both_10
Definition CbmKresEta.h:541
TH1D * multi_EMT_Eta_two_Both_5
Definition CbmKresEta.h:606
TH2D * WrongEta_pt_vs_rap_est_all_Outside
Definition CbmKresEta.h:419
std::vector< std::vector< TVector3 > > EMT_eta_gg_pair_momenta_Target
Definition CbmKresEta.h:135
vector< TH1 * > fHistoList_Eta_all_Target
Definition CbmKresEta.h:359
TH1D * GammasRapidity_fromEta_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:251
TH1D * PlaneAngles_first_wrongpairs_oneInRich_Both
Definition CbmKresEta.h:336
TH1D * PlaneAngles_first_fromEta_zeroInRich_Outside
Definition CbmKresEta.h:271
TH1D * EMT_Eta_InvMass_two_Both
Definition CbmKresEta.h:582
std::vector< std::vector< TVector3 > > Gammas_onetwo_Both
Definition CbmKresEta.h:188
TH1D * AngleBetweenGammas_wrongEta_before_cuts_all_Target
Definition CbmKresEta.h:361
TH1D * multi_EMT_Eta_onetwo_Both_11
Definition CbmKresEta.h:559
TH1D * TrueEta_InvMass_after_cuts_two_Both
Definition CbmKresEta.h:578
TH1D * GammasRapidity_fromEta_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:237
vector< CbmRichRing * > VRichRing_minus_Outside
Definition CbmKresEta.h:100
TH1D * multi_EMT_Eta_onetwo_Both_15
Definition CbmKresEta.h:563
TH1D * multi_InvMass_Eta_two_Both_4
Definition CbmKresEta.h:588
TH1D * Particle_pull_X_Outside
Definition CbmKresEta.h:211
TH1D * multi_InvMass_Eta_all_Both_7
Definition CbmKresEta.h:485
void Mixing_Outside()
TH1D * multi_EMT_Eta_all_Both_9
Definition CbmKresEta.h:504
std::vector< int > EMT_eta_gg_Event_Target
Definition CbmKresEta.h:134
TH1D * multi_InvMass_Eta_two_Both_9
Definition CbmKresEta.h:593
TH1D * GammasMomentum_wrongpairs_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:333
std::vector< std::vector< int > > Gammas_MCIndex_onetwo_Target
Definition CbmKresEta.h:160
TH1D * multi_EMT_Eta_two_Both_2
Definition CbmKresEta.h:603
TH1D * multi_InvMass_Eta_two_Both_6
Definition CbmKresEta.h:590
TH1D * Particles_PDG_two_Target
Definition CbmKresEta.h:393
vector< TH1 * > fHistoList_rap_pt_Eta_all_Both
Definition CbmKresEta.h:478
std::vector< int > frefId
Definition CbmKresEta.h:129
TH1D * GammasMomentum_fromEta_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:313
void Mixing_Target()
TH1D * GammasMomentum_fromEta_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:282
TH1D * GammasRapidity_fromEta_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:223
TH1D * InvMass_GammaCombinations_all_Target
Definition CbmKresEta.h:362
TH1D * Particles_PDG_all_Target
Definition CbmKresEta.h:363
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Outside
Definition CbmKresEta.h:182
TClonesArray * fMcTracks
Definition CbmKresEta.h:78
TH1D * GammasOpenAngle_fromEta_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:326
TH1D * GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:221
TH1D * GammasMomentum_wrongpairs_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:256
TH1D * GammasInvMass_wrongpairs_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:286
TH1D * AngleBetweenGammas_trueEta_before_cuts_onetwo_Target
Definition CbmKresEta.h:375
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:273
TH1D * Particles_PDG_all_Both
Definition CbmKresEta.h:467
TH1D * multi_EMT_Eta_all_Both_2
Definition CbmKresEta.h:497
std::vector< int > EMT_eta_gg_Event_Both
Definition CbmKresEta.h:144
TH1D * Particle_pull_px_Target
Definition CbmKresEta.h:205
TH1D * AngleBetweenGammas_trueEta_before_cuts_two_Target
Definition CbmKresEta.h:390
TH1D * multi_EMT_Eta_two_Both_9
Definition CbmKresEta.h:610
TH1D * PlaneAngles_first_fromEta_oneInRich_Outside
Definition CbmKresEta.h:285
TH2D * TrueEta_pt_vs_rap_est_all_Outside
Definition CbmKresEta.h:418
TH2D * TrueEta_pt_vs_rap_est_two_Both
Definition CbmKresEta.h:576
TH1D * multi_EMT_Eta_all_Both_15
Definition CbmKresEta.h:510
TH1D * Particles_PDG_two_Outside
Definition CbmKresEta.h:445
vector< CbmRichRing * > VRichRing_plus_Target
Definition CbmKresEta.h:125
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:318
TH1D * multi_EMT_Eta_onetwo_Both_9
Definition CbmKresEta.h:557
TH1D * WrongEta_InvMass_after_cuts_all_Both
Definition CbmKresEta.h:473
TH1D * Particle_pull_pz_Target
Definition CbmKresEta.h:207
vector< CbmStsTrack * > VStsTrack_minus_Target
Definition CbmKresEta.h:113
std::vector< std::vector< int > > Gammas_MCIndex_onetwo_Both
Definition CbmKresEta.h:196
TH1D * GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:229
TH1D * InvMass_GammaCombinations_onetwo_Outside
Definition CbmKresEta.h:429
vector< TH1 * > fHistoList_Eta_cuts_Outside
Definition CbmKresEta.h:263
TH1D * multi_InvMass_Eta_onetwo_Both_15
Definition CbmKresEta.h:546
TH1D * GammasMomentum_fromEta_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:296
TH2D * WrongEta_pt_vs_rap_two_Outside
Definition CbmKresEta.h:447
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Outside
Definition CbmKresEta.h:180
TClonesArray * fStsTracks
Definition CbmKresEta.h:80
TH1D * PlaneAngles_first_wrongpairs_zeroInRich_Both
Definition CbmKresEta.h:322
TH1D * AngleBetweenGammas_trueEta_before_cuts_all_Target
Definition CbmKresEta.h:360
std::vector< std::vector< int > > Gammas_stsIndex_all_Both
Definition CbmKresEta.h:190
TH1D * multi_EMT_Eta_two_Both_17
Definition CbmKresEta.h:618
TH2D * TrueEta_pt_vs_rap_est_two_Outside
Definition CbmKresEta.h:448
vector< TH1 * > fHistoList_rap_pt_Eta_two_Target
Definition CbmKresEta.h:406
TH1D * PlaneAngles_last_fromEta_afterCuts_oneInRich_Both
Definition CbmKresEta.h:337
std::vector< std::vector< int > > Gammas_stsIndex_two_Target
Definition CbmKresEta.h:155
TH1D * PlaneAngles_last_fromEta_zeroInRich_Target
Definition CbmKresEta.h:224
TH1D * GammasMomentum_wrongpairs_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:288
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Target
Definition CbmKresEta.h:163
TH1D * InvMass_GammaCombinations_all_Both
Definition CbmKresEta.h:466
TH2D * TrueEta_pt_vs_rap_est_onetwo_Target
Definition CbmKresEta.h:381
std::vector< int > VStsIndex_minus_Target
Definition CbmKresEta.h:116
TH1D * multi_InvMass_Eta_all_Both_5
Definition CbmKresEta.h:483
TH1D * GammasInvMass_wrongpairs_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:300
TH1D * PlaneAngles_first_wrongpairs_oneInRich_Outside
Definition CbmKresEta.h:291
TH1D * multi_EMT_Eta_onetwo_Both_6
Definition CbmKresEta.h:554
TH1D * GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:274
TH1D * PlaneAngles_last_fromEta_twoInRich_Both
Definition CbmKresEta.h:347
TH1D * Particles_PDG_onetwo_Both
Definition CbmKresEta.h:520
TH1D * WrongEta_InvMass_after_cuts_two_Both
Definition CbmKresEta.h:579
TH1D * multi_EMT_Eta_onetwo_Both_10
Definition CbmKresEta.h:558
std::vector< int > VRings_minus_Outside
Definition CbmKresEta.h:98
vector< TH1 * > fHistoList_rap_pt_Eta_onetwo_Both
Definition CbmKresEta.h:531
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Both
Definition CbmKresEta.h:192
TH1D * GammasOpenAngle_fromEta_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:249
TH1D * multi_EMT_Eta_all_Both_7
Definition CbmKresEta.h:502
TH1D * GammasInvMass_fromEta_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:266
TH1D * PlaneAngles_first_fromEta_oneInRich_Both
Definition CbmKresEta.h:330
vector< TH1 * > fHistoList_Eta_two_Outside
Definition CbmKresEta.h:441
std::vector< int > VStsIndex_plus_Outside
Definition CbmKresEta.h:106
TH1D * WrongEta_AngleBetweenGammas_after_cuts_all_Outside
Definition CbmKresEta.h:423
TH1D * InvMass_GammaCombinations_two_Outside
Definition CbmKresEta.h:444
TH1D * AngleBetweenGammas_wrongEta_before_cuts_two_Target
Definition CbmKresEta.h:391
vector< TH1 * > fHistoList_Eta_cuts_Target
Definition CbmKresEta.h:217
TH2D * WrongEta_pt_vs_rap_onetwo_Target
Definition CbmKresEta.h:380
TH1D * AngleBetweenGammas_trueEta_before_cuts_two_Outside
Definition CbmKresEta.h:442
TH2D * TrueEta_pt_vs_rap_est_two_Target
Definition CbmKresEta.h:396
TH1D * GammasRapidity_wrongpairs_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:289
TH1D * WrongEta_InvMass_after_cuts_onetwo_Target
Definition CbmKresEta.h:384
TH1D * multi_EMT_Eta_all_Both_17
Definition CbmKresEta.h:512
TH1D * multi_InvMass_Eta_all_Both_12
Definition CbmKresEta.h:490
TH1D * WrongEta_InvMass_after_cuts_all_Outside
Definition CbmKresEta.h:421
std::vector< std::vector< TVector3 > > Gammas_two_Outside
Definition CbmKresEta.h:169
TH1D * multi_EMT_Eta_two_Both_8
Definition CbmKresEta.h:609
std::vector< std::vector< int > > Gammas_stsIndex_two_Outside
Definition CbmKresEta.h:173
TH1D * GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:228
TH1D * multi_EMT_Eta_onetwo_Both_2
Definition CbmKresEta.h:550
std::vector< std::vector< TVector3 > > Gammas_onetwo_Outside
Definition CbmKresEta.h:170
TH2D * TrueEta_pt_vs_rap_all_Target
Definition CbmKresEta.h:364
TH1D * PlaneAngles_first_fromEta_zeroInRich_Both
Definition CbmKresEta.h:316
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:255
std::vector< int > fMCId
Definition CbmKresEta.h:130
TH1D * PlaneAngles_last_wrongpairs_twoInRich_Both
Definition CbmKresEta.h:353
vector< TH1 * > fHistoList_Eta_onetwo_Target
Definition CbmKresEta.h:374
TH1D * multi_InvMass_Eta_two_Both_12
Definition CbmKresEta.h:596
TH1D * EMT_Eta_InvMass_onetwo_Both
Definition CbmKresEta.h:529
TH1D * multi_InvMass_Eta_two_Both_7
Definition CbmKresEta.h:591
TClonesArray * fArrayStsHit
Definition CbmKresEta.h:87
TH1D * TrueEta_InvMass_after_cuts_onetwo_Outside
Definition CbmKresEta.h:435
TH1D * GammasInvMass_wrongpairs_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:331
void Exec(int fEventNumEta, double OpeningAngleCut, double GammaInvMassCut, int RealPID)
TH1D * multi_InvMass_Eta_two_Both_1
Definition CbmKresEta.h:585
vector< TH1 * > fHistoList_Eta_two_Target
Definition CbmKresEta.h:389
std::vector< int > VRings_plus_Outside
Definition CbmKresEta.h:105
TH1D * multi_EMT_Eta_all_Both_16
Definition CbmKresEta.h:511
TH2D * TrueEta_pt_vs_rap_onetwo_Target
Definition CbmKresEta.h:379
TH1D * TrueEta_AngleBetweenGammas_after_cuts_onetwo_Outside
Definition CbmKresEta.h:437
TH1D * PlaneAngles_last_fromEta_zeroInRich_Both
Definition CbmKresEta.h:315
void FindGammasOutside(int EventNumEta, double AngleCut, double InvMassCut, int RealPID, vector< CbmMCTrack * > MCtracks_minus_Outside, vector< CbmMCTrack * > MCtracks_plus_Outside, vector< CbmStsTrack * > StsTrack_minus_Outside, vector< CbmStsTrack * > StsTrack_plus_Outside, std::vector< int > Rings_minus_Outside, std::vector< int > Rings_plus_Outside, std::vector< int > stsIndex_minus_Outside, std::vector< int > stsIndex_plus_Outside, vector< CbmRichRing * > richRing_minus_Outside, vector< CbmRichRing * > richRing_plus_Outside, vector< Int_t > MCIndex_minus_Outside, vector< Int_t > MCIndex_plus_Outside)
void FindEta(TString mod, TString position, vector< vector< TVector3 > > Gammas, vector< vector< int > > StsIndex, vector< vector< int > > MCIndex, vector< vector< CbmMCTrack * > > GammasMC, vector< TH1 * > gg, vector< TH1 * > rap_pt_separation)
TH1D * multi_EMT_Eta_two_Both_12
Definition CbmKresEta.h:613
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Both
Definition CbmKresEta.h:199
TH1D * TrueEta_InvMass_after_cuts_all_Outside
Definition CbmKresEta.h:420
TH1D * InvMass_GammaCombinations_two_Target
Definition CbmKresEta.h:392
TH1D * multi_InvMass_Eta_two_Both_5
Definition CbmKresEta.h:589
TH1D * PlaneAngles_last_fromEta_zeroInRich_Outside
Definition CbmKresEta.h:270
TH1D * TrueEta_AngleBetweenGammas_after_cuts_two_Outside
Definition CbmKresEta.h:452
TH1D * GammasMomentum_wrongpairs_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:351
TH1D * multi_EMT_Eta_two_Both_4
Definition CbmKresEta.h:605
double CalculatePlaneAngle_last(CbmStsTrack *Sts_1, CbmStsTrack *Sts_2)
TH1D * multi_EMT_Eta_all_Both_3
Definition CbmKresEta.h:498
TH1D * multi_EMT_Eta_onetwo_Both_7
Definition CbmKresEta.h:555
int FindInRich(int richInd, int stsMcTrackId)
TH1D * multi_EMT_Eta_all_Both_14
Definition CbmKresEta.h:509
TH1D * multi_InvMass_Eta_onetwo_Both_12
Definition CbmKresEta.h:543
TH1D * PlaneAngles_first_fromEta_twoInRich_Target
Definition CbmKresEta.h:253
TH1D * multi_EMT_Eta_onetwo_Both_16
Definition CbmKresEta.h:564
TH2D * TrueEta_pt_vs_rap_two_Both
Definition CbmKresEta.h:574
TH2D * WrongEta_pt_vs_rap_all_Both
Definition CbmKresEta.h:469
TH1D * WrongEta_AngleBetweenGammas_after_cuts_two_Target
Definition CbmKresEta.h:401
TH1D * GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:275
TH1D * multi_InvMass_Eta_onetwo_Both_7
Definition CbmKresEta.h:538
TH1D * multi_InvMass_Eta_all_Both_13
Definition CbmKresEta.h:491
void InitHistograms()
TH1D * multi_InvMass_Eta_all_Both_4
Definition CbmKresEta.h:482
vector< CbmMCTrack * > VMCtracks_plus_Outside
Definition CbmKresEta.h:104
TH1D * multi_EMT_Eta_two_Both_14
Definition CbmKresEta.h:615
TH1D * multi_InvMass_Eta_two_Both_8
Definition CbmKresEta.h:592
TH1D * GammasOpenAngle_fromEta_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:295
TH1D * multi_InvMass_Eta_onetwo_Both_4
Definition CbmKresEta.h:535
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Both
Definition CbmKresEta.h:198
TH1D * GammasMomentum_wrongpairs_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:242
TH1D * multi_InvMass_Eta_two_Both_15
Definition CbmKresEta.h:599
CbmRichRingFitterEllipseTau * fTauFit
Definition CbmKresEta.h:92
TH1D * PlaneAngles_first_fromEta_zeroInRich_Target
Definition CbmKresEta.h:225
TH1D * PlaneAngles_last_wrongpairs_zeroInRich_Target
Definition CbmKresEta.h:230
TH2D * WrongEta_pt_vs_rap_all_Outside
Definition CbmKresEta.h:417
TH1D * multi_InvMass_Eta_all_Both_16
Definition CbmKresEta.h:494
TH1D * multi_InvMass_Eta_onetwo_Both_14
Definition CbmKresEta.h:545
vector< CbmMCTrack * > VMCtracks_minus_Outside
Definition CbmKresEta.h:97
std::vector< int > EMT_eta_gg_NofRings_Target
Definition CbmKresEta.h:136
TH1D * GammasRapidity_wrongpairs_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:352
vector< CbmStsTrack * > VStsTrack_minus_Outside
Definition CbmKresEta.h:96
std::vector< std::vector< TVector3 > > Gammas_onetwo_Target
Definition CbmKresEta.h:152
TH2D * WrongEta_pt_vs_rap_est_onetwo_Both
Definition CbmKresEta.h:524
vector< CbmRichRing * > VRichRing_plus_Outside
Definition CbmKresEta.h:107
TH1D * multi_EMT_Eta_all_Both_4
Definition CbmKresEta.h:499
vector< Int_t > VMCIndex_plus_Outside
Definition CbmKresEta.h:108
vector< CbmRichRing * > VRichRing_minus_Target
Definition CbmKresEta.h:117
vector< TH1 * > fHistoList_Eta_cuts_oneInRich_Outside
Definition CbmKresEta.h:279
TH2D * TrueEta_pt_vs_rap_all_Outside
Definition CbmKresEta.h:416
vector< TH1 * > fHistoList_Eta_cuts_oneInRich_Target
Definition CbmKresEta.h:233
TH1D * PlaneAngles_first_fromEta_oneInRich_Target
Definition CbmKresEta.h:239
TH1D * GammasInvMass_fromEta_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:311
std::vector< std::vector< int > > Gammas_MCIndex_two_Outside
Definition CbmKresEta.h:177
TH1D * multi_InvMass_Eta_two_Both_10
Definition CbmKresEta.h:594
TH1D * WrongEta_AngleBetweenGammas_after_cuts_all_Target
Definition CbmKresEta.h:371
TH1D * multi_InvMass_Eta_two_Both_13
Definition CbmKresEta.h:597
TH1D * PlaneAngles_last_fromEta_twoInRich_Target
Definition CbmKresEta.h:252
TH1D * multi_EMT_Eta_two_Both_1
Definition CbmKresEta.h:602
TClonesArray * fRichRings
Definition CbmKresEta.h:83
TH1D * multi_InvMass_Eta_all_Both_8
Definition CbmKresEta.h:486
TH1D * Particle_pull_py_Outside
Definition CbmKresEta.h:209
std::vector< int > VRings_plus_Target
Definition CbmKresEta.h:123
TH1D * AngleBetweenGammas_trueEta_before_cuts_all_Both
Definition CbmKresEta.h:464
TH1D * AngleBetweenGammas_wrongEta_before_cuts_onetwo_Outside
Definition CbmKresEta.h:428
TH1D * AngleBetweenGammas_trueEta_before_cuts_onetwo_Both
Definition CbmKresEta.h:517
TH1D * PlaneAngles_first_wrongpairs_oneInRich_Target
Definition CbmKresEta.h:245
TH1D * AngleBetweenGammas_trueEta_before_cuts_all_Outside
Definition CbmKresEta.h:412
vector< TVector3 > VMomenta_plus_Target
Definition CbmKresEta.h:122
TH1D * multi_InvMass_Eta_all_Both_10
Definition CbmKresEta.h:488
std::vector< std::vector< TVector3 > > Gammas_all_Both
Definition CbmKresEta.h:186
TH1D * GammasRapidity_wrongpairs_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:303
vector< TH1 * > fHistoList_Eta_cuts_oneInRich_Both
Definition CbmKresEta.h:324
TH2D * TrueEta_pt_vs_rap_est_onetwo_Both
Definition CbmKresEta.h:523
TH1D * multi_EMT_Eta_onetwo_Both_17
Definition CbmKresEta.h:565
vector< TH1 * > fHistoList_Eta_two_Both
Definition CbmKresEta.h:569
TH1D * PlaneAngles_last_fromEta_oneInRich_Target
Definition CbmKresEta.h:238
virtual ~CbmKresEta()
TH1D * PlaneAngles_first_wrongpairs_twoInRich_Both
Definition CbmKresEta.h:354
TClonesArray * fArrayMvdHit
Definition CbmKresEta.h:86
TH1D * GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:317
std::vector< std::vector< TVector3 > > Gammas_all_Target
Definition CbmKresEta.h:150
TH1D * multi_InvMass_Eta_all_Both_9
Definition CbmKresEta.h:487
TH1D * Particle_pull_pz_Outside
Definition CbmKresEta.h:210
TH1D * GammasInvMass_fromEta_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:220
TH1D * TrueEta_AngleBetweenGammas_after_cuts_onetwo_Both
Definition CbmKresEta.h:527
TH1D * GammasOpenAngle_fromEta_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:235
TH1D * TrueEta_AngleBetweenGammas_after_cuts_all_Both
Definition CbmKresEta.h:474
TH1D * GammasMomentum_fromEta_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:345
std::vector< std::vector< int > > Gammas_MCIndex_all_Target
Definition CbmKresEta.h:158
TH2D * WrongEta_pt_vs_rap_est_all_Both
Definition CbmKresEta.h:471
TH1D * multi_InvMass_Eta_all_Both_14
Definition CbmKresEta.h:492
vector< TH1 * > fHistoList_Eta_cuts_twoInRich_Outside
Definition CbmKresEta.h:293
TH1D * WrongEta_InvMass_after_cuts_all_Target
Definition CbmKresEta.h:369
vector< TH1 * > fHistoList_rap_pt_Eta_all_Target
Definition CbmKresEta.h:404
TH1D * PlaneAngles_last_wrongpairs_zeroInRich_Outside
Definition CbmKresEta.h:276
TH1D * Particle_pull_Y_Outside
Definition CbmKresEta.h:212
TH1D * WrongEta_AngleBetweenGammas_after_cuts_onetwo_Both
Definition CbmKresEta.h:528
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:241
TH1D * EMT_Eta_InvMass_onetwo_Target
Definition CbmKresEta.h:387
TH1D * multi_InvMass_Eta_all_Both_2
Definition CbmKresEta.h:480
TH1D * TrueEta_AngleBetweenGammas_after_cuts_two_Both
Definition CbmKresEta.h:580
TH1D * TrueEta_AngleBetweenGammas_after_cuts_two_Target
Definition CbmKresEta.h:400
TH1D * multi_EMT_Eta_all_Both_10
Definition CbmKresEta.h:505
TClonesArray * fRichProjections
Definition CbmKresEta.h:82
TH1D * GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:319
TH1D * AngleBetweenGammas_wrongEta_before_cuts_two_Outside
Definition CbmKresEta.h:443
TH1D * TrueEta_InvMass_after_cuts_two_Outside
Definition CbmKresEta.h:450
TH1D * TrueEta_AngleBetweenGammas_after_cuts_onetwo_Target
Definition CbmKresEta.h:385
TH1D * multi_InvMass_Eta_all_Both_11
Definition CbmKresEta.h:489
TH1D * EMT_Eta_InvMass_two_Target
Definition CbmKresEta.h:402
TH1D * multi_InvMass_Eta_two_Both_14
Definition CbmKresEta.h:598
TH2D * TrueEta_pt_vs_rap_onetwo_Outside
Definition CbmKresEta.h:431
void FindGammasTarget(int EventNumEta, double AngleCut, double InvMassCut, int RealPID, vector< CbmMCTrack * > MCtracks_minus, vector< CbmMCTrack * > MCtracks_plus, vector< CbmStsTrack * > StsTrack_minus, vector< CbmStsTrack * > StsTrack_plus, vector< TVector3 > Momenta_minus, vector< TVector3 > Momenta_plus, std::vector< int > Rings_minus, std::vector< int > Rings_plus, std::vector< int > stsIndex_minus, std::vector< int > stsIndex_plus, vector< CbmRichRing * > richRing_minus, vector< CbmRichRing * > richRing_plus, vector< Int_t > MCIndex_minus, vector< Int_t > MCIndex_plus)
TH1D * GammasInvMass_fromEta_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:280
TH1D * multi_EMT_Eta_all_Both_6
Definition CbmKresEta.h:501
vector< CbmMCTrack * > VMCtracks_plus_Target
Definition CbmKresEta.h:120
TH1D * PlaneAngles_last_wrongpairs_afterCuts_oneInRich_Both
Definition CbmKresEta.h:339
TH1D * GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:226
TH1D * multi_EMT_Eta_onetwo_Both_12
Definition CbmKresEta.h:560
TH1D * GammasInvMass_fromEta_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:234
TClonesArray * fStsTrackMatches
Definition CbmKresEta.h:81
TH1D * Particles_PDG_onetwo_Target
Definition CbmKresEta.h:378
double CalculatePlaneAngle_first(CbmStsTrack *Sts_1, CbmStsTrack *Sts_2)
std::vector< std::vector< int > > Gammas_MCIndex_all_Both
Definition CbmKresEta.h:194
TH1D * multi_EMT_Eta_two_Both_16
Definition CbmKresEta.h:617
TH2D * TrueEta_pt_vs_rap_est_all_Both
Definition CbmKresEta.h:470
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:301
vector< Int_t > VMCIndex_plus_Target
Definition CbmKresEta.h:126
TH2D * WrongEta_pt_vs_rap_est_onetwo_Outside
Definition CbmKresEta.h:434
TH1D * multi_InvMass_Eta_two_Both_11
Definition CbmKresEta.h:595
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:227
TH1D * InvMass_GammaCombinations_onetwo_Both
Definition CbmKresEta.h:519
TH1D * PlaneAngles_last_fromEta_oneInRich_Outside
Definition CbmKresEta.h:284
TH1D * Particle_pull_Z_Outside
Definition CbmKresEta.h:213
vector< TH1 * > fHistoList_Eta_all_Outside
Definition CbmKresEta.h:411
TH1D * GammasInvMass_fromEta_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:248
TH1D * multi_EMT_Eta_two_Both_6
Definition CbmKresEta.h:607
std::vector< int > EMT_eta_gg_NofRings_Both
Definition CbmKresEta.h:146
int CheckIfElectron(CbmRichRing *ring, double momentum)
CbmKFVertex fKFVertex
Definition CbmKresEta.h:90
TH1D * AngleBetweenGammas_trueEta_before_cuts_onetwo_Outside
Definition CbmKresEta.h:427
TH1D * multi_EMT_Eta_onetwo_Both_13
Definition CbmKresEta.h:561
TH1D * InvMass_GammaCombinations_onetwo_Target
Definition CbmKresEta.h:377
TH1D * multi_InvMass_Eta_onetwo_Both_11
Definition CbmKresEta.h:542
TH1D * multi_EMT_Eta_onetwo_Both_4
Definition CbmKresEta.h:552
TH1D * GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:320
TH1D * PlaneAngles_first_wrongpairs_zeroInRich_Outside
Definition CbmKresEta.h:277
TH1D * multi_InvMass_Eta_onetwo_Both_8
Definition CbmKresEta.h:539
TH2D * WrongEta_pt_vs_rap_two_Target
Definition CbmKresEta.h:395
TH1D * multi_InvMass_Eta_onetwo_Both_1
Definition CbmKresEta.h:532
vector< CbmMCTrack * > VMCtracks_minus_Target
Definition CbmKresEta.h:112
TH1D * EMT_Eta_InvMass_all_Outside
Definition CbmKresEta.h:424
TH1D * GammasInvMass_wrongpairs_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:254
TH2D * WrongEta_pt_vs_rap_est_two_Target
Definition CbmKresEta.h:397
TH1D * PlaneAngles_last_fromEta_oneInRich_Both
Definition CbmKresEta.h:329
TH2D * WrongEta_pt_vs_rap_est_two_Both
Definition CbmKresEta.h:577
TH1D * multi_EMT_Eta_all_Both_12
Definition CbmKresEta.h:507
vector< Int_t > VMCIndex_minus_Target
Definition CbmKresEta.h:118
vector< TH1 * > fHistoList_Eta_onetwo_Outside
Definition CbmKresEta.h:426
std::vector< std::vector< TVector3 > > Gammas_two_Both
Definition CbmKresEta.h:187
std::vector< int > EMT_eta_gg_Event_Outside
Definition CbmKresEta.h:139
TH1D * multi_EMT_Eta_all_Both_8
Definition CbmKresEta.h:503
TH1D * GammasMomentum_fromEta_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:268
TH1D * Particle_pull_px_Outside
Definition CbmKresEta.h:208
TH1D * EMT_Eta_InvMass_all_Target
Definition CbmKresEta.h:372
TH1D * WrongEta_InvMass_after_cuts_two_Outside
Definition CbmKresEta.h:451
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:287
TH1D * multi_InvMass_Eta_all_Both_17
Definition CbmKresEta.h:495
TH1D * multi_EMT_Eta_onetwo_Both_14
Definition CbmKresEta.h:562
CbmVertex * fPrimVertex
Definition CbmKresEta.h:89
TH1D * TrueEta_InvMass_after_cuts_onetwo_Both
Definition CbmKresEta.h:525
TH1D * GammasMomentum_fromEta_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:250
TH1D * PlaneAngles_last_wrongpairs_twoInRich_Target
Definition CbmKresEta.h:258
TH1D * GammasRapidity_fromEta_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:328
vector< Int_t > VMCIndex_minus_Outside
Definition CbmKresEta.h:101
TH1D * Particles_PDG_onetwo_Outside
Definition CbmKresEta.h:430
vector< TH1 * > fHistoList_Eta_cuts_zeroInRich_Target
Definition CbmKresEta.h:219
std::vector< std::vector< int > > Gammas_MCIndex_two_Target
Definition CbmKresEta.h:159
TH1D * multi_InvMass_Eta_onetwo_Both_3
Definition CbmKresEta.h:534
TH1D * GammasOpenAngle_fromEta_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:281
TH2D * WrongEta_pt_vs_rap_onetwo_Outside
Definition CbmKresEta.h:432
TH1D * GammasRapidity_fromEta_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:314
TH1D * PlaneAngles_last_wrongpairs_oneInRich_Target
Definition CbmKresEta.h:244
TH1D * AngleBetweenGammas_wrongEta_before_cuts_onetwo_Target
Definition CbmKresEta.h:376
TH2D * TrueEta_pt_vs_rap_est_all_Target
Definition CbmKresEta.h:366
TH1D * multi_InvMass_Eta_all_Both_1
Definition CbmKresEta.h:479
TH1D * GammasRapidity_fromEta_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:297
vector< TH1 * > fHistoList_Eta_onetwo_Both
Definition CbmKresEta.h:516
TH1D * PlaneAngles_first_fromEta_afterCuts_oneInRich_Both
Definition CbmKresEta.h:338
TH1D * multi_EMT_Eta_all_Both_13
Definition CbmKresEta.h:508
TH1D * Particles_PDG_two_Both
Definition CbmKresEta.h:573
TH1D * PlaneAngles_first_fromEta_twoInRich_Outside
Definition CbmKresEta.h:299
TH1D * PlaneAngles_first_wrongpairs_afterCuts_oneInRich_Both
Definition CbmKresEta.h:340
TH1D * GammasRapidity_wrongpairs_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:243
std::vector< std::vector< TVector3 > > EMT_eta_gg_pair_momenta_Outside
Definition CbmKresEta.h:140
TH1D * TrueEta_AngleBetweenGammas_after_cuts_all_Target
Definition CbmKresEta.h:370
vector< TH1 * > fHistoList_rap_pt_Eta_all_Outside
Definition CbmKresEta.h:456
vector< TVector3 > VMomenta_minus_Target
Definition CbmKresEta.h:114
TH1D * GammasInvMass_fromEta_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:343
TH1D * GammasRapidity_fromEta_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:269
vector< CbmStsTrack * > VStsTrack_plus_Target
Definition CbmKresEta.h:121
TH1D * GammasRapidity_wrongpairs_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:257
TH1D * PlaneAngles_last_wrongpairs_twoInRich_Outside
Definition CbmKresEta.h:304
std::vector< CbmMCTrack * > fMCtracks
Definition CbmKresEta.h:131
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Outside
Definition CbmKresEta.h:181
vector< TH1 * > fHistoList_Eta_cuts_zeroInRich_Both
Definition CbmKresEta.h:310
TH1D * Particle_pull_py_Target
Definition CbmKresEta.h:206
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Both
Definition CbmKresEta.h:200
TH1D * WrongEta_InvMass_after_cuts_two_Target
Definition CbmKresEta.h:399
TH1D * GammasInvMass_fromEta_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:325
TH1D * PlaneAngles_first_wrongpairs_twoInRich_Target
Definition CbmKresEta.h:259
TH1D * multi_InvMass_Eta_all_Both_3
Definition CbmKresEta.h:481
TH2D * TrueEta_pt_vs_rap_onetwo_Both
Definition CbmKresEta.h:521
vector< TH1 * > fHistoList_rap_pt_Eta_two_Both
Definition CbmKresEta.h:584
vector< CbmStsTrack * > VStsTrack_plus_Outside
Definition CbmKresEta.h:103
std::vector< int > EMT_eta_gg_NofRings_Outside
Definition CbmKresEta.h:141
TH1D * InvMass_GammaCombinations_two_Both
Definition CbmKresEta.h:572
TH1D * WrongEta_AngleBetweenGammas_after_cuts_two_Both
Definition CbmKresEta.h:581
TH1D * multi_InvMass_Eta_all_Both_15
Definition CbmKresEta.h:493
TH1D * WrongEta_AngleBetweenGammas_after_cuts_onetwo_Outside
Definition CbmKresEta.h:438
std::vector< std::vector< TVector3 > > EMT_eta_gg_pair_momenta_Both
Definition CbmKresEta.h:145
TH1D * InvMass_GammaCombinations_all_Outside
Definition CbmKresEta.h:414
TH1D * multi_EMT_Eta_two_Both_7
Definition CbmKresEta.h:608
vector< TH1 * > fHistoList_Eta_cuts_Both
Definition CbmKresEta.h:308
TH1D * multi_InvMass_Eta_onetwo_Both_5
Definition CbmKresEta.h:536
std::vector< std::vector< int > > Gammas_stsIndex_all_Outside
Definition CbmKresEta.h:172
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:332
TH1D * multi_InvMass_Eta_all_Both_6
Definition CbmKresEta.h:484
void SaveOutsideTracks(CbmMCTrack *mcTrack1, CbmStsTrack *stsTrack, double charge, int stsInd, int richInd, int stsMcTrackId, CbmRichRing *RING)
TH1D * TrueEta_InvMass_after_cuts_onetwo_Target
Definition CbmKresEta.h:383
vector< TH1 * > fHistoList_Eta_cuts_twoInRich_Target
Definition CbmKresEta.h:247
TH1D * TrueEta_AngleBetweenGammas_after_cuts_all_Outside
Definition CbmKresEta.h:422
TH1D * multi_InvMass_Eta_onetwo_Both_17
Definition CbmKresEta.h:548
TH1D * multi_EMT_Eta_onetwo_Both_1
Definition CbmKresEta.h:549
TH1D * WrongEta_InvMass_after_cuts_onetwo_Outside
Definition CbmKresEta.h:436
TH1D * TrueEta_InvMass_after_cuts_all_Target
Definition CbmKresEta.h:368
TH1D * WrongEta_InvMass_after_cuts_onetwo_Both
Definition CbmKresEta.h:526
TH1D * AngleBetweenGammas_wrongEta_before_cuts_all_Both
Definition CbmKresEta.h:465
TH1D * multi_EMT_Eta_all_Both_5
Definition CbmKresEta.h:500
vector< TH1 * > fHistoList_Eta_all_Both
Definition CbmKresEta.h:463
TH1D * multi_EMT_Eta_onetwo_Both_3
Definition CbmKresEta.h:551
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Outside
Definition CbmKresEta.h:174
vector< TH1 * > fHistoList_rap_pt_Eta_onetwo_Target
Definition CbmKresEta.h:405
TH1D * GammasRapidity_wrongpairs_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:334
TH1D * multi_InvMass_Eta_onetwo_Both_16
Definition CbmKresEta.h:547
TH1D * WrongEta_AngleBetweenGammas_after_cuts_two_Outside
Definition CbmKresEta.h:453
std::vector< std::vector< int > > Gammas_stsIndex_two_Both
Definition CbmKresEta.h:191
TH1D * AngleBetweenGammas_wrongEta_before_cuts_two_Both
Definition CbmKresEta.h:571
std::vector< std::vector< int > > Gammas_MCIndex_two_Both
Definition CbmKresEta.h:195
TH1D * GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:272
TH1D * multi_EMT_Eta_two_Both_3
Definition CbmKresEta.h:604
TH1D * Particles_PDG_all_Outside
Definition CbmKresEta.h:415
TH1D * GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:312
TH2D * WrongEta_pt_vs_rap_est_two_Outside
Definition CbmKresEta.h:449
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:350
TH1D * GammasRapidity_fromEta_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:346
std::vector< std::vector< TVector3 > > Gammas_all_Outside
Definition CbmKresEta.h:168
TH2D * TrueEta_pt_vs_rap_two_Outside
Definition CbmKresEta.h:446
TH1D * PlaneAngles_first_wrongpairs_zeroInRich_Target
Definition CbmKresEta.h:231
TH2D * WrongEta_pt_vs_rap_onetwo_Both
Definition CbmKresEta.h:522
vector< TVector3 > frefmomenta
Definition CbmKresEta.h:128
TH1D * multi_InvMass_Eta_onetwo_Both_9
Definition CbmKresEta.h:540
vector< TH1 * > fHistoList_rap_pt_Eta_two_Outside
Definition CbmKresEta.h:458
std::vector< int > VRings_minus_Target
Definition CbmKresEta.h:115
std::vector< std::vector< TVector3 > > Gammas_two_Target
Definition CbmKresEta.h:151
vector< TH1 * > fHistoList_Eta_cuts_zeroInRich_Outside
Definition CbmKresEta.h:265
TH1D * multi_EMT_Eta_two_Both_10
Definition CbmKresEta.h:611
TH1D * EMT_Eta_InvMass_onetwo_Outside
Definition CbmKresEta.h:439
TH1D * WrongEta_AngleBetweenGammas_after_cuts_all_Both
Definition CbmKresEta.h:475
TH1D * multi_InvMass_Eta_onetwo_Both_2
Definition CbmKresEta.h:533
TH2D * WrongEta_pt_vs_rap_two_Both
Definition CbmKresEta.h:575
TH1D * multi_EMT_Eta_two_Both_15
Definition CbmKresEta.h:616
void Mixing_Both(vector< TH1 * > rap_pt_separation_all, vector< TH1 * > rap_pt_separation_onetwo, vector< TH1 * > rap_pt_separation_two)
std::vector< std::vector< int > > Gammas_stsIndex_all_Target
Definition CbmKresEta.h:154
TH1D * PlaneAngles_first_fromEta_twoInRich_Both
Definition CbmKresEta.h:348
TH1D * GammasMomentum_wrongpairs_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:302
TH1D * GammasRapidity_fromEta_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:283
TH1D * multi_EMT_Eta_all_Both_1
Definition CbmKresEta.h:496
void SaveTargetTracks(CbmMCTrack *mcTrack1, CbmStsTrack *stsTrack, TVector3 refmom, double charge, int stsInd, int richInd, int stsMcTrackId, CbmRichRing *RING)
std::vector< int > VStsIndex_minus_Outside
Definition CbmKresEta.h:99
TH1D * PlaneAngles_first_wrongpairs_twoInRich_Outside
Definition CbmKresEta.h:305
TH1D * EMT_Eta_InvMass_two_Outside
Definition CbmKresEta.h:454
TH2D * TrueEta_pt_vs_rap_est_onetwo_Outside
Definition CbmKresEta.h:433
TH2D * TrueEta_pt_vs_rap_two_Target
Definition CbmKresEta.h:394
void FindGammasBoth()
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Target
Definition CbmKresEta.h:156
TH2D * WrongEta_pt_vs_rap_all_Target
Definition CbmKresEta.h:365
TH1D * AngleBetweenGammas_wrongEta_before_cuts_all_Outside
Definition CbmKresEta.h:413
vector< TH1 * > fHistoList_Eta_cuts_twoInRich_Both
Definition CbmKresEta.h:342
std::vector< std::vector< int > > Gammas_MCIndex_onetwo_Outside
Definition CbmKresEta.h:178
TH1D * multi_InvMass_Eta_onetwo_Both_13
Definition CbmKresEta.h:544
std::vector< int > VStsIndex_plus_Target
Definition CbmKresEta.h:124
vector< TH1 * > fHistoList_Eta
Definition CbmKresEta.h:204
TH2D * TrueEta_pt_vs_rap_all_Both
Definition CbmKresEta.h:468
TH1D * PlaneAngles_last_wrongpairs_zeroInRich_Both
Definition CbmKresEta.h:321
TClonesArray * fGlobalTracks
Definition CbmKresEta.h:79
TH1D * PlaneAngles_last_wrongpairs_oneInRich_Outside
Definition CbmKresEta.h:290
TH1D * multi_InvMass_Eta_two_Both_16
Definition CbmKresEta.h:600
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Target
Definition CbmKresEta.h:164
TH1D * TrueEta_InvMass_after_cuts_all_Both
Definition CbmKresEta.h:472
TH1D * GammasInvMass_wrongpairs_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:240
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Target
Definition CbmKresEta.h:162
TH1D * PlaneAngles_last_wrongpairs_oneInRich_Both
Definition CbmKresEta.h:335
vector< TH1 * > fHistoList_rap_pt_Eta_onetwo_Outside
Definition CbmKresEta.h:457
TH1D * GammasMomentum_fromEta_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:222
TH1D * GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:267
TH1D * GammasOpenAngle_fromEta_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:344
TH2D * WrongEta_pt_vs_rap_est_all_Target
Definition CbmKresEta.h:367
TClonesArray * fRichRingMatches
Definition CbmKresEta.h:84
TH1D * multi_EMT_Eta_two_Both_13
Definition CbmKresEta.h:614
TH1D * multi_InvMass_Eta_two_Both_17
Definition CbmKresEta.h:601
TH1D * multi_InvMass_Eta_two_Both_2
Definition CbmKresEta.h:586
TH1D * WrongEta_AngleBetweenGammas_after_cuts_onetwo_Target
Definition CbmKresEta.h:386
TH1D * AngleBetweenGammas_trueEta_before_cuts_two_Both
Definition CbmKresEta.h:570
TH1D * TrueEta_InvMass_after_cuts_two_Target
Definition CbmKresEta.h:398
TH1D * EMT_Eta_InvMass_all_Both
Definition CbmKresEta.h:476
TH1D * GammasInvMass_fromEta_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:294
TH1D * multi_EMT_Eta_two_Both_11
Definition CbmKresEta.h:612
TH2D * WrongEta_pt_vs_rap_est_onetwo_Target
Definition CbmKresEta.h:382
TH1D * PlaneAngles_last_fromEta_twoInRich_Outside
Definition CbmKresEta.h:298
TH1D * multi_EMT_Eta_onetwo_Both_8
Definition CbmKresEta.h:556
TH1D * multi_InvMass_Eta_two_Both_3
Definition CbmKresEta.h:587
TH1D * GammasInvMass_wrongpairs_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:349
TH1D * GammasMomentum_fromEta_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:236
TH1D * GammasMomentum_fromEta_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:327
TH1D * AngleBetweenGammas_wrongEta_before_cuts_onetwo_Both
Definition CbmKresEta.h:518
TH1D * multi_EMT_Eta_all_Both_11
Definition CbmKresEta.h:506
TH1D * multi_InvMass_Eta_onetwo_Both_6
Definition CbmKresEta.h:537
static Double_t CalculateOpeningAngleBetweenGammas_Reco(TVector3 electron1, TVector3 electron2, TVector3 electron3, TVector3 electron4)
static TVector3 FitToVertexAndGetChi(CbmStsTrack *stsTrack, double x, double y, double z, double &chi)
static Double_t CalculateOpeningAngle_Reco(TVector3 electron1, TVector3 electron2)
static LmvmKinePar CalculateKinematicParams_4particles(const TVector3 part1, const TVector3 part2, const TVector3 part3, const TVector3 part4)
static double Invmass_2particles_RECO(const TVector3 part1, const TVector3 part2)
static LmvmKinePar CalculateKinematicParamsReco(const TVector3 electron1, const TVector3 electron2)
double GetPz() const
Definition CbmMCTrack.h:72
double GetPx() const
Definition CbmMCTrack.h:70
double GetStartZ() const
Definition CbmMCTrack.h:75
int32_t GetMotherId() const
Definition CbmMCTrack.h:69
double GetStartX() const
Definition CbmMCTrack.h:73
int32_t GetPdgCode() const
Definition CbmMCTrack.h:68
double GetStartY() const
Definition CbmMCTrack.h:74
double GetPy() const
Definition CbmMCTrack.h:71
int32_t GetNofLinks() const
Definition CbmMatch.h:42
const CbmLink & GetMatchedLink() const
Definition CbmMatch.h:41
double GetY() const
Definition CbmPixelHit.h:74
double GetX() const
Definition CbmPixelHit.h:73
Here the ring is fitted with Taubin algorithm from A. Ayriyan, G. Ososkov, N. Chernov.
virtual void DoFit(CbmRichRingLight *ring)
Inherited from CbmRichRingFitterBase.
float GetAaxis() const
float GetBaxis() const
void AddHit(CbmRichHitLight hit)
Add new hit to the ring.
uint32_t GetHit(int32_t i) const
Definition CbmRichRing.h:39
int32_t GetNofHits() const
Definition CbmRichRing.h:37
data class for a reconstructed 3-d hit in the STS
Definition CbmStsHit.h:35
int32_t GetNofMvdHits() const
Definition CbmStsTrack.h:87
int32_t GetMvdHitIndex(int32_t iHit) const
Definition CbmStsTrack.h:75
int32_t GetStsHitIndex(int32_t iHit) const
int32_t GetNofStsHits() const
Definition CbmStsTrack.h:93
const FairTrackParam * GetParamFirst() const
Definition CbmTrack.h:68
Double_t fMinv
Definition LmvmKinePar.h:22
Double_t fPt
Definition LmvmKinePar.h:20
Double_t fMomentumMag
Definition LmvmKinePar.h:19
Double_t fRapidity
Definition LmvmKinePar.h:21
Hash for CbmL1LinkKey.