CbmRoot
Loading...
Searching...
No Matches
CbmKresConversionManualmbias4.cxx
Go to the documentation of this file.
1/* Copyright (C) 2017-2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
2 SPDX-License-Identifier: GPL-3.0-only
3 Authors: Ievgenii Kres, Florian Uhlig [committer] */
4
24
25#include "CbmGlobalTrack.h"
27#include "CbmKresConversionBG.h"
28#include "CbmKresFunctions.h"
29#include "CbmMCTrack.h"
30#include "CbmMvdHit.h"
31#include "CbmRichHit.h"
32#include "CbmRichRing.h"
34#include "CbmRichRingLight.h"
35#include "CbmStsHit.h"
36#include "CbmStsTrack.h"
37#include "CbmTrackMatchNew.h"
38
39#include "FairRootManager.h"
40
41#include "TDirectory.h"
42
43#include <iostream>
44
45#include "KFParticle.h"
46
47
48using namespace std;
49
51 : fMcTracks(nullptr)
52 , fGlobalTracks(nullptr)
53 , fStsTracks(nullptr)
54 , fStsTrackMatches(nullptr)
55 , fRichProjections(nullptr)
56 , fRichRings(nullptr)
57 , fRichRingMatches(nullptr)
58 , fRichHits(nullptr)
59 , fArrayMvdHit(nullptr)
60 , fArrayStsHit(nullptr)
61 , fArrayCentrality(nullptr)
62 , fPrimVertex(nullptr)
63 , fKFVertex()
64 , fTauFit(nullptr)
65 , fAnaBG(nullptr)
66 , VStsTrack_minus_Outside_mbias4()
67 , VMCtracks_minus_Outside_mbias4()
68 , VRings_minus_Outside_mbias4()
69 , VStsIndex_minus_Outside_mbias4()
70 , VRichRing_minus_Outside_mbias4()
71 , VStsTrack_plus_Outside_mbias4()
72 , VMCtracks_plus_Outside_mbias4()
73 , VRings_plus_Outside_mbias4()
74 , VStsIndex_plus_Outside_mbias4()
75 , VRichRing_plus_Outside_mbias4()
76 , VMCtracks_minus_Target_mbias4()
77 , VStsTrack_minus_Target_mbias4()
78 , VMomenta_minus_Target_mbias4()
79 , VRings_minus_Target_mbias4()
80 , VStsIndex_minus_Target_mbias4()
81 , VRichRing_minus_Target_mbias4()
82 , VMCtracks_plus_Target_mbias4()
83 , VStsTrack_plus_Target_mbias4()
84 , VMomenta_plus_Target_mbias4()
85 , VRings_plus_Target_mbias4()
86 , VStsIndex_plus_Target_mbias4()
87 , VRichRing_plus_Target_mbias4()
88 , frefmomenta()
89 , frefId()
90 , fMCtracks()
91 , EMT_man_Event_Target_mbias4()
92 , EMT_man_pair_momenta_Target_mbias4()
93 , EMT_man_NofRings_Target_mbias4()
94 , EMT_man_Event_Outside_mbias4()
95 , EMT_man_pair_momenta_Outside_mbias4()
96 , EMT_man_NofRings_Outside_mbias4()
97 , EMT_man_Event_Both_mbias4()
98 , EMT_man_pair_momenta_Both_mbias4()
99 , EMT_man_NofRings_Both_mbias4()
100 , Gammas_all_Target_mbias4()
101 , Gammas_zero_Target_mbias4()
102 , Gammas_one_Target_mbias4()
103 , Gammas_two_Target_mbias4()
104 , Gammas_onetwo_Target_mbias4()
105 , Gammas_stsIndex_all_Target_mbias4()
106 , Gammas_stsIndex_zero_Target_mbias4()
107 , Gammas_stsIndex_one_Target_mbias4()
108 , Gammas_stsIndex_two_Target_mbias4()
109 , Gammas_stsIndex_onetwo_Target_mbias4()
110 , Gammas_MC_all_Target_mbias4()
111 , Gammas_MC_zero_Target_mbias4()
112 , Gammas_MC_one_Target_mbias4()
113 , Gammas_MC_two_Target_mbias4()
114 , Gammas_MC_onetwo_Target_mbias4()
115 , Gammas_all_Outside_mbias4()
116 , Gammas_zero_Outside_mbias4()
117 , Gammas_one_Outside_mbias4()
118 , Gammas_two_Outside_mbias4()
119 , Gammas_onetwo_Outside_mbias4()
120 , Gammas_stsIndex_all_Outside_mbias4()
121 , Gammas_stsIndex_zero_Outside_mbias4()
122 , Gammas_stsIndex_one_Outside_mbias4()
123 , Gammas_stsIndex_two_Outside_mbias4()
124 , Gammas_stsIndex_onetwo_Outside_mbias4()
125 , Gammas_MC_all_Outside_mbias4()
126 , Gammas_MC_zero_Outside_mbias4()
127 , Gammas_MC_one_Outside_mbias4()
128 , Gammas_MC_two_Outside_mbias4()
129 , Gammas_MC_onetwo_Outside_mbias4()
130 , Gammas_all_Both_mbias4()
131 , Gammas_zero_Both_mbias4()
132 , Gammas_one_Both_mbias4()
133 , Gammas_two_Both_mbias4()
134 , Gammas_onetwo_Both_mbias4()
135 , Gammas_stsIndex_all_Both_mbias4()
136 , Gammas_stsIndex_zero_Both_mbias4()
137 , Gammas_stsIndex_one_Both_mbias4()
138 , Gammas_stsIndex_two_Both_mbias4()
139 , Gammas_stsIndex_onetwo_Both_mbias4()
140 , Gammas_MC_all_Both_mbias4()
141 , Gammas_MC_zero_Both_mbias4()
142 , Gammas_MC_one_Both_mbias4()
143 , Gammas_MC_two_Both_mbias4()
144 , Gammas_MC_onetwo_Both_mbias4()
145 , fHistoList_man_cuts_Both_mbias4()
146 , InvMass_vs_OA_candidates_Both_mbias4(nullptr)
147 , InvMass_vs_OA_fromPi0_Both_mbias4(nullptr)
148 , GammasInvMass_candidates_Both_mbias4(nullptr)
149 , GammasOA_candidates_Both_mbias4(nullptr)
150 , GammasInvMass_fromPi0_Both_mbias4(nullptr)
151 , GammasOA_fromPi0_Both_mbias4(nullptr)
152 , PlaneAngles_last_candidates_Both_mbias4(nullptr)
153 , PlaneAngles_last_fromPi0_Both_mbias4(nullptr)
154 , PlaneAngles_first_candidates_Both_mbias4(nullptr)
155 , PlaneAngles_first_fromPi0_Both_mbias4(nullptr)
156 , fHistoList_man_cuts_Target_mbias4()
157 , InvMass_vs_OA_candidates_Target_mbias4(nullptr)
158 , InvMass_vs_OA_fromPi0_Target_mbias4(nullptr)
159 , GammasInvMass_candidates_Target_mbias4(nullptr)
160 , GammasOA_candidates_Target_mbias4(nullptr)
161 , GammasInvMass_fromPi0_Target_mbias4(nullptr)
162 , GammasOA_fromPi0_Target_mbias4(nullptr)
163 , PlaneAngles_last_candidates_Target_mbias4(nullptr)
164 , PlaneAngles_last_fromPi0_Target_mbias4(nullptr)
165 , PlaneAngles_first_candidates_Target_mbias4(nullptr)
166 , PlaneAngles_first_fromPi0_Target_mbias4(nullptr)
167 , fHistoList_man_cuts_Outside_mbias4()
168 , InvMass_vs_OA_candidates_Outside_mbias4(nullptr)
169 , InvMass_vs_OA_fromPi0_Outside_mbias4(nullptr)
170 , GammasInvMass_candidates_Outside_mbias4(nullptr)
171 , GammasOA_candidates_Outside_mbias4(nullptr)
172 , GammasInvMass_fromPi0_Outside_mbias4(nullptr)
173 , GammasOA_fromPi0_Outside_mbias4(nullptr)
174 , PlaneAngles_last_candidates_Outside_mbias4(nullptr)
175 , PlaneAngles_last_fromPi0_Outside_mbias4(nullptr)
176 , PlaneAngles_first_candidates_Outside_mbias4(nullptr)
177 , PlaneAngles_first_fromPi0_Outside_mbias4(nullptr)
178 , fHistoList_man_all_Target_mbias4()
179 , GammaInvMassReco_all_Target_mbias4(nullptr)
180 , GammaOpeningAngleReco_all_Target_mbias4(nullptr)
181 , Pdg_all_Target_mbias4(nullptr)
182 , P_reco_all_Target_mbias4(nullptr)
183 , Pt_reco_all_Target_mbias4(nullptr)
184 , Pi0InvMassReco_all_Target_mbias4(nullptr)
185 , EMT_InvMass_all_Target_mbias4(nullptr)
186 , Pi0_pt_vs_rap_all_Target_mbias4(nullptr)
187 , Pi0_pt_vs_rap_est_all_Target_mbias4(nullptr)
188 , fHistoList_man_zero_Target_mbias4()
189 , GammaInvMassReco_zero_Target_mbias4(nullptr)
190 , GammaOpeningAngleReco_zero_Target_mbias4(nullptr)
191 , Pdg_zero_Target_mbias4(nullptr)
192 , P_reco_zero_Target_mbias4(nullptr)
193 , Pt_reco_zero_Target_mbias4(nullptr)
194 , Pi0InvMassReco_zero_Target_mbias4(nullptr)
195 , EMT_InvMass_zero_Target_mbias4(nullptr)
196 , Pi0_pt_vs_rap_zero_Target_mbias4(nullptr)
197 , Pi0_pt_vs_rap_est_zero_Target_mbias4(nullptr)
198 , fHistoList_man_one_Target_mbias4()
199 , GammaInvMassReco_one_Target_mbias4(nullptr)
200 , GammaOpeningAngleReco_one_Target_mbias4(nullptr)
201 , Pdg_one_Target_mbias4(nullptr)
202 , P_reco_one_Target_mbias4(nullptr)
203 , Pt_reco_one_Target_mbias4(nullptr)
204 , Pi0InvMassReco_one_Target_mbias4(nullptr)
205 , EMT_InvMass_one_Target_mbias4(nullptr)
206 , Pi0_pt_vs_rap_one_Target_mbias4(nullptr)
207 , Pi0_pt_vs_rap_est_one_Target_mbias4(nullptr)
208 , fHistoList_man_two_Target_mbias4()
209 , GammaInvMassReco_two_Target_mbias4(nullptr)
210 , GammaOpeningAngleReco_two_Target_mbias4(nullptr)
211 , Pdg_two_Target_mbias4(nullptr)
212 , P_reco_two_Target_mbias4(nullptr)
213 , Pt_reco_two_Target_mbias4(nullptr)
214 , Pi0InvMassReco_two_Target_mbias4(nullptr)
215 , EMT_InvMass_two_Target_mbias4(nullptr)
216 , Pi0_pt_vs_rap_two_Target_mbias4(nullptr)
217 , Pi0_pt_vs_rap_est_two_Target_mbias4(nullptr)
218 , fHistoList_man_onetwo_Target_mbias4()
219 , GammaInvMassReco_onetwo_Target_mbias4(nullptr)
220 , GammaOpeningAngleReco_onetwo_Target_mbias4(nullptr)
221 , Pdg_onetwo_Target_mbias4(nullptr)
222 , P_reco_onetwo_Target_mbias4(nullptr)
223 , Pt_reco_onetwo_Target_mbias4(nullptr)
224 , Pi0InvMassReco_onetwo_Target_mbias4(nullptr)
225 , EMT_InvMass_onetwo_Target_mbias4(nullptr)
226 , Pi0_pt_vs_rap_onetwo_Target_mbias4(nullptr)
227 , Pi0_pt_vs_rap_est_onetwo_Target_mbias4(nullptr)
228 , fHistoList_man_all_Outside_mbias4()
229 , GammaInvMassReco_all_Outside_mbias4(nullptr)
230 , GammaOpeningAngleReco_all_Outside_mbias4(nullptr)
231 , Pdg_all_Outside_mbias4(nullptr)
232 , P_reco_all_Outside_mbias4(nullptr)
233 , Pt_reco_all_Outside_mbias4(nullptr)
234 , Pi0InvMassReco_all_Outside_mbias4(nullptr)
235 , EMT_InvMass_all_Outside_mbias4(nullptr)
236 , Pi0_pt_vs_rap_all_Outside_mbias4(nullptr)
237 , Pi0_pt_vs_rap_est_all_Outside_mbias4(nullptr)
238 , fHistoList_man_zero_Outside_mbias4()
239 , GammaInvMassReco_zero_Outside_mbias4(nullptr)
240 , GammaOpeningAngleReco_zero_Outside_mbias4(nullptr)
241 , Pdg_zero_Outside_mbias4(nullptr)
242 , P_reco_zero_Outside_mbias4(nullptr)
243 , Pt_reco_zero_Outside_mbias4(nullptr)
244 , Pi0InvMassReco_zero_Outside_mbias4(nullptr)
245 , EMT_InvMass_zero_Outside_mbias4(nullptr)
246 , Pi0_pt_vs_rap_zero_Outside_mbias4(nullptr)
247 , Pi0_pt_vs_rap_est_zero_Outside_mbias4(nullptr)
248 , fHistoList_man_one_Outside_mbias4()
249 , GammaInvMassReco_one_Outside_mbias4(nullptr)
250 , GammaOpeningAngleReco_one_Outside_mbias4(nullptr)
251 , Pdg_one_Outside_mbias4(nullptr)
252 , P_reco_one_Outside_mbias4(nullptr)
253 , Pt_reco_one_Outside_mbias4(nullptr)
254 , Pi0InvMassReco_one_Outside_mbias4(nullptr)
255 , EMT_InvMass_one_Outside_mbias4(nullptr)
256 , Pi0_pt_vs_rap_one_Outside_mbias4(nullptr)
257 , Pi0_pt_vs_rap_est_one_Outside_mbias4(nullptr)
258 , fHistoList_man_two_Outside_mbias4()
259 , GammaInvMassReco_two_Outside_mbias4(nullptr)
260 , GammaOpeningAngleReco_two_Outside_mbias4(nullptr)
261 , Pdg_two_Outside_mbias4(nullptr)
262 , P_reco_two_Outside_mbias4(nullptr)
263 , Pt_reco_two_Outside_mbias4(nullptr)
264 , Pi0InvMassReco_two_Outside_mbias4(nullptr)
265 , EMT_InvMass_two_Outside_mbias4(nullptr)
266 , Pi0_pt_vs_rap_two_Outside_mbias4(nullptr)
267 , Pi0_pt_vs_rap_est_two_Outside_mbias4(nullptr)
268 , fHistoList_man_onetwo_Outside_mbias4()
269 , GammaInvMassReco_onetwo_Outside_mbias4(nullptr)
270 , GammaOpeningAngleReco_onetwo_Outside_mbias4(nullptr)
271 , Pdg_onetwo_Outside_mbias4(nullptr)
272 , P_reco_onetwo_Outside_mbias4(nullptr)
273 , Pt_reco_onetwo_Outside_mbias4(nullptr)
274 , Pi0InvMassReco_onetwo_Outside_mbias4(nullptr)
275 , EMT_InvMass_onetwo_Outside_mbias4(nullptr)
276 , Pi0_pt_vs_rap_onetwo_Outside_mbias4(nullptr)
277 , Pi0_pt_vs_rap_est_onetwo_Outside_mbias4(nullptr)
278 , fHistoList_man_all_Both_mbias4()
279 , GammaInvMassReco_all_Both_mbias4(nullptr)
280 , GammaOpeningAngleReco_all_Both_mbias4(nullptr)
281 , Pdg_all_Both_mbias4(nullptr)
282 , P_reco_all_Both_mbias4(nullptr)
283 , Pt_reco_all_Both_mbias4(nullptr)
284 , Pi0InvMassReco_all_Both_mbias4(nullptr)
285 , EMT_InvMass_all_Both_mbias4(nullptr)
286 , Pi0_pt_vs_rap_all_Both_mbias4(nullptr)
287 , Pi0_pt_vs_rap_est_all_Both_mbias4(nullptr)
288 , fHistoList_man_zero_Both_mbias4()
289 , GammaInvMassReco_zero_Both_mbias4(nullptr)
290 , GammaOpeningAngleReco_zero_Both_mbias4(nullptr)
291 , Pdg_zero_Both_mbias4(nullptr)
292 , P_reco_zero_Both_mbias4(nullptr)
293 , Pt_reco_zero_Both_mbias4(nullptr)
294 , Pi0InvMassReco_zero_Both_mbias4(nullptr)
295 , EMT_InvMass_zero_Both_mbias4(nullptr)
296 , Pi0_pt_vs_rap_zero_Both_mbias4(nullptr)
297 , Pi0_pt_vs_rap_est_zero_Both_mbias4(nullptr)
298 , fHistoList_man_one_Both_mbias4()
299 , GammaInvMassReco_one_Both_mbias4(nullptr)
300 , GammaOpeningAngleReco_one_Both_mbias4(nullptr)
301 , Pdg_one_Both_mbias4(nullptr)
302 , P_reco_one_Both_mbias4(nullptr)
303 , Pt_reco_one_Both_mbias4(nullptr)
304 , Pi0InvMassReco_one_Both_mbias4(nullptr)
305 , EMT_InvMass_one_Both_mbias4(nullptr)
306 , Pi0_pt_vs_rap_one_Both_mbias4(nullptr)
307 , Pi0_pt_vs_rap_est_one_Both_mbias4(nullptr)
308 , fHistoList_man_two_Both_mbias4()
309 , GammaInvMassReco_two_Both_mbias4(nullptr)
310 , GammaOpeningAngleReco_two_Both_mbias4(nullptr)
311 , Pdg_two_Both_mbias4(nullptr)
312 , P_reco_two_Both_mbias4(nullptr)
313 , Pt_reco_two_Both_mbias4(nullptr)
314 , Pi0InvMassReco_two_Both_mbias4(nullptr)
315 , EMT_InvMass_two_Both_mbias4(nullptr)
316 , Pi0_pt_vs_rap_two_Both_mbias4(nullptr)
317 , Pi0_pt_vs_rap_est_two_Both_mbias4(nullptr)
318 , fHistoList_man_onetwo_Both_mbias4()
319 , GammaInvMassReco_onetwo_Both_mbias4(nullptr)
320 , GammaOpeningAngleReco_onetwo_Both_mbias4(nullptr)
321 , Pdg_onetwo_Both_mbias4(nullptr)
322 , P_reco_onetwo_Both_mbias4(nullptr)
323 , Pt_reco_onetwo_Both_mbias4(nullptr)
324 , Pi0InvMassReco_onetwo_Both_mbias4(nullptr)
325 , EMT_InvMass_onetwo_Both_mbias4(nullptr)
326 , Pi0_pt_vs_rap_onetwo_Both_mbias4(nullptr)
327 , Pi0_pt_vs_rap_est_onetwo_Both_mbias4(nullptr)
328 , fHistoList_man_Both_mbias4()
329 , Pdg_vs_Distance_mbias4(nullptr)
330 , P_vs_Distance_mbias4(nullptr)
331 , fHistoList_multiplicity_man_Target_mbias4()
332 , MultiplicityGamma_all_Target_mbias4(nullptr)
333 , MultiplicityGamma_zero_Target_mbias4(nullptr)
334 , MultiplicityGamma_one_Target_mbias4(nullptr)
335 , MultiplicityGamma_two_Target_mbias4(nullptr)
336 , MultiplicityGamma_onetwo_Target_mbias4(nullptr)
337 , MultiplicityChargedParticles_all_Target_mbias4(nullptr)
338 , MultiplicityChargedParticles_zero_Target_mbias4(nullptr)
339 , MultiplicityChargedParticles_one_Target_mbias4(nullptr)
340 , MultiplicityChargedParticles_two_Target_mbias4(nullptr)
341 , MultiplicityChargedParticles_onetwo_Target_mbias4(nullptr)
342 , fHistoList_multiplicity_man_Outside_mbias4()
343 , MultiplicityGamma_all_Outside_mbias4(nullptr)
344 , MultiplicityGamma_zero_Outside_mbias4(nullptr)
345 , MultiplicityGamma_one_Outside_mbias4(nullptr)
346 , MultiplicityGamma_two_Outside_mbias4(nullptr)
347 , MultiplicityGamma_onetwo_Outside_mbias4(nullptr)
348 , MultiplicityChargedParticles_all_Outside_mbias4(nullptr)
349 , MultiplicityChargedParticles_zero_Outside_mbias4(nullptr)
350 , MultiplicityChargedParticles_one_Outside_mbias4(nullptr)
351 , MultiplicityChargedParticles_two_Outside_mbias4(nullptr)
352 , MultiplicityChargedParticles_onetwo_Outside_mbias4(nullptr)
353 , fHistoList_multiplicity_man_Both_mbias4()
354 , MultiplicityGamma_all_Both_mbias4(nullptr)
355 , MultiplicityGamma_zero_Both_mbias4(nullptr)
356 , MultiplicityGamma_one_Both_mbias4(nullptr)
357 , MultiplicityGamma_two_Both_mbias4(nullptr)
358 , MultiplicityGamma_onetwo_Both_mbias4(nullptr)
359 , MultiplicityChargedParticles_all_Both_mbias4(nullptr)
360 , MultiplicityChargedParticles_zero_Both_mbias4(nullptr)
361 , MultiplicityChargedParticles_one_Both_mbias4(nullptr)
362 , MultiplicityChargedParticles_two_Both_mbias4(nullptr)
363 , MultiplicityChargedParticles_onetwo_Both_mbias4(nullptr)
364 , fHistoList_rap_vs_pt_InM_mbias4()
365 , rap_vs_Pt_InM_1_mbias4(nullptr)
366 , rap_vs_Pt_InM_2_mbias4(nullptr)
367 , rap_vs_Pt_InM_3_mbias4(nullptr)
368 , rap_vs_Pt_InM_4_mbias4(nullptr)
369 , rap_vs_Pt_InM_5_mbias4(nullptr)
370 , rap_vs_Pt_InM_6_mbias4(nullptr)
371 , rap_vs_Pt_InM_7_mbias4(nullptr)
372 , rap_vs_Pt_InM_8_mbias4(nullptr)
373 , rap_vs_Pt_InM_9_mbias4(nullptr)
374 , rap_vs_Pt_InM_10_mbias4(nullptr)
375 , rap_vs_Pt_InM_11_mbias4(nullptr)
376 , rap_vs_Pt_InM_12_mbias4(nullptr)
377 , rap_vs_Pt_InM_13_mbias4(nullptr)
378 , rap_vs_Pt_InM_14_mbias4(nullptr)
379 , rap_vs_Pt_InM_15_mbias4(nullptr)
380 , rap_vs_Pt_InM_16_mbias4(nullptr)
381 , rap_vs_Pt_InM_17_mbias4(nullptr)
382 , rap_vs_Pt_InM_18_mbias4(nullptr)
383 , rap_vs_Pt_InM_19_mbias4(nullptr)
384 , rap_vs_Pt_InM_20_mbias4(nullptr)
385 , rap_vs_Pt_InM_21_mbias4(nullptr)
386 , rap_vs_Pt_InM_22_mbias4(nullptr)
387 , rap_vs_Pt_InM_23_mbias4(nullptr)
388 , rap_vs_Pt_InM_24_mbias4(nullptr)
389 , rap_vs_Pt_InM_25_mbias4(nullptr)
390 , rap_vs_Pt_InM_26_mbias4(nullptr)
391 , rap_vs_Pt_InM_27_mbias4(nullptr)
392 , rap_vs_Pt_InM_28_mbias4(nullptr)
393 , rap_vs_Pt_InM_29_mbias4(nullptr)
394 , rap_vs_Pt_InM_30_mbias4(nullptr)
395 , rap_vs_Pt_InM_mixing_1_mbias4(nullptr)
396 , rap_vs_Pt_InM_mixing_2_mbias4(nullptr)
397 , rap_vs_Pt_InM_mixing_3_mbias4(nullptr)
398 , rap_vs_Pt_InM_mixing_4_mbias4(nullptr)
399 , rap_vs_Pt_InM_mixing_5_mbias4(nullptr)
400 , rap_vs_Pt_InM_mixing_6_mbias4(nullptr)
401 , rap_vs_Pt_InM_mixing_7_mbias4(nullptr)
402 , rap_vs_Pt_InM_mixing_8_mbias4(nullptr)
403 , rap_vs_Pt_InM_mixing_9_mbias4(nullptr)
404 , rap_vs_Pt_InM_mixing_10_mbias4(nullptr)
405 , rap_vs_Pt_InM_mixing_11_mbias4(nullptr)
406 , rap_vs_Pt_InM_mixing_12_mbias4(nullptr)
407 , rap_vs_Pt_InM_mixing_13_mbias4(nullptr)
408 , rap_vs_Pt_InM_mixing_14_mbias4(nullptr)
409 , rap_vs_Pt_InM_mixing_15_mbias4(nullptr)
410 , rap_vs_Pt_InM_mixing_16_mbias4(nullptr)
411 , rap_vs_Pt_InM_mixing_17_mbias4(nullptr)
412 , rap_vs_Pt_InM_mixing_18_mbias4(nullptr)
413 , rap_vs_Pt_InM_mixing_19_mbias4(nullptr)
414 , rap_vs_Pt_InM_mixing_20_mbias4(nullptr)
415 , rap_vs_Pt_InM_mixing_21_mbias4(nullptr)
416 , rap_vs_Pt_InM_mixing_22_mbias4(nullptr)
417 , rap_vs_Pt_InM_mixing_23_mbias4(nullptr)
418 , rap_vs_Pt_InM_mixing_24_mbias4(nullptr)
419 , rap_vs_Pt_InM_mixing_25_mbias4(nullptr)
420 , rap_vs_Pt_InM_mixing_26_mbias4(nullptr)
421 , rap_vs_Pt_InM_mixing_27_mbias4(nullptr)
422 , rap_vs_Pt_InM_mixing_28_mbias4(nullptr)
423 , rap_vs_Pt_InM_mixing_29_mbias4(nullptr)
424 , rap_vs_Pt_InM_mixing_30_mbias4(nullptr)
425 , fHistoList_bg_InM_all_Target_mbias4()
426 , fHistoList_bg_InM_zero_Target_mbias4()
427 , fHistoList_bg_InM_one_Target_mbias4()
428 , fHistoList_bg_InM_two_Target_mbias4()
429 , fHistoList_bg_InM_onetwo_Target_mbias4()
430 , fHistoList_bg_InM_all_Outside_mbias4()
431 , fHistoList_bg_InM_zero_Outside_mbias4()
432 , fHistoList_bg_InM_one_Outside_mbias4()
433 , fHistoList_bg_InM_two_Outside_mbias4()
434 , fHistoList_bg_InM_onetwo_Outside_mbias4()
435 , fHistoList_bg_InM_all_Both_mbias4()
436 , BG1_InM_all_Both_mbias4(nullptr)
437 , BG2_InM_all_Both_mbias4(nullptr)
438 , BG3_InM_all_Both_mbias4(nullptr)
439 , BG4_InM_all_Both_mbias4(nullptr)
440 , BG5_InM_all_Both_mbias4(nullptr)
441 , BG6_InM_all_Both_mbias4(nullptr)
442 , BG7_InM_all_Both_mbias4(nullptr)
443 , BG8_InM_all_Both_mbias4(nullptr)
444 , BG9_InM_all_Both_mbias4(nullptr)
445 , BG10_InM_all_Both_mbias4(nullptr)
446 , PdgCase8_InM_all_Both_mbias4(nullptr)
447 , PdgCase8mothers_InM_all_Both_mbias4(nullptr)
448 , sameMIDcase8_InM_all_Both_mbias4(nullptr)
449 , sameGRIDcase8_InM_all_Both_mbias4(nullptr)
450 , Case1ZYPos_InM_all_Both_mbias4(nullptr)
451 , sameMIDcase8_mothedPDG_InM_all_Both_mbias4(nullptr)
452 , PdgCase8NonEComeFromTarget_mbias4_InM_all_Both_mbias4(nullptr)
453 , PdgCase8NonE_NOT_FromTarget_mbias4_InM_all_Both_mbias4(nullptr)
454 , PdgCase8motherNonE_InM_all_Both_mbias4(nullptr)
455 , Case8ElFromDalitz_InM_all_Both_mbias4(nullptr)
456 , Case8NonElFrom_pn_InM_all_Both_mbias4(nullptr)
457 , Case8NonElFrom_eta_InM_all_Both_mbias4(nullptr)
458 , Case8NonElFrom_kaon_InM_all_Both_mbias4(nullptr)
459 , sameMIDcase8NonEPdg_InM_all_Both_mbias4(nullptr)
460 , sameMIDcase8NonEMotherPdg_InM_all_Both_mbias4(nullptr)
461 , sameMIDcase8NonEMotherIM_InM_all_Both_mbias4(nullptr)
462 , sameMIDcase8NonEPdgFromTarget_mbias4_InM_all_Both_mbias4(nullptr)
463 , sameMIDcase8NonEComeFromTarget_mbias4IM_InM_all_Both_mbias4(nullptr)
464 , sameMIDcase8NonEComeFromTarget_mbias4P_InM_all_Both_mbias4(nullptr)
465 , sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_all_Both_mbias4(nullptr)
466 , fHistoList_bg_InM_zero_Both_mbias4()
467 , BG1_InM_zero_Both_mbias4(nullptr)
468 , BG2_InM_zero_Both_mbias4(nullptr)
469 , BG3_InM_zero_Both_mbias4(nullptr)
470 , BG4_InM_zero_Both_mbias4(nullptr)
471 , BG5_InM_zero_Both_mbias4(nullptr)
472 , BG6_InM_zero_Both_mbias4(nullptr)
473 , BG7_InM_zero_Both_mbias4(nullptr)
474 , BG8_InM_zero_Both_mbias4(nullptr)
475 , BG9_InM_zero_Both_mbias4(nullptr)
476 , BG10_InM_zero_Both_mbias4(nullptr)
477 , PdgCase8_InM_zero_Both_mbias4(nullptr)
478 , PdgCase8mothers_InM_zero_Both_mbias4(nullptr)
479 , sameMIDcase8_InM_zero_Both_mbias4(nullptr)
480 , sameGRIDcase8_InM_zero_Both_mbias4(nullptr)
481 , Case1ZYPos_InM_zero_Both_mbias4(nullptr)
482 , sameMIDcase8_mothedPDG_InM_zero_Both_mbias4(nullptr)
483 , PdgCase8NonEComeFromTarget_mbias4_InM_zero_Both_mbias4(nullptr)
484 , PdgCase8NonE_NOT_FromTarget_mbias4_InM_zero_Both_mbias4(nullptr)
485 , PdgCase8motherNonE_InM_zero_Both_mbias4(nullptr)
486 , Case8ElFromDalitz_InM_zero_Both_mbias4(nullptr)
487 , Case8NonElFrom_pn_InM_zero_Both_mbias4(nullptr)
488 , Case8NonElFrom_eta_InM_zero_Both_mbias4(nullptr)
489 , Case8NonElFrom_kaon_InM_zero_Both_mbias4(nullptr)
490 , sameMIDcase8NonEPdg_InM_zero_Both_mbias4(nullptr)
491 , sameMIDcase8NonEMotherPdg_InM_zero_Both_mbias4(nullptr)
492 , sameMIDcase8NonEMotherIM_InM_zero_Both_mbias4(nullptr)
493 , sameMIDcase8NonEPdgFromTarget_mbias4_InM_zero_Both_mbias4(nullptr)
494 , sameMIDcase8NonEComeFromTarget_mbias4IM_InM_zero_Both_mbias4(nullptr)
495 , sameMIDcase8NonEComeFromTarget_mbias4P_InM_zero_Both_mbias4(nullptr)
496 , sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_zero_Both_mbias4(nullptr)
497 , BG1_InM_one_Both_mbias4(nullptr)
498 , BG2_InM_one_Both_mbias4(nullptr)
499 , BG3_InM_one_Both_mbias4(nullptr)
500 , BG4_InM_one_Both_mbias4(nullptr)
501 , BG5_InM_one_Both_mbias4(nullptr)
502 , BG6_InM_one_Both_mbias4(nullptr)
503 , BG7_InM_one_Both_mbias4(nullptr)
504 , BG8_InM_one_Both_mbias4(nullptr)
505 , BG9_InM_one_Both_mbias4(nullptr)
506 , BG10_InM_one_Both_mbias4(nullptr)
507 , PdgCase8_InM_one_Both_mbias4(nullptr)
508 , PdgCase8mothers_InM_one_Both_mbias4(nullptr)
509 , sameMIDcase8_InM_one_Both_mbias4(nullptr)
510 , sameGRIDcase8_InM_one_Both_mbias4(nullptr)
511 , Case1ZYPos_InM_one_Both_mbias4(nullptr)
512 , sameMIDcase8_mothedPDG_InM_one_Both_mbias4(nullptr)
513 , PdgCase8NonEComeFromTarget_mbias4_InM_one_Both_mbias4(nullptr)
514 , PdgCase8NonE_NOT_FromTarget_mbias4_InM_one_Both_mbias4(nullptr)
515 , PdgCase8motherNonE_InM_one_Both_mbias4(nullptr)
516 , Case8ElFromDalitz_InM_one_Both_mbias4(nullptr)
517 , Case8NonElFrom_pn_InM_one_Both_mbias4(nullptr)
518 , Case8NonElFrom_eta_InM_one_Both_mbias4(nullptr)
519 , Case8NonElFrom_kaon_InM_one_Both_mbias4(nullptr)
520 , sameMIDcase8NonEPdg_InM_one_Both_mbias4(nullptr)
521 , sameMIDcase8NonEMotherPdg_InM_one_Both_mbias4(nullptr)
522 , sameMIDcase8NonEMotherIM_InM_one_Both_mbias4(nullptr)
523 , sameMIDcase8NonEPdgFromTarget_mbias4_InM_one_Both_mbias4(nullptr)
524 , sameMIDcase8NonEComeFromTarget_mbias4IM_InM_one_Both_mbias4(nullptr)
525 , sameMIDcase8NonEComeFromTarget_mbias4P_InM_one_Both_mbias4(nullptr)
526 , sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_one_Both_mbias4(nullptr)
527 , fHistoList_bg_InM_two_Both_mbias4()
528 , BG1_InM_two_Both_mbias4(nullptr)
529 , BG2_InM_two_Both_mbias4(nullptr)
530 , BG3_InM_two_Both_mbias4(nullptr)
531 , BG4_InM_two_Both_mbias4(nullptr)
532 , BG5_InM_two_Both_mbias4(nullptr)
533 , BG6_InM_two_Both_mbias4(nullptr)
534 , BG7_InM_two_Both_mbias4(nullptr)
535 , BG8_InM_two_Both_mbias4(nullptr)
536 , BG9_InM_two_Both_mbias4(nullptr)
537 , BG10_InM_two_Both_mbias4(nullptr)
538 , PdgCase8_InM_two_Both_mbias4(nullptr)
539 , PdgCase8mothers_InM_two_Both_mbias4(nullptr)
540 , sameMIDcase8_InM_two_Both_mbias4(nullptr)
541 , sameGRIDcase8_InM_two_Both_mbias4(nullptr)
542 , Case1ZYPos_InM_two_Both_mbias4(nullptr)
543 , sameMIDcase8_mothedPDG_InM_two_Both_mbias4(nullptr)
544 , PdgCase8NonEComeFromTarget_mbias4_InM_two_Both_mbias4(nullptr)
545 , PdgCase8NonE_NOT_FromTarget_mbias4_InM_two_Both_mbias4(nullptr)
546 , PdgCase8motherNonE_InM_two_Both_mbias4(nullptr)
547 , Case8ElFromDalitz_InM_two_Both_mbias4(nullptr)
548 , Case8NonElFrom_pn_InM_two_Both_mbias4(nullptr)
549 , Case8NonElFrom_eta_InM_two_Both_mbias4(nullptr)
550 , Case8NonElFrom_kaon_InM_two_Both_mbias4(nullptr)
551 , sameMIDcase8NonEPdg_InM_two_Both_mbias4(nullptr)
552 , sameMIDcase8NonEMotherPdg_InM_two_Both_mbias4(nullptr)
553 , sameMIDcase8NonEMotherIM_InM_two_Both_mbias4(nullptr)
554 , sameMIDcase8NonEPdgFromTarget_mbias4_InM_two_Both_mbias4(nullptr)
555 , sameMIDcase8NonEComeFromTarget_mbias4IM_InM_two_Both_mbias4(nullptr)
556 , sameMIDcase8NonEComeFromTarget_mbias4P_InM_two_Both_mbias4(nullptr)
557 , sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_two_Both_mbias4(nullptr)
558 , fHistoList_bg_InM_onetwo_Both_mbias4()
559 , BG1_InM_onetwo_Both_mbias4(nullptr)
560 , BG2_InM_onetwo_Both_mbias4(nullptr)
561 , BG3_InM_onetwo_Both_mbias4(nullptr)
562 , BG4_InM_onetwo_Both_mbias4(nullptr)
563 , BG5_InM_onetwo_Both_mbias4(nullptr)
564 , BG6_InM_onetwo_Both_mbias4(nullptr)
565 , BG7_InM_onetwo_Both_mbias4(nullptr)
566 , BG8_InM_onetwo_Both_mbias4(nullptr)
567 , BG9_InM_onetwo_Both_mbias4(nullptr)
568 , BG10_InM_onetwo_Both_mbias4(nullptr)
569 , PdgCase8_InM_onetwo_Both_mbias4(nullptr)
570 , PdgCase8mothers_InM_onetwo_Both_mbias4(nullptr)
571 , sameMIDcase8_InM_onetwo_Both_mbias4(nullptr)
572 , sameGRIDcase8_InM_onetwo_Both_mbias4(nullptr)
573 , Case1ZYPos_InM_onetwo_Both_mbias4(nullptr)
574 , sameMIDcase8_mothedPDG_InM_onetwo_Both_mbias4(nullptr)
575 , PdgCase8NonEComeFromTarget_mbias4_InM_onetwo_Both_mbias4(nullptr)
576 , PdgCase8NonE_NOT_FromTarget_mbias4_InM_onetwo_Both_mbias4(nullptr)
577 , PdgCase8motherNonE_InM_onetwo_Both_mbias4(nullptr)
578 , Case8ElFromDalitz_InM_onetwo_Both_mbias4(nullptr)
579 , Case8NonElFrom_pn_InM_onetwo_Both_mbias4(nullptr)
580 , Case8NonElFrom_eta_InM_onetwo_Both_mbias4(nullptr)
581 , Case8NonElFrom_kaon_InM_onetwo_Both_mbias4(nullptr)
582 , sameMIDcase8NonEPdg_InM_onetwo_Both_mbias4(nullptr)
583 , sameMIDcase8NonEMotherPdg_InM_onetwo_Both_mbias4(nullptr)
584 , sameMIDcase8NonEMotherIM_InM_onetwo_Both_mbias4(nullptr)
585 , sameMIDcase8NonEPdgFromTarget_mbias4_InM_onetwo_Both_mbias4(nullptr)
586 , sameMIDcase8NonEComeFromTarget_mbias4IM_InM_onetwo_Both_mbias4(nullptr)
587 , sameMIDcase8NonEComeFromTarget_mbias4P_InM_onetwo_Both_mbias4(nullptr)
588 , sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_onetwo_Both_mbias4(nullptr)
589{
590}
591
593
595{
596 FairRootManager* ioman = FairRootManager::Instance();
597 if (nullptr == ioman) { Fatal("CbmKresConversionManualmbias4::Init", "RootManager not instantised!"); }
598
599 fMcTracks = (TClonesArray*) ioman->GetObject("MCTrack");
600 if (nullptr == fMcTracks) { Fatal("CbmKresConversionManualmbias4::Init", "No MCTrack array!"); }
601
602 fPrimVertex = dynamic_cast<CbmVertex*>(ioman->GetObject("PrimaryVertex."));
603 if (nullptr == fPrimVertex) { fPrimVertex = dynamic_cast<CbmVertex*>(ioman->GetObject("PrimaryVertex")); }
604 if (nullptr == fPrimVertex) { LOG(fatal) << "CbmKresConversionManualmbias4::Init No PrimaryVertex array!"; }
605
606 fGlobalTracks = (TClonesArray*) ioman->GetObject("GlobalTrack");
607 if (nullptr == fGlobalTracks) { Fatal("CbmKresConversionManualmbias4::Init", "No GlobalTrack array!"); }
608
609 fStsTracks = (TClonesArray*) ioman->GetObject("StsTrack");
610 if (nullptr == fStsTracks) { Fatal("CbmKresConversionManualmbias4::Init", "No StsTrack array!"); }
611
612 fStsTrackMatches = (TClonesArray*) ioman->GetObject("StsTrackMatch");
613 if (nullptr == fStsTrackMatches) { Fatal("CbmKresConversionManualmbias4::Init", "No StsTrackMatch array!"); }
614
615 fRichProjections = (TClonesArray*) ioman->GetObject("RichProjection");
616 if (nullptr == fRichProjections) { Fatal("CbmKresConversionManualmbias4::Init", "No RichProjection array!"); }
617
618 fRichRings = (TClonesArray*) ioman->GetObject("RichRing");
619 if (nullptr == fRichRings) { Fatal("CbmKresConversionManualmbias4::Init", "No RichRing array!"); }
620
621 fRichRingMatches = (TClonesArray*) ioman->GetObject("RichRingMatch");
622 if (nullptr == fRichRingMatches) { Fatal("CbmKresConversionManualmbias4::Init", "No RichRingMatch array!"); }
623
624 fRichHits = (TClonesArray*) ioman->GetObject("RichHit");
625 if (nullptr == fRichHits) { Fatal("CbmKresConversionManualmbias4::Init", "No RichHit array!"); }
626
627 fArrayMvdHit = (TClonesArray*) ioman->GetObject("MvdHit");
628 if (nullptr == fArrayMvdHit) { Fatal("CbmKresConversionManualmbias4::Init", "No MvdHit array!"); }
629
630 fArrayStsHit = (TClonesArray*) ioman->GetObject("StsHit");
631 if (nullptr == fArrayStsHit) { Fatal("CbmKresConversionManualmbias4::Init", "No StsHit array!"); }
632
633 fArrayCentrality = (FairMCEventHeader*) ioman->GetObject("MCEventHeader.");
634 if (nullptr == fArrayCentrality) { Fatal("CbmAnaConversion2::Init", "No fArrayCentrality array!"); }
635
636
638
640
642 fAnaBG->Init();
643}
644
645void CbmKresConversionManualmbias4::Exec(int fEventNumMan, double OpeningAngleCut, double GammaInvMassCut, int RealPID)
646{
647 cout << "CbmKresConversionManualmbias4, event No. " << fEventNumMan << endl;
648
649 double centrality_mbias4 = fArrayCentrality->GetB();
650
651 if (fEventNumMan == 1000) {
666 Mixing_Both();
670
675
680 }
681 if (centrality_mbias4 < 10) return;
682
683
684 if (fPrimVertex != nullptr) { fKFVertex = CbmKFVertex(*fPrimVertex); }
685 else {
686 return;
687 }
688
689
735
736
747
760
761
762 Int_t ngTracks = fGlobalTracks->GetEntriesFast();
763 for (Int_t i = 0; i < ngTracks; i++) {
764 CbmGlobalTrack* gTrack = (CbmGlobalTrack*) fGlobalTracks->At(i);
765 if (nullptr == gTrack) continue;
766 Int_t stsInd = gTrack->GetStsTrackIndex();
767 Int_t richInd = gTrack->GetRichRingIndex();
768
769 if (stsInd < 0) continue;
770 CbmStsTrack* stsTrack = (CbmStsTrack*) fStsTracks->At(stsInd);
771 if (stsTrack == nullptr) continue;
772 CbmTrackMatchNew* stsMatch = (CbmTrackMatchNew*) fStsTrackMatches->At(stsInd);
773 if (stsMatch == nullptr) continue;
774 if (stsMatch->GetNofLinks() <= 0) continue;
775 Int_t stsMcTrackId = stsMatch->GetMatchedLink().GetIndex();
776 if (stsMcTrackId < 0) continue;
777 CbmMCTrack* mcTrack = (CbmMCTrack*) fMcTracks->At(stsMcTrackId);
778 if (mcTrack == nullptr) continue;
779
780
781 FairTrackParam* proj = (FairTrackParam*) fRichProjections->At(i);
782 if (richInd < 0 && proj->GetX() > -115 && proj->GetX() < 115
783 && ((proj->GetY() < -120 && proj->GetY() > -200) || (proj->GetY() > 120 && proj->GetY() < 200)))
784 continue;
785 CbmRichRing* Ring = nullptr;
786 if (richInd > -1) {
787 Ring = static_cast<CbmRichRing*>(fRichRings->At(richInd));
788 //Pdg_vs_Distance_mbias4->Fill(TMath::Abs(mcTrack->GetPdgCode()), Ring->GetDistance());
789 //if (TMath::Abs(mcTrack->GetPdgCode()) == 11) P_vs_Distance_mbias4->Fill(mcTrack->GetP(), Ring->GetDistance());
790 }
791
792
793 // Doing fit with Fit To primary Vertex and calculate chi2 to primary vertex
794 double chi2 = 0;
796 fKFVertex.GetRefZ(), chi2);
797 const FairTrackParam* track_par = stsTrack->GetParamFirst();
798 double charge = track_par->GetQp();
799
800
801 if (chi2 != chi2) continue;
802 if (chi2 == 0) continue;
803
804 if (chi2 > 3) { SaveOutsideTracks(mcTrack, stsTrack, charge, stsInd, richInd, stsMcTrackId, Ring); }
805 if (chi2 > 3) continue;
806
807 SaveTargetTracks(mcTrack, stsTrack, Momentum, charge, stsInd, richInd, stsMcTrackId, Ring);
808 }
809
810
811 FindGammasTarget(fEventNumMan, OpeningAngleCut, GammaInvMassCut, RealPID, VMCtracks_minus_Target_mbias4,
816
817 FindGammasOutside(fEventNumMan, OpeningAngleCut, GammaInvMassCut, RealPID, VMCtracks_minus_Outside_mbias4,
821
823
824
845
866
887}
888
889
891 int stsInd, int richInd, int stsMcTrackId, CbmRichRing* RING)
892{
893 int InRich = FindInRich(richInd, stsMcTrackId);
894 if (charge < 0) {
895 VMCtracks_minus_Outside_mbias4.push_back(mcTrack1);
896 VStsTrack_minus_Outside_mbias4.push_back(stsTrack);
897 VRings_minus_Outside_mbias4.push_back(InRich);
898 VStsIndex_minus_Outside_mbias4.push_back(stsInd);
899 VRichRing_minus_Outside_mbias4.push_back(RING);
900 }
901 if (charge > 0) {
902 VMCtracks_plus_Outside_mbias4.push_back(mcTrack1);
903 VStsTrack_plus_Outside_mbias4.push_back(stsTrack);
904 VRings_plus_Outside_mbias4.push_back(InRich);
905 VStsIndex_plus_Outside_mbias4.push_back(stsInd);
906 VRichRing_plus_Outside_mbias4.push_back(RING);
907 }
908}
909
911 double charge, int stsInd, int richInd, int stsMcTrackId,
912 CbmRichRing* RING)
913{
914 int InRich = FindInRich(richInd, stsMcTrackId);
915 if (charge < 0) {
916 VMCtracks_minus_Target_mbias4.push_back(mcTrack1);
917 VStsTrack_minus_Target_mbias4.push_back(stsTrack);
918 VMomenta_minus_Target_mbias4.push_back(refmom);
919 VRings_minus_Target_mbias4.push_back(InRich);
920 VStsIndex_minus_Target_mbias4.push_back(stsInd);
921 VRichRing_minus_Target_mbias4.push_back(RING);
922 }
923 if (charge > 0) {
924 VMCtracks_plus_Target_mbias4.push_back(mcTrack1);
925 VStsTrack_plus_Target_mbias4.push_back(stsTrack);
926 VMomenta_plus_Target_mbias4.push_back(refmom);
927 VRings_plus_Target_mbias4.push_back(InRich);
928 VStsIndex_plus_Target_mbias4.push_back(stsInd);
929 VRichRing_plus_Target_mbias4.push_back(RING);
930 }
931}
932
933
935 int EventNumMan, double AngleCut, double InvMassCut, int RealPID, vector<CbmMCTrack*> MCtracks_minus,
936 vector<CbmMCTrack*> MCtracks_plus, vector<CbmStsTrack*> StsTrack_minus, vector<CbmStsTrack*> StsTrack_plus,
937 vector<TVector3> Momenta_minus, vector<TVector3> Momenta_plus, std::vector<int> Rings_minus,
938 std::vector<int> Rings_plus, std::vector<int> stsIndex_minus, std::vector<int> stsIndex_plus,
939 vector<CbmRichRing*> richRing_minus, vector<CbmRichRing*> richRing_plus)
940{
941 for (size_t i = 0; i < Momenta_minus.size(); i++) {
942 for (size_t j = 0; j < Momenta_plus.size(); j++) {
943
944 CbmStsTrack* sts1 = StsTrack_minus.at(i);
945 CbmStsTrack* sts2 = StsTrack_plus.at(j);
946
947 TVector3 part1 = Momenta_minus[i];
948 TVector3 part2 = Momenta_plus[j];
949 CbmMCTrack* part1MC = MCtracks_minus[i];
950 CbmMCTrack* part2MC = MCtracks_plus[j];
951
952 int richcheck_0 = 0;
953 int richcheck_1 = 0;
954 if (RealPID == 1) {
955 // Real PID
956 richcheck_0 = CheckIfElectron(richRing_minus[i], Momenta_minus[i].Mag());
957 richcheck_1 = CheckIfElectron(richRing_plus[j], Momenta_plus[j].Mag());
958 }
959 else {
960 // MC PID
961 richcheck_0 = Rings_minus[i];
962 richcheck_1 = Rings_plus[j];
963 }
964 int richcheck = richcheck_0 + richcheck_1;
965
966 Double_t InvmassReco = CbmKresFunctions::Invmass_2particles_RECO(part1, part2);
967 Double_t OpeningAngle = CbmKresFunctions::CalculateOpeningAngle_Reco(part1, part2);
968
969 InvMass_vs_OA_candidates_Target_mbias4->Fill(InvmassReco, OpeningAngle);
970 InvMass_vs_OA_candidates_Both_mbias4->Fill(InvmassReco, OpeningAngle);
972 GammasInvMass_candidates_Both_mbias4->Fill(InvmassReco);
973 GammasOA_candidates_Target_mbias4->Fill(OpeningAngle);
974 GammasOA_candidates_Both_mbias4->Fill(OpeningAngle);
975
976 double PlaneAngle_last = CalculatePlaneAngle_last(sts1, sts2);
977 double PlaneAngle_first = CalculatePlaneAngle_first(sts1, sts2);
978
979 //if (PlaneAngle_last > 20) continue;
980
981 if (part1MC->GetMotherId() == part2MC->GetMotherId() && part1MC->GetMotherId() != -1) {
982 CbmMCTrack* mcTrackmama = (CbmMCTrack*) fMcTracks->At(part1MC->GetMotherId());
983 if (mcTrackmama->GetMotherId() != -1 && mcTrackmama->GetPdgCode() == 22) {
984 CbmMCTrack* mcTrackgrmama = (CbmMCTrack*) fMcTracks->At(mcTrackmama->GetMotherId());
985 if (mcTrackgrmama->GetPdgCode() == 111) {
986 GammasInvMass_fromPi0_Target_mbias4->Fill(InvmassReco);
987 GammasInvMass_fromPi0_Both_mbias4->Fill(InvmassReco);
988 GammasOA_fromPi0_Target_mbias4->Fill(OpeningAngle);
989 GammasOA_fromPi0_Both_mbias4->Fill(OpeningAngle);
990 InvMass_vs_OA_fromPi0_Target_mbias4->Fill(InvmassReco, OpeningAngle);
991 InvMass_vs_OA_fromPi0_Both_mbias4->Fill(InvmassReco, OpeningAngle);
992 PlaneAngles_last_fromPi0_Target_mbias4->Fill(PlaneAngle_last);
993 PlaneAngles_last_fromPi0_Both_mbias4->Fill(PlaneAngle_last);
994 PlaneAngles_first_fromPi0_Target_mbias4->Fill(PlaneAngle_first);
995 PlaneAngles_first_fromPi0_Both_mbias4->Fill(PlaneAngle_first);
996 }
997 }
998 }
999
1000 PlaneAngles_last_candidates_Target_mbias4->Fill(PlaneAngle_last);
1001 PlaneAngles_last_candidates_Both_mbias4->Fill(PlaneAngle_last);
1002 PlaneAngles_first_candidates_Target_mbias4->Fill(PlaneAngle_first);
1003 PlaneAngles_first_candidates_Both_mbias4->Fill(PlaneAngle_first);
1004
1005 // cuts
1006 if (TMath::Abs(OpeningAngle) > AngleCut) continue;
1007 if (TMath::Abs(InvmassReco) > InvMassCut) continue;
1008
1009 frefmomenta.clear();
1010 frefmomenta.push_back(part1);
1011 frefmomenta.push_back(part2);
1012 frefId.clear();
1013 frefId.push_back(stsIndex_minus[i]);
1014 frefId.push_back(stsIndex_plus[j]);
1015 fMCtracks.clear();
1016 fMCtracks.push_back(part1MC);
1017 fMCtracks.push_back(part2MC);
1018
1019 // for event mixing
1020 EMT_man_Event_Target_mbias4.push_back(EventNumMan);
1022 EMT_man_NofRings_Target_mbias4.push_back(richcheck);
1023
1024
1025 // everything (RICH == 0, RICH == 1, RICH == 2) together
1026 if (richcheck == 0 || richcheck == 1 || richcheck == 2) {
1030 GammaInvMassReco_all_Target_mbias4->Fill(InvmassReco);
1031 GammaOpeningAngleReco_all_Target_mbias4->Fill(OpeningAngle);
1032 Pdg_all_Target_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1033 Pdg_all_Target_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1034 P_reco_all_Target_mbias4->Fill(part1.Mag());
1035 P_reco_all_Target_mbias4->Fill(part2.Mag());
1036 Pt_reco_all_Target_mbias4->Fill(part1.Perp());
1037 Pt_reco_all_Target_mbias4->Fill(part2.Perp());
1038 GammaInvMassReco_all_Both_mbias4->Fill(InvmassReco);
1039 GammaOpeningAngleReco_all_Both_mbias4->Fill(OpeningAngle);
1040 Pdg_all_Both_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1041 Pdg_all_Both_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1042 P_reco_all_Both_mbias4->Fill(part1.Mag());
1043 P_reco_all_Both_mbias4->Fill(part2.Mag());
1044 Pt_reco_all_Both_mbias4->Fill(part1.Perp());
1045 Pt_reco_all_Both_mbias4->Fill(part2.Perp());
1046 }
1047
1048 // only cases, when RICH == 0
1049 if (richcheck == 0) {
1053 GammaInvMassReco_zero_Target_mbias4->Fill(InvmassReco);
1054 GammaOpeningAngleReco_zero_Target_mbias4->Fill(OpeningAngle);
1055 Pdg_zero_Target_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1056 Pdg_zero_Target_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1057 P_reco_zero_Target_mbias4->Fill(part1.Mag());
1058 P_reco_zero_Target_mbias4->Fill(part2.Mag());
1059 Pt_reco_zero_Target_mbias4->Fill(part1.Perp());
1060 Pt_reco_zero_Target_mbias4->Fill(part2.Perp());
1061 GammaInvMassReco_zero_Both_mbias4->Fill(InvmassReco);
1062 GammaOpeningAngleReco_zero_Both_mbias4->Fill(OpeningAngle);
1063 Pdg_zero_Both_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1064 Pdg_zero_Both_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1065 P_reco_zero_Both_mbias4->Fill(part1.Mag());
1066 P_reco_zero_Both_mbias4->Fill(part2.Mag());
1067 Pt_reco_zero_Both_mbias4->Fill(part1.Perp());
1068 Pt_reco_zero_Both_mbias4->Fill(part2.Perp());
1069 }
1070
1071 // only cases, when RICH == 1
1072 if (richcheck == 1) {
1076 GammaInvMassReco_one_Target_mbias4->Fill(InvmassReco);
1077 GammaOpeningAngleReco_one_Target_mbias4->Fill(OpeningAngle);
1078 Pdg_one_Target_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1079 Pdg_one_Target_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1080 P_reco_one_Target_mbias4->Fill(part1.Mag());
1081 P_reco_one_Target_mbias4->Fill(part2.Mag());
1082 Pt_reco_one_Target_mbias4->Fill(part1.Perp());
1083 Pt_reco_one_Target_mbias4->Fill(part2.Perp());
1084 GammaInvMassReco_one_Both_mbias4->Fill(InvmassReco);
1085 GammaOpeningAngleReco_one_Both_mbias4->Fill(OpeningAngle);
1086 Pdg_one_Both_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1087 Pdg_one_Both_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1088 P_reco_one_Both_mbias4->Fill(part1.Mag());
1089 P_reco_one_Both_mbias4->Fill(part2.Mag());
1090 Pt_reco_one_Both_mbias4->Fill(part1.Perp());
1091 Pt_reco_one_Both_mbias4->Fill(part2.Perp());
1092 }
1093
1094 // only cases, when RICH == 2
1095 if (richcheck == 2) {
1099 GammaInvMassReco_two_Target_mbias4->Fill(InvmassReco);
1100 GammaOpeningAngleReco_two_Target_mbias4->Fill(OpeningAngle);
1101 Pdg_two_Target_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1102 Pdg_two_Target_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1103 P_reco_two_Target_mbias4->Fill(part1.Mag());
1104 P_reco_two_Target_mbias4->Fill(part2.Mag());
1105 Pt_reco_two_Target_mbias4->Fill(part1.Perp());
1106 Pt_reco_two_Target_mbias4->Fill(part2.Perp());
1107 GammaInvMassReco_two_Both_mbias4->Fill(InvmassReco);
1108 GammaOpeningAngleReco_two_Both_mbias4->Fill(OpeningAngle);
1109 Pdg_two_Both_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1110 Pdg_two_Both_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1111 P_reco_two_Both_mbias4->Fill(part1.Mag());
1112 P_reco_two_Both_mbias4->Fill(part2.Mag());
1113 Pt_reco_two_Both_mbias4->Fill(part1.Perp());
1114 Pt_reco_two_Both_mbias4->Fill(part2.Perp());
1115 }
1116
1117 // cases, when RICH == 1 or RICH == 2 together
1118 if (richcheck == 1 || richcheck == 2) {
1122 GammaInvMassReco_onetwo_Target_mbias4->Fill(InvmassReco);
1124 Pdg_onetwo_Target_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1125 Pdg_onetwo_Target_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1126 P_reco_onetwo_Target_mbias4->Fill(part1.Mag());
1127 P_reco_onetwo_Target_mbias4->Fill(part2.Mag());
1128 Pt_reco_onetwo_Target_mbias4->Fill(part1.Perp());
1129 Pt_reco_onetwo_Target_mbias4->Fill(part2.Perp());
1130 GammaInvMassReco_onetwo_Both_mbias4->Fill(InvmassReco);
1131 GammaOpeningAngleReco_onetwo_Both_mbias4->Fill(OpeningAngle);
1132 Pdg_onetwo_Both_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1133 Pdg_onetwo_Both_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1134 P_reco_onetwo_Both_mbias4->Fill(part1.Mag());
1135 P_reco_onetwo_Both_mbias4->Fill(part2.Mag());
1136 Pt_reco_onetwo_Both_mbias4->Fill(part1.Perp());
1137 Pt_reco_onetwo_Both_mbias4->Fill(part2.Perp());
1138 }
1139 }
1140 }
1141
1142 //cout << "number of gammas Target with 0-2 electron identified in RICH = " << Gammas_all_Target_mbias4.size() << endl;
1143 //cout << "number of gammas Target with 0 electron identified in RICH = " << Gammas_zero_Target_mbias4.size() << endl;
1144 //cout << "number of gammas Target with 1 electron identified in RICH = " << Gammas_one_Target_mbias4.size() << endl;
1145 //cout << "number of gammas Target with 2 electron identified in RICH = " << Gammas_two_Target_mbias4.size() << endl;
1146 //cout << "number of gammas Target with 1-2 electron identified in RICH = " << Gammas_onetwo_Target_mbias4.size() << endl;
1147}
1148
1149
1151 int EventNumMan, double AngleCut, double InvMassCut, int RealPID, vector<CbmMCTrack*> MCtracks_minus_Outside_mbias4,
1152 vector<CbmMCTrack*> MCtracks_plus_Outside_mbias4, vector<CbmStsTrack*> StsTrack_minus_Outside_mbias4,
1153 vector<CbmStsTrack*> StsTrack_plus_Outside_mbias4, std::vector<int> Rings_minus_Outside_mbias4,
1154 std::vector<int> Rings_plus_Outside_mbias4, std::vector<int> stsIndex_minus_Outside_mbias4,
1155 std::vector<int> stsIndex_plus_Outside_mbias4, vector<CbmRichRing*> richRing_minus_Outside_mbias4,
1156 vector<CbmRichRing*> richRing_plus_Outside_mbias4)
1157{
1158 for (size_t i = 0; i < StsTrack_minus_Outside_mbias4.size(); i++) {
1159 for (size_t j = 0; j < StsTrack_plus_Outside_mbias4.size(); j++) {
1160
1161 CbmStsTrack* part1STS = StsTrack_minus_Outside_mbias4[i];
1162 CbmStsTrack* part2STS = StsTrack_plus_Outside_mbias4[j];
1163 CbmMCTrack* part1MC = MCtracks_minus_Outside_mbias4[i];
1164 CbmMCTrack* part2MC = MCtracks_plus_Outside_mbias4[j];
1165
1166 KFParticle electron;
1168 KFParticle positron;
1169 CbmKFParticleInterface::SetKFParticleFromStsTrack(part2STS, &positron, -11);
1170 const KFParticle* daughters[2] = {&electron, &positron};
1171 KFParticle intersection;
1172 intersection.Construct(daughters, 2);
1173
1174 if (intersection.GetZ() > 75 || intersection.GetZ() < -5) continue; // kick weird intersections
1175
1176
1177 TVector3 part1 =
1178 CbmKresFunctions::FitToVertex(part1STS, intersection.GetX(), intersection.GetY(), intersection.GetZ());
1179 TVector3 part2 =
1180 CbmKresFunctions::FitToVertex(part2STS, intersection.GetX(), intersection.GetY(), intersection.GetZ());
1181
1182
1183 //cout << "=================" << endl;
1184 //cout << "FitToVertex part1: px = " << part1.X() << "; \t py = " << part1.Y() << "; \t pz = " << part1.Z() << endl;
1185 //cout << "FitToVertex part2: px = " << part2.X() << "; \t py = " << part2.Y() << "; \t pz = " << part2.Z() << endl;
1186 //cout << "KFParticle part1: px = " << electron.KFParticleBase::Px() << "; \t py = " << electron.KFParticleBase::Py() << "; \t pz = " << electron.KFParticleBase::Pz() << endl;
1187 //cout << "KFParticle part2: px = " << positron.KFParticleBase::Px() << "; \t py = " << positron.KFParticleBase::Py() << "; \t pz = " << positron.KFParticleBase::Pz() << endl;
1188 //TVector3 part1(electron.KFParticleBase::Px(), electron.KFParticleBase::Py(), electron.KFParticleBase::Pz());
1189 //TVector3 part2(positron.KFParticleBase::Px(), positron.KFParticleBase::Py(), positron.KFParticleBase::Pz());
1190
1191
1192 int richcheck = 0;
1193 int richcheck_0 = 0;
1194 int richcheck_1 = 0;
1195 if (RealPID == 1) {
1196 // Real PID
1197 richcheck_0 = CheckIfElectron(richRing_minus_Outside_mbias4[i], part1.Mag());
1198 richcheck_1 = CheckIfElectron(richRing_plus_Outside_mbias4[j], part2.Mag());
1199 richcheck = richcheck_0 + richcheck_1;
1200 }
1201 else {
1202 // MC PID
1203 richcheck = Rings_minus_Outside_mbias4[i] + Rings_plus_Outside_mbias4[j];
1204 }
1205
1206
1207 Double_t InvmassReco = CbmKresFunctions::Invmass_2particles_RECO(part1, part2);
1208 Double_t OpeningAngle = CbmKresFunctions::CalculateOpeningAngle_Reco(part1, part2);
1209
1210 InvMass_vs_OA_candidates_Outside_mbias4->Fill(InvmassReco, OpeningAngle);
1211 InvMass_vs_OA_candidates_Both_mbias4->Fill(InvmassReco, OpeningAngle);
1212 GammasInvMass_candidates_Outside_mbias4->Fill(InvmassReco);
1213 GammasInvMass_candidates_Both_mbias4->Fill(InvmassReco);
1214 GammasOA_candidates_Outside_mbias4->Fill(OpeningAngle);
1215 GammasOA_candidates_Both_mbias4->Fill(OpeningAngle);
1216
1217 double PlaneAngle_last = CalculatePlaneAngle_last(part1STS, part2STS);
1218 double PlaneAngle_first = CalculatePlaneAngle_first(part1STS, part2STS);
1219
1220 // if (PlaneAngle_last > 20) continue;
1221
1222 if (part1MC->GetMotherId() == part2MC->GetMotherId() && part1MC->GetMotherId() != -1) {
1223 CbmMCTrack* mcTrackmama = (CbmMCTrack*) fMcTracks->At(part1MC->GetMotherId());
1224 if (mcTrackmama->GetMotherId() != -1 && mcTrackmama->GetPdgCode() == 22) {
1225 CbmMCTrack* mcTrackgrmama = (CbmMCTrack*) fMcTracks->At(mcTrackmama->GetMotherId());
1226 if (mcTrackgrmama->GetPdgCode() == 111) {
1227 GammasInvMass_fromPi0_Outside_mbias4->Fill(InvmassReco);
1228 GammasInvMass_fromPi0_Both_mbias4->Fill(InvmassReco);
1229 GammasOA_fromPi0_Outside_mbias4->Fill(OpeningAngle);
1230 GammasOA_fromPi0_Both_mbias4->Fill(OpeningAngle);
1231 InvMass_vs_OA_fromPi0_Outside_mbias4->Fill(InvmassReco, OpeningAngle);
1232 InvMass_vs_OA_fromPi0_Both_mbias4->Fill(InvmassReco, OpeningAngle);
1233 PlaneAngles_last_fromPi0_Outside_mbias4->Fill(PlaneAngle_last);
1234 PlaneAngles_last_fromPi0_Both_mbias4->Fill(PlaneAngle_last);
1235 PlaneAngles_first_fromPi0_Outside_mbias4->Fill(PlaneAngle_first);
1236 PlaneAngles_first_fromPi0_Both_mbias4->Fill(PlaneAngle_first);
1237 }
1238 }
1239 }
1240
1241 PlaneAngles_last_candidates_Outside_mbias4->Fill(PlaneAngle_last);
1242 PlaneAngles_last_candidates_Both_mbias4->Fill(PlaneAngle_last);
1243 PlaneAngles_first_candidates_Outside_mbias4->Fill(PlaneAngle_first);
1244 PlaneAngles_first_candidates_Both_mbias4->Fill(PlaneAngle_first);
1245
1246 // cuts
1247 if (TMath::Abs(OpeningAngle) > AngleCut) continue;
1248 if (TMath::Abs(InvmassReco) > InvMassCut) continue;
1249
1250 frefmomenta.clear();
1251 frefmomenta.push_back(part1);
1252 frefmomenta.push_back(part2);
1253
1254 frefId.clear();
1255 frefId.push_back(stsIndex_minus_Outside_mbias4[i]);
1256 frefId.push_back(stsIndex_plus_Outside_mbias4[j]);
1257 fMCtracks.clear();
1258 fMCtracks.push_back(part1MC);
1259 fMCtracks.push_back(part2MC);
1260
1261
1262 // for event mixing
1263 EMT_man_Event_Outside_mbias4.push_back(EventNumMan);
1265 EMT_man_NofRings_Outside_mbias4.push_back(richcheck);
1266
1267
1268 // everything (RICH == 0, RICH == 1, RICH == 2) together
1269 if (richcheck == 0 || richcheck == 1 || richcheck == 2) {
1273 GammaInvMassReco_all_Outside_mbias4->Fill(InvmassReco);
1274 GammaOpeningAngleReco_all_Outside_mbias4->Fill(OpeningAngle);
1275 Pdg_all_Outside_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1276 Pdg_all_Outside_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1277 P_reco_all_Outside_mbias4->Fill(part1.Mag());
1278 P_reco_all_Outside_mbias4->Fill(part2.Mag());
1279 Pt_reco_all_Outside_mbias4->Fill(part1.Perp());
1280 Pt_reco_all_Outside_mbias4->Fill(part2.Perp());
1281 GammaInvMassReco_all_Both_mbias4->Fill(InvmassReco);
1282 GammaOpeningAngleReco_all_Both_mbias4->Fill(OpeningAngle);
1283 Pdg_all_Both_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1284 Pdg_all_Both_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1285 P_reco_all_Both_mbias4->Fill(part1.Mag());
1286 P_reco_all_Both_mbias4->Fill(part2.Mag());
1287 Pt_reco_all_Both_mbias4->Fill(part1.Perp());
1288 Pt_reco_all_Both_mbias4->Fill(part2.Perp());
1289 }
1290
1291 // only cases, when RICH == 0
1292 if (richcheck == 0) {
1296 GammaInvMassReco_zero_Outside_mbias4->Fill(InvmassReco);
1298 Pdg_zero_Outside_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1299 Pdg_zero_Outside_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1300 P_reco_zero_Outside_mbias4->Fill(part1.Mag());
1301 P_reco_zero_Outside_mbias4->Fill(part2.Mag());
1302 Pt_reco_zero_Outside_mbias4->Fill(part1.Perp());
1303 Pt_reco_zero_Outside_mbias4->Fill(part2.Perp());
1304 GammaInvMassReco_zero_Both_mbias4->Fill(InvmassReco);
1305 GammaOpeningAngleReco_zero_Both_mbias4->Fill(OpeningAngle);
1306 Pdg_zero_Both_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1307 Pdg_zero_Both_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1308 P_reco_zero_Both_mbias4->Fill(part1.Mag());
1309 P_reco_zero_Both_mbias4->Fill(part2.Mag());
1310 Pt_reco_zero_Both_mbias4->Fill(part1.Perp());
1311 Pt_reco_zero_Both_mbias4->Fill(part2.Perp());
1312 }
1313
1314 // only cases, when RICH == 1
1315 if (richcheck == 1) {
1319 GammaInvMassReco_one_Outside_mbias4->Fill(InvmassReco);
1320 GammaOpeningAngleReco_one_Outside_mbias4->Fill(OpeningAngle);
1321 Pdg_one_Outside_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1322 Pdg_one_Outside_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1323 P_reco_one_Outside_mbias4->Fill(part1.Mag());
1324 P_reco_one_Outside_mbias4->Fill(part2.Mag());
1325 Pt_reco_one_Outside_mbias4->Fill(part1.Perp());
1326 Pt_reco_one_Outside_mbias4->Fill(part2.Perp());
1327 GammaInvMassReco_one_Both_mbias4->Fill(InvmassReco);
1328 GammaOpeningAngleReco_one_Both_mbias4->Fill(OpeningAngle);
1329 Pdg_one_Both_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1330 Pdg_one_Both_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1331 P_reco_one_Both_mbias4->Fill(part1.Mag());
1332 P_reco_one_Both_mbias4->Fill(part2.Mag());
1333 Pt_reco_one_Both_mbias4->Fill(part1.Perp());
1334 Pt_reco_one_Both_mbias4->Fill(part2.Perp());
1335 }
1336
1337 // only cases, when RICH == 2
1338 if (richcheck == 2) {
1342 GammaInvMassReco_two_Outside_mbias4->Fill(InvmassReco);
1343 GammaOpeningAngleReco_two_Outside_mbias4->Fill(OpeningAngle);
1344 Pdg_two_Outside_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1345 Pdg_two_Outside_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1346 P_reco_two_Outside_mbias4->Fill(part1.Mag());
1347 P_reco_two_Outside_mbias4->Fill(part2.Mag());
1348 Pt_reco_two_Outside_mbias4->Fill(part1.Perp());
1349 Pt_reco_two_Outside_mbias4->Fill(part2.Perp());
1350 GammaInvMassReco_two_Both_mbias4->Fill(InvmassReco);
1351 GammaOpeningAngleReco_two_Both_mbias4->Fill(OpeningAngle);
1352 Pdg_two_Both_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1353 Pdg_two_Both_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1354 P_reco_two_Both_mbias4->Fill(part1.Mag());
1355 P_reco_two_Both_mbias4->Fill(part2.Mag());
1356 Pt_reco_two_Both_mbias4->Fill(part1.Perp());
1357 Pt_reco_two_Both_mbias4->Fill(part2.Perp());
1358 }
1359
1360 // cases, when RICH == 1 or RICH == 2 together
1361 if (richcheck == 1 || richcheck == 2) {
1365 GammaInvMassReco_onetwo_Outside_mbias4->Fill(InvmassReco);
1367 Pdg_onetwo_Outside_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1368 Pdg_onetwo_Outside_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1369 P_reco_onetwo_Outside_mbias4->Fill(part1.Mag());
1370 P_reco_onetwo_Outside_mbias4->Fill(part2.Mag());
1371 Pt_reco_onetwo_Outside_mbias4->Fill(part1.Perp());
1372 Pt_reco_onetwo_Outside_mbias4->Fill(part2.Perp());
1373 GammaInvMassReco_onetwo_Both_mbias4->Fill(InvmassReco);
1374 GammaOpeningAngleReco_onetwo_Both_mbias4->Fill(OpeningAngle);
1375 Pdg_onetwo_Both_mbias4->Fill(TMath::Abs(part1MC->GetPdgCode()));
1376 Pdg_onetwo_Both_mbias4->Fill(TMath::Abs(part2MC->GetPdgCode()));
1377 P_reco_onetwo_Both_mbias4->Fill(part1.Mag());
1378 P_reco_onetwo_Both_mbias4->Fill(part2.Mag());
1379 Pt_reco_onetwo_Both_mbias4->Fill(part1.Perp());
1380 Pt_reco_onetwo_Both_mbias4->Fill(part2.Perp());
1381 }
1382 }
1383 }
1384
1385 //cout << "number of gammas Outside with 0-2 electron identified in RICH = " << Gammas_all_Outside_mbias4.size() << endl;
1386 //cout << "number of gammas Outside with 0 electron identified in RICH = " << Gammas_zero_Outside_mbias4.size() << endl;
1387 //cout << "number of gammas Outside with 1 electron identified in RICH = " << Gammas_one_Outside_mbias4.size() << endl;
1388 //cout << "number of gammas Outside with 2 electron identified in RICH = " << Gammas_two_Outside_mbias4.size() << endl;
1389 //cout << "number of gammas Outside with 1-2 electron identified in RICH = " << Gammas_onetwo_Outside_mbias4.size() << endl;
1390}
1391
1392
1394{
1415
1446
1467
1468 //cout << "number of gammas Both with 0-2 electron identified in RICH = " << Gammas_all_Both_mbias4.size() << endl;
1469 //cout << "number of gammas Both with 0 electron identified in RICH = " << Gammas_zero_Both_mbias4.size() << endl;
1470 //cout << "number of gammas Both with 1 electron identified in RICH = " << Gammas_one_Both_mbias4.size() << endl;
1471 //cout << "number of gammas Both with 2 electron identified in RICH = " << Gammas_two_Both_mbias4.size() << endl;
1472 //cout << "number of gammas Both with 1-2 electron identified in RICH = " << Gammas_onetwo_Both_mbias4.size() << endl;
1473}
1474
1475
1476void CbmKresConversionManualmbias4::FindPi0(TString mod, TString position, vector<vector<TVector3>> Gammas,
1477 vector<vector<int>> StsIndex, vector<vector<CbmMCTrack*>> GammasMC,
1478 TH1D* Pi0InvMassReco, TH2D* Pi0_pt_vs_rap, TH2D* Pi0_pt_vs_rap_est,
1479 TH2D* MultiplicityGamma, TH2D* MultiplicityChargedParticles,
1480 vector<TH1*> BGCases)
1481{
1482 // combine all gamma in pi0 --> calculate inv mass for gammas and pi0 // not the case, when one particle used twice for different gammas
1483 if (Gammas.size() < 2) return; // min 2 gammas to form pi0 are required
1484 for (size_t gamma1 = 0; gamma1 < Gammas.size() - 1; gamma1++) {
1485 for (size_t gamma2 = gamma1 + 1; gamma2 < Gammas.size(); gamma2++) {
1486 // 4 reconstructed particles from gammas
1487 TVector3 e11 = Gammas[gamma1][0];
1488 TVector3 e12 = Gammas[gamma1][1];
1489 TVector3 e21 = Gammas[gamma2][0];
1490 TVector3 e22 = Gammas[gamma2][1];
1491
1492 // MC true data for this particles
1493 CbmMCTrack* mcTrack1 = GammasMC[gamma1][0];
1494 CbmMCTrack* mcTrack2 = GammasMC[gamma1][1];
1495 CbmMCTrack* mcTrack3 = GammasMC[gamma2][0];
1496 CbmMCTrack* mcTrack4 = GammasMC[gamma2][1];
1497
1498 if (StsIndex[gamma1][0] == StsIndex[gamma2][0] || StsIndex[gamma1][0] == StsIndex[gamma2][1]
1499 || StsIndex[gamma1][1] == StsIndex[gamma2][0] || StsIndex[gamma1][1] == StsIndex[gamma2][1])
1500 continue; // particles not used twice --> different
1501
1503
1504 //Double_t OA1 = CbmKresFunctions::CalculateOpeningAngle_Reco(e11, e12);
1505 //Double_t OA2 = CbmKresFunctions::CalculateOpeningAngle_Reco(e21, e22);
1506 //if (params.fRapidity > 2.8 && params.fRapidity <= 3.6 && params.fPt > 0.0 && params.fPt <= 0.4 && OA1 > 1 && OA2 > 1) continue;
1507
1508 Pi0InvMassReco->Fill(params.fMinv);
1509 MultiplicityGamma->Fill(Gammas.size(), params.fMinv);
1510 MultiplicityChargedParticles->Fill(fGlobalTracks->GetEntriesFast(), params.fMinv);
1511
1512 // separate by rap and Pt only for onetwo
1513 if (mod == "OneTwo" && position == "Both") {
1514 if (params.fRapidity > 1.2 && params.fRapidity <= 1.6) {
1515 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_1_mbias4->Fill(params.fMinv);
1516 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_2_mbias4->Fill(params.fMinv);
1517 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_3_mbias4->Fill(params.fMinv);
1518 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_4_mbias4->Fill(params.fMinv);
1519 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_5_mbias4->Fill(params.fMinv);
1520 }
1521 if (params.fRapidity > 1.6 && params.fRapidity <= 2.0) {
1522 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_6_mbias4->Fill(params.fMinv);
1523 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_7_mbias4->Fill(params.fMinv);
1524 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_8_mbias4->Fill(params.fMinv);
1525 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_9_mbias4->Fill(params.fMinv);
1526 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_10_mbias4->Fill(params.fMinv);
1527 }
1528 if (params.fRapidity > 2.0 && params.fRapidity <= 2.4) {
1529 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_11_mbias4->Fill(params.fMinv);
1530 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_12_mbias4->Fill(params.fMinv);
1531 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_13_mbias4->Fill(params.fMinv);
1532 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_14_mbias4->Fill(params.fMinv);
1533 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_15_mbias4->Fill(params.fMinv);
1534 }
1535 if (params.fRapidity > 2.4 && params.fRapidity <= 2.8) {
1536 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_16_mbias4->Fill(params.fMinv);
1537 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_17_mbias4->Fill(params.fMinv);
1538 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_18_mbias4->Fill(params.fMinv);
1539 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_19_mbias4->Fill(params.fMinv);
1540 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_20_mbias4->Fill(params.fMinv);
1541 }
1542 if (params.fRapidity > 2.8 && params.fRapidity <= 3.2) {
1543 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_21_mbias4->Fill(params.fMinv);
1544 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_22_mbias4->Fill(params.fMinv);
1545 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_23_mbias4->Fill(params.fMinv);
1546 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_24_mbias4->Fill(params.fMinv);
1547 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_25_mbias4->Fill(params.fMinv);
1548 }
1549 if (params.fRapidity > 3.2 && params.fRapidity <= 3.6) {
1550 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_26_mbias4->Fill(params.fMinv);
1551 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_27_mbias4->Fill(params.fMinv);
1552 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_28_mbias4->Fill(params.fMinv);
1553 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_29_mbias4->Fill(params.fMinv);
1554 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_30_mbias4->Fill(params.fMinv);
1555 }
1556 }
1557
1558
1559 if (position == "Both") fAnaBG->Exec(mcTrack1, mcTrack2, mcTrack3, mcTrack4, params.fMinv, BGCases);
1560
1561
1562 // fill histos with rapidity and Pt for correctly reconstructed pi0(dalitz) and pi0(g+g)
1563 if (nullptr == mcTrack1 || nullptr == mcTrack2 || nullptr == mcTrack3 || nullptr == mcTrack4) continue;
1564 if (TMath::Abs(mcTrack1->GetPdgCode()) != 11 || TMath::Abs(mcTrack2->GetPdgCode()) != 11
1565 || TMath::Abs(mcTrack3->GetPdgCode()) != 11 || TMath::Abs(mcTrack4->GetPdgCode()) != 11)
1566 continue;
1567 if (mcTrack1->GetPdgCode() + mcTrack2->GetPdgCode() != 0) continue;
1568 if (mcTrack3->GetPdgCode() + mcTrack4->GetPdgCode() != 0) continue;
1569 int motherId1 = mcTrack1->GetMotherId();
1570 int motherId2 = mcTrack2->GetMotherId();
1571 int motherId3 = mcTrack3->GetMotherId();
1572 int motherId4 = mcTrack4->GetMotherId();
1573 if (motherId1 == -1 || motherId2 == -1 || motherId3 == -1 || motherId4 == -1) continue;
1574 if (motherId1 != motherId2 || motherId3 != motherId4) continue;
1575 CbmMCTrack* mother1 = (CbmMCTrack*) fMcTracks->At(motherId1);
1576 CbmMCTrack* mother2 = (CbmMCTrack*) fMcTracks->At(motherId2);
1577 CbmMCTrack* mother3 = (CbmMCTrack*) fMcTracks->At(motherId3);
1578 CbmMCTrack* mother4 = (CbmMCTrack*) fMcTracks->At(motherId4);
1579 if (nullptr == mother1 || nullptr == mother2 || nullptr == mother3 || nullptr == mother4) continue;
1580 int mcMotherPdg1 = mother1->GetPdgCode();
1581 int mcMotherPdg2 = mother2->GetPdgCode();
1582 int mcMotherPdg3 = mother3->GetPdgCode();
1583 int mcMotherPdg4 = mother4->GetPdgCode();
1584 int grandmotherId1 = mother1->GetMotherId();
1585 int grandmotherId2 = mother2->GetMotherId();
1586 int grandmotherId3 = mother3->GetMotherId();
1587 int grandmotherId4 = mother4->GetMotherId();
1588
1589 if (mcMotherPdg1 == 22 && mcMotherPdg2 == 22 && mcMotherPdg3 == 111 && mcMotherPdg4 == 111) {
1590 if (grandmotherId1 != motherId3) continue;
1591 Pi0_pt_vs_rap->Fill(params.fRapidity, params.fPt);
1592 Pi0_pt_vs_rap_est->Fill(params.fRapidity, params.fPt);
1593 }
1594
1595 if (mcMotherPdg1 == 111 && mcMotherPdg2 == 111 && mcMotherPdg3 == 22 && mcMotherPdg4 == 22) {
1596 if (grandmotherId3 != motherId1) continue;
1597 Pi0_pt_vs_rap->Fill(params.fRapidity, params.fPt);
1598 Pi0_pt_vs_rap_est->Fill(params.fRapidity, params.fPt);
1599 }
1600
1601 if (mcMotherPdg1 == 22 && mcMotherPdg2 == 22 && mcMotherPdg3 == 22 && mcMotherPdg4 == 22) {
1602 if (grandmotherId1 != grandmotherId2 || grandmotherId3 != grandmotherId4 || grandmotherId1 != grandmotherId3)
1603 continue;
1604 if (grandmotherId1 == -1) continue;
1605 CbmMCTrack* grandmother1 = (CbmMCTrack*) fMcTracks->At(grandmotherId1);
1606 if (nullptr == grandmother1) continue;
1607 int mcGrandMotherPdg1 = grandmother1->GetPdgCode();
1608 if (mcGrandMotherPdg1 != 111) continue;
1609 Pi0_pt_vs_rap->Fill(params.fRapidity, params.fPt);
1610 Pi0_pt_vs_rap_est->Fill(params.fRapidity, params.fPt);
1611 }
1612 }
1613 }
1614}
1615
1616
1618// TARGET
1619{
1620 Int_t nof_Target = EMT_man_Event_Target_mbias4.size();
1621 cout << "Mixing in Manual(Target_mbias4) - nof entries " << nof_Target << endl;
1622 for (Int_t a = 0; a < nof_Target - 1; a++) {
1623 for (Int_t b = a + 1; b < nof_Target; b++) {
1625 continue; // to make sure that the photons are from two different events
1626 TVector3 e11 = EMT_man_pair_momenta_Target_mbias4[a][0];
1627 TVector3 e12 = EMT_man_pair_momenta_Target_mbias4[a][1];
1628 TVector3 e21 = EMT_man_pair_momenta_Target_mbias4[b][0];
1629 TVector3 e22 = EMT_man_pair_momenta_Target_mbias4[b][1];
1631
1642 }
1643 }
1644}
1645
1646
1648// OUTSIDE
1649{
1650 Int_t nof_Outside = EMT_man_Event_Outside_mbias4.size();
1651 cout << "Mixing in Manual(Outside_mbias4) - nof entries " << nof_Outside << endl;
1652 for (Int_t a = 0; a < nof_Outside - 1; a++) {
1653 for (Int_t b = a + 1; b < nof_Outside; b++) {
1655 continue; // to make sure that the photons are from two different events
1656 TVector3 e11 = EMT_man_pair_momenta_Outside_mbias4[a][0];
1657 TVector3 e12 = EMT_man_pair_momenta_Outside_mbias4[a][1];
1658 TVector3 e21 = EMT_man_pair_momenta_Outside_mbias4[b][0];
1659 TVector3 e22 = EMT_man_pair_momenta_Outside_mbias4[b][1];
1661
1672 }
1673 }
1674}
1675
1676
1678// BOTH
1679{
1680 Int_t nof_Both = EMT_man_Event_Both_mbias4.size();
1681 cout << "Mixing in Manual(Both_mbias4) - nof entries " << nof_Both << endl;
1682 for (Int_t a = 0; a < nof_Both - 1; a++) {
1683 for (Int_t b = a + 1; b < nof_Both; b++) {
1685 continue; // to make sure that the photons are from two different events
1686 TVector3 e11 = EMT_man_pair_momenta_Both_mbias4[a][0];
1687 TVector3 e12 = EMT_man_pair_momenta_Both_mbias4[a][1];
1688 TVector3 e21 = EMT_man_pair_momenta_Both_mbias4[b][0];
1689 TVector3 e22 = EMT_man_pair_momenta_Both_mbias4[b][1];
1691
1692 EMT_InvMass_all_Both_mbias4->Fill(params.fMinv);
1694 EMT_InvMass_zero_Both_mbias4->Fill(params.fMinv);
1696 EMT_InvMass_one_Both_mbias4->Fill(params.fMinv);
1698 EMT_InvMass_two_Both_mbias4->Fill(params.fMinv);
1702
1703 // separate by rap and Pt only for onetwo
1706 if (params.fRapidity > 1.2 && params.fRapidity <= 1.6) {
1707 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_mixing_1_mbias4->Fill(params.fMinv);
1708 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_mixing_2_mbias4->Fill(params.fMinv);
1709 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_mixing_3_mbias4->Fill(params.fMinv);
1710 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_mixing_4_mbias4->Fill(params.fMinv);
1711 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_mixing_5_mbias4->Fill(params.fMinv);
1712 }
1713 if (params.fRapidity > 1.6 && params.fRapidity <= 2.0) {
1714 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_mixing_6_mbias4->Fill(params.fMinv);
1715 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_mixing_7_mbias4->Fill(params.fMinv);
1716 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_mixing_8_mbias4->Fill(params.fMinv);
1717 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_mixing_9_mbias4->Fill(params.fMinv);
1718 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_mixing_10_mbias4->Fill(params.fMinv);
1719 }
1720 if (params.fRapidity > 2.0 && params.fRapidity <= 2.4) {
1721 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_mixing_11_mbias4->Fill(params.fMinv);
1722 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_mixing_12_mbias4->Fill(params.fMinv);
1723 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_mixing_13_mbias4->Fill(params.fMinv);
1724 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_mixing_14_mbias4->Fill(params.fMinv);
1725 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_mixing_15_mbias4->Fill(params.fMinv);
1726 }
1727 if (params.fRapidity > 2.4 && params.fRapidity <= 2.8) {
1728 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_mixing_16_mbias4->Fill(params.fMinv);
1729 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_mixing_17_mbias4->Fill(params.fMinv);
1730 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_mixing_18_mbias4->Fill(params.fMinv);
1731 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_mixing_19_mbias4->Fill(params.fMinv);
1732 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_mixing_20_mbias4->Fill(params.fMinv);
1733 }
1734 if (params.fRapidity > 2.8 && params.fRapidity <= 3.2) {
1735 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_mixing_21_mbias4->Fill(params.fMinv);
1736 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_mixing_22_mbias4->Fill(params.fMinv);
1737 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_mixing_23_mbias4->Fill(params.fMinv);
1738 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_mixing_24_mbias4->Fill(params.fMinv);
1739 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_mixing_25_mbias4->Fill(params.fMinv);
1740 }
1741 if (params.fRapidity > 3.2 && params.fRapidity <= 3.6) {
1742 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_mixing_26_mbias4->Fill(params.fMinv);
1743 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_mixing_27_mbias4->Fill(params.fMinv);
1744 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_mixing_28_mbias4->Fill(params.fMinv);
1745 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_mixing_29_mbias4->Fill(params.fMinv);
1746 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_mixing_30_mbias4->Fill(params.fMinv);
1747 }
1748 }
1749 }
1750 }
1751}
1752
1753
1754int CbmKresConversionManualmbias4::FindInRich(int richInd, int stsMcTrackId)
1755{
1756 int RingsInRich = 0;
1757 if (richInd > 0) {
1758 CbmTrackMatchNew* richMatch = (CbmTrackMatchNew*) fRichRingMatches->At(richInd);
1759 if (richMatch != nullptr && richMatch->GetNofLinks() > 0) {
1760 int richMcTrackId = richMatch->GetMatchedLink().GetIndex();
1761 if (richMcTrackId > 0) {
1762 if (stsMcTrackId == richMcTrackId) { // check that global track was matched correctly for STS and RICH together
1763 CbmMCTrack* mcTrack2 = (CbmMCTrack*) fMcTracks->At(richMcTrackId);
1764 if (mcTrack2 != nullptr) {
1765 int pdgRICH = mcTrack2->GetPdgCode();
1766 if (TMath::Abs(pdgRICH) == 11) RingsInRich++;
1767 }
1768 }
1769 }
1770 }
1771 }
1772 return RingsInRich;
1773}
1774
1775
1777{
1778 int identified = 0;
1779
1780 if (nullptr != ring) {
1781 CbmRichRingLight ringHit;
1782 int nofHits = ring->GetNofHits();
1783 for (int i = 0; i < nofHits; i++) {
1784 int hitInd = ring->GetHit(i);
1785 CbmRichHit* hit = (CbmRichHit*) fRichHits->At(hitInd);
1786 if (nullptr == hit) continue;
1787 CbmRichHitLight hl(hit->GetX(), hit->GetY());
1788 ringHit.AddHit(hl);
1789 }
1790 fTauFit->DoFit(&ringHit);
1791 if (ringHit.GetAaxis() > 4 && ringHit.GetAaxis() < 6 && ringHit.GetBaxis() > 4 && ringHit.GetBaxis() < 6
1792 && momentum > 0.2 && momentum < 4.)
1793 identified++;
1794 //if (ring->GetDistance() < 2 && ringHit.GetAaxis() > 4 && ringHit.GetAaxis() < 6 && ringHit.GetBaxis() > 4 && ringHit.GetBaxis() < 6 && momentum > 0.2 && momentum < 4.) identified ++;
1795 }
1796
1797 return identified;
1798}
1799
1800
1802{
1803 double FinalAngle = 400;
1804 Int_t hits1sts = Sts_1->GetNofStsHits();
1805 Int_t hits2sts = Sts_2->GetNofStsHits();
1806 Int_t hits1mvd = Sts_1->GetNofMvdHits();
1807 Int_t hits2mvd = Sts_2->GetNofMvdHits();
1808
1809 double Xpart1 = 0;
1810 double Ypart1 = 0;
1811 double Zpart1 = 0;
1812 double Xpart2 = 0;
1813 double Ypart2 = 0;
1814 double Zpart2 = 0;
1815 if (hits1sts > 0) {
1816 Int_t stsHitIndex1 = Sts_1->GetStsHitIndex(hits1sts - 1);
1817 CbmStsHit* stsHit1 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex1);
1818 Xpart1 = stsHit1->GetX();
1819 Ypart1 = stsHit1->GetY();
1820 Zpart1 = stsHit1->GetZ();
1821 }
1822 else {
1823 Int_t mvdHitIndex1 = Sts_1->GetMvdHitIndex(hits1mvd - 1);
1824 CbmMvdHit* mvdHit1 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex1);
1825 Xpart1 = mvdHit1->GetX();
1826 Ypart1 = mvdHit1->GetY();
1827 Zpart1 = mvdHit1->GetZ();
1828 }
1829
1830 if (hits2sts > 0) {
1831 Int_t stsHitIndex2 = Sts_2->GetStsHitIndex(hits2sts - 1);
1832 CbmStsHit* stsHit2 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex2);
1833 Xpart2 = stsHit2->GetX();
1834 Ypart2 = stsHit2->GetY();
1835 Zpart2 = stsHit2->GetZ();
1836 }
1837 else {
1838 Int_t mvdHitIndex2 = Sts_2->GetMvdHitIndex(hits2mvd - 1);
1839 CbmMvdHit* mvdHit2 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex2);
1840 Xpart2 = mvdHit2->GetX();
1841 Ypart2 = mvdHit2->GetY();
1842 Zpart2 = mvdHit2->GetZ();
1843 }
1844
1845 // check difference in 2 cm, because of two slices of every STS and Mvd stations
1846 if (TMath::Abs(Zpart1 - Zpart2) > 2
1847 && Zpart1 > Zpart2) { // if last hits are in different stations --> try to find the latest common station
1848 for (int i = hits1sts - 2; i > -1; i--) { // start from second last station
1849 Int_t stsHitIndex = Sts_1->GetStsHitIndex(i);
1850 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
1851 Xpart1 = stsHit->GetX();
1852 Ypart1 = stsHit->GetY();
1853 Zpart1 = stsHit->GetZ();
1854 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1855 }
1856 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 > Zpart2) {
1857 for (int i = hits1mvd - 2; i > -1; i--) {
1858 Int_t mvdHitIndex = Sts_1->GetMvdHitIndex(i);
1859 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
1860 Xpart1 = mvdHit->GetX();
1861 Ypart1 = mvdHit->GetY();
1862 Zpart1 = mvdHit->GetZ();
1863 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1864 }
1865 }
1866 }
1867
1868 if (TMath::Abs(Zpart1 - Zpart2) > 2
1869 && Zpart1 < Zpart2) { // if last hits are in different stations --> try to find the latest common station
1870 for (int i = hits2sts - 2; i > -1; i--) { // start from second last station
1871 Int_t stsHitIndex = Sts_2->GetStsHitIndex(i);
1872 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
1873 Xpart2 = stsHit->GetX();
1874 Ypart2 = stsHit->GetY();
1875 Zpart2 = stsHit->GetZ();
1876 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1877 }
1878 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 < Zpart2) {
1879 for (int i = hits2mvd - 2; i > -1; i--) {
1880 Int_t mvdHitIndex = Sts_2->GetMvdHitIndex(i);
1881 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
1882 Xpart2 = mvdHit->GetX();
1883 Ypart2 = mvdHit->GetY();
1884 Zpart2 = mvdHit->GetZ();
1885 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1886 }
1887 }
1888 }
1889
1890 // calculate angle if we found common station
1891 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
1892 FinalAngle = TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 / TMath::Pi());
1893 }
1894
1895 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
1896}
1897
1899{
1900 double FinalAngle = 400;
1901 Int_t hits1sts = Sts_1->GetNofStsHits();
1902 Int_t hits2sts = Sts_2->GetNofStsHits();
1903 Int_t hits1mvd = Sts_1->GetNofMvdHits();
1904 Int_t hits2mvd = Sts_2->GetNofMvdHits();
1905
1906 double Xpart1 = 0;
1907 double Ypart1 = 0;
1908 double Zpart1 = 0;
1909 double Xpart2 = 0;
1910 double Ypart2 = 0;
1911 double Zpart2 = 0;
1912 if (hits1mvd > 0) {
1913 Int_t mvdHitIndex1 = Sts_1->GetMvdHitIndex(0);
1914 CbmMvdHit* mvdHit1 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex1);
1915 Xpart1 = mvdHit1->GetX();
1916 Ypart1 = mvdHit1->GetY();
1917 Zpart1 = mvdHit1->GetZ();
1918 }
1919 else {
1920 Int_t stsHitIndex1 = Sts_1->GetStsHitIndex(0);
1921 CbmStsHit* stsHit1 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex1);
1922 Xpart1 = stsHit1->GetX();
1923 Ypart1 = stsHit1->GetY();
1924 Zpart1 = stsHit1->GetZ();
1925 }
1926
1927 if (hits2mvd > 0) {
1928 Int_t mvdHitIndex2 = Sts_2->GetMvdHitIndex(0);
1929 CbmMvdHit* mvdHit2 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex2);
1930 Xpart2 = mvdHit2->GetX();
1931 Ypart2 = mvdHit2->GetY();
1932 Zpart2 = mvdHit2->GetZ();
1933 }
1934 else {
1935 Int_t stsHitIndex2 = Sts_2->GetStsHitIndex(0);
1936 CbmStsHit* stsHit2 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex2);
1937 Xpart2 = stsHit2->GetX();
1938 Ypart2 = stsHit2->GetY();
1939 Zpart2 = stsHit2->GetZ();
1940 }
1941
1942 // check difference in 2 cm, because of two slices of every STS and Mvd stations
1943 if (TMath::Abs(Zpart1 - Zpart2) > 2
1944 && Zpart1 < Zpart2) { // if first hits are in different stations --> try to find the earliest common station
1945 for (int i = 1; i < hits1mvd; i++) { // start from second hit
1946 Int_t mvdHitIndex = Sts_1->GetMvdHitIndex(i);
1947 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
1948 Xpart1 = mvdHit->GetX();
1949 Ypart1 = mvdHit->GetY();
1950 Zpart1 = mvdHit->GetZ();
1951 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1952 }
1953 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 < Zpart2) {
1954 for (int i = 0; i < hits1sts; i++) {
1955 Int_t stsHitIndex = Sts_1->GetStsHitIndex(i);
1956 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
1957 Xpart1 = stsHit->GetX();
1958 Ypart1 = stsHit->GetY();
1959 Zpart1 = stsHit->GetZ();
1960 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1961 }
1962 }
1963 }
1964
1965 if (TMath::Abs(Zpart1 - Zpart2) > 2
1966 && Zpart1 > Zpart2) { // if first hits are in different stations --> try to find the earliest common station
1967 for (int i = 1; i < hits2mvd; i++) { // start from second hit
1968 Int_t mvdHitIndex = Sts_2->GetMvdHitIndex(i);
1969 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
1970 Xpart2 = mvdHit->GetX();
1971 Ypart2 = mvdHit->GetY();
1972 Zpart2 = mvdHit->GetZ();
1973 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1974 }
1975 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 > Zpart2) {
1976 for (int i = 0; i < hits2sts; i++) {
1977 Int_t stsHitIndex = Sts_2->GetStsHitIndex(i);
1978 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
1979 Xpart2 = stsHit->GetX();
1980 Ypart2 = stsHit->GetY();
1981 Zpart2 = stsHit->GetZ();
1982 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1983 }
1984 }
1985 }
1986
1987 // calculate angle if we found common station
1988 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
1989 FinalAngle = TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 / TMath::Pi());
1990 }
1991
1992 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
1993}
1994
1995
1997{
1998 gDirectory->mkdir("Target_mbias4");
1999 gDirectory->cd("Target_mbias4");
2000 gDirectory->mkdir("CheckCuts_Target_mbias4");
2001 gDirectory->cd("CheckCuts_Target_mbias4");
2002 for (UInt_t i = 0; i < fHistoList_man_cuts_Target_mbias4.size(); i++) {
2004 }
2005 gDirectory->cd("..");
2006 gDirectory->mkdir("multiplicity");
2007 gDirectory->cd("multiplicity");
2008 for (UInt_t i = 0; i < fHistoList_multiplicity_man_Target_mbias4.size(); i++) {
2010 }
2011 gDirectory->cd("..");
2012 gDirectory->mkdir("all");
2013 gDirectory->cd("all");
2014 for (UInt_t i = 0; i < fHistoList_man_all_Target_mbias4.size(); i++) {
2016 }
2017 gDirectory->cd("..");
2018 gDirectory->mkdir("zero");
2019 gDirectory->cd("zero");
2020 for (UInt_t i = 0; i < fHistoList_man_zero_Target_mbias4.size(); i++) {
2022 }
2023 gDirectory->cd("..");
2024 gDirectory->mkdir("one");
2025 gDirectory->cd("one");
2026 for (UInt_t i = 0; i < fHistoList_man_one_Target_mbias4.size(); i++) {
2028 }
2029 gDirectory->cd("..");
2030 gDirectory->mkdir("two");
2031 gDirectory->cd("two");
2032 for (UInt_t i = 0; i < fHistoList_man_two_Target_mbias4.size(); i++) {
2034 }
2035 gDirectory->cd("..");
2036 gDirectory->mkdir("onetwo");
2037 gDirectory->cd("onetwo");
2038 for (UInt_t i = 0; i < fHistoList_man_onetwo_Target_mbias4.size(); i++) {
2040 }
2041 gDirectory->cd("..");
2042
2043 gDirectory->cd("..");
2044
2045
2046 gDirectory->mkdir("Outside_mbias4");
2047 gDirectory->cd("Outside_mbias4");
2048 gDirectory->mkdir("CheckCuts_Outside_mbias4");
2049 gDirectory->cd("CheckCuts_Outside_mbias4");
2050 for (UInt_t i = 0; i < fHistoList_man_cuts_Outside_mbias4.size(); i++) {
2052 }
2053 gDirectory->cd("..");
2054 gDirectory->mkdir("multiplicity");
2055 gDirectory->cd("multiplicity");
2056 for (UInt_t i = 0; i < fHistoList_multiplicity_man_Outside_mbias4.size(); i++) {
2058 }
2059 gDirectory->cd("..");
2060 gDirectory->mkdir("all");
2061 gDirectory->cd("all");
2062 for (UInt_t i = 0; i < fHistoList_man_all_Outside_mbias4.size(); i++) {
2064 }
2065 gDirectory->cd("..");
2066 gDirectory->mkdir("zero");
2067 gDirectory->cd("zero");
2068 for (UInt_t i = 0; i < fHistoList_man_zero_Outside_mbias4.size(); i++) {
2070 }
2071 gDirectory->cd("..");
2072 gDirectory->mkdir("one");
2073 gDirectory->cd("one");
2074 for (UInt_t i = 0; i < fHistoList_man_one_Outside_mbias4.size(); i++) {
2076 }
2077 gDirectory->cd("..");
2078 gDirectory->mkdir("two");
2079 gDirectory->cd("two");
2080 for (UInt_t i = 0; i < fHistoList_man_two_Outside_mbias4.size(); i++) {
2082 }
2083 gDirectory->cd("..");
2084 gDirectory->mkdir("onetwo");
2085 gDirectory->cd("onetwo");
2086 for (UInt_t i = 0; i < fHistoList_man_onetwo_Outside_mbias4.size(); i++) {
2088 }
2089 gDirectory->cd("..");
2090 gDirectory->cd("..");
2091
2092
2093 gDirectory->mkdir("Both_mbias4");
2094 gDirectory->cd("Both_mbias4");
2095 gDirectory->mkdir("CheckCuts_Both_mbias4");
2096 gDirectory->cd("CheckCuts_Both_mbias4");
2097 for (UInt_t i = 0; i < fHistoList_man_cuts_Both_mbias4.size(); i++) {
2099 }
2100 gDirectory->cd("..");
2101 gDirectory->mkdir("multiplicity");
2102 gDirectory->cd("multiplicity");
2103 for (UInt_t i = 0; i < fHistoList_multiplicity_man_Both_mbias4.size(); i++) {
2105 }
2106 gDirectory->cd("..");
2107 gDirectory->mkdir("rap_vs_Pt_onetwo_Both_mbias4");
2108 gDirectory->cd("rap_vs_Pt_onetwo_Both_mbias4");
2109 for (UInt_t i = 0; i < fHistoList_rap_vs_pt_InM_mbias4.size(); i++) {
2111 }
2112 gDirectory->cd("..");
2113 gDirectory->mkdir("BG_Both_mbias4");
2114 gDirectory->cd("BG_Both_mbias4");
2115 gDirectory->mkdir("BG_Both_mbias4_all");
2116 gDirectory->cd("BG_Both_mbias4_all");
2117 for (UInt_t i = 0; i < fHistoList_bg_InM_all_Both_mbias4.size(); i++) {
2119 }
2120 gDirectory->cd("..");
2121 gDirectory->mkdir("BG_Both_mbias4_zero");
2122 gDirectory->cd("BG_Both_mbias4_zero");
2123 for (UInt_t i = 0; i < fHistoList_bg_InM_zero_Both_mbias4.size(); i++) {
2125 }
2126 gDirectory->cd("..");
2127 gDirectory->mkdir("BG_Both_mbias4_one");
2128 gDirectory->cd("BG_Both_mbias4_one");
2129 for (UInt_t i = 0; i < fHistoList_bg_InM_one_Both_mbias4.size(); i++) {
2131 }
2132 gDirectory->cd("..");
2133 gDirectory->mkdir("BG_Both_mbias4_two");
2134 gDirectory->cd("BG_Both_mbias4_two");
2135 for (UInt_t i = 0; i < fHistoList_bg_InM_two_Both_mbias4.size(); i++) {
2137 }
2138 gDirectory->cd("..");
2139 gDirectory->mkdir("BG_Both_mbias4_onetwo");
2140 gDirectory->cd("BG_Both_mbias4_onetwo");
2141 for (UInt_t i = 0; i < fHistoList_bg_InM_onetwo_Both_mbias4.size(); i++) {
2143 }
2144 gDirectory->cd("..");
2145 gDirectory->cd("..");
2146 gDirectory->mkdir("all");
2147 gDirectory->cd("all");
2148 for (UInt_t i = 0; i < fHistoList_man_all_Both_mbias4.size(); i++) {
2150 }
2151 gDirectory->cd("..");
2152 gDirectory->mkdir("zero");
2153 gDirectory->cd("zero");
2154 for (UInt_t i = 0; i < fHistoList_man_zero_Both_mbias4.size(); i++) {
2156 }
2157 gDirectory->cd("..");
2158 gDirectory->mkdir("one");
2159 gDirectory->cd("one");
2160 for (UInt_t i = 0; i < fHistoList_man_one_Both_mbias4.size(); i++) {
2162 }
2163 gDirectory->cd("..");
2164 gDirectory->mkdir("two");
2165 gDirectory->cd("two");
2166 for (UInt_t i = 0; i < fHistoList_man_two_Both_mbias4.size(); i++) {
2168 }
2169 gDirectory->cd("..");
2170 gDirectory->mkdir("onetwo");
2171 gDirectory->cd("onetwo");
2172 for (UInt_t i = 0; i < fHistoList_man_onetwo_Both_mbias4.size(); i++) {
2174 }
2175 gDirectory->cd("..");
2176 for (UInt_t i = 0; i < fHistoList_man_Both_mbias4.size(); i++) {
2177 fHistoList_man_Both_mbias4[i]->Write();
2178 }
2179 gDirectory->cd("..");
2180}
2181
2183{
2185 InvMass_vs_OA_candidates_Both_mbias4 = new TH2D("InvMass_vs_OA_candidates_Both_mbias4",
2186 "InvMass_vs_OA_candidates_Both_mbias4; invariant mass in "
2187 "GeV/c^{2}; opening angle in degree",
2188 500, 0, 0.5, 500, 0, 50);
2190 InvMass_vs_OA_fromPi0_Both_mbias4 = new TH2D("InvMass_vs_OA_fromPi0_Both_mbias4",
2191 "InvMass_vs_OA_fromPi0_Both_mbias4; invariant mass in GeV/c^{2}; "
2192 "opening angle in degree",
2193 500, 0, 0.5, 500, 0, 50);
2196 new TH1D("GammasInvMass_candidates_Both_mbias4",
2197 "GammasInvMass_candidates_Both_mbias4; invariant mass in GeV/c^{2};#", 510, -0.01, 0.5);
2200 "GammasOA_candidates_Both_mbias4", "GammasOA_candidates_Both_mbias4; opening angle in degree;#", 300, -0.1, 29.9);
2203 new TH1D("GammasInvMass_fromPi0_Both_mbias4", "GammasInvMass_fromPi0_Both_mbias4; invariant mass in GeV/c^{2};#",
2204 510, -0.01, 0.5);
2206 GammasOA_fromPi0_Both_mbias4 = new TH1D("GammasOA_fromPi0_Both_mbias4",
2207 "GammasOA_fromPi0_Both_mbias4; opening angle in degree;#", 300, -0.1, 29.9);
2210 new TH1D("PlaneAngles_last_candidates_Both_mbias4",
2211 "PlaneAngles_last_candidates_Both_mbias4; #theta angle in degree;#", 720, -1., 179.);
2214 new TH1D("PlaneAngles_last_fromPi0_Both_mbias4", "PlaneAngles_last_fromPi0_Both_mbias4; #theta angle in degree;#",
2215 720, -1., 179.);
2218 new TH1D("PlaneAngles_first_candidates_Both_mbias4",
2219 "PlaneAngles_first_candidates_Both_mbias4; #theta angle in degree;#", 720, -1., 179.);
2222 new TH1D("PlaneAngles_first_fromPi0_Both_mbias4", "PlaneAngles_first_fromPi0_Both_mbias4; #theta angle in degree;#",
2223 720, -1., 179.);
2225
2226
2227 InvMass_vs_OA_candidates_Target_mbias4 = new TH2D("InvMass_vs_OA_candidates_Target_mbias4",
2228 "InvMass_vs_OA_candidates_Target_mbias4; invariant mass in "
2229 "GeV/c^{2}; opening angle in degree",
2230 500, 0, 0.5, 500, 0, 50);
2232 InvMass_vs_OA_fromPi0_Target_mbias4 = new TH2D("InvMass_vs_OA_fromPi0_Target_mbias4",
2233 "InvMass_vs_OA_fromPi0_Target_mbias4; invariant mass in "
2234 "GeV/c^{2}; opening angle in degree",
2235 500, 0, 0.5, 500, 0, 50);
2238 new TH1D("GammasInvMass_candidates_Target_mbias4",
2239 "GammasInvMass_candidates_Target_mbias4; invariant mass in GeV/c^{2};#", 510, -0.01, 0.5);
2242 new TH1D("GammasOA_candidates_Target_mbias4", "GammasOA_candidates_Target_mbias4; opening angle in degree;#", 300,
2243 -0.1, 29.9);
2246 new TH1D("GammasInvMass_fromPi0_Target_mbias4",
2247 "GammasInvMass_fromPi0_Target_mbias4; invariant mass in GeV/c^{2};#", 510, -0.01, 0.5);
2250 "GammasOA_fromPi0_Target_mbias4", "GammasOA_fromPi0_Target_mbias4; opening angle in degree;#", 300, -0.1, 29.9);
2253 new TH1D("PlaneAngles_last_candidates_Target_mbias4",
2254 "PlaneAngles_last_candidates_Target_mbias4; #theta angle in degree;#", 720, -1., 179.);
2257 new TH1D("PlaneAngles_last_fromPi0_Target_mbias4",
2258 "PlaneAngles_last_fromPi0_Target_mbias4; #theta angle in degree;#", 720, -1., 179.);
2261 new TH1D("PlaneAngles_first_candidates_Target_mbias4",
2262 "PlaneAngles_first_candidates_Target_mbias4; #theta angle in degree;#", 720, -1., 179.);
2265 new TH1D("PlaneAngles_first_fromPi0_Target_mbias4",
2266 "PlaneAngles_first_fromPi0_Target_mbias4; #theta angle in degree;#", 720, -1., 179.);
2268
2269
2270 InvMass_vs_OA_candidates_Outside_mbias4 = new TH2D("InvMass_vs_OA_candidates_Outside_mbias4",
2271 "InvMass_vs_OA_candidates_Outside_mbias4; invariant mass in "
2272 "GeV/c^{2}; opening angle in degree",
2273 500, 0, 0.5, 500, 0, 50);
2275 InvMass_vs_OA_fromPi0_Outside_mbias4 = new TH2D("InvMass_vs_OA_fromPi0_Outside_mbias4",
2276 "InvMass_vs_OA_fromPi0_Outside_mbias4; invariant mass in "
2277 "GeV/c^{2}; opening angle in degree",
2278 500, 0, 0.5, 500, 0, 50);
2281 new TH1D("GammasInvMass_candidates_Outside_mbias4",
2282 "GammasInvMass_candidates_Outside_mbias4; invariant mass in GeV/c^{2};#", 510, -0.01, 0.5);
2285 new TH1D("GammasOA_candidates_Outside_mbias4", "GammasOA_candidates_Outside_mbias4; opening angle in degree;#", 300,
2286 -0.1, 29.9);
2289 new TH1D("GammasInvMass_fromPi0_Outside_mbias4",
2290 "GammasInvMass_fromPi0_Outside_mbias4; invariant mass in GeV/c^{2};#", 510, -0.01, 0.5);
2293 "GammasOA_fromPi0_Outside_mbias4", "GammasOA_fromPi0_Outside_mbias4; opening angle in degree;#", 300, -0.1, 29.9);
2296 new TH1D("PlaneAngles_last_candidates_Outside_mbias4",
2297 "PlaneAngles_last_candidates_Outside_mbias4; #theta angle in degree;#", 720, -1., 179.);
2300 new TH1D("PlaneAngles_last_fromPi0_Outside_mbias4",
2301 "PlaneAngles_last_fromPi0_Outside_mbias4; #theta angle in degree;#", 720, -1., 179.);
2304 new TH1D("PlaneAngles_first_candidates_Outside_mbias4",
2305 "PlaneAngles_first_candidates_Outside_mbias4; #theta angle in degree;#", 720, -1., 179.);
2308 new TH1D("PlaneAngles_first_fromPi0_Outside_mbias4",
2309 "PlaneAngles_first_fromPi0_Outside_mbias4; #theta angle in degree;#", 720, -1., 179.);
2311
2312
2313 // Target_mbias4 => all
2315 new TH1D("GammaInvMassReco_all_Target_mbias4", "GammaInvMassReco_all_Target_mbias4; invariant mass in GeV/c^{2};#",
2316 110, -0.01, 0.1);
2319 new TH1D("GammaOpeningAngleReco_all_Target_mbias4", "GammaOpeningAngleReco_all_Target_mbias4; angle [deg];#", 200,
2320 -0.1, 19.9);
2322 Pdg_all_Target_mbias4 = new TH1D("Pdg_all_Target_mbias4", "Pdg_all_Target_mbias4; Id;#", 800, 0, 400);
2325 new TH1D("P_reco_all_Target_mbias4", "P_reco_all_Target_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2328 new TH1D("Pt_reco_all_Target_mbias4", "Pt_reco_all_Target_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2331 new TH1D("Pi0InvMassReco_all_Target_mbias4", "Pi0InvMassReco_all_Target_mbias4; invariant mass in GeV/c^{2};#",
2332 1000, 0, 2.0);
2335 "EMT_InvMass_all_Target_mbias4", "EMT_InvMass_all_Target_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2338 new TH2D("Pi0_pt_vs_rap_all_Target_mbias4", "Pi0_pt_vs_rap_all_Target_mbias4; rapidity y; p_{t} in GeV/c ", 90, -2.,
2339 7., 60, -1., 5.);
2342 new TH2D("Pi0_pt_vs_rap_est_all_Target_mbias4", "Pi0_pt_vs_rap_est_all_Target_mbias4; rapidity y; p_{t} in GeV/c ",
2343 10, 0., 4., 10, 0., 4.);
2345
2346
2347 // Target_mbias4 => zero
2349 new TH1D("GammaInvMassReco_zero_Target_mbias4",
2350 "GammaInvMassReco_zero_Target_mbias4; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2353 new TH1D("GammaOpeningAngleReco_zero_Target_mbias4", "GammaOpeningAngleReco_zero_Target_mbias4; angle [deg];#", 200,
2354 -0.1, 19.9);
2356 Pdg_zero_Target_mbias4 = new TH1D("Pdg_zero_Target_mbias4", "Pdg_zero_Target_mbias4; Id;#", 800, 0, 400);
2359 new TH1D("P_reco_zero_Target_mbias4", "P_reco_zero_Target_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2362 new TH1D("Pt_reco_zero_Target_mbias4", "Pt_reco_zero_Target_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2365 new TH1D("Pi0InvMassReco_zero_Target_mbias4", "Pi0InvMassReco_zero_Target_mbias4; invariant mass in GeV/c^{2};#",
2366 1000, 0, 2.0);
2369 "EMT_InvMass_zero_Target_mbias4", "EMT_InvMass_zero_Target_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2372 new TH2D("Pi0_pt_vs_rap_zero_Target_mbias4", "Pi0_pt_vs_rap_zero_Target_mbias4; rapidity y; p_{t} in GeV/c ", 90,
2373 -2., 7., 60, -1., 5.);
2376 new TH2D("Pi0_pt_vs_rap_est_zero_Target_mbias4",
2377 "Pi0_pt_vs_rap_est_zero_Target_mbias4; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2379
2380
2381 // Target_mbias4 => one
2383 new TH1D("GammaInvMassReco_one_Target_mbias4", "GammaInvMassReco_one_Target_mbias4; invariant mass in GeV/c^{2};#",
2384 110, -0.01, 0.1);
2387 new TH1D("GammaOpeningAngleReco_one_Target_mbias4", "GammaOpeningAngleReco_one_Target_mbias4; angle [deg];#", 200,
2388 -0.1, 19.9);
2390 Pdg_one_Target_mbias4 = new TH1D("Pdg_one_Target_mbias4", "Pdg_one_Target_mbias4; Id;#", 800, 0, 400);
2393 new TH1D("P_reco_one_Target_mbias4", "P_reco_one_Target_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2396 new TH1D("Pt_reco_one_Target_mbias4", "Pt_reco_one_Target_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2399 new TH1D("Pi0InvMassReco_one_Target_mbias4", "Pi0InvMassReco_one_Target_mbias4; invariant mass in GeV/c^{2};#",
2400 1000, 0, 2.0);
2403 "EMT_InvMass_one_Target_mbias4", "EMT_InvMass_one_Target_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2406 new TH2D("Pi0_pt_vs_rap_one_Target_mbias4", "Pi0_pt_vs_rap_one_Target_mbias4; rapidity y; p_{t} in GeV/c ", 90, -2.,
2407 7., 60, -1., 5.);
2410 new TH2D("Pi0_pt_vs_rap_est_one_Target_mbias4", "Pi0_pt_vs_rap_est_one_Target_mbias4; rapidity y; p_{t} in GeV/c ",
2411 10, 0., 4., 10, 0., 4.);
2413
2414
2415 // Target_mbias4 => two
2417 new TH1D("GammaInvMassReco_two_Target_mbias4", "GammaInvMassReco_two_Target_mbias4; invariant mass in GeV/c^{2};#",
2418 110, -0.01, 0.1);
2421 new TH1D("GammaOpeningAngleReco_two_Target_mbias4", "GammaOpeningAngleReco_two_Target_mbias4; angle [deg];#", 200,
2422 -0.1, 19.9);
2424 Pdg_two_Target_mbias4 = new TH1D("Pdg_two_Target_mbias4", "Pdg_two_Target_mbias4; Id;#", 800, 0, 400);
2427 new TH1D("P_reco_two_Target_mbias4", "P_reco_two_Target_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2430 new TH1D("Pt_reco_two_Target_mbias4", "Pt_reco_two_Target_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2433 new TH1D("Pi0InvMassReco_two_Target_mbias4", "Pi0InvMassReco_two_Target_mbias4; invariant mass in GeV/c^{2};#",
2434 1000, 0, 2.0);
2437 "EMT_InvMass_two_Target_mbias4", "EMT_InvMass_two_Target_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2440 new TH2D("Pi0_pt_vs_rap_two_Target_mbias4", "Pi0_pt_vs_rap_two_Target_mbias4; rapidity y; p_{t} in GeV/c ", 90, -2.,
2441 7., 60, -1., 5.);
2444 new TH2D("Pi0_pt_vs_rap_est_two_Target_mbias4", "Pi0_pt_vs_rap_est_two_Target_mbias4; rapidity y; p_{t} in GeV/c ",
2445 10, 0., 4., 10, 0., 4.);
2447
2448
2449 // Target_mbias4 => onetwo
2451 new TH1D("GammaInvMassReco_onetwo_Target_mbias4",
2452 "GammaInvMassReco_onetwo_Target_mbias4; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2455 new TH1D("GammaOpeningAngleReco_onetwo_Target_mbias4", "GammaOpeningAngleReco_onetwo_Target_mbias4; angle [deg];#",
2456 200, -0.1, 19.9);
2458 Pdg_onetwo_Target_mbias4 = new TH1D("Pdg_onetwo_Target_mbias4", "Pdg_onetwo_Target_mbias4; Id;#", 800, 0, 400);
2461 new TH1D("P_reco_onetwo_Target_mbias4", "P_reco_onetwo_Target_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2464 new TH1D("Pt_reco_onetwo_Target_mbias4", "Pt_reco_onetwo_Target_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2467 new TH1D("Pi0InvMassReco_onetwo_Target_mbias4",
2468 "Pi0InvMassReco_onetwo_Target_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2471 new TH1D("EMT_InvMass_onetwo_Target_mbias4", "EMT_InvMass_onetwo_Target_mbias4; invariant mass in GeV/c^{2};#",
2472 1000, 0, 2.0);
2475 new TH2D("Pi0_pt_vs_rap_onetwo_Target_mbias4", "Pi0_pt_vs_rap_onetwo_Target_mbias4; rapidity y; p_{t} in GeV/c ",
2476 90, -2., 7., 60, -1., 5.);
2479 new TH2D("Pi0_pt_vs_rap_est_onetwo_Target_mbias4",
2480 "Pi0_pt_vs_rap_est_onetwo_Target_mbias4; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2482
2483
2484 // Outside_mbias4 => all
2486 new TH1D("GammaInvMassReco_all_Outside_mbias4",
2487 "GammaInvMassReco_all_Outside_mbias4; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2490 new TH1D("GammaOpeningAngleReco_all_Outside_mbias4", "GammaOpeningAngleReco_all_Outside_mbias4; angle [deg];#", 200,
2491 -0.1, 19.9);
2493 Pdg_all_Outside_mbias4 = new TH1D("Pdg_all_Outside_mbias4", "Pdg_all_Outside_mbias4; Id;#", 800, 0, 400);
2496 new TH1D("P_reco_all_Outside_mbias4", "P_reco_all_Outside_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2499 new TH1D("Pt_reco_all_Outside_mbias4", "Pt_reco_all_Outside_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2502 new TH1D("Pi0InvMassReco_all_Outside_mbias4", "Pi0InvMassReco_all_Outside_mbias4; invariant mass in GeV/c^{2};#",
2503 1000, 0, 2.0);
2506 "EMT_InvMass_all_Outside_mbias4", "EMT_InvMass_all_Outside_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2509 new TH2D("Pi0_pt_vs_rap_all_Outside_mbias4", "Pi0_pt_vs_rap_all_Outside_mbias4; rapidity y; p_{t} in GeV/c ", 90,
2510 -2., 7., 60, -1., 5.);
2513 new TH2D("Pi0_pt_vs_rap_est_all_Outside_mbias4",
2514 "Pi0_pt_vs_rap_est_all_Outside_mbias4; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2516
2517
2518 // Outside_mbias4 => zero
2520 new TH1D("GammaInvMassReco_zero_Outside_mbias4",
2521 "GammaInvMassReco_zero_Outside_mbias4; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2524 new TH1D("GammaOpeningAngleReco_zero_Outside_mbias4", "GammaOpeningAngleReco_zero_Outside_mbias4; angle [deg];#",
2525 200, -0.1, 19.9);
2527 Pdg_zero_Outside_mbias4 = new TH1D("Pdg_zero_Outside_mbias4", "Pdg_zero_Outside_mbias4; Id;#", 800, 0, 400);
2530 new TH1D("P_reco_zero_Outside_mbias4", "P_reco_zero_Outside_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2533 new TH1D("Pt_reco_zero_Outside_mbias4", "Pt_reco_zero_Outside_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2536 new TH1D("Pi0InvMassReco_zero_Outside_mbias4", "Pi0InvMassReco_zero_Outside_mbias4; invariant mass in GeV/c^{2};#",
2537 1000, 0, 2.0);
2540 "EMT_InvMass_zero_Outside_mbias4", "EMT_InvMass_zero_Outside_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2543 new TH2D("Pi0_pt_vs_rap_zero_Outside_mbias4", "Pi0_pt_vs_rap_zero_Outside_mbias4; rapidity y; p_{t} in GeV/c ", 90,
2544 -2., 7., 60, -1., 5.);
2547 new TH2D("Pi0_pt_vs_rap_est_zero_Outside_mbias4",
2548 "Pi0_pt_vs_rap_est_zero_Outside_mbias4; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2550
2551
2552 // Outside_mbias4 => one
2554 new TH1D("GammaInvMassReco_one_Outside_mbias4",
2555 "GammaInvMassReco_one_Outside_mbias4; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2558 new TH1D("GammaOpeningAngleReco_one_Outside_mbias4", "GammaOpeningAngleReco_one_Outside_mbias4; angle [deg];#", 200,
2559 -0.1, 19.9);
2561 Pdg_one_Outside_mbias4 = new TH1D("Pdg_one_Outside_mbias4", "Pdg_one_Outside_mbias4; Id;#", 800, 0, 400);
2564 new TH1D("P_reco_one_Outside_mbias4", "P_reco_one_Outside_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2567 new TH1D("Pt_reco_one_Outside_mbias4", "Pt_reco_one_Outside_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2570 new TH1D("Pi0InvMassReco_one_Outside_mbias4", "Pi0InvMassReco_one_Outside_mbias4; invariant mass in GeV/c^{2};#",
2571 1000, 0, 2.0);
2574 "EMT_InvMass_one_Outside_mbias4", "EMT_InvMass_one_Outside_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2577 new TH2D("Pi0_pt_vs_rap_one_Outside_mbias4", "Pi0_pt_vs_rap_one_Outside_mbias4; rapidity y; p_{t} in GeV/c ", 90,
2578 -2., 7., 60, -1., 5.);
2581 new TH2D("Pi0_pt_vs_rap_est_one_Outside_mbias4",
2582 "Pi0_pt_vs_rap_est_one_Outside_mbias4; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2584
2585
2586 // Outside_mbias4 => two
2588 new TH1D("GammaInvMassReco_two_Outside_mbias4",
2589 "GammaInvMassReco_two_Outside_mbias4; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2592 new TH1D("GammaOpeningAngleReco_two_Outside_mbias4", "GammaOpeningAngleReco_two_Outside_mbias4; angle [deg];#", 200,
2593 -0.1, 19.9);
2595 Pdg_two_Outside_mbias4 = new TH1D("Pdg_two_Outside_mbias4", "Pdg_two_Outside_mbias4; Id;#", 800, 0, 400);
2598 new TH1D("P_reco_two_Outside_mbias4", "P_reco_two_Outside_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2601 new TH1D("Pt_reco_two_Outside_mbias4", "Pt_reco_two_Outside_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2604 new TH1D("Pi0InvMassReco_two_Outside_mbias4", "Pi0InvMassReco_two_Outside_mbias4; invariant mass in GeV/c^{2};#",
2605 1000, 0, 2.0);
2608 "EMT_InvMass_two_Outside_mbias4", "EMT_InvMass_two_Outside_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2611 new TH2D("Pi0_pt_vs_rap_two_Outside_mbias4", "Pi0_pt_vs_rap_two_Outside_mbias4; rapidity y; p_{t} in GeV/c ", 90,
2612 -2., 7., 60, -1., 5.);
2615 new TH2D("Pi0_pt_vs_rap_est_two_Outside_mbias4",
2616 "Pi0_pt_vs_rap_est_two_Outside_mbias4; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2618
2619
2620 // Outside_mbias4 => onetwo
2622 new TH1D("GammaInvMassReco_onetwo_Outside_mbias4",
2623 "GammaInvMassReco_onetwo_Outside_mbias4; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2626 new TH1D("GammaOpeningAngleReco_onetwo_Outside_mbias4",
2627 "GammaOpeningAngleReco_onetwo_Outside_mbias4; angle [deg];#", 200, -0.1, 19.9);
2629 Pdg_onetwo_Outside_mbias4 = new TH1D("Pdg_onetwo_Outside_mbias4", "Pdg_onetwo_Outside_mbias4; Id;#", 800, 0, 400);
2632 new TH1D("P_reco_onetwo_Outside_mbias4", "P_reco_onetwo_Outside_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2635 new TH1D("Pt_reco_onetwo_Outside_mbias4", "Pt_reco_onetwo_Outside_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2638 new TH1D("Pi0InvMassReco_onetwo_Outside_mbias4",
2639 "Pi0InvMassReco_onetwo_Outside_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2642 new TH1D("EMT_InvMass_onetwo_Outside_mbias4", "EMT_InvMass_onetwo_Outside_mbias4; invariant mass in GeV/c^{2};#",
2643 1000, 0, 2.0);
2646 new TH2D("Pi0_pt_vs_rap_onetwo_Outside_mbias4", "Pi0_pt_vs_rap_onetwo_Outside_mbias4; rapidity y; p_{t} in GeV/c ",
2647 90, -2., 7., 60, -1., 5.);
2650 new TH2D("Pi0_pt_vs_rap_est_onetwo_Outside_mbias4",
2651 "Pi0_pt_vs_rap_est_onetwo_Outside_mbias4; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2653
2654
2655 // Both_mbias4 => all
2657 new TH1D("GammaInvMassReco_all_Both_mbias4", "GammaInvMassReco_all_Both_mbias4; invariant mass in GeV/c^{2};#", 110,
2658 -0.01, 0.1);
2661 "GammaOpeningAngleReco_all_Both_mbias4", "GammaOpeningAngleReco_all_Both_mbias4; angle [deg];#", 200, -0.1, 19.9);
2663 Pdg_all_Both_mbias4 = new TH1D("Pdg_all_Both_mbias4", "Pdg_all_Both_mbias4; Id;#", 800, 0, 400);
2665 P_reco_all_Both_mbias4 = new TH1D("P_reco_all_Both_mbias4", "P_reco_all_Both_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2668 new TH1D("Pt_reco_all_Both_mbias4", "Pt_reco_all_Both_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2671 "Pi0InvMassReco_all_Both_mbias4", "Pi0InvMassReco_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2674 new TH1D("EMT_InvMass_all_Both_mbias4", "EMT_InvMass_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2677 new TH2D("Pi0_pt_vs_rap_all_Both_mbias4", "Pi0_pt_vs_rap_all_Both_mbias4; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2678 60, -1., 5.);
2681 new TH2D("Pi0_pt_vs_rap_est_all_Both_mbias4", "Pi0_pt_vs_rap_est_all_Both_mbias4; rapidity y; p_{t} in GeV/c ", 10,
2682 0., 4., 10, 0., 4.);
2684
2685
2686 // Both_mbias4 => zero
2688 new TH1D("GammaInvMassReco_zero_Both_mbias4", "GammaInvMassReco_zero_Both_mbias4; invariant mass in GeV/c^{2};#",
2689 110, -0.01, 0.1);
2692 "GammaOpeningAngleReco_zero_Both_mbias4", "GammaOpeningAngleReco_zero_Both_mbias4; angle [deg];#", 200, -0.1, 19.9);
2694 Pdg_zero_Both_mbias4 = new TH1D("Pdg_zero_Both_mbias4", "Pdg_zero_Both_mbias4; Id;#", 800, 0, 400);
2696 P_reco_zero_Both_mbias4 = new TH1D("P_reco_zero_Both_mbias4", "P_reco_zero_Both_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2699 new TH1D("Pt_reco_zero_Both_mbias4", "Pt_reco_zero_Both_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2702 "Pi0InvMassReco_zero_Both_mbias4", "Pi0InvMassReco_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2704 EMT_InvMass_zero_Both_mbias4 = new TH1D("EMT_InvMass_zero_Both_mbias4",
2705 "EMT_InvMass_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2708 new TH2D("Pi0_pt_vs_rap_zero_Both_mbias4", "Pi0_pt_vs_rap_zero_Both_mbias4; rapidity y; p_{t} in GeV/c ", 90, -2.,
2709 7., 60, -1., 5.);
2712 new TH2D("Pi0_pt_vs_rap_est_zero_Both_mbias4", "Pi0_pt_vs_rap_est_zero_Both_mbias4; rapidity y; p_{t} in GeV/c ",
2713 10, 0., 4., 10, 0., 4.);
2715
2716
2717 // Both_mbias4 => one
2719 new TH1D("GammaInvMassReco_one_Both_mbias4", "GammaInvMassReco_one_Both_mbias4; invariant mass in GeV/c^{2};#", 110,
2720 -0.01, 0.1);
2723 "GammaOpeningAngleReco_one_Both_mbias4", "GammaOpeningAngleReco_one_Both_mbias4; angle [deg];#", 200, -0.1, 19.9);
2725 Pdg_one_Both_mbias4 = new TH1D("Pdg_one_Both_mbias4", "Pdg_one_Both_mbias4; Id;#", 800, 0, 400);
2727 P_reco_one_Both_mbias4 = new TH1D("P_reco_one_Both_mbias4", "P_reco_one_Both_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2730 new TH1D("Pt_reco_one_Both_mbias4", "Pt_reco_one_Both_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2733 "Pi0InvMassReco_one_Both_mbias4", "Pi0InvMassReco_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2736 new TH1D("EMT_InvMass_one_Both_mbias4", "EMT_InvMass_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2739 new TH2D("Pi0_pt_vs_rap_one_Both_mbias4", "Pi0_pt_vs_rap_one_Both_mbias4; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2740 60, -1., 5.);
2743 new TH2D("Pi0_pt_vs_rap_est_one_Both_mbias4", "Pi0_pt_vs_rap_est_one_Both_mbias4; rapidity y; p_{t} in GeV/c ", 10,
2744 0., 4., 10, 0., 4.);
2746
2747
2748 // Both_mbias4 => two
2750 new TH1D("GammaInvMassReco_two_Both_mbias4", "GammaInvMassReco_two_Both_mbias4; invariant mass in GeV/c^{2};#", 110,
2751 -0.01, 0.1);
2754 "GammaOpeningAngleReco_two_Both_mbias4", "GammaOpeningAngleReco_two_Both_mbias4; angle [deg];#", 200, -0.1, 19.9);
2756 Pdg_two_Both_mbias4 = new TH1D("Pdg_two_Both_mbias4", "Pdg_two_Both_mbias4; Id;#", 800, 0, 400);
2758 P_reco_two_Both_mbias4 = new TH1D("P_reco_two_Both_mbias4", "P_reco_two_Both_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2761 new TH1D("Pt_reco_two_Both_mbias4", "Pt_reco_two_Both_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2764 "Pi0InvMassReco_two_Both_mbias4", "Pi0InvMassReco_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2767 new TH1D("EMT_InvMass_two_Both_mbias4", "EMT_InvMass_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2770 new TH2D("Pi0_pt_vs_rap_two_Both_mbias4", "Pi0_pt_vs_rap_two_Both_mbias4; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2771 60, -1., 5.);
2774 new TH2D("Pi0_pt_vs_rap_est_two_Both_mbias4", "Pi0_pt_vs_rap_est_two_Both_mbias4; rapidity y; p_{t} in GeV/c ", 10,
2775 0., 4., 10, 0., 4.);
2777
2778
2779 // Both_mbias4 => onetwo
2781 new TH1D("GammaInvMassReco_onetwo_Both_mbias4",
2782 "GammaInvMassReco_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2785 new TH1D("GammaOpeningAngleReco_onetwo_Both_mbias4", "GammaOpeningAngleReco_onetwo_Both_mbias4; angle [deg];#", 200,
2786 -0.1, 19.9);
2788 Pdg_onetwo_Both_mbias4 = new TH1D("Pdg_onetwo_Both_mbias4", "Pdg_onetwo_Both_mbias4; Id;#", 800, 0, 400);
2791 new TH1D("P_reco_onetwo_Both_mbias4", "P_reco_onetwo_Both_mbias4; P in GeV/c^{2};#", 600, 0, 6);
2794 new TH1D("Pt_reco_onetwo_Both_mbias4", "Pt_reco_onetwo_Both_mbias4; P_{t} in GeV/c^{2};#", 300, 0, 3);
2797 new TH1D("Pi0InvMassReco_onetwo_Both_mbias4", "Pi0InvMassReco_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#",
2798 1000, 0, 2.0);
2801 "EMT_InvMass_onetwo_Both_mbias4", "EMT_InvMass_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2804 new TH2D("Pi0_pt_vs_rap_onetwo_Both_mbias4", "Pi0_pt_vs_rap_onetwo_Both_mbias4; rapidity y; p_{t} in GeV/c ", 90,
2805 -2., 7., 60, -1., 5.);
2808 new TH2D("Pi0_pt_vs_rap_est_onetwo_Both_mbias4",
2809 "Pi0_pt_vs_rap_est_onetwo_Both_mbias4; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2811
2812
2813 // Both_mbias4
2815 new TH2D("Pdg_vs_Distance_mbias4", "Pdg_vs_Distance_mbias4; pdg; distance in cm", 2500, 0, 2499, 500, 0, 50);
2817 P_vs_Distance_mbias4 = new TH2D("P_vs_Distance_mbias4",
2818 "Distance between projected track and center of the ring (for e+ "
2819 "and e-); P in GeV/c^{2}; distance in cm",
2820 300, 0, 3, 300, 0, 15);
2822
2823
2824 // Multiplicity Target_mbias4
2825 MultiplicityGamma_all_Target_mbias4 = new TH2D("MultiplicityGamma_all_Target_mbias4",
2826 "MultiplicityGamma_all_Target_mbias4; Nof gammas in event; "
2827 "invariant mass in GeV/c^{2};#",
2828 400, 0, 30, 1000, 0, 2.0);
2830 MultiplicityGamma_zero_Target_mbias4 = new TH2D("MultiplicityGamma_zero_Target_mbias4",
2831 "MultiplicityGamma_zero_Target_mbias4; Nof gammas in event; "
2832 "invariant mass in GeV/c^{2};#",
2833 400, 0, 30, 1000, 0, 2.0);
2835 MultiplicityGamma_one_Target_mbias4 = new TH2D("MultiplicityGamma_one_Target_mbias4",
2836 "MultiplicityGamma_one_Target_mbias4; Nof gammas in event; "
2837 "invariant mass in GeV/c^{2};#",
2838 400, 0, 30, 1000, 0, 2.0);
2840 MultiplicityGamma_two_Target_mbias4 = new TH2D("MultiplicityGamma_two_Target_mbias4",
2841 "MultiplicityGamma_two_Target_mbias4; Nof gammas in event; "
2842 "invariant mass in GeV/c^{2};#",
2843 400, 0, 30, 1000, 0, 2.0);
2845 MultiplicityGamma_onetwo_Target_mbias4 = new TH2D("MultiplicityGamma_onetwo_Target_mbias4",
2846 "MultiplicityGamma_onetwo_Target_mbias4; Nof gammas in event; "
2847 "invariant mass in GeV/c^{2};#",
2848 400, 0, 30, 1000, 0, 2.0);
2850
2852 new TH2D("MultiplicityChargedParticles_all_Target_mbias4",
2853 "MultiplicityChargedParticles_all_Target_mbias4; Nof charged "
2854 "particles in event; invariant mass in GeV/c^{2};#",
2855 1000, 0, 1000, 1000, 0, 2.0);
2858 new TH2D("MultiplicityChargedParticles_zero_Target_mbias4",
2859 "MultiplicityChargedParticles_zero_Target_mbias4; Nof charged "
2860 "particles in event; invariant mass in GeV/c^{2};#",
2861 1000, 0, 1000, 1000, 0, 2.0);
2864 new TH2D("MultiplicityChargedParticles_one_Target_mbias4",
2865 "MultiplicityChargedParticles_one_Target_mbias4; Nof charged "
2866 "particles in event; invariant mass in GeV/c^{2};#",
2867 1000, 0, 1000, 1000, 0, 2.0);
2870 new TH2D("MultiplicityChargedParticles_two_Target_mbias4",
2871 "MultiplicityChargedParticles_two_Target_mbias4; Nof charged "
2872 "particles in event; invariant mass in GeV/c^{2};#",
2873 1000, 0, 1000, 1000, 0, 2.0);
2876 new TH2D("MultiplicityChargedParticles_onetwo_Target_mbias4",
2877 "MultiplicityChargedParticles_onetwo_Target_mbias4; Nof charged "
2878 "particles in event; invariant mass in GeV/c^{2};#",
2879 1000, 0, 1000, 1000, 0, 2.0);
2881
2882 // Multiplicity Outside_mbias4
2883 MultiplicityGamma_all_Outside_mbias4 = new TH2D("MultiplicityGamma_all_Outside_mbias4",
2884 "MultiplicityGamma_all_Outside_mbias4; Nof gammas in event; "
2885 "invariant mass in GeV/c^{2};#",
2886 400, 0, 30, 1000, 0, 2.0);
2888 MultiplicityGamma_zero_Outside_mbias4 = new TH2D("MultiplicityGamma_zero_Outside_mbias4",
2889 "MultiplicityGamma_zero_Outside_mbias4; Nof gammas in event; "
2890 "invariant mass in GeV/c^{2};#",
2891 400, 0, 30, 1000, 0, 2.0);
2893 MultiplicityGamma_one_Outside_mbias4 = new TH2D("MultiplicityGamma_one_Outside_mbias4",
2894 "MultiplicityGamma_one_Outside_mbias4; Nof gammas in event; "
2895 "invariant mass in GeV/c^{2};#",
2896 400, 0, 30, 1000, 0, 2.0);
2898 MultiplicityGamma_two_Outside_mbias4 = new TH2D("MultiplicityGamma_two_Outside_mbias4",
2899 "MultiplicityGamma_two_Outside_mbias4; Nof gammas in event; "
2900 "invariant mass in GeV/c^{2};#",
2901 400, 0, 30, 1000, 0, 2.0);
2903 MultiplicityGamma_onetwo_Outside_mbias4 = new TH2D("MultiplicityGamma_onetwo_Outside_mbias4",
2904 "MultiplicityGamma_onetwo_Outside_mbias4; Nof gammas in event; "
2905 "invariant mass in GeV/c^{2};#",
2906 400, 0, 30, 1000, 0, 2.0);
2908
2910 new TH2D("MultiplicityChargedParticles_all_Outside_mbias4",
2911 "MultiplicityChargedParticles_all_Outside_mbias4; Nof charged "
2912 "particles in event; invariant mass in GeV/c^{2};#",
2913 1000, 0, 1000, 1000, 0, 2.0);
2916 new TH2D("MultiplicityChargedParticles_zero_Outside_mbias4",
2917 "MultiplicityChargedParticles_zero_Outside_mbias4; Nof charged "
2918 "particles in event; invariant mass in GeV/c^{2};#",
2919 1000, 0, 1000, 1000, 0, 2.0);
2922 new TH2D("MultiplicityChargedParticles_one_Outside_mbias4",
2923 "MultiplicityChargedParticles_one_Outside_mbias4; Nof charged "
2924 "particles in event; invariant mass in GeV/c^{2};#",
2925 1000, 0, 1000, 1000, 0, 2.0);
2928 new TH2D("MultiplicityChargedParticles_two_Outside_mbias4",
2929 "MultiplicityChargedParticles_two_Outside_mbias4; Nof charged "
2930 "particles in event; invariant mass in GeV/c^{2};#",
2931 1000, 0, 1000, 1000, 0, 2.0);
2934 new TH2D("MultiplicityChargedParticles_onetwo_Outside_mbias4",
2935 "MultiplicityChargedParticles_onetwo_Outside_mbias4; Nof charged "
2936 "particles in event; invariant mass in GeV/c^{2};#",
2937 1000, 0, 1000, 1000, 0, 2.0);
2939
2940
2941 // Multiplicity Both_mbias4
2942 MultiplicityGamma_all_Both_mbias4 = new TH2D("MultiplicityGamma_all_Both_mbias4",
2943 "MultiplicityGamma_all_Both_mbias4; Nof gammas in event; "
2944 "invariant mass in GeV/c^{2};#",
2945 400, 0, 30, 1000, 0, 2.0);
2947 MultiplicityGamma_zero_Both_mbias4 = new TH2D("MultiplicityGamma_zero_Both_mbias4",
2948 "MultiplicityGamma_zero_Both_mbias4; Nof gammas in event; "
2949 "invariant mass in GeV/c^{2};#",
2950 400, 0, 30, 1000, 0, 2.0);
2952 MultiplicityGamma_one_Both_mbias4 = new TH2D("MultiplicityGamma_one_Both_mbias4",
2953 "MultiplicityGamma_one_Both_mbias4; Nof gammas in event; "
2954 "invariant mass in GeV/c^{2};#",
2955 400, 0, 30, 1000, 0, 2.0);
2957 MultiplicityGamma_two_Both_mbias4 = new TH2D("MultiplicityGamma_two_Both_mbias4",
2958 "MultiplicityGamma_two_Both_mbias4; Nof gammas in event; "
2959 "invariant mass in GeV/c^{2};#",
2960 400, 0, 30, 1000, 0, 2.0);
2962 MultiplicityGamma_onetwo_Both_mbias4 = new TH2D("MultiplicityGamma_onetwo_Both_mbias4",
2963 "MultiplicityGamma_onetwo_Both_mbias4; Nof gammas in event; "
2964 "invariant mass in GeV/c^{2};#",
2965 400, 0, 30, 1000, 0, 2.0);
2967
2968 MultiplicityChargedParticles_all_Both_mbias4 = new TH2D("MultiplicityChargedParticles_all_Both_mbias4",
2969 "MultiplicityChargedParticles_all_Both_mbias4; Nof charged "
2970 "particles in event; invariant mass in GeV/c^{2};#",
2971 1000, 0, 1000, 1000, 0, 2.0);
2973 MultiplicityChargedParticles_zero_Both_mbias4 = new TH2D("MultiplicityChargedParticles_zero_Both_mbias4",
2974 "MultiplicityChargedParticles_zero_Both_mbias4; Nof charged "
2975 "particles in event; invariant mass in GeV/c^{2};#",
2976 1000, 0, 1000, 1000, 0, 2.0);
2978 MultiplicityChargedParticles_one_Both_mbias4 = new TH2D("MultiplicityChargedParticles_one_Both_mbias4",
2979 "MultiplicityChargedParticles_one_Both_mbias4; Nof charged "
2980 "particles in event; invariant mass in GeV/c^{2};#",
2981 1000, 0, 1000, 1000, 0, 2.0);
2983 MultiplicityChargedParticles_two_Both_mbias4 = new TH2D("MultiplicityChargedParticles_two_Both_mbias4",
2984 "MultiplicityChargedParticles_two_Both_mbias4; Nof charged "
2985 "particles in event; invariant mass in GeV/c^{2};#",
2986 1000, 0, 1000, 1000, 0, 2.0);
2989 new TH2D("MultiplicityChargedParticles_onetwo_Both_mbias4",
2990 "MultiplicityChargedParticles_onetwo_Both_mbias4; Nof charged "
2991 "particles in event; invariant mass in GeV/c^{2};#",
2992 1000, 0, 1000, 1000, 0, 2.0);
2994
2995
2996 // rap_vs_Pt for "OneTwo" and "Both_mbias4"
2998 rap_vs_Pt_InM_1_mbias4 = new TH1D("rap_vs_Pt_InM_1_mbias4",
2999 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3000 "mass in GeV/c^{2};#",
3001 1000, 0, 2.0);
3003 rap_vs_Pt_InM_2_mbias4 = new TH1D("rap_vs_Pt_InM_2_mbias4",
3004 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3005 "mass in GeV/c^{2};#",
3006 1000, 0, 2.0);
3008 rap_vs_Pt_InM_3_mbias4 = new TH1D("rap_vs_Pt_InM_3_mbias4",
3009 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3010 "mass in GeV/c^{2};#",
3011 1000, 0, 2.0);
3013 rap_vs_Pt_InM_4_mbias4 = new TH1D("rap_vs_Pt_InM_4_mbias4",
3014 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3015 "mass in GeV/c^{2};#",
3016 1000, 0, 2.0);
3018 rap_vs_Pt_InM_5_mbias4 = new TH1D("rap_vs_Pt_InM_5_mbias4",
3019 "rapidity = (1.2-1.6) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3020 "mass in GeV/c^{2};#",
3021 1000, 0, 2.0);
3023 rap_vs_Pt_InM_6_mbias4 = new TH1D("rap_vs_Pt_InM_6_mbias4",
3024 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3025 "mass in GeV/c^{2};#",
3026 1000, 0, 2.0);
3028 rap_vs_Pt_InM_7_mbias4 = new TH1D("rap_vs_Pt_InM_7_mbias4",
3029 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3030 "mass in GeV/c^{2};#",
3031 1000, 0, 2.0);
3033 rap_vs_Pt_InM_8_mbias4 = new TH1D("rap_vs_Pt_InM_8_mbias4",
3034 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3035 "mass in GeV/c^{2};#",
3036 1000, 0, 2.0);
3038 rap_vs_Pt_InM_9_mbias4 = new TH1D("rap_vs_Pt_InM_9_mbias4",
3039 "rapidity = (1.6-2.0) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3040 "mass in GeV/c^{2};#",
3041 1000, 0, 2.0);
3043 rap_vs_Pt_InM_10_mbias4 = new TH1D("rap_vs_Pt_InM_10_mbias4",
3044 "rapidity = (1.6-2.0) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3045 "mass in GeV/c^{2};#",
3046 1000, 0, 2.0);
3048 rap_vs_Pt_InM_11_mbias4 = new TH1D("rap_vs_Pt_InM_11_mbias4",
3049 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3050 "mass in GeV/c^{2};#",
3051 1000, 0, 2.0);
3053 rap_vs_Pt_InM_12_mbias4 = new TH1D("rap_vs_Pt_InM_12_mbias4",
3054 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3055 "mass in GeV/c^{2};#",
3056 1000, 0, 2.0);
3058 rap_vs_Pt_InM_13_mbias4 = new TH1D("rap_vs_Pt_InM_13_mbias4",
3059 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3060 "mass in GeV/c^{2};#",
3061 1000, 0, 2.0);
3063 rap_vs_Pt_InM_14_mbias4 = new TH1D("rap_vs_Pt_InM_14_mbias4",
3064 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3065 "mass in GeV/c^{2};#",
3066 1000, 0, 2.0);
3068 rap_vs_Pt_InM_15_mbias4 = new TH1D("rap_vs_Pt_InM_15_mbias4",
3069 "rapidity = (2.0-2.4) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3070 "mass in GeV/c^{2};#",
3071 1000, 0, 2.0);
3073 rap_vs_Pt_InM_16_mbias4 = new TH1D("rap_vs_Pt_InM_16_mbias4",
3074 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3075 "mass in GeV/c^{2};#",
3076 1000, 0, 2.0);
3078 rap_vs_Pt_InM_17_mbias4 = new TH1D("rap_vs_Pt_InM_17_mbias4",
3079 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3080 "mass in GeV/c^{2};#",
3081 1000, 0, 2.0);
3083 rap_vs_Pt_InM_18_mbias4 = new TH1D("rap_vs_Pt_InM_18_mbias4",
3084 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3085 "mass in GeV/c^{2};#",
3086 1000, 0, 2.0);
3088 rap_vs_Pt_InM_19_mbias4 = new TH1D("rap_vs_Pt_InM_19_mbias4",
3089 "rapidity = (2.4-2.8) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3090 "mass in GeV/c^{2};#",
3091 1000, 0, 2.0);
3093 rap_vs_Pt_InM_20_mbias4 = new TH1D("rap_vs_Pt_InM_20_mbias4",
3094 "rapidity = (2.4-2.8) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3095 "mass in GeV/c^{2};#",
3096 1000, 0, 2.0);
3098 rap_vs_Pt_InM_21_mbias4 = new TH1D("rap_vs_Pt_InM_21_mbias4",
3099 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3100 "mass in GeV/c^{2};#",
3101 1000, 0, 2.0);
3103 rap_vs_Pt_InM_22_mbias4 = new TH1D("rap_vs_Pt_InM_22_mbias4",
3104 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3105 "mass in GeV/c^{2};#",
3106 1000, 0, 2.0);
3108 rap_vs_Pt_InM_23_mbias4 = new TH1D("rap_vs_Pt_InM_23_mbias4",
3109 "rapidity = (2.8-3.2) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3110 "mass in GeV/c^{2};#",
3111 1000, 0, 2.0);
3113 rap_vs_Pt_InM_24_mbias4 = new TH1D("rap_vs_Pt_InM_24_mbias4",
3114 "rapidity = (2.8-3.2) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3115 "mass in GeV/c^{2};#",
3116 1000, 0, 2.0);
3118 rap_vs_Pt_InM_25_mbias4 = new TH1D("rap_vs_Pt_InM_25_mbias4",
3119 "rapidity = (2.8-3.2) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3120 "mass in GeV/c^{2};#",
3121 1000, 0, 2.0);
3123 rap_vs_Pt_InM_26_mbias4 = new TH1D("rap_vs_Pt_InM_26_mbias4",
3124 "rapidity = (3.2-3.6) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3125 "mass in GeV/c^{2};#",
3126 1000, 0, 2.0);
3128 rap_vs_Pt_InM_27_mbias4 = new TH1D("rap_vs_Pt_InM_27_mbias4",
3129 "rapidity = (3.2-3.6) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3130 "mass in GeV/c^{2};#",
3131 1000, 0, 2.0);
3133 rap_vs_Pt_InM_28_mbias4 = new TH1D("rap_vs_Pt_InM_28_mbias4",
3134 "rapidity = (3.2-3.6) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3135 "mass in GeV/c^{2};#",
3136 1000, 0, 2.0);
3138 rap_vs_Pt_InM_29_mbias4 = new TH1D("rap_vs_Pt_InM_29_mbias4",
3139 "rapidity = (3.2-3.6) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3140 "mass in GeV/c^{2};#",
3141 1000, 0, 2.0);
3143 rap_vs_Pt_InM_30_mbias4 = new TH1D("rap_vs_Pt_InM_30_mbias4",
3144 "rapidity = (3.2-3.6) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3145 "mass in GeV/c^{2};#",
3146 1000, 0, 2.0);
3148
3149 rap_vs_Pt_InM_mixing_1_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_1_mbias4",
3150 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3151 "mass in GeV/c^{2};#",
3152 1000, 0, 2.0);
3154 rap_vs_Pt_InM_mixing_2_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_2_mbias4",
3155 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3156 "mass in GeV/c^{2};#",
3157 1000, 0, 2.0);
3159 rap_vs_Pt_InM_mixing_3_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_3_mbias4",
3160 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3161 "mass in GeV/c^{2};#",
3162 1000, 0, 2.0);
3164 rap_vs_Pt_InM_mixing_4_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_4_mbias4",
3165 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3166 "mass in GeV/c^{2};#",
3167 1000, 0, 2.0);
3169 rap_vs_Pt_InM_mixing_5_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_5_mbias4",
3170 "rapidity = (1.2-1.6) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3171 "mass in GeV/c^{2};#",
3172 1000, 0, 2.0);
3174 rap_vs_Pt_InM_mixing_6_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_6_mbias4",
3175 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3176 "mass in GeV/c^{2};#",
3177 1000, 0, 2.0);
3179 rap_vs_Pt_InM_mixing_7_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_7_mbias4",
3180 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3181 "mass in GeV/c^{2};#",
3182 1000, 0, 2.0);
3184 rap_vs_Pt_InM_mixing_8_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_8_mbias4",
3185 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3186 "mass in GeV/c^{2};#",
3187 1000, 0, 2.0);
3189 rap_vs_Pt_InM_mixing_9_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_9_mbias4",
3190 "rapidity = (1.6-2.0) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3191 "mass in GeV/c^{2};#",
3192 1000, 0, 2.0);
3194 rap_vs_Pt_InM_mixing_10_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_10_mbias4",
3195 "rapidity = (1.6-2.0) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3196 "mass in GeV/c^{2};#",
3197 1000, 0, 2.0);
3199 rap_vs_Pt_InM_mixing_11_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_11_mbias4",
3200 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3201 "mass in GeV/c^{2};#",
3202 1000, 0, 2.0);
3204 rap_vs_Pt_InM_mixing_12_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_12_mbias4",
3205 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3206 "mass in GeV/c^{2};#",
3207 1000, 0, 2.0);
3209 rap_vs_Pt_InM_mixing_13_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_13_mbias4",
3210 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3211 "mass in GeV/c^{2};#",
3212 1000, 0, 2.0);
3214 rap_vs_Pt_InM_mixing_14_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_14_mbias4",
3215 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3216 "mass in GeV/c^{2};#",
3217 1000, 0, 2.0);
3219 rap_vs_Pt_InM_mixing_15_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_15_mbias4",
3220 "rapidity = (2.0-2.4) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3221 "mass in GeV/c^{2};#",
3222 1000, 0, 2.0);
3224 rap_vs_Pt_InM_mixing_16_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_16_mbias4",
3225 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3226 "mass in GeV/c^{2};#",
3227 1000, 0, 2.0);
3229 rap_vs_Pt_InM_mixing_17_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_17_mbias4",
3230 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3231 "mass in GeV/c^{2};#",
3232 1000, 0, 2.0);
3234 rap_vs_Pt_InM_mixing_18_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_18_mbias4",
3235 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3236 "mass in GeV/c^{2};#",
3237 1000, 0, 2.0);
3239 rap_vs_Pt_InM_mixing_19_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_19_mbias4",
3240 "rapidity = (2.4-2.8) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3241 "mass in GeV/c^{2};#",
3242 1000, 0, 2.0);
3244 rap_vs_Pt_InM_mixing_20_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_20_mbias4",
3245 "rapidity = (2.4-2.8) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3246 "mass in GeV/c^{2};#",
3247 1000, 0, 2.0);
3249 rap_vs_Pt_InM_mixing_21_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_21_mbias4",
3250 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3251 "mass in GeV/c^{2};#",
3252 1000, 0, 2.0);
3254 rap_vs_Pt_InM_mixing_22_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_22_mbias4",
3255 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3256 "mass in GeV/c^{2};#",
3257 1000, 0, 2.0);
3259 rap_vs_Pt_InM_mixing_23_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_23_mbias4",
3260 "rapidity = (2.8-3.2) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3261 "mass in GeV/c^{2};#",
3262 1000, 0, 2.0);
3264 rap_vs_Pt_InM_mixing_24_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_24_mbias4",
3265 "rapidity = (2.8-3.2) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3266 "mass in GeV/c^{2};#",
3267 1000, 0, 2.0);
3269 rap_vs_Pt_InM_mixing_25_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_25_mbias4",
3270 "rapidity = (2.8-3.2) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3271 "mass in GeV/c^{2};#",
3272 1000, 0, 2.0);
3274 rap_vs_Pt_InM_mixing_26_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_26_mbias4",
3275 "rapidity = (3.2-3.6) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3276 "mass in GeV/c^{2};#",
3277 1000, 0, 2.0);
3279 rap_vs_Pt_InM_mixing_27_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_27_mbias4",
3280 "rapidity = (3.2-3.6) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3281 "mass in GeV/c^{2};#",
3282 1000, 0, 2.0);
3284 rap_vs_Pt_InM_mixing_28_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_28_mbias4",
3285 "rapidity = (3.2-3.6) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3286 "mass in GeV/c^{2};#",
3287 1000, 0, 2.0);
3289 rap_vs_Pt_InM_mixing_29_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_29_mbias4",
3290 "rapidity = (3.2-3.6) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3291 "mass in GeV/c^{2};#",
3292 1000, 0, 2.0);
3294 rap_vs_Pt_InM_mixing_30_mbias4 = new TH1D("rap_vs_Pt_InM_mixing_30_mbias4",
3295 "rapidity = (3.2-3.6) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3296 "mass in GeV/c^{2};#",
3297 1000, 0, 2.0);
3299
3300
3301 // BG cases
3302 //Both_mbias4 all
3304 new TH1D("BG1_InM_all_Both_mbias4", "BG1_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3307 new TH1D("BG2_InM_all_Both_mbias4", "BG2_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3310 new TH1D("BG3_InM_all_Both_mbias4", "BG3_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3313 new TH1D("BG4_InM_all_Both_mbias4", "BG4_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3316 new TH1D("BG5_InM_all_Both_mbias4", "BG5_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3319 new TH1D("BG6_InM_all_Both_mbias4", "BG6_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3322 new TH1D("BG7_InM_all_Both_mbias4", "BG7_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3325 new TH1D("BG8_InM_all_Both_mbias4", "BG8_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3328 new TH1D("BG9_InM_all_Both_mbias4", "BG9_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3331 new TH1D("BG10_InM_all_Both_mbias4", "BG10_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3334 new TH1D("PdgCase8_InM_all_Both_mbias4", "PdgCase8_InM_all_Both_mbias4; Id ;#", 5000, -2500, 2500);
3337 new TH1D("PdgCase8mothers_InM_all_Both_mbias4", "PdgCase8mothers_InM_all_Both_mbias4; Id ;#", 5000, -2500, 2500);
3340 new TH1D("sameMIDcase8_InM_all_Both_mbias4", "sameMIDcase8_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#",
3341 1000, 0, 2.0);
3344 new TH1D("sameGRIDcase8_InM_all_Both_mbias4", "sameGRIDcase8_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#",
3345 1000, 0, 2.0);
3347 Case1ZYPos_InM_all_Both_mbias4 = new TH2D("Case1ZYPos_InM_all_Both_mbias4",
3348 "Case1ZYPos_InM_all_Both_mbias4; z[cm]; y[cm]", 400, -1, 200, 200, -50, 50);
3351 new TH1D("sameMIDcase8_mothedPDG_InM_all_Both_mbias4", "sameMIDcase8_mothedPDG_InM_all_Both_mbias4; Id ;#", 5000,
3352 -2500, 2500);
3355 new TH1D("PdgCase8NonEComeFromTarget_mbias4_InM_all_Both_mbias4",
3356 "PdgCase8NonEComeFromTarget_mbias4_InM_all_Both_mbias4; Id ;#", 5000, -2500, 2500);
3359 new TH1D("PdgCase8NonE_NOT_FromTarget_mbias4_InM_all_Both_mbias4",
3360 "PdgCase8NonE_NOT_FromTarget_mbias4_InM_all_Both_mbias4; Id ;#", 5000, -2500, 2500);
3362 PdgCase8motherNonE_InM_all_Both_mbias4 = new TH1D("PdgCase8motherNonE_InM_all_Both_mbias4",
3363 "PdgCase8motherNonE_InM_all_Both_mbias4; Id ;#", 5000, -2500, 2500);
3366 new TH1D("Case8ElFromDalitz_InM_all_Both_mbias4",
3367 "Case8ElFromDalitz_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3370 new TH1D("Case8NonElFrom_pn_InM_all_Both_mbias4",
3371 "Case8NonElFrom_pn_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3374 new TH1D("Case8NonElFrom_eta_InM_all_Both_mbias4",
3375 "Case8NonElFrom_eta_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3378 new TH1D("Case8NonElFrom_kaon_InM_all_Both_mbias4",
3379 "Case8NonElFrom_kaon_InM_all_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3382 "sameMIDcase8NonEPdg_InM_all_Both_mbias4", "sameMIDcase8NonEPdg_InM_all_Both_mbias4; Id ;#", 5000, -2500, 2500);
3385 new TH1D("sameMIDcase8NonEMotherPdg_InM_all_Both_mbias4", "sameMIDcase8NonEMotherPdg_InM_all_Both_mbias4; Id ;#",
3386 5000, -2500, 2500);
3389 new TH1D("sameMIDcase8NonEMotherIM_InM_all_Both_mbias4",
3390 "sameMIDcase8NonEMotherIM_InM_all_Both_mbias4; invariant mass in "
3391 "GeV/c^{2};#",
3392 1000, 0, 2.0);
3395 new TH1D("sameMIDcase8NonEPdgFromTarget_mbias4_InM_all_Both_mbias4",
3396 "sameMIDcase8NonEPdgFromTarget_mbias4_InM_all_Both_mbias4; Id ;#", 5000, -2500, 2500);
3399 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4IM_InM_all_Both_mbias4",
3400 "sameMIDcase8NonEComeFromTarget_mbias4IM_InM_all_Both_mbias4; "
3401 "invariant mass in GeV/c^{2};#",
3402 1000, 0, 2.0);
3405 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4P_InM_all_Both_mbias4",
3406 "sameMIDcase8NonEComeFromTarget_mbias4P_InM_all_Both_mbias4; P in "
3407 "GeV/c^{2} ;#",
3408 200, 0, 10);
3411 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_all_Both_mbias4",
3412 "sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_all_Both_mbias4; Pt "
3413 "in GeV/c^{2} ;#",
3414 200, 0, 10);
3416 //Both_mbias4 zero
3418 new TH1D("BG1_InM_zero_Both_mbias4", "BG1_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3421 new TH1D("BG2_InM_zero_Both_mbias4", "BG2_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3424 new TH1D("BG3_InM_zero_Both_mbias4", "BG3_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3427 new TH1D("BG4_InM_zero_Both_mbias4", "BG4_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3430 new TH1D("BG5_InM_zero_Both_mbias4", "BG5_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3433 new TH1D("BG6_InM_zero_Both_mbias4", "BG6_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3436 new TH1D("BG7_InM_zero_Both_mbias4", "BG7_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3439 new TH1D("BG8_InM_zero_Both_mbias4", "BG8_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3442 new TH1D("BG9_InM_zero_Both_mbias4", "BG9_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3445 new TH1D("BG10_InM_zero_Both_mbias4", "BG10_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3448 new TH1D("PdgCase8_InM_zero_Both_mbias4", "PdgCase8_InM_zero_Both_mbias4; Id ;#", 5000, -2500, 2500);
3451 new TH1D("PdgCase8mothers_InM_zero_Both_mbias4", "PdgCase8mothers_InM_zero_Both_mbias4; Id ;#", 5000, -2500, 2500);
3454 new TH1D("sameMIDcase8_InM_zero_Both_mbias4", "sameMIDcase8_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#",
3455 1000, 0, 2.0);
3458 new TH1D("sameGRIDcase8_InM_zero_Both_mbias4", "sameGRIDcase8_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#",
3459 1000, 0, 2.0);
3462 "Case1ZYPos_InM_zero_Both_mbias4", "Case1ZYPos_InM_zero_Both_mbias4; z[cm]; y[cm]", 400, -1, 200, 200, -50, 50);
3465 new TH1D("sameMIDcase8_mothedPDG_InM_zero_Both_mbias4", "sameMIDcase8_mothedPDG_InM_zero_Both_mbias4; Id ;#", 5000,
3466 -2500, 2500);
3469 new TH1D("PdgCase8NonEComeFromTarget_mbias4_InM_zero_Both_mbias4",
3470 "PdgCase8NonEComeFromTarget_mbias4_InM_zero_Both_mbias4; Id ;#", 5000, -2500, 2500);
3473 new TH1D("PdgCase8NonE_NOT_FromTarget_mbias4_InM_zero_Both_mbias4",
3474 "PdgCase8NonE_NOT_FromTarget_mbias4_InM_zero_Both_mbias4; Id ;#", 5000, -2500, 2500);
3477 "PdgCase8motherNonE_InM_zero_Both_mbias4", "PdgCase8motherNonE_InM_zero_Both_mbias4; Id ;#", 5000, -2500, 2500);
3480 new TH1D("Case8ElFromDalitz_InM_zero_Both_mbias4",
3481 "Case8ElFromDalitz_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3484 new TH1D("Case8NonElFrom_pn_InM_zero_Both_mbias4",
3485 "Case8NonElFrom_pn_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3488 new TH1D("Case8NonElFrom_eta_InM_zero_Both_mbias4",
3489 "Case8NonElFrom_eta_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3492 new TH1D("Case8NonElFrom_kaon_InM_zero_Both_mbias4",
3493 "Case8NonElFrom_kaon_InM_zero_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3496 "sameMIDcase8NonEPdg_InM_zero_Both_mbias4", "sameMIDcase8NonEPdg_InM_zero_Both_mbias4; Id ;#", 5000, -2500, 2500);
3499 new TH1D("sameMIDcase8NonEMotherPdg_InM_zero_Both_mbias4", "sameMIDcase8NonEMotherPdg_InM_zero_Both_mbias4; Id ;#",
3500 5000, -2500, 2500);
3503 new TH1D("sameMIDcase8NonEMotherIM_InM_zero_Both_mbias4",
3504 "sameMIDcase8NonEMotherIM_InM_zero_Both_mbias4; invariant mass in "
3505 "GeV/c^{2};#",
3506 1000, 0, 2.0);
3509 new TH1D("sameMIDcase8NonEPdgFromTarget_mbias4_InM_zero_Both_mbias4",
3510 "sameMIDcase8NonEPdgFromTarget_mbias4_InM_zero_Both_mbias4; Id ;#", 5000, -2500, 2500);
3513 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4IM_InM_zero_Both_mbias4",
3514 "sameMIDcase8NonEComeFromTarget_mbias4IM_InM_zero_Both_mbias4; "
3515 "invariant mass in GeV/c^{2};#",
3516 1000, 0, 2.0);
3519 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4P_InM_zero_Both_mbias4",
3520 "sameMIDcase8NonEComeFromTarget_mbias4P_InM_zero_Both_mbias4; P "
3521 "in GeV/c^{2} ;#",
3522 200, 0, 10);
3525 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_zero_Both_mbias4",
3526 "sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_zero_Both_mbias4; Pt "
3527 "in GeV/c^{2} ;#",
3528 200, 0, 10);
3530 //Both_mbias4 one
3532 new TH1D("BG1_InM_one_Both_mbias4", "BG1_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3535 new TH1D("BG2_InM_one_Both_mbias4", "BG2_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3538 new TH1D("BG3_InM_one_Both_mbias4", "BG3_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3541 new TH1D("BG4_InM_one_Both_mbias4", "BG4_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3544 new TH1D("BG5_InM_one_Both_mbias4", "BG5_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3547 new TH1D("BG6_InM_one_Both_mbias4", "BG6_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3550 new TH1D("BG7_InM_one_Both_mbias4", "BG7_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3553 new TH1D("BG8_InM_one_Both_mbias4", "BG8_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3556 new TH1D("BG9_InM_one_Both_mbias4", "BG9_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3559 new TH1D("BG10_InM_one_Both_mbias4", "BG10_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3562 new TH1D("PdgCase8_InM_one_Both_mbias4", "PdgCase8_InM_one_Both_mbias4; Id ;#", 5000, -2500, 2500);
3565 new TH1D("PdgCase8mothers_InM_one_Both_mbias4", "PdgCase8mothers_InM_one_Both_mbias4; Id ;#", 5000, -2500, 2500);
3568 new TH1D("sameMIDcase8_InM_one_Both_mbias4", "sameMIDcase8_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#",
3569 1000, 0, 2.0);
3572 new TH1D("sameGRIDcase8_InM_one_Both_mbias4", "sameGRIDcase8_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#",
3573 1000, 0, 2.0);
3575 Case1ZYPos_InM_one_Both_mbias4 = new TH2D("Case1ZYPos_InM_one_Both_mbias4",
3576 "Case1ZYPos_InM_one_Both_mbias4; z[cm]; y[cm]", 400, -1, 200, 200, -50, 50);
3579 new TH1D("sameMIDcase8_mothedPDG_InM_one_Both_mbias4", "sameMIDcase8_mothedPDG_InM_one_Both_mbias4; Id ;#", 5000,
3580 -2500, 2500);
3583 new TH1D("PdgCase8NonEComeFromTarget_mbias4_InM_one_Both_mbias4",
3584 "PdgCase8NonEComeFromTarget_mbias4_InM_one_Both_mbias4; Id ;#", 5000, -2500, 2500);
3587 new TH1D("PdgCase8NonE_NOT_FromTarget_mbias4_InM_one_Both_mbias4",
3588 "PdgCase8NonE_NOT_FromTarget_mbias4_InM_one_Both_mbias4; Id ;#", 5000, -2500, 2500);
3590 PdgCase8motherNonE_InM_one_Both_mbias4 = new TH1D("PdgCase8motherNonE_InM_one_Both_mbias4",
3591 "PdgCase8motherNonE_InM_one_Both_mbias4; Id ;#", 5000, -2500, 2500);
3594 new TH1D("Case8ElFromDalitz_InM_one_Both_mbias4",
3595 "Case8ElFromDalitz_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3598 new TH1D("Case8NonElFrom_pn_InM_one_Both_mbias4",
3599 "Case8NonElFrom_pn_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3602 new TH1D("Case8NonElFrom_eta_InM_one_Both_mbias4",
3603 "Case8NonElFrom_eta_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3606 new TH1D("Case8NonElFrom_kaon_InM_one_Both_mbias4",
3607 "Case8NonElFrom_kaon_InM_one_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3610 "sameMIDcase8NonEPdg_InM_one_Both_mbias4", "sameMIDcase8NonEPdg_InM_one_Both_mbias4; Id ;#", 5000, -2500, 2500);
3613 new TH1D("sameMIDcase8NonEMotherPdg_InM_one_Both_mbias4", "sameMIDcase8NonEMotherPdg_InM_one_Both_mbias4; Id ;#",
3614 5000, -2500, 2500);
3617 new TH1D("sameMIDcase8NonEMotherIM_InM_one_Both_mbias4",
3618 "sameMIDcase8NonEMotherIM_InM_one_Both_mbias4; invariant mass in "
3619 "GeV/c^{2};#",
3620 1000, 0, 2.0);
3623 new TH1D("sameMIDcase8NonEPdgFromTarget_mbias4_InM_one_Both_mbias4",
3624 "sameMIDcase8NonEPdgFromTarget_mbias4_InM_one_Both_mbias4; Id ;#", 5000, -2500, 2500);
3627 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4IM_InM_one_Both_mbias4",
3628 "sameMIDcase8NonEComeFromTarget_mbias4IM_InM_one_Both_mbias4; "
3629 "invariant mass in GeV/c^{2};#",
3630 1000, 0, 2.0);
3633 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4P_InM_one_Both_mbias4",
3634 "sameMIDcase8NonEComeFromTarget_mbias4P_InM_one_Both_mbias4; P in "
3635 "GeV/c^{2} ;#",
3636 200, 0, 10);
3639 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_one_Both_mbias4",
3640 "sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_one_Both_mbias4; Pt "
3641 "in GeV/c^{2} ;#",
3642 200, 0, 10);
3644 //Both_mbias4 two
3646 new TH1D("BG1_InM_two_Both_mbias4", "BG1_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3649 new TH1D("BG2_InM_two_Both_mbias4", "BG2_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3652 new TH1D("BG3_InM_two_Both_mbias4", "BG3_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3655 new TH1D("BG4_InM_two_Both_mbias4", "BG4_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3658 new TH1D("BG5_InM_two_Both_mbias4", "BG5_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3661 new TH1D("BG6_InM_two_Both_mbias4", "BG6_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3664 new TH1D("BG7_InM_two_Both_mbias4", "BG7_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3667 new TH1D("BG8_InM_two_Both_mbias4", "BG8_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3670 new TH1D("BG9_InM_two_Both_mbias4", "BG9_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3673 new TH1D("BG10_InM_two_Both_mbias4", "BG10_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3676 new TH1D("PdgCase8_InM_two_Both_mbias4", "PdgCase8_InM_two_Both_mbias4; Id ;#", 5000, -2500, 2500);
3679 new TH1D("PdgCase8mothers_InM_two_Both_mbias4", "PdgCase8mothers_InM_two_Both_mbias4; Id ;#", 5000, -2500, 2500);
3682 new TH1D("sameMIDcase8_InM_two_Both_mbias4", "sameMIDcase8_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#",
3683 1000, 0, 2.0);
3686 new TH1D("sameGRIDcase8_InM_two_Both_mbias4", "sameGRIDcase8_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#",
3687 1000, 0, 2.0);
3689 Case1ZYPos_InM_two_Both_mbias4 = new TH2D("Case1ZYPos_InM_two_Both_mbias4",
3690 "Case1ZYPos_InM_two_Both_mbias4; z[cm]; y[cm]", 400, -1, 200, 200, -50, 50);
3693 new TH1D("sameMIDcase8_mothedPDG_InM_two_Both_mbias4", "sameMIDcase8_mothedPDG_InM_two_Both_mbias4; Id ;#", 5000,
3694 -2500, 2500);
3697 new TH1D("PdgCase8NonEComeFromTarget_mbias4_InM_two_Both_mbias4",
3698 "PdgCase8NonEComeFromTarget_mbias4_InM_two_Both_mbias4; Id ;#", 5000, -2500, 2500);
3701 new TH1D("PdgCase8NonE_NOT_FromTarget_mbias4_InM_two_Both_mbias4",
3702 "PdgCase8NonE_NOT_FromTarget_mbias4_InM_two_Both_mbias4; Id ;#", 5000, -2500, 2500);
3704 PdgCase8motherNonE_InM_two_Both_mbias4 = new TH1D("PdgCase8motherNonE_InM_two_Both_mbias4",
3705 "PdgCase8motherNonE_InM_two_Both_mbias4; Id ;#", 5000, -2500, 2500);
3708 new TH1D("Case8ElFromDalitz_InM_two_Both_mbias4",
3709 "Case8ElFromDalitz_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3712 new TH1D("Case8NonElFrom_pn_InM_two_Both_mbias4",
3713 "Case8NonElFrom_pn_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3716 new TH1D("Case8NonElFrom_eta_InM_two_Both_mbias4",
3717 "Case8NonElFrom_eta_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3720 new TH1D("Case8NonElFrom_kaon_InM_two_Both_mbias4",
3721 "Case8NonElFrom_kaon_InM_two_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3724 "sameMIDcase8NonEPdg_InM_two_Both_mbias4", "sameMIDcase8NonEPdg_InM_two_Both_mbias4; Id ;#", 5000, -2500, 2500);
3727 new TH1D("sameMIDcase8NonEMotherPdg_InM_two_Both_mbias4", "sameMIDcase8NonEMotherPdg_InM_two_Both_mbias4; Id ;#",
3728 5000, -2500, 2500);
3731 new TH1D("sameMIDcase8NonEMotherIM_InM_two_Both_mbias4",
3732 "sameMIDcase8NonEMotherIM_InM_two_Both_mbias4; invariant mass in "
3733 "GeV/c^{2};#",
3734 1000, 0, 2.0);
3737 new TH1D("sameMIDcase8NonEPdgFromTarget_mbias4_InM_two_Both_mbias4",
3738 "sameMIDcase8NonEPdgFromTarget_mbias4_InM_two_Both_mbias4; Id ;#", 5000, -2500, 2500);
3741 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4IM_InM_two_Both_mbias4",
3742 "sameMIDcase8NonEComeFromTarget_mbias4IM_InM_two_Both_mbias4; "
3743 "invariant mass in GeV/c^{2};#",
3744 1000, 0, 2.0);
3747 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4P_InM_two_Both_mbias4",
3748 "sameMIDcase8NonEComeFromTarget_mbias4P_InM_two_Both_mbias4; P in "
3749 "GeV/c^{2} ;#",
3750 200, 0, 10);
3753 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_two_Both_mbias4",
3754 "sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_two_Both_mbias4; Pt "
3755 "in GeV/c^{2} ;#",
3756 200, 0, 10);
3758 //Both_mbias4 onetwo
3760 new TH1D("BG1_InM_onetwo_Both_mbias4", "BG1_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3763 new TH1D("BG2_InM_onetwo_Both_mbias4", "BG2_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3766 new TH1D("BG3_InM_onetwo_Both_mbias4", "BG3_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3769 new TH1D("BG4_InM_onetwo_Both_mbias4", "BG4_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3772 new TH1D("BG5_InM_onetwo_Both_mbias4", "BG5_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3775 new TH1D("BG6_InM_onetwo_Both_mbias4", "BG6_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3778 new TH1D("BG7_InM_onetwo_Both_mbias4", "BG7_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3781 new TH1D("BG8_InM_onetwo_Both_mbias4", "BG8_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3784 new TH1D("BG9_InM_onetwo_Both_mbias4", "BG9_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3787 new TH1D("BG10_InM_onetwo_Both_mbias4", "BG10_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3790 new TH1D("PdgCase8_InM_onetwo_Both_mbias4", "PdgCase8_InM_onetwo_Both_mbias4; Id ;#", 5000, -2500, 2500);
3792 PdgCase8mothers_InM_onetwo_Both_mbias4 = new TH1D("PdgCase8mothers_InM_onetwo_Both_mbias4",
3793 "PdgCase8mothers_InM_onetwo_Both_mbias4; Id ;#", 5000, -2500, 2500);
3796 new TH1D("sameMIDcase8_InM_onetwo_Both_mbias4",
3797 "sameMIDcase8_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3800 new TH1D("sameGRIDcase8_InM_onetwo_Both_mbias4",
3801 "sameGRIDcase8_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3804 "Case1ZYPos_InM_onetwo_Both_mbias4", "Case1ZYPos_InM_onetwo_Both_mbias4; z[cm]; y[cm]", 400, -1, 200, 200, -50, 50);
3807 new TH1D("sameMIDcase8_mothedPDG_InM_onetwo_Both_mbias4", "sameMIDcase8_mothedPDG_InM_onetwo_Both_mbias4; Id ;#",
3808 5000, -2500, 2500);
3811 new TH1D("PdgCase8NonEComeFromTarget_mbias4_InM_onetwo_Both_mbias4",
3812 "PdgCase8NonEComeFromTarget_mbias4_InM_onetwo_Both_mbias4; Id ;#", 5000, -2500, 2500);
3815 new TH1D("PdgCase8NonE_NOT_FromTarget_mbias4_InM_onetwo_Both_mbias4",
3816 "PdgCase8NonE_NOT_FromTarget_mbias4_InM_onetwo_Both_mbias4; Id ;#", 5000, -2500, 2500);
3819 "PdgCase8motherNonE_InM_onetwo_Both_mbias4", "PdgCase8motherNonE_InM_onetwo_Both_mbias4; Id ;#", 5000, -2500, 2500);
3822 new TH1D("Case8ElFromDalitz_InM_onetwo_Both_mbias4",
3823 "Case8ElFromDalitz_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3826 new TH1D("Case8NonElFrom_pn_InM_onetwo_Both_mbias4",
3827 "Case8NonElFrom_pn_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3830 new TH1D("Case8NonElFrom_eta_InM_onetwo_Both_mbias4",
3831 "Case8NonElFrom_eta_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3834 new TH1D("Case8NonElFrom_kaon_InM_onetwo_Both_mbias4",
3835 "Case8NonElFrom_kaon_InM_onetwo_Both_mbias4; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3838 new TH1D("sameMIDcase8NonEPdg_InM_onetwo_Both_mbias4", "sameMIDcase8NonEPdg_InM_onetwo_Both_mbias4; Id ;#", 5000,
3839 -2500, 2500);
3842 new TH1D("sameMIDcase8NonEMotherPdg_InM_onetwo_Both_mbias4",
3843 "sameMIDcase8NonEMotherPdg_InM_onetwo_Both_mbias4; Id ;#", 5000, -2500, 2500);
3846 new TH1D("sameMIDcase8NonEMotherIM_InM_onetwo_Both_mbias4",
3847 "sameMIDcase8NonEMotherIM_InM_onetwo_Both_mbias4; invariant mass "
3848 "in GeV/c^{2};#",
3849 1000, 0, 2.0);
3852 new TH1D("sameMIDcase8NonEPdgFromTarget_mbias4_InM_onetwo_Both_mbias4",
3853 "sameMIDcase8NonEPdgFromTarget_mbias4_InM_onetwo_Both_mbias4; Id ;#", 5000, -2500, 2500);
3856 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4IM_InM_onetwo_Both_mbias4",
3857 "sameMIDcase8NonEComeFromTarget_mbias4IM_InM_onetwo_Both_mbias4; "
3858 "invariant mass in GeV/c^{2};#",
3859 1000, 0, 2.0);
3862 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4P_InM_onetwo_Both_mbias4",
3863 "sameMIDcase8NonEComeFromTarget_mbias4P_InM_onetwo_Both_mbias4; P "
3864 "in GeV/c^{2} ;#",
3865 200, 0, 10);
3868 new TH1D("sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_onetwo_Both_mbias4",
3869 "sameMIDcase8NonEComeFromTarget_mbias4Pt_InM_onetwo_Both_mbias4; "
3870 "Pt in GeV/c^{2} ;#",
3871 200, 0, 10);
3873}
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
void Exec(CbmMCTrack *mctrack1, CbmMCTrack *mctrack2, CbmMCTrack *mctrack3, CbmMCTrack *mctrack4, Double_t invmassRecoPi0, vector< TH1 * > BGCases)
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_zero_Outside_mbias4
vector< CbmRichRing * > VRichRing_plus_Outside_mbias4
void FindGammasOutside(int EventNumMan, double AngleCut, double InvMassCut, int RealPID, vector< CbmMCTrack * > MCtracks_minus_Outside_mbias4, vector< CbmMCTrack * > MCtracks_plus_Outside_mbias4, vector< CbmStsTrack * > StsTrack_minus_Outside_mbias4, vector< CbmStsTrack * > StsTrack_plus_Outside_mbias4, std::vector< int > Rings_minus_Outside_mbias4, std::vector< int > Rings_plus_Outside_mbias4, std::vector< int > stsIndex_minus_Outside_mbias4, std::vector< int > stsIndex_plus_Outside_mbias4, vector< CbmRichRing * > richRing_minus_Outside_mbias4, vector< CbmRichRing * > richRing_plus_Outside_mbias4)
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Outside_mbias4
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Both_mbias4
std::vector< std::vector< TVector3 > > Gammas_two_Target_mbias4
std::vector< std::vector< TVector3 > > Gammas_all_Both_mbias4
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_zero_Both_mbias4
std::vector< std::vector< TVector3 > > Gammas_one_Both_mbias4
CbmRichRingFitterEllipseTau * fTauFit
void FindGammasTarget(int EventNumMan, 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)
std::vector< std::vector< int > > Gammas_stsIndex_zero_Both_mbias4
std::vector< std::vector< TVector3 > > Gammas_all_Target_mbias4
int FindInRich(int richInd, int stsMcTrackId)
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Target_mbias4
double CalculatePlaneAngle_first(CbmStsTrack *Sts_1, CbmStsTrack *Sts_2)
std::vector< std::vector< int > > Gammas_stsIndex_zero_Target_mbias4
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Both_mbias4
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Both_mbias4
std::vector< std::vector< TVector3 > > Gammas_one_Target_mbias4
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Outside_mbias4
std::vector< std::vector< int > > Gammas_stsIndex_zero_Outside_mbias4
int CheckIfElectron(CbmRichRing *ring, double momentum)
std::vector< std::vector< int > > Gammas_stsIndex_one_Target_mbias4
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_one_Both_mbias4
void FindPi0(TString mod, TString position, vector< vector< TVector3 > > Gammas, vector< vector< int > > StsIndex, vector< vector< CbmMCTrack * > > GammasMC, TH1D *Pi0InvMassReco, TH2D *Pi0_pt_vs_rap, TH2D *Pi0_pt_vs_rap_est, TH2D *MultiplicityGamma, TH2D *MultiplicityChargedParticles, vector< TH1 * > BGCases)
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Both_mbias4
std::vector< std::vector< TVector3 > > EMT_man_pair_momenta_Target_mbias4
std::vector< std::vector< TVector3 > > Gammas_two_Outside_mbias4
vector< CbmRichRing * > VRichRing_plus_Target_mbias4
std::vector< std::vector< TVector3 > > Gammas_onetwo_Target_mbias4
std::vector< std::vector< TVector3 > > Gammas_one_Outside_mbias4
std::vector< std::vector< TVector3 > > EMT_man_pair_momenta_Both_mbias4
void SaveOutsideTracks(CbmMCTrack *mcTrack1, CbmStsTrack *stsTrack, double charge, int stsInd, int richInd, int stsMcTrackId, CbmRichRing *RING)
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Outside_mbias4
vector< CbmStsTrack * > VStsTrack_minus_Outside_mbias4
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Target_mbias4
std::vector< std::vector< TVector3 > > Gammas_two_Both_mbias4
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_zero_Target_mbias4
std::vector< std::vector< int > > Gammas_stsIndex_two_Both_mbias4
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Outside_mbias4
void Exec(int fEventNumMan, double OpeningAngleCut, double GammaInvMassCut, int RealPID)
std::vector< std::vector< int > > Gammas_stsIndex_two_Outside_mbias4
vector< CbmStsTrack * > VStsTrack_plus_Outside_mbias4
std::vector< std::vector< TVector3 > > Gammas_zero_Outside_mbias4
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Target_mbias4
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Target_mbias4
std::vector< std::vector< int > > Gammas_stsIndex_one_Outside_mbias4
std::vector< std::vector< TVector3 > > Gammas_onetwo_Outside_mbias4
std::vector< std::vector< int > > Gammas_stsIndex_two_Target_mbias4
std::vector< std::vector< TVector3 > > Gammas_all_Outside_mbias4
std::vector< std::vector< int > > Gammas_stsIndex_all_Outside_mbias4
std::vector< std::vector< TVector3 > > Gammas_zero_Both_mbias4
std::vector< std::vector< TVector3 > > Gammas_zero_Target_mbias4
vector< CbmStsTrack * > VStsTrack_minus_Target_mbias4
double CalculatePlaneAngle_last(CbmStsTrack *Sts_1, CbmStsTrack *Sts_2)
vector< CbmRichRing * > VRichRing_minus_Target_mbias4
void SaveTargetTracks(CbmMCTrack *mcTrack1, CbmStsTrack *stsTrack, TVector3 refmom, double charge, int stsInd, int richInd, int stsMcTrackId, CbmRichRing *RING)
std::vector< std::vector< TVector3 > > EMT_man_pair_momenta_Outside_mbias4
vector< CbmMCTrack * > VMCtracks_minus_Outside_mbias4
std::vector< std::vector< int > > Gammas_stsIndex_all_Target_mbias4
std::vector< std::vector< int > > Gammas_stsIndex_one_Both_mbias4
vector< CbmRichRing * > VRichRing_minus_Outside_mbias4
std::vector< std::vector< TVector3 > > Gammas_onetwo_Both_mbias4
vector< CbmStsTrack * > VStsTrack_plus_Target_mbias4
std::vector< std::vector< int > > Gammas_stsIndex_all_Both_mbias4
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_one_Target_mbias4
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_one_Outside_mbias4
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 TVector3 FitToVertex(CbmStsTrack *stsTrack, double x, double y, double z)
int32_t GetMotherId() const
Definition CbmMCTrack.h:69
int32_t GetPdgCode() const
Definition CbmMCTrack.h:68
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 fRapidity
Definition LmvmKinePar.h:21
Hash for CbmL1LinkKey.