28 #include "CbmRichUtil.h"
32 #include "FairRootManager.h"
33 #include "KFParticle.h"
41 : fTrainPhotons(nullptr)
43 , fAnnPhotonsSelection(nullptr)
46 , fGlobalTracks(nullptr)
48 , fStsTrackMatches(nullptr)
49 , fRichProjections(nullptr)
51 , fRichRingMatches(nullptr)
53 , fArrayMvdHit(nullptr)
54 , fArrayStsHit(nullptr)
55 , fPrimVertex(nullptr)
58 , VStsTrack_minus_Outside()
59 , VMCtracks_minus_Outside()
60 , VRings_minus_Outside()
61 , VStsIndex_minus_Outside()
62 , VRichRing_minus_Outside()
63 , VStsTrack_plus_Outside()
64 , VMCtracks_plus_Outside()
65 , VRings_plus_Outside()
66 , VStsIndex_plus_Outside()
67 , VRichRing_plus_Outside()
68 , VMCtracks_minus_Target()
69 , VStsTrack_minus_Target()
70 , VMomenta_minus_Target()
71 , VRings_minus_Target()
72 , VStsIndex_minus_Target()
73 , VRichRing_minus_Target()
74 , VMCtracks_plus_Target()
75 , VStsTrack_plus_Target()
76 , VMomenta_plus_Target()
77 , VRings_plus_Target()
78 , VStsIndex_plus_Target()
79 , VRichRing_plus_Target()
83 , DP_LK_EMT_momenta_minus_Target()
84 , DP_LK_EMT_NofRings_minus_Target()
85 , DP_LK_EMT_STS_minus_Target()
86 , DP_LK_EMT_STS_minus_index_Target()
87 , DP_LK_EMT_momenta_plus_Target()
88 , DP_LK_EMT_NofRings_plus_Target()
89 , DP_LK_EMT_STS_plus_Target()
90 , DP_LK_EMT_STS_plus_index_Target()
91 , DP_LK_EMT_NofRings_minus_Outside()
92 , DP_LK_EMT_STS_minus_Outside()
93 , DP_LK_EMT_STS_minus_index_Outside()
94 , DP_LK_EMT_NofRings_plus_Outside()
95 , DP_LK_EMT_STS_plus_Outside()
96 , DP_LK_EMT_STS_plus_index_Outside()
97 , DP_EMT_Event_minus_Target()
98 , DP_EMT_momenta_minus_Target()
99 , DP_EMT_NofRings_minus_Target()
100 , DP_EMT_Hits_minus_Target()
101 , DP_EMT_Event_plus_Target()
102 , DP_EMT_momenta_plus_Target()
103 , DP_EMT_NofRings_plus_Target()
104 , DP_EMT_Hits_plus_Target()
105 , DP_EMT_Event_minus_Outside()
106 , DP_EMT_STS_minus_Outside()
107 , DP_EMT_KFTrack_minus_Outside()
108 , DP_EMT_NofRings_minus_Outside()
109 , DP_EMT_Hits_minus_Outside()
110 , DP_EMT_Event_plus_Outside()
111 , DP_EMT_STS_plus_Outside()
112 , DP_EMT_KFTrack_plus_Outside()
113 , DP_EMT_NofRings_plus_Outside()
114 , DP_EMT_Hits_plus_Outside()
115 , Gammas_all_Target()
116 , Gammas_zero_Target()
117 , Gammas_one_Target()
118 , Gammas_two_Target()
119 , Gammas_onetwo_Target()
120 , Gammas_stsIndex_all_Target()
121 , Gammas_stsIndex_zero_Target()
122 , Gammas_stsIndex_one_Target()
123 , Gammas_stsIndex_two_Target()
124 , Gammas_stsIndex_onetwo_Target()
125 , Gammas_MC_all_Target()
126 , Gammas_MC_zero_Target()
127 , Gammas_MC_one_Target()
128 , Gammas_MC_two_Target()
129 , Gammas_MC_onetwo_Target()
130 , Gammas_all_Outside()
131 , Gammas_zero_Outside()
132 , Gammas_one_Outside()
133 , Gammas_two_Outside()
134 , Gammas_onetwo_Outside()
135 , Gammas_stsIndex_all_Outside()
136 , Gammas_stsIndex_zero_Outside()
137 , Gammas_stsIndex_one_Outside()
138 , Gammas_stsIndex_two_Outside()
139 , Gammas_stsIndex_onetwo_Outside()
140 , Gammas_MC_all_Outside()
141 , Gammas_MC_zero_Outside()
142 , Gammas_MC_one_Outside()
143 , Gammas_MC_two_Outside()
144 , Gammas_MC_onetwo_Outside()
149 , Gammas_onetwo_Both()
150 , Gammas_stsIndex_all_Both()
151 , Gammas_stsIndex_zero_Both()
152 , Gammas_stsIndex_one_Both()
153 , Gammas_stsIndex_two_Both()
154 , Gammas_stsIndex_onetwo_Both()
155 , Gammas_MC_all_Both()
156 , Gammas_MC_zero_Both()
157 , Gammas_MC_one_Both()
158 , Gammas_MC_two_Both()
159 , Gammas_MC_onetwo_Both()
160 , fHistoList_dp_Target()
161 , Mother_PDG_Target(nullptr)
162 , GrandMother_PDG_Target(nullptr)
163 , fHistoList_dp_Outside()
164 , Mother_PDG_Outside(nullptr)
165 , GrandMother_PDG_Outside(nullptr)
166 , fHistoList_dp_Both()
167 , Pdg_vs_Distance_for_dp(nullptr)
168 , P_vs_Distance_for_dp(nullptr)
169 , DP_AnnTruePairs(nullptr)
170 , DP_AnnFalsePairs(nullptr)
171 , DP_AnnTruePairs_AfterCuts(nullptr)
172 , DP_AnnFalsePairs_AfterCuts(nullptr)
173 , fHistoList_dp_cuts_Both()
174 , DP_candidates_InvMass_vs_OA_Both(nullptr)
175 , DP_InvMass_vs_OA_Both(nullptr)
176 , DP_candidates_InvMass_Both(nullptr)
177 , DP_InvMass_Both(nullptr)
178 , DP_candidates_OA_Both(nullptr)
179 , DP_OA_Both(nullptr)
180 , DP_candidates_PlaneAngles_last_Both(nullptr)
181 , DP_PlaneAngles_last_Both(nullptr)
182 , DP_candidates_PlaneAngles_first_Both(nullptr)
183 , DP_PlaneAngles_first_Both(nullptr)
184 , fHistoList_dp_cuts_Target()
185 , DP_candidates_InvMass_vs_OA_Target(nullptr)
186 , DP_InvMass_vs_OA_Target(nullptr)
187 , DP_candidates_InvMass_Target(nullptr)
188 , DP_InvMass_Target(nullptr)
189 , DP_candidates_OA_Target(nullptr)
190 , DP_OA_Target(nullptr)
191 , DP_candidates_PlaneAngles_last_Target(nullptr)
192 , DP_PlaneAngles_last_Target(nullptr)
193 , DP_candidates_PlaneAngles_first_Target(nullptr)
194 , DP_PlaneAngles_first_Target(nullptr)
195 , fHistoList_dp_cuts_Outside()
196 , DP_candidates_InvMass_vs_OA_Outside(nullptr)
197 , DP_InvMass_vs_OA_Outside(nullptr)
198 , DP_candidates_InvMass_Outside(nullptr)
199 , DP_InvMass_Outside(nullptr)
200 , DP_candidates_OA_Outside(nullptr)
201 , DP_OA_Outside(nullptr)
202 , DP_candidates_PlaneAngles_last_Outside(nullptr)
203 , DP_PlaneAngles_last_Outside(nullptr)
204 , DP_candidates_PlaneAngles_first_Outside(nullptr)
205 , DP_PlaneAngles_first_Outside(nullptr)
206 , fHistoList_dp_all_Target()
207 , DP_InvMassReco_all_Target(nullptr)
208 , DP_OpeningAngleReco_all_Target(nullptr)
209 , DP_Pdg_all_Target(nullptr)
210 , DP_P_reco_all_Target(nullptr)
211 , DP_Pt_reco_all_Target(nullptr)
212 , Ph_fromTarget_Pt_reco_all_Target(nullptr)
213 , Ph_fromPions_Pt_reco_all_Target(nullptr)
214 , Ph_fromEtas_Pt_reco_all_Target(nullptr)
215 , Ph_fromDalitz_Pt_reco_all_Target(nullptr)
216 , fHistoList_dp_zero_Target()
217 , DP_InvMassReco_zero_Target(nullptr)
218 , DP_OpeningAngleReco_zero_Target(nullptr)
219 , DP_Pdg_zero_Target(nullptr)
220 , DP_P_reco_zero_Target(nullptr)
221 , DP_Pt_reco_zero_Target(nullptr)
222 , Ph_fromTarget_Pt_reco_zero_Target(nullptr)
223 , Ph_fromPions_Pt_reco_zero_Target(nullptr)
224 , Ph_fromEtas_Pt_reco_zero_Target(nullptr)
225 , Ph_fromDalitz_Pt_reco_zero_Target(nullptr)
226 , fHistoList_dp_one_Target()
227 , DP_InvMassReco_one_Target(nullptr)
228 , DP_OpeningAngleReco_one_Target(nullptr)
229 , DP_Pdg_one_Target(nullptr)
230 , DP_P_reco_one_Target(nullptr)
231 , DP_Pt_reco_one_Target(nullptr)
232 , Ph_fromTarget_Pt_reco_one_Target(nullptr)
233 , Ph_fromPions_Pt_reco_one_Target(nullptr)
234 , Ph_fromEtas_Pt_reco_one_Target(nullptr)
235 , Ph_fromDalitz_Pt_reco_one_Target(nullptr)
236 , fHistoList_dp_two_Target()
237 , DP_InvMassReco_two_Target(nullptr)
238 , DP_OpeningAngleReco_two_Target(nullptr)
239 , DP_Pdg_two_Target(nullptr)
240 , DP_P_reco_two_Target(nullptr)
241 , DP_Pt_reco_two_Target(nullptr)
242 , Ph_fromTarget_Pt_reco_two_Target(nullptr)
243 , Ph_fromPions_Pt_reco_two_Target(nullptr)
244 , Ph_fromEtas_Pt_reco_two_Target(nullptr)
245 , Ph_fromDalitz_Pt_reco_two_Target(nullptr)
246 , Ph_fromXi_Pt_reco_two_Target(nullptr)
247 , Ph_fromOther_Pt_reco_two_Target(nullptr)
248 , Ph_twoFromTarget_Pt_reco_two_Target(nullptr)
249 , Ph_fromCombinatorial_Pt_reco_two_Target(nullptr)
250 , Ph_fromConversion_Pt_reco_two_Target(nullptr)
251 , twoFromTarget_PDG_two_Target(nullptr)
252 , fromCombinatorial_PDG_two_Target(nullptr)
253 , CombinatorialGrMotherPdg_two_Target(nullptr)
254 , CombinatorialMotherPdg_two_Target(nullptr)
255 , Electrons_two_Target(nullptr)
256 , Pions_two_Target(nullptr)
257 , PionElectron_two_Target(nullptr)
258 , elsePionOrElectron_two_Target(nullptr)
259 , DalitzAndConversion_Pt_two_Target(nullptr)
260 , DoubleConversion_Pt_two_Target(nullptr)
261 , fromFireball_P_two_Target(nullptr)
262 , twoFromTarget_P_two_Target(nullptr)
263 , fromCombinatorial_electron_P_two_Target(nullptr)
264 , fromCombinatorial_NOTelectron_P_two_Target(nullptr)
265 , fHistoList_dp_onetwo_Target()
266 , DP_InvMassReco_onetwo_Target(nullptr)
267 , DP_OpeningAngleReco_onetwo_Target(nullptr)
268 , DP_Pdg_onetwo_Target(nullptr)
269 , DP_P_reco_onetwo_Target(nullptr)
270 , DP_Pt_reco_onetwo_Target(nullptr)
271 , Ph_fromTarget_Pt_reco_onetwo_Target(nullptr)
272 , Ph_fromPions_Pt_reco_onetwo_Target(nullptr)
273 , Ph_fromEtas_Pt_reco_onetwo_Target(nullptr)
274 , Ph_fromDalitz_Pt_reco_onetwo_Target(nullptr)
275 , Ph_fromXi_Pt_reco_onetwo_Target(nullptr)
276 , Ph_fromOther_Pt_reco_onetwo_Target(nullptr)
277 , Ph_twoFromTarget_Pt_reco_onetwo_Target(nullptr)
278 , Ph_fromCombinatorial_Pt_reco_onetwo_Target(nullptr)
279 , Ph_fromConversion_Pt_reco_onetwo_Target(nullptr)
280 , twoFromTarget_PDG_onetwo_Target(nullptr)
281 , fromCombinatorial_PDG_onetwo_Target(nullptr)
282 , CombinatorialGrMotherPdg_onetwo_Target(nullptr)
283 , CombinatorialMotherPdg_onetwo_Target(nullptr)
284 , Electrons_onetwo_Target(nullptr)
285 , Pions_onetwo_Target(nullptr)
286 , PionElectron_onetwo_Target(nullptr)
287 , elsePionOrElectron_onetwo_Target(nullptr)
288 , DalitzAndConversion_Pt_onetwo_Target(nullptr)
289 , DoubleConversion_Pt_onetwo_Target(nullptr)
290 , fromFireball_P_onetwo_Target(nullptr)
291 , twoFromTarget_P_onetwo_Target(nullptr)
292 , fromCombinatorial_electron_P_onetwo_Target(nullptr)
293 , fromCombinatorial_NOTelectron_P_onetwo_Target(nullptr)
294 , fHistoList_dp_all_Outside()
295 , DP_InvMassReco_all_Outside(nullptr)
296 , DP_OpeningAngleReco_all_Outside(nullptr)
297 , DP_Pdg_all_Outside(nullptr)
298 , DP_P_reco_all_Outside(nullptr)
299 , DP_Pt_reco_all_Outside(nullptr)
300 , Ph_fromTarget_Pt_reco_all_Outside(nullptr)
301 , Ph_fromPions_Pt_reco_all_Outside(nullptr)
302 , Ph_fromEtas_Pt_reco_all_Outside(nullptr)
303 , Ph_fromDalitz_Pt_reco_all_Outside(nullptr)
304 , fHistoList_dp_zero_Outside()
305 , DP_InvMassReco_zero_Outside(nullptr)
306 , DP_OpeningAngleReco_zero_Outside(nullptr)
307 , DP_Pdg_zero_Outside(nullptr)
308 , DP_P_reco_zero_Outside(nullptr)
309 , DP_Pt_reco_zero_Outside(nullptr)
310 , Ph_fromTarget_Pt_reco_zero_Outside(nullptr)
311 , Ph_fromPions_Pt_reco_zero_Outside(nullptr)
312 , Ph_fromEtas_Pt_reco_zero_Outside(nullptr)
313 , Ph_fromDalitz_Pt_reco_zero_Outside(nullptr)
314 , fHistoList_dp_one_Outside()
315 , DP_InvMassReco_one_Outside(nullptr)
316 , DP_OpeningAngleReco_one_Outside(nullptr)
317 , DP_Pdg_one_Outside(nullptr)
318 , DP_P_reco_one_Outside(nullptr)
319 , DP_Pt_reco_one_Outside(nullptr)
320 , Ph_fromTarget_Pt_reco_one_Outside(nullptr)
321 , Ph_fromPions_Pt_reco_one_Outside(nullptr)
322 , Ph_fromEtas_Pt_reco_one_Outside(nullptr)
323 , Ph_fromDalitz_Pt_reco_one_Outside(nullptr)
324 , fHistoList_dp_two_Outside()
325 , DP_InvMassReco_two_Outside(nullptr)
326 , DP_OpeningAngleReco_two_Outside(nullptr)
327 , DP_Pdg_two_Outside(nullptr)
328 , DP_P_reco_two_Outside(nullptr)
329 , DP_Pt_reco_two_Outside(nullptr)
330 , Ph_fromTarget_Pt_reco_two_Outside(nullptr)
331 , Ph_fromPions_Pt_reco_two_Outside(nullptr)
332 , Ph_fromEtas_Pt_reco_two_Outside(nullptr)
333 , Ph_fromDalitz_Pt_reco_two_Outside(nullptr)
334 , Ph_fromXi_Pt_reco_two_Outside(nullptr)
335 , Ph_fromOther_Pt_reco_two_Outside(nullptr)
336 , Ph_twoFromTarget_Pt_reco_two_Outside(nullptr)
337 , Ph_fromCombinatorial_Pt_reco_two_Outside(nullptr)
338 , Ph_fromConversion_Pt_reco_two_Outside(nullptr)
339 , twoFromTarget_PDG_two_Outside(nullptr)
340 , fromCombinatorial_PDG_two_Outside(nullptr)
341 , CombinatorialGrMotherPdg_two_Outside(nullptr)
342 , CombinatorialMotherPdg_two_Outside(nullptr)
343 , Electrons_two_Outside(nullptr)
344 , Pions_two_Outside(nullptr)
345 , PionElectron_two_Outside(nullptr)
346 , elsePionOrElectron_two_Outside(nullptr)
347 , DalitzAndConversion_Pt_two_Outside(nullptr)
348 , DoubleConversion_Pt_two_Outside(nullptr)
349 , fromFireball_P_two_Outside(nullptr)
350 , twoFromTarget_P_two_Outside(nullptr)
351 , fromCombinatorial_electron_P_two_Outside(nullptr)
352 , fromCombinatorial_NOTelectron_P_two_Outside(nullptr)
353 , fHistoList_dp_onetwo_Outside()
354 , DP_InvMassReco_onetwo_Outside(nullptr)
355 , DP_OpeningAngleReco_onetwo_Outside(nullptr)
356 , DP_Pdg_onetwo_Outside(nullptr)
357 , DP_P_reco_onetwo_Outside(nullptr)
358 , DP_Pt_reco_onetwo_Outside(nullptr)
359 , Ph_fromTarget_Pt_reco_onetwo_Outside(nullptr)
360 , Ph_fromPions_Pt_reco_onetwo_Outside(nullptr)
361 , Ph_fromEtas_Pt_reco_onetwo_Outside(nullptr)
362 , Ph_fromDalitz_Pt_reco_onetwo_Outside(nullptr)
363 , Ph_fromXi_Pt_reco_onetwo_Outside(nullptr)
364 , Ph_fromOther_Pt_reco_onetwo_Outside(nullptr)
365 , Ph_twoFromTarget_Pt_reco_onetwo_Outside(nullptr)
366 , Ph_fromCombinatorial_Pt_reco_onetwo_Outside(nullptr)
367 , Ph_fromConversion_Pt_reco_onetwo_Outside(nullptr)
368 , twoFromTarget_PDG_onetwo_Outside(nullptr)
369 , fromCombinatorial_PDG_onetwo_Outside(nullptr)
370 , CombinatorialGrMotherPdg_onetwo_Outside(nullptr)
371 , CombinatorialMotherPdg_onetwo_Outside(nullptr)
372 , Electrons_onetwo_Outside(nullptr)
373 , Pions_onetwo_Outside(nullptr)
374 , PionElectron_onetwo_Outside(nullptr)
375 , elsePionOrElectron_onetwo_Outside(nullptr)
376 , DalitzAndConversion_Pt_onetwo_Outside(nullptr)
377 , DoubleConversion_Pt_onetwo_Outside(nullptr)
378 , fromFireball_P_onetwo_Outside(nullptr)
379 , twoFromTarget_P_onetwo_Outside(nullptr)
380 , fromCombinatorial_electron_P_onetwo_Outside(nullptr)
381 , fromCombinatorial_NOTelectron_P_onetwo_Outside(nullptr)
382 , fHistoList_dp_all_Both()
383 , DP_InvMassReco_all_Both(nullptr)
384 , DP_OpeningAngleReco_all_Both(nullptr)
385 , DP_Pdg_all_Both(nullptr)
386 , DP_P_reco_all_Both(nullptr)
387 , DP_Pt_reco_all_Both(nullptr)
388 , Ph_fromTarget_Pt_reco_all_Both(nullptr)
389 , Ph_fromPions_Pt_reco_all_Both(nullptr)
390 , Ph_fromEtas_Pt_reco_all_Both(nullptr)
391 , Ph_fromDalitz_Pt_reco_all_Both(nullptr)
392 , Ph_pt_vs_rap_est_all_Both(nullptr)
393 , fHistoList_dp_zero_Both()
394 , DP_InvMassReco_zero_Both(nullptr)
395 , DP_OpeningAngleReco_zero_Both(nullptr)
396 , DP_Pdg_zero_Both(nullptr)
397 , DP_P_reco_zero_Both(nullptr)
398 , DP_Pt_reco_zero_Both(nullptr)
399 , Ph_fromTarget_Pt_reco_zero_Both(nullptr)
400 , Ph_fromPions_Pt_reco_zero_Both(nullptr)
401 , Ph_fromEtas_Pt_reco_zero_Both(nullptr)
402 , Ph_fromDalitz_Pt_reco_zero_Both(nullptr)
403 , Ph_pt_vs_rap_est_zero_Both(nullptr)
404 , fHistoList_dp_one_Both()
405 , DP_InvMassReco_one_Both(nullptr)
406 , DP_OpeningAngleReco_one_Both(nullptr)
407 , DP_Pdg_one_Both(nullptr)
408 , DP_P_reco_one_Both(nullptr)
409 , DP_Pt_reco_one_Both(nullptr)
410 , Ph_fromTarget_Pt_reco_one_Both(nullptr)
411 , Ph_fromPions_Pt_reco_one_Both(nullptr)
412 , Ph_fromEtas_Pt_reco_one_Both(nullptr)
413 , Ph_fromDalitz_Pt_reco_one_Both(nullptr)
414 , Ph_pt_vs_rap_est_one_Both(nullptr)
415 , fHistoList_dp_two_Both()
416 , DP_InvMassReco_two_Both(nullptr)
417 , DP_OpeningAngleReco_two_Both(nullptr)
418 , DP_Pdg_two_Both(nullptr)
419 , DP_P_reco_two_Both(nullptr)
420 , DP_Pt_reco_two_Both(nullptr)
421 , Ph_fromTarget_Pt_reco_two_Both(nullptr)
422 , Ph_fromPions_Pt_reco_two_Both(nullptr)
423 , Ph_fromEtas_Pt_reco_two_Both(nullptr)
424 , Ph_fromDalitz_Pt_reco_two_Both(nullptr)
425 , Ph_twoFromTarget_Pt_reco_two_Both(nullptr)
426 , Ph_fromCombinatorial_Pt_reco_two_Both(nullptr)
427 , Ph_fromConversion_Pt_reco_two_Both(nullptr)
428 , Ph_pt_vs_rap_est_two_Both(nullptr)
429 , fHistoList_dp_onetwo_Both()
430 , DP_InvMassReco_onetwo_Both(nullptr)
431 , DP_OpeningAngleReco_onetwo_Both(nullptr)
432 , DP_Pdg_onetwo_Both(nullptr)
433 , DP_P_reco_onetwo_Both(nullptr)
434 , DP_Pt_reco_onetwo_Both(nullptr)
435 , Ph_fromTarget_Pt_reco_onetwo_Both(nullptr)
436 , Ph_fromPions_Pt_reco_onetwo_Both(nullptr)
437 , Ph_fromEtas_Pt_reco_onetwo_Both(nullptr)
438 , Ph_fromDalitz_Pt_reco_onetwo_Both(nullptr)
439 , Ph_twoFromTarget_Pt_reco_onetwo_Both(nullptr)
440 , Ph_fromCombinatorial_Pt_reco_onetwo_Both(nullptr)
441 , Ph_fromConversion_Pt_reco_onetwo_Both(nullptr)
442 , Ph_pt_vs_rap_est_onetwo_Both(nullptr)
443 , DP_EMT_Pt_all_Target(nullptr)
444 , DP_EMT_Pt_zero_Target(nullptr)
445 , DP_EMT_Pt_one_Target(nullptr)
446 , DP_EMT_Pt_two_Target(nullptr)
447 , DP_EMT_Pt_onetwo_Target(nullptr)
448 , DP_EMT_Pt_all_Outside(nullptr)
449 , DP_EMT_Pt_zero_Outside(nullptr)
450 , DP_EMT_Pt_one_Outside(nullptr)
451 , DP_EMT_Pt_two_Outside(nullptr)
452 , DP_EMT_Pt_onetwo_Outside(nullptr)
453 , DP_EMT_Pt_all_Both(nullptr)
454 , DP_EMT_Pt_zero_Both(nullptr)
455 , DP_EMT_Pt_one_Both(nullptr)
456 , DP_EMT_Pt_two_Both(nullptr)
457 , DP_EMT_Pt_onetwo_Both(nullptr)
458 , DP_LK_EMT_neg_Pt_all_Target(nullptr)
459 , DP_LK_EMT_neg_Pt_zero_Target(nullptr)
460 , DP_LK_EMT_neg_Pt_one_Target(nullptr)
461 , DP_LK_EMT_neg_Pt_two_Target(nullptr)
462 , DP_LK_EMT_neg_Pt_onetwo_Target(nullptr)
463 , DP_LK_EMT_pos_Pt_all_Target(nullptr)
464 , DP_LK_EMT_pos_Pt_zero_Target(nullptr)
465 , DP_LK_EMT_pos_Pt_one_Target(nullptr)
466 , DP_LK_EMT_pos_Pt_two_Target(nullptr)
467 , DP_LK_EMT_pos_Pt_onetwo_Target(nullptr)
468 , DP_LK_EMT_Pt_all_Target(nullptr)
469 , DP_LK_EMT_Pt_zero_Target(nullptr)
470 , DP_LK_EMT_Pt_one_Target(nullptr)
471 , DP_LK_EMT_Pt_two_Target(nullptr)
472 , DP_LK_EMT_Pt_onetwo_Target(nullptr)
473 , DP_LK_EMT_neg_Pt_all_Outside(nullptr)
474 , DP_LK_EMT_neg_Pt_zero_Outside(nullptr)
475 , DP_LK_EMT_neg_Pt_one_Outside(nullptr)
476 , DP_LK_EMT_neg_Pt_two_Outside(nullptr)
477 , DP_LK_EMT_neg_Pt_onetwo_Outside(nullptr)
478 , DP_LK_EMT_pos_Pt_all_Outside(nullptr)
479 , DP_LK_EMT_pos_Pt_zero_Outside(nullptr)
480 , DP_LK_EMT_pos_Pt_one_Outside(nullptr)
481 , DP_LK_EMT_pos_Pt_two_Outside(nullptr)
482 , DP_LK_EMT_pos_Pt_onetwo_Outside(nullptr)
483 , DP_LK_EMT_Pt_all_Outside(nullptr)
484 , DP_LK_EMT_Pt_zero_Outside(nullptr)
485 , DP_LK_EMT_Pt_one_Outside(nullptr)
486 , DP_LK_EMT_Pt_two_Outside(nullptr)
487 , DP_LK_EMT_Pt_onetwo_Outside(nullptr) {}
492 FairRootManager* ioman = FairRootManager::Instance();
493 if (
nullptr == ioman) {
494 Fatal(
"CbmKresConversionPhotons::Init",
"RootManager not instantised!");
497 fMcTracks = (TClonesArray*) ioman->GetObject(
"MCTrack");
499 Fatal(
"CbmKresConversionPhotons::Init",
"No MCTrack array!");
507 LOG(fatal) <<
"CbmKresConversionPhotons::Init No PrimaryVertex array!";
510 fGlobalTracks = (TClonesArray*) ioman->GetObject(
"GlobalTrack");
512 Fatal(
"CbmKresConversionPhotons::Init",
"No GlobalTrack array!");
515 fStsTracks = (TClonesArray*) ioman->GetObject(
"StsTrack");
517 Fatal(
"CbmKresConversionPhotons::Init",
"No StsTrack array!");
522 Fatal(
"CbmKresConversionPhotons::Init",
"No StsTrackMatch array!");
527 Fatal(
"CbmKresConversionPhotons::Init",
"No RichProjection array!");
530 fRichRings = (TClonesArray*) ioman->GetObject(
"RichRing");
532 Fatal(
"CbmKresConversionPhotons::Init",
"No RichRing array!");
537 Fatal(
"CbmKresConversionPhotons::Init",
"No RichRingMatch array!");
540 fRichHits = (TClonesArray*) ioman->GetObject(
"RichHit");
542 Fatal(
"CbmKresConversionPhotons::Init",
"No RichHit array!");
545 fArrayMvdHit = (TClonesArray*) ioman->GetObject(
"MvdHit");
547 Fatal(
"CbmKresConversionPhotons::Init",
"No MvdHit array!");
550 fArrayStsHit = (TClonesArray*) ioman->GetObject(
"StsHit");
552 Fatal(
"CbmKresConversionPhotons::Init",
"No StsHit array!");
568 cout <<
"\t *** CbmKresTrainAnnDirectPhotons ==> is "
578 cout <<
"\t *** CbmKresSelectAnnPhotons ==> is "
586 double OpeningAngleCut,
587 double GammaInvMassCut,
589 cout <<
"CbmKresConversionPhotons, event No. " << fEventNumDP << endl;
594 Fatal(
"CbmKresConversionPhotons::Exec",
"No PrimaryVertex array!");
671 for (Int_t
i = 0;
i < ngTracks;
i++) {
673 if (
nullptr == gTrack)
continue;
677 if (stsInd < 0)
continue;
679 if (stsTrack ==
nullptr)
continue;
682 if (stsMatch ==
nullptr)
continue;
685 if (stsMcTrackId < 0)
continue;
687 if (mcTrack ==
nullptr)
continue;
691 if (richInd < 0 && proj->GetX() > -115 && proj->GetX() < 115
692 && ((proj->GetY() < -120 && proj->GetY() > -200)
693 || (proj->GetY() > 120 && proj->GetY() < 200)))
714 double charge = track_par->GetQp();
717 if (chi2 != chi2)
continue;
718 if (chi2 == 0)
continue;
722 mcTrack, stsTrack, charge, stsInd, richInd, stsMcTrackId, Ring);
724 if (chi2 > 3)
continue;
727 mcTrack, stsTrack, Momentum, charge, stsInd, richInd, stsMcTrackId, Ring);
767 int numformix = 1000;
768 if (fEventNumDP % numformix == 0) {
781 if (fEventNumDP % numformix == 0) {
826 int InRich =
FindInRich(richInd, stsMcTrackId);
851 int InRich =
FindInRich(richInd, stsMcTrackId);
876 vector<CbmMCTrack*> MCtracks_minus,
877 vector<CbmMCTrack*> MCtracks_plus,
878 vector<CbmStsTrack*> StsTrack_minus,
879 vector<CbmStsTrack*> StsTrack_plus,
880 vector<TVector3> Momenta_minus,
881 vector<TVector3> Momenta_plus,
882 std::vector<int> Rings_minus,
883 std::vector<int> Rings_plus,
884 std::vector<int> stsIndex_minus,
885 std::vector<int> stsIndex_plus,
886 vector<CbmRichRing*> richRing_minus,
887 vector<CbmRichRing*> richRing_plus) {
888 for (
size_t i = 0;
i < Momenta_minus.size();
i++) {
889 for (
size_t j = 0; j < Momenta_plus.size(); j++) {
894 TVector3 part1 = Momenta_minus[
i];
895 TVector3 part2 = Momenta_plus[j];
898 int sts1_index = stsIndex_minus[
i];
899 int sts2_index = stsIndex_plus[j];
907 richcheck_1 =
CheckIfElectron(richRing_plus[j], Momenta_plus[j].Mag());
910 richcheck_0 = Rings_minus[
i];
911 richcheck_1 = Rings_plus[j];
913 int richcheck = richcheck_0 + richcheck_1;
915 Double_t InvmassReco =
917 Double_t OpeningAngle =
967 double AnnValue = 999;
987 if (TMath::Abs(OpeningAngle) > AngleCut)
continue;
988 if (TMath::Abs(InvmassReco) > InvMassCut)
continue;
993 int fromFireball = 0;
997 int fromCombinatorial = 0;
998 int fromConversion = 0;
999 int twoFromTarget = 0;
1010 if (mcTrackGrmama->
GetPdgCode() == 111) fromPions = 1;
1011 if (mcTrackGrmama->
GetPdgCode() == 221) fromEta = 1;
1012 if (mcTrackGrmama->
GetPdgCode() == 3212) fromXi = 1;
1016 if (mcTrackmama->
GetPdgCode() == 111) fromDalitz = 1;
1017 if (mcTrackmama->
GetPdgCode() != 111) fromOther = 1;
1023 fromCombinatorial = 1;
1052 frefId.push_back(stsIndex_minus[
i]);
1053 frefId.push_back(stsIndex_plus[j]);
1060 std::vector<TVector3> minusTrack =
SaveAllHits(sts1);
1061 std::vector<TVector3> plusTrack =
SaveAllHits(sts2);
1086 if (richcheck == 0 || richcheck == 1 || richcheck == 2) {
1102 if (fromFireball == 1) {
1107 if (fromPions == 1) {
1115 if (fromDalitz == 1) {
1122 if (richcheck == 0) {
1138 if (fromFireball == 1) {
1143 if (fromPions == 1) {
1151 if (fromDalitz == 1) {
1158 if (richcheck == 1) {
1174 if (fromFireball == 1) {
1179 if (fromPions == 1) {
1187 if (fromDalitz == 1) {
1194 if (richcheck == 2) {
1210 if (fromFireball == 1) {
1217 if (fromPions == 1) {
1225 if (fromDalitz == 1) {
1230 if (fromOther == 1) {
1233 if (fromCombinatorial == 1) {
1273 && TMath::Abs(part2MC->
GetPdgCode()) == 11) {
1275 }
else if (TMath::Abs(part1MC->
GetPdgCode()) == 211
1276 && TMath::Abs(part2MC->
GetPdgCode()) == 211) {
1278 }
else if ((TMath::Abs(part1MC->
GetPdgCode()) == 211
1281 && TMath::Abs(part2MC->
GetPdgCode()) == 211)) {
1289 && TMath::Abs(part2MC->
GetPdgCode()) == 11) {
1297 (mama1PDG == 22 && mama2PDG == 111
1299 || (mama1PDG == 111 && mama2PDG == 22
1306 mama1PDG == 22 && mama2PDG == 22
1321 if (TMath::Abs(part1MC->
GetPdgCode() == 11)) {
1326 if (TMath::Abs(part2MC->
GetPdgCode() == 11)) {
1332 if (fromConversion == 1) {
1336 if (twoFromTarget == 1) {
1347 if (richcheck == 1 || richcheck == 2) {
1363 if (fromFireball == 1) {
1370 if (fromPions == 1) {
1378 if (fromDalitz == 1) {
1383 if (fromOther == 1) {
1386 if (fromCombinatorial == 1) {
1426 && TMath::Abs(part2MC->
GetPdgCode()) == 11) {
1428 }
else if (TMath::Abs(part1MC->
GetPdgCode()) == 211
1429 && TMath::Abs(part2MC->
GetPdgCode()) == 211) {
1431 }
else if ((TMath::Abs(part1MC->
GetPdgCode()) == 211
1434 && TMath::Abs(part2MC->
GetPdgCode()) == 211)) {
1442 && TMath::Abs(part2MC->
GetPdgCode()) == 11) {
1450 (mama1PDG == 22 && mama2PDG == 111
1452 || (mama1PDG == 111 && mama2PDG == 22
1459 mama1PDG == 22 && mama2PDG == 22
1474 if (TMath::Abs(part1MC->
GetPdgCode() == 11)) {
1479 if (TMath::Abs(part2MC->
GetPdgCode() == 11)) {
1485 if (fromConversion == 1) {
1489 if (twoFromTarget == 1) {
1510 vector<CbmMCTrack*> MCtracks_minus_Outside,
1511 vector<CbmMCTrack*> MCtracks_plus_Outside,
1512 vector<CbmStsTrack*> StsTrack_minus_Outside,
1513 vector<CbmStsTrack*> StsTrack_plus_Outside,
1514 std::vector<int> Rings_minus_Outside,
1515 std::vector<int> Rings_plus_Outside,
1516 std::vector<int> stsIndex_minus_Outside,
1517 std::vector<int> stsIndex_plus_Outside,
1518 vector<CbmRichRing*> richRing_minus_Outside,
1519 vector<CbmRichRing*> richRing_plus_Outside) {
1520 for (
size_t i = 0;
i < StsTrack_minus_Outside.size();
i++) {
1521 for (
size_t j = 0; j < StsTrack_plus_Outside.size(); j++) {
1526 CbmMCTrack* part2MC = MCtracks_plus_Outside[j];
1527 int sts1_index = stsIndex_minus_Outside[
i];
1530 KFParticle electron;
1532 part1STS, &electron, 11);
1533 KFParticle positron;
1535 part2STS, &positron, -11);
1536 const KFParticle* daughters[2] = {&electron, &positron};
1537 KFParticle intersection;
1538 intersection.Construct(daughters, 2);
1540 if (intersection.GetZ() > 75 || intersection.GetZ() < -5)
1545 intersection.GetX(),
1546 intersection.GetY(),
1547 intersection.GetZ());
1549 intersection.GetX(),
1550 intersection.GetY(),
1551 intersection.GetZ());
1557 int richcheck_0 = 0;
1558 int richcheck_1 = 0;
1563 richcheck = richcheck_0 + richcheck_1;
1566 richcheck = Rings_minus_Outside[
i] + Rings_plus_Outside[j];
1570 Double_t InvmassReco =
1572 Double_t OpeningAngle =
1616 intersection.GetZ(),
1622 double AnnValue = 999;
1627 intersection.GetZ(),
1640 if (TMath::Abs(OpeningAngle) > AngleCut)
continue;
1641 if (TMath::Abs(InvmassReco) > InvMassCut)
continue;
1647 int fromFireball = 0;
1651 int fromCombinatorial = 0;
1652 int fromConversion = 0;
1653 int twoFromTarget = 0;
1664 if (mcTrackGrmama->
GetPdgCode() == 111) fromPions = 1;
1665 if (mcTrackGrmama->
GetPdgCode() == 221) fromEta = 1;
1666 if (mcTrackGrmama->
GetPdgCode() == 3212) fromXi = 1;
1670 if (mcTrackmama->
GetPdgCode() == 111) fromDalitz = 1;
1671 if (mcTrackmama->
GetPdgCode() != 111) fromOther = 1;
1677 fromCombinatorial = 1;
1707 frefId.push_back(stsIndex_minus_Outside[
i]);
1708 frefId.push_back(stsIndex_plus_Outside[j]);
1715 std::vector<TVector3> minusTrack =
SaveAllHits(part1STS);
1716 std::vector<TVector3> plusTrack =
SaveAllHits(part2STS);
1742 if (richcheck == 0 || richcheck == 1 || richcheck == 2) {
1758 if (fromFireball == 1) {
1763 if (fromPions == 1) {
1771 if (fromDalitz == 1) {
1778 if (richcheck == 0) {
1794 if (fromFireball == 1) {
1799 if (fromPions == 1) {
1807 if (fromDalitz == 1) {
1814 if (richcheck == 1) {
1830 if (fromFireball == 1) {
1835 if (fromPions == 1) {
1843 if (fromDalitz == 1) {
1850 if (richcheck == 2) {
1866 if (fromFireball == 1) {
1873 if (fromPions == 1) {
1881 if (fromDalitz == 1) {
1886 if (fromOther == 1) {
1889 if (fromCombinatorial == 1) {
1929 && TMath::Abs(part2MC->
GetPdgCode()) == 11) {
1931 }
else if (TMath::Abs(part1MC->
GetPdgCode()) == 211
1932 && TMath::Abs(part2MC->
GetPdgCode()) == 211) {
1934 }
else if ((TMath::Abs(part1MC->
GetPdgCode()) == 211
1937 && TMath::Abs(part2MC->
GetPdgCode()) == 211)) {
1945 && TMath::Abs(part2MC->
GetPdgCode()) == 11) {
1953 (mama1PDG == 22 && mama2PDG == 111
1955 || (mama1PDG == 111 && mama2PDG == 22
1962 mama1PDG == 22 && mama2PDG == 22
1977 if (TMath::Abs(part1MC->
GetPdgCode() == 11)) {
1982 if (TMath::Abs(part2MC->
GetPdgCode() == 11)) {
1988 if (fromConversion == 1) {
1992 if (twoFromTarget == 1) {
2005 if (richcheck == 1 || richcheck == 2) {
2021 if (fromFireball == 1) {
2028 if (fromPions == 1) {
2036 if (fromDalitz == 1) {
2041 if (fromOther == 1) {
2044 if (fromCombinatorial == 1) {
2084 && TMath::Abs(part2MC->
GetPdgCode()) == 11) {
2086 }
else if (TMath::Abs(part1MC->
GetPdgCode()) == 211
2087 && TMath::Abs(part2MC->
GetPdgCode()) == 211) {
2089 }
else if ((TMath::Abs(part1MC->
GetPdgCode()) == 211
2092 && TMath::Abs(part2MC->
GetPdgCode()) == 211)) {
2100 && TMath::Abs(part2MC->
GetPdgCode()) == 11) {
2108 (mama1PDG == 22 && mama2PDG == 111
2110 || (mama1PDG == 111 && mama2PDG == 22
2117 mama1PDG == 22 && mama2PDG == 22
2132 if (TMath::Abs(part1MC->
GetPdgCode() == 11)) {
2137 if (TMath::Abs(part2MC->
GetPdgCode() == 11)) {
2143 if (fromConversion == 1) {
2147 if (twoFromTarget == 1) {
2256 int RingsInRich = 0;
2260 if (richMatch !=
nullptr && richMatch->
GetNofLinks() > 0) {
2262 if (richMcTrackId > 0) {
2267 if (mcTrack2 !=
nullptr) {
2269 if (TMath::Abs(pdgRICH) == 11) RingsInRich++;
2283 if (
nullptr != ring) {
2286 for (
int i = 0;
i < nofHits;
i++) {
2289 if (
nullptr == hit)
continue;
2304 std::vector<TVector3>
2306 std::vector<TVector3> AllHitsOfTrack;
2307 AllHitsOfTrack.clear();
2312 for (
int i = 0;
i < hitsMVD;
i++) {
2317 AllHitsOfTrack.push_back(position);
2319 for (
int i = 0;
i < hitsSTS;
i++) {
2324 AllHitsOfTrack.push_back(position);
2327 return AllHitsOfTrack;
2331 std::vector<TVector3> track_1,
2332 std::vector<TVector3> track_2) {
2333 double FinalAngle = 400;
2335 int hits_1 = track_1.size();
2336 int hits_2 = track_2.size();
2345 TVector3 last1 = track_1[hits_1 - 1];
2346 TVector3 last2 = track_2[hits_2 - 1];
2350 TMath::Abs(last1.Z() - last2.Z()) > 2
2354 for (
int i = hits_1 - 2;
i > -1;
i--) {
2355 if (TMath::Abs(last1.Z() - last2.Z()) < 2)
continue;
2364 TMath::Abs(last1.Z() - last2.Z()) > 2
2368 for (
int i = hits_2 - 2;
i > -1;
i--) {
2369 if (TMath::Abs(last1.Z() - last2.Z()) < 2)
continue;
2378 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
2380 TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 /
TMath::Pi());
2383 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
2389 double FinalAngle = 400;
2404 Xpart1 = stsHit1->
GetX();
2405 Ypart1 = stsHit1->
GetY();
2406 Zpart1 = stsHit1->
GetZ();
2410 Xpart1 = mvdHit1->
GetX();
2411 Ypart1 = mvdHit1->
GetY();
2412 Zpart1 = mvdHit1->
GetZ();
2418 Xpart2 = stsHit2->
GetX();
2419 Ypart2 = stsHit2->
GetY();
2420 Zpart2 = stsHit2->
GetZ();
2424 Xpart2 = mvdHit2->
GetX();
2425 Ypart2 = mvdHit2->
GetY();
2426 Zpart2 = mvdHit2->
GetZ();
2431 TMath::Abs(Zpart1 - Zpart2) > 2
2434 for (
int i = hits1sts - 2;
i > -1;
i--) {
2435 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
2438 Xpart1 = stsHit->
GetX();
2439 Ypart1 = stsHit->
GetY();
2440 Zpart1 = stsHit->
GetZ();
2442 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 > Zpart2) {
2443 for (
int i = hits1mvd - 2;
i > -1;
i--) {
2444 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
2447 Xpart1 = mvdHit->
GetX();
2448 Ypart1 = mvdHit->
GetY();
2449 Zpart1 = mvdHit->
GetZ();
2455 TMath::Abs(Zpart1 - Zpart2) > 2
2458 for (
int i = hits2sts - 2;
i > -1;
i--) {
2459 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
2462 Xpart2 = stsHit->
GetX();
2463 Ypart2 = stsHit->
GetY();
2464 Zpart2 = stsHit->
GetZ();
2466 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 < Zpart2) {
2467 for (
int i = hits2mvd - 2;
i > -1;
i--) {
2468 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
2471 Xpart2 = mvdHit->
GetX();
2472 Ypart2 = mvdHit->
GetY();
2473 Zpart2 = mvdHit->
GetZ();
2479 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
2481 TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 /
TMath::Pi());
2484 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
2489 double FinalAngle = 400;
2504 Xpart1 = mvdHit1->
GetX();
2505 Ypart1 = mvdHit1->
GetY();
2506 Zpart1 = mvdHit1->
GetZ();
2510 Xpart1 = stsHit1->
GetX();
2511 Ypart1 = stsHit1->
GetY();
2512 Zpart1 = stsHit1->
GetZ();
2518 Xpart2 = mvdHit2->
GetX();
2519 Ypart2 = mvdHit2->
GetY();
2520 Zpart2 = mvdHit2->
GetZ();
2524 Xpart2 = stsHit2->
GetX();
2525 Ypart2 = stsHit2->
GetY();
2526 Zpart2 = stsHit2->
GetZ();
2531 TMath::Abs(Zpart1 - Zpart2) > 2
2534 for (
int i = 1;
i < hits1mvd;
i++) {
2535 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
2538 Xpart1 = mvdHit->
GetX();
2539 Ypart1 = mvdHit->
GetY();
2540 Zpart1 = mvdHit->
GetZ();
2542 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 < Zpart2) {
2543 for (
int i = 0;
i < hits1sts;
i++) {
2544 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
2547 Xpart1 = stsHit->
GetX();
2548 Ypart1 = stsHit->
GetY();
2549 Zpart1 = stsHit->
GetZ();
2555 TMath::Abs(Zpart1 - Zpart2) > 2
2558 for (
int i = 1;
i < hits2mvd;
i++) {
2559 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
2562 Xpart2 = mvdHit->
GetX();
2563 Ypart2 = mvdHit->
GetY();
2564 Zpart2 = mvdHit->
GetZ();
2566 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 > Zpart2) {
2567 for (
int i = 0;
i < hits2sts;
i++) {
2568 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
2571 Xpart2 = stsHit->
GetX();
2572 Ypart2 = stsHit->
GetY();
2573 Zpart2 = stsHit->
GetZ();
2579 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
2581 TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 /
TMath::Pi());
2584 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
2592 for (
int a = 0; a < nof_minus - 1; a++) {
2593 for (
int b = a + 1; b < nof_minus; b++) {
2603 double OpeningAngle =
2616 if (TMath::Abs(OpeningAngle) > AngleCut)
continue;
2617 if (TMath::Abs(InvmassReco) > InvMassCut)
continue;
2630 if ((rings_amount == 1 || rings_amount == 2)) {
2637 if ((rings_amount == 1 || rings_amount == 2))
2644 for (
int a = 0; a < nof_plus - 1; a++) {
2645 for (
int b = a + 1; b < nof_plus; b++) {
2655 double OpeningAngle =
2668 if (TMath::Abs(OpeningAngle) > AngleCut)
continue;
2669 if (TMath::Abs(InvmassReco) > InvMassCut)
continue;
2682 if ((rings_amount == 1 || rings_amount == 2)) {
2689 if ((rings_amount == 1 || rings_amount == 2))
2700 for (
int a = 0; a < nof_minus - 1; a++) {
2701 for (
int b = a + 1; b < nof_minus; b++) {
2708 KFParticle electron;
2710 KFParticle positron;
2712 const KFParticle* daughters[2] = {&electron, &positron};
2713 KFParticle intersection;
2714 intersection.Construct(daughters, 2);
2716 if (intersection.GetZ() > 75 || intersection.GetZ() < -5)
2720 sts1, intersection.GetX(), intersection.GetY(), intersection.GetZ());
2722 sts2, intersection.GetX(), intersection.GetY(), intersection.GetZ());
2725 double OpeningAngle =
2738 if (TMath::Abs(OpeningAngle) > AngleCut)
continue;
2739 if (TMath::Abs(InvmassReco) > InvMassCut)
continue;
2752 if ((rings_amount == 1 || rings_amount == 2))
2758 if ((rings_amount == 1 || rings_amount == 2))
2765 for (
int a = 0; a < nof_plus - 1; a++) {
2766 for (
int b = a + 1; b < nof_plus; b++) {
2773 KFParticle electron;
2775 KFParticle positron;
2777 const KFParticle* daughters[2] = {&electron, &positron};
2778 KFParticle intersection;
2779 intersection.Construct(daughters, 2);
2781 if (intersection.GetZ() > 75 || intersection.GetZ() < -5)
2785 sts1, intersection.GetX(), intersection.GetY(), intersection.GetZ());
2787 sts2, intersection.GetX(), intersection.GetY(), intersection.GetZ());
2790 double OpeningAngle =
2803 if (TMath::Abs(OpeningAngle) > AngleCut)
continue;
2804 if (TMath::Abs(InvmassReco) > InvMassCut)
continue;
2817 if ((rings_amount == 1 || rings_amount == 2))
2823 if ((rings_amount == 1 || rings_amount == 2))
2835 cout <<
"Mixing for direct photons in Target - nof entries " << nof_Target
2848 double OpeningAngle =
2861 if (TMath::Abs(OpeningAngle) > AngleCut)
continue;
2862 if (TMath::Abs(InvmassReco) > InvMassCut)
continue;
2872 if (rings_amount == 0) {
2876 if (rings_amount == 1) {
2880 if (rings_amount == 2) {
2884 if ((rings_amount == 1 || rings_amount == 2)) {
2898 cout <<
"Mixing for direct photons in Outside - nof entries " << nof_Outside
2909 const KFParticle* daughters[2] = {&electron, &positron};
2910 KFParticle intersection;
2911 intersection.Construct(daughters, 2);
2913 if (intersection.GetZ() > 75 || intersection.GetZ() < -5)
2919 TVector3 e1(electron.KFParticleBase::Px(),
2920 electron.KFParticleBase::Py(),
2921 electron.KFParticleBase::Pz());
2922 TVector3 e2(positron.KFParticleBase::Px(),
2923 positron.KFParticleBase::Py(),
2924 positron.KFParticleBase::Pz());
2930 double OpeningAngle =
2938 intersection.GetZ(),
2943 if (TMath::Abs(OpeningAngle) > AngleCut)
continue;
2944 if (TMath::Abs(InvmassReco) > InvMassCut)
continue;
2954 if (rings_amount == 0) {
2958 if (rings_amount == 1) {
2962 if (rings_amount == 2) {
2966 if ((rings_amount == 1 || rings_amount == 2)) {
2976 gDirectory->mkdir(
"direct photons");
2977 gDirectory->cd(
"direct photons");
2980 gDirectory->mkdir(
"Target");
2981 gDirectory->cd(
"Target");
2982 gDirectory->mkdir(
"CheckCuts_Target");
2983 gDirectory->cd(
"CheckCuts_Target");
2987 gDirectory->cd(
"..");
2988 gDirectory->mkdir(
"all");
2989 gDirectory->cd(
"all");
2993 gDirectory->cd(
"..");
2994 gDirectory->mkdir(
"zero");
2995 gDirectory->cd(
"zero");
2999 gDirectory->cd(
"..");
3000 gDirectory->mkdir(
"one");
3001 gDirectory->cd(
"one");
3005 gDirectory->cd(
"..");
3006 gDirectory->mkdir(
"two");
3007 gDirectory->cd(
"two");
3011 gDirectory->cd(
"..");
3012 gDirectory->mkdir(
"onetwo");
3013 gDirectory->cd(
"onetwo");
3017 gDirectory->cd(
"..");
3021 gDirectory->cd(
"..");
3024 gDirectory->mkdir(
"Outside");
3025 gDirectory->cd(
"Outside");
3026 gDirectory->mkdir(
"CheckCuts_Outside");
3027 gDirectory->cd(
"CheckCuts_Outside");
3031 gDirectory->cd(
"..");
3032 gDirectory->mkdir(
"all");
3033 gDirectory->cd(
"all");
3037 gDirectory->cd(
"..");
3038 gDirectory->mkdir(
"zero");
3039 gDirectory->cd(
"zero");
3043 gDirectory->cd(
"..");
3044 gDirectory->mkdir(
"one");
3045 gDirectory->cd(
"one");
3049 gDirectory->cd(
"..");
3050 gDirectory->mkdir(
"two");
3051 gDirectory->cd(
"two");
3055 gDirectory->cd(
"..");
3056 gDirectory->mkdir(
"onetwo");
3057 gDirectory->cd(
"onetwo");
3061 gDirectory->cd(
"..");
3065 gDirectory->cd(
"..");
3068 gDirectory->mkdir(
"Both");
3069 gDirectory->cd(
"Both");
3070 gDirectory->mkdir(
"CheckCuts_Both");
3071 gDirectory->cd(
"CheckCuts_Both");
3075 gDirectory->cd(
"..");
3076 gDirectory->mkdir(
"all");
3077 gDirectory->cd(
"all");
3081 gDirectory->cd(
"..");
3082 gDirectory->mkdir(
"zero");
3083 gDirectory->cd(
"zero");
3087 gDirectory->cd(
"..");
3088 gDirectory->mkdir(
"one");
3089 gDirectory->cd(
"one");
3093 gDirectory->cd(
"..");
3094 gDirectory->mkdir(
"two");
3095 gDirectory->cd(
"two");
3099 gDirectory->cd(
"..");
3100 gDirectory->mkdir(
"onetwo");
3101 gDirectory->cd(
"onetwo");
3105 gDirectory->cd(
"..");
3109 gDirectory->cd(
"..");
3112 gDirectory->cd(
"..");
3119 new TH1D(
"Mother_PDG_Target",
"Mother_PDG_Target; Id;#", 1000, -10, 490);
3122 "GrandMother_PDG_Target",
"GrandMother_PDG_Target; Id;#", 1000, -10, 490);
3127 new TH1D(
"Mother_PDG_Outside",
"Mother_PDG_Outside; Id;#", 1000, -10, 490);
3130 "GrandMother_PDG_Outside",
"GrandMother_PDG_Outside; Id;#", 1000, -10, 490);
3135 new TH2D(
"Pdg_vs_Distance_for_dp",
3136 "Pdg_vs_Distance_for_dp; pdg; distance in cm",
3145 new TH2D(
"P_vs_Distance_for_dp",
3146 "Distance between projected track and center of the ring (for e+ "
3147 "and e-); P in GeV/c^{2}; distance in cm",
3156 "DP_AnnTruePairs",
"DP_AnnTruePairs; Ann value ;#", 100, -1.2, 1.2);
3159 "DP_AnnFalsePairs",
"DP_AnnFalsePairs; Ann value ;#", 100, -1.2, 1.2);
3162 new TH1D(
"DP_AnnTruePairs_AfterCuts",
3163 "DP_AnnTruePairs_AfterCuts; Ann value ;#",
3169 new TH1D(
"DP_AnnFalsePairs_AfterCuts",
3170 "DP_AnnFalsePairs_AfterCuts; Ann value ;#",
3180 new TH2D(
"DP_candidates_InvMass_vs_OA_Both",
3181 "DP_candidates_InvMass_vs_OA_Both; invariant mass in GeV/c^{2}; "
3182 "opening angle in degree",
3191 "DP_InvMass_vs_OA_Both; invariant mass in "
3192 "GeV/c^{2}; opening angle in degree",
3201 new TH1D(
"DP_candidates_InvMass_Both",
3202 "DP_candidates_InvMass_Both; invariant mass in GeV/c^{2};#",
3208 "DP_InvMass_Both; invariant mass in GeV/c^{2};#",
3214 new TH1D(
"DP_candidates_OA_Both",
3215 "DP_candidates_OA_Both; opening angle in degree;#",
3221 "DP_OA_Both",
"DP_OA_Both; opening angle in degree;#", 300, -0.1, 29.9);
3224 new TH1D(
"DP_candidates_PlaneAngles_last_Both",
3225 "DP_candidates_PlaneAngles_last_Both; #theta angle in degree;#",
3231 new TH1D(
"DP_PlaneAngles_last_Both",
3232 "DP_PlaneAngles_last_Both; #theta angle in degree;#",
3238 new TH1D(
"DP_candidates_PlaneAngles_first_Both",
3239 "DP_candidates_PlaneAngles_first_Both; #theta angle in degree;#",
3245 new TH1D(
"DP_PlaneAngles_first_Both",
3246 "DP_PlaneAngles_first_Both; #theta angle in degree;#",
3254 new TH2D(
"DP_candidates_InvMass_vs_OA_Target",
3255 "DP_candidates_InvMass_vs_OA_Target; invariant mass in GeV/c^{2}; "
3256 "opening angle in degree",
3265 "DP_InvMass_vs_OA_Target; invariant mass "
3266 "in GeV/c^{2}; opening angle in degree",
3275 new TH1D(
"DP_candidates_InvMass_Target",
3276 "DP_candidates_InvMass_Target; invariant mass in GeV/c^{2};#",
3282 new TH1D(
"DP_InvMass_Target",
3283 "DP_InvMass_Target; invariant mass in GeV/c^{2};#",
3289 new TH1D(
"DP_candidates_OA_Target",
3290 "DP_candidates_OA_Target; opening angle in degree;#",
3296 "DP_OA_Target",
"DP_OA_Target; opening angle in degree;#", 300, -0.1, 29.9);
3299 new TH1D(
"DP_candidates_PlaneAngles_last_Target",
3300 "DP_candidates_PlaneAngles_last_Target; #theta angle in degree;#",
3306 new TH1D(
"DP_PlaneAngles_last_Target",
3307 "DP_PlaneAngles_last_Target; #theta angle in degree;#",
3313 new TH1D(
"DP_candidates_PlaneAngles_first_Target",
3314 "DP_candidates_PlaneAngles_first_Target; #theta angle in degree;#",
3320 new TH1D(
"DP_PlaneAngles_first_Target",
3321 "DP_PlaneAngles_first_Target; #theta angle in degree;#",
3329 new TH2D(
"DP_candidates_InvMass_vs_OA_Outside",
3330 "DP_candidates_InvMass_vs_OA_Outside; invariant mass in "
3331 "GeV/c^{2}; opening angle in degree",
3340 new TH2D(
"DP_InvMass_vs_OA_Outside",
3341 "DP_InvMass_vs_OA_Outside; invariant mass in GeV/c^{2}; opening "
3351 new TH1D(
"DP_candidates_InvMass_Outside",
3352 "DP_candidates_InvMass_Outside; invariant mass in GeV/c^{2};#",
3358 new TH1D(
"DP_InvMass_Outside",
3359 "DP_InvMass_Outside; invariant mass in GeV/c^{2};#",
3365 new TH1D(
"DP_candidates_OA_Outside",
3366 "DP_candidates_OA_Outside; opening angle in degree;#",
3372 "DP_OA_Outside; opening angle in degree;#",
3378 new TH1D(
"DP_candidates_PlaneAngles_last_Outside",
3379 "DP_candidates_PlaneAngles_last_Outside; #theta angle in degree;#",
3385 new TH1D(
"DP_PlaneAngles_last_Outside",
3386 "DP_PlaneAngles_last_Outside; #theta angle in degree;#",
3392 "DP_candidates_PlaneAngles_first_Outside",
3393 "DP_candidates_PlaneAngles_first_Outside; #theta angle in degree;#",
3399 new TH1D(
"DP_PlaneAngles_first_Outside",
3400 "DP_PlaneAngles_first_Outside; #theta angle in degree;#",
3409 new TH1D(
"DP_InvMassReco_all_Target",
3410 "DP_InvMassReco_all_Target; invariant mass in GeV/c^{2};#",
3416 new TH1D(
"DP_OpeningAngleReco_all_Target",
3417 "DP_OpeningAngleReco_all_Target; angle [deg];#",
3423 new TH1D(
"DP_Pdg_all_Target",
"DP_Pdg_all_Target; Id;#", 1000, -10, 490);
3426 "DP_P_reco_all_Target; P in GeV/c^{2};#",
3432 new TH1D(
"DP_Pt_reco_all_Target",
3433 "DP_Pt_reco_all_Target; P_{t} in GeV/c^{2};#",
3439 new TH1D(
"Ph_fromTarget_Pt_reco_all_Target",
3440 "Ph_fromTarget_Pt_reco_all_Target; P_{t} in GeV/c^{2};#",
3446 new TH1D(
"Ph_fromPions_Pt_reco_all_Target",
3447 "Ph_fromPions_Pt_reco_all_Target; P_{t} in GeV/c^{2};#",
3453 new TH1D(
"Ph_fromEtas_Pt_reco_all_Target",
3454 "Ph_fromEtas_Pt_reco_all_Target; P_{t} in GeV/c^{2};#",
3460 new TH1D(
"Ph_fromDalitz_Pt_reco_all_Target",
3461 "Ph_fromDalitz_Pt_reco_all_Target; P_{t} in GeV/c^{2};#",
3469 new TH1D(
"DP_InvMassReco_zero_Target",
3470 "DP_InvMassReco_zero_Target; invariant mass in GeV/c^{2};#",
3476 new TH1D(
"DP_OpeningAngleReco_zero_Target",
3477 "DP_OpeningAngleReco_zero_Target; angle [deg];#",
3483 new TH1D(
"DP_Pdg_zero_Target",
"DP_Pdg_zero_Target; Id;#", 1000, -10, 490);
3486 "DP_P_reco_zero_Target; P in GeV/c^{2};#",
3492 new TH1D(
"DP_Pt_reco_zero_Target",
3493 "DP_Pt_reco_zero_Target; P_{t} in GeV/c^{2};#",
3499 new TH1D(
"Ph_fromTarget_Pt_reco_zero_Target",
3500 "Ph_fromTarget_Pt_reco_zero_Target; P_{t} in GeV/c^{2};#",
3506 new TH1D(
"Ph_fromPions_Pt_reco_zero_Target",
3507 "Ph_fromPions_Pt_reco_zero_Target; P_{t} in GeV/c^{2};#",
3513 new TH1D(
"Ph_fromEtas_Pt_reco_zero_Target",
3514 "Ph_fromEtas_Pt_reco_zero_Target; P_{t} in GeV/c^{2};#",
3520 new TH1D(
"Ph_fromDalitz_Pt_reco_zero_Target",
3521 "Ph_fromDalitz_Pt_reco_zero_Target; P_{t} in GeV/c^{2};#",
3529 new TH1D(
"DP_InvMassReco_one_Target",
3530 "DP_InvMassReco_one_Target; invariant mass in GeV/c^{2};#",
3536 new TH1D(
"DP_OpeningAngleReco_one_Target",
3537 "DP_OpeningAngleReco_one_Target; angle [deg];#",
3543 new TH1D(
"DP_Pdg_one_Target",
"DP_Pdg_one_Target; Id;#", 1000, -10, 490);
3546 "DP_P_reco_one_Target; P in GeV/c^{2};#",
3552 new TH1D(
"DP_Pt_reco_one_Target",
3553 "DP_Pt_reco_one_Target; P_{t} in GeV/c^{2};#",
3559 new TH1D(
"Ph_fromTarget_Pt_reco_one_Target",
3560 "Ph_fromTarget_Pt_reco_one_Target; P_{t} in GeV/c^{2};#",
3566 new TH1D(
"Ph_fromPions_Pt_reco_one_Target",
3567 "Ph_fromPions_Pt_reco_one_Target; P_{t} in GeV/c^{2};#",
3573 new TH1D(
"Ph_fromEtas_Pt_reco_one_Target",
3574 "Ph_fromEtas_Pt_reco_one_Target; P_{t} in GeV/c^{2};#",
3580 new TH1D(
"Ph_fromDalitz_Pt_reco_one_Target",
3581 "Ph_fromDalitz_Pt_reco_one_Target; P_{t} in GeV/c^{2};#",
3589 new TH1D(
"DP_InvMassReco_two_Target",
3590 "DP_InvMassReco_two_Target; invariant mass in GeV/c^{2};#",
3596 new TH1D(
"DP_OpeningAngleReco_two_Target",
3597 "DP_OpeningAngleReco_two_Target; angle [deg];#",
3603 new TH1D(
"DP_Pdg_two_Target",
"DP_Pdg_two_Target; Id;#", 1000, -10, 490);
3606 "DP_P_reco_two_Target; P in GeV/c^{2};#",
3612 new TH1D(
"DP_Pt_reco_two_Target",
3613 "DP_Pt_reco_two_Target; P_{t} in GeV/c^{2};#",
3619 new TH1D(
"Ph_fromTarget_Pt_reco_two_Target",
3620 "Ph_fromTarget_Pt_reco_two_Target; P_{t} in GeV/c^{2};#",
3626 new TH1D(
"Ph_fromPions_Pt_reco_two_Target",
3627 "Ph_fromPions_Pt_reco_two_Target; P_{t} in GeV/c^{2};#",
3633 new TH1D(
"Ph_fromEtas_Pt_reco_two_Target",
3634 "Ph_fromEtas_Pt_reco_two_Target; P_{t} in GeV/c^{2};#",
3640 new TH1D(
"Ph_fromXi_Pt_reco_two_Target",
3641 "Ph_fromXi_Pt_reco_two_Target; P_{t} in GeV/c^{2};#",
3647 new TH1D(
"Ph_fromDalitz_Pt_reco_two_Target",
3648 "Ph_fromDalitz_Pt_reco_two_Target; P_{t} in GeV/c^{2};#",
3654 new TH1D(
"Ph_fromOther_Pt_reco_two_Target",
3655 "Ph_fromOther_Pt_reco_two_Target; P_{t} in GeV/c^{2};#",
3661 new TH1D(
"Ph_twoFromTarget_Pt_reco_two_Target",
3662 "Ph_twoFromTarget_Pt_reco_two_Target; P_{t} in GeV/c^{2};#",
3668 new TH1D(
"Ph_fromCombinatorial_Pt_reco_two_Target",
3669 "Ph_fromCombinatorial_Pt_reco_two_Target; P_{t} in GeV/c^{2};#",
3675 new TH1D(
"Ph_fromConversion_Pt_reco_two_Target",
3676 "Ph_fromConversion_Pt_reco_two_Target; P_{t} in GeV/c^{2};#",
3683 "twoFromTarget_PDG_two_Target; Id;#",
3689 new TH1D(
"fromCombinatorial_PDG_two_Target",
3690 "fromCombinatorial_PDG_two_Target; Id;#",
3696 new TH1D(
"CombinatorialMotherPdg_two_Target",
3697 "CombinatorialMotherPdg_two_Target; Id;#",
3703 new TH1D(
"CombinatorialGrMotherPdg_two_Target",
3704 "CombinatorialGrMotherPdg_two_Target; Id;#",
3710 "Electrons_two_Target; P_{t} in GeV/c^{2};#",
3716 "Pions_two_Target",
"Pions_two_Target; P_{t} in GeV/c^{2};#", 30, 0, 3);
3719 new TH1D(
"PionElectron_two_Target",
3720 "PionElectron_two_Target; P_{t} in GeV/c^{2};#",
3726 new TH1D(
"elsePionOrElectron_two_Target",
3727 "elsePionOrElectron_two_Target; P_{t} in GeV/c^{2};#",
3733 new TH1D(
"DalitzAndConversion_Pt_two_Target",
3734 "DalitzAndConversion_Pt_two_Target; P_{t} in GeV/c^{2};#",
3740 new TH1D(
"DoubleConversion_Pt_two_Target",
3741 "DoubleConversion_Pt_two_Target; P_{t} in GeV/c^{2};#",
3747 new TH1D(
"fromFireball_P_two_Target",
3748 "fromFireball_P_two_Target; P_{t} in GeV/c^{2};#",
3754 new TH1D(
"twoFromTarget_P_two_Target",
3755 "twoFromTarget_P_two_Target; P_{t} in GeV/c^{2};#",
3761 new TH1D(
"fromCombinatorial_electron_P_two_Target",
3762 "fromCombinatorial_electron_P_two_Target; P_{t} in GeV/c^{2};#",
3768 new TH1D(
"fromCombinatorial_NOTelectron_P_two_Target",
3769 "fromCombinatorial_NOTelectron_P_two_Target; P_{t} in GeV/c^{2};#",
3778 new TH1D(
"DP_InvMassReco_onetwo_Target",
3779 "DP_InvMassReco_onetwo_Target; invariant mass in GeV/c^{2};#",
3785 new TH1D(
"DP_OpeningAngleReco_onetwo_Target",
3786 "DP_OpeningAngleReco_onetwo_Target; angle [deg];#",
3792 "DP_Pdg_onetwo_Target",
"DP_Pdg_onetwo_Target; Id;#", 1000, -10, 490);
3795 new TH1D(
"DP_P_reco_onetwo_Target",
3796 "DP_P_reco_onetwo_Target; P in GeV/c^{2};#",
3802 new TH1D(
"DP_Pt_reco_onetwo_Target",
3803 "DP_Pt_reco_onetwo_Target; P_{t} in GeV/c^{2};#",
3809 new TH1D(
"Ph_fromTarget_Pt_reco_onetwo_Target",
3810 "Ph_fromTarget_Pt_reco_onetwo_Target; P_{t} in GeV/c^{2};#",
3816 new TH1D(
"Ph_fromPions_Pt_reco_onetwo_Target",
3817 "Ph_fromPions_Pt_reco_onetwo_Target; P_{t} in GeV/c^{2};#",
3823 new TH1D(
"Ph_fromEtas_Pt_reco_onetwo_Target",
3824 "Ph_fromEtas_Pt_reco_onetwo_Target; P_{t} in GeV/c^{2};#",
3830 new TH1D(
"Ph_fromXi_Pt_reco_onetwo_Target",
3831 "Ph_fromXi_Pt_reco_onetwo_Target; P_{t} in GeV/c^{2};#",
3837 new TH1D(
"Ph_fromDalitz_Pt_reco_onetwo_Target",
3838 "Ph_fromDalitz_Pt_reco_onetwo_Target; P_{t} in GeV/c^{2};#",
3844 new TH1D(
"Ph_fromOther_Pt_reco_onetwo_Target",
3845 "Ph_fromOther_Pt_reco_onetwo_Target; P_{t} in GeV/c^{2};#",
3851 new TH1D(
"Ph_twoFromTarget_Pt_reco_onetwo_Target",
3852 "Ph_twoFromTarget_Pt_reco_onetwo_Target; P_{t} in GeV/c^{2};#",
3858 new TH1D(
"Ph_fromCombinatorial_Pt_reco_onetwo_Target",
3859 "Ph_fromCombinatorial_Pt_reco_onetwo_Target; P_{t} in GeV/c^{2};#",
3866 new TH1D(
"Ph_fromConversion_Pt_reco_onetwo_Target",
3867 "Ph_fromConversion_Pt_reco_onetwo_Target; P_{t} in GeV/c^{2};#",
3875 new TH1D(
"twoFromTarget_PDG_onetwo_Target",
3876 "twoFromTarget_PDG_onetwo_Target; Id;#",
3882 new TH1D(
"fromCombinatorial_PDG_onetwo_Target",
3883 "fromCombinatorial_PDG_onetwo_Target; Id;#",
3889 new TH1D(
"CombinatorialMotherPdg_onetwo_Target",
3890 "CombinatorialMotherPdg_onetwo_Target; Id;#",
3896 new TH1D(
"CombinatorialGrMotherPdg_onetwo_Target",
3897 "CombinatorialGrMotherPdg_onetwo_Target; Id;#",
3903 new TH1D(
"Electrons_onetwo_Target",
3904 "Electrons_onetwo_Target; P_{t} in GeV/c^{2};#",
3910 "Pions_onetwo_Target; P_{t} in GeV/c^{2};#",
3916 new TH1D(
"PionElectron_onetwo_Target",
3917 "PionElectron_onetwo_Target; P_{t} in GeV/c^{2};#",
3923 new TH1D(
"elsePionOrElectron_onetwo_Target",
3924 "elsePionOrElectron_onetwo_Target; P_{t} in GeV/c^{2};#",
3930 new TH1D(
"DalitzAndConversion_Pt_onetwo_Target",
3931 "DalitzAndConversion_Pt_onetwo_Target; P_{t} in GeV/c^{2};#",
3937 new TH1D(
"DoubleConversion_Pt_onetwo_Target",
3938 "DoubleConversion_Pt_onetwo_Target; P_{t} in GeV/c^{2};#",
3944 new TH1D(
"fromFireball_P_onetwo_Target",
3945 "fromFireball_P_onetwo_Target; P_{t} in GeV/c^{2};#",
3951 new TH1D(
"twoFromTarget_P_onetwo_Target",
3952 "twoFromTarget_P_onetwo_Target; P_{t} in GeV/c^{2};#",
3958 new TH1D(
"fromCombinatorial_electron_P_onetwo_Target",
3959 "fromCombinatorial_electron_P_onetwo_Target; P_{t} in GeV/c^{2};#",
3966 "fromCombinatorial_NOTelectron_P_onetwo_Target",
3967 "fromCombinatorial_NOTelectron_P_onetwo_Target; P_{t} in GeV/c^{2};#",
3977 new TH1D(
"DP_InvMassReco_all_Outside",
3978 "DP_InvMassReco_all_Outside; invariant mass in GeV/c^{2};#",
3984 new TH1D(
"DP_OpeningAngleReco_all_Outside",
3985 "DP_OpeningAngleReco_all_Outside; angle [deg];#",
3991 new TH1D(
"DP_Pdg_all_Outside",
"DP_Pdg_all_Outside; Id;#", 1000, -10, 490);
3994 "DP_P_reco_all_Outside; P in GeV/c^{2};#",
4000 new TH1D(
"DP_Pt_reco_all_Outside",
4001 "DP_Pt_reco_all_Outside; P_{t} in GeV/c^{2};#",
4007 new TH1D(
"Ph_fromTarget_Pt_reco_all_Outside",
4008 "Ph_fromTarget_Pt_reco_all_Outside; P_{t} in GeV/c^{2};#",
4014 new TH1D(
"Ph_fromPions_Pt_reco_all_Outside",
4015 "Ph_fromPions_Pt_reco_all_Outside; P_{t} in GeV/c^{2};#",
4021 new TH1D(
"Ph_fromEtas_Pt_reco_all_Outside",
4022 "Ph_fromEtas_Pt_reco_all_Outside; P_{t} in GeV/c^{2};#",
4028 new TH1D(
"Ph_fromDalitz_Pt_reco_all_Outside",
4029 "Ph_fromDalitz_Pt_reco_all_Outside; P_{t} in GeV/c^{2};#",
4037 new TH1D(
"DP_InvMassReco_zero_Outside",
4038 "DP_InvMassReco_zero_Outside; invariant mass in GeV/c^{2};#",
4044 new TH1D(
"DP_OpeningAngleReco_zero_Outside",
4045 "DP_OpeningAngleReco_zero_Outside; angle [deg];#",
4051 "DP_Pdg_zero_Outside",
"DP_Pdg_zero_Outside; Id;#", 1000, -10, 490);
4054 "DP_P_reco_zero_Outside; P in GeV/c^{2};#",
4060 new TH1D(
"DP_Pt_reco_zero_Outside",
4061 "DP_Pt_reco_zero_Outside; P_{t} in GeV/c^{2};#",
4067 new TH1D(
"Ph_fromTarget_Pt_reco_zero_Outside",
4068 "Ph_fromTarget_Pt_reco_zero_Outside; P_{t} in GeV/c^{2};#",
4074 new TH1D(
"Ph_fromPions_Pt_reco_zero_Outside",
4075 "Ph_fromPions_Pt_reco_zero_Outside; P_{t} in GeV/c^{2};#",
4081 new TH1D(
"Ph_fromEtas_Pt_reco_zero_Outside",
4082 "Ph_fromEtas_Pt_reco_zero_Outside; P_{t} in GeV/c^{2};#",
4088 new TH1D(
"Ph_fromDalitz_Pt_reco_zero_Outside",
4089 "Ph_fromDalitz_Pt_reco_zero_Outside; P_{t} in GeV/c^{2};#",
4097 new TH1D(
"DP_InvMassReco_one_Outside",
4098 "DP_InvMassReco_one_Outside; invariant mass in GeV/c^{2};#",
4104 new TH1D(
"DP_OpeningAngleReco_one_Outside",
4105 "DP_OpeningAngleReco_one_Outside; angle [deg];#",
4111 new TH1D(
"DP_Pdg_one_Outside",
"DP_Pdg_one_Outside; Id;#", 1000, -10, 490);
4114 "DP_P_reco_one_Outside; P in GeV/c^{2};#",
4120 new TH1D(
"DP_Pt_reco_one_Outside",
4121 "DP_Pt_reco_one_Outside; P_{t} in GeV/c^{2};#",
4127 new TH1D(
"Ph_fromTarget_Pt_reco_one_Outside",
4128 "Ph_fromTarget_Pt_reco_one_Outside; P_{t} in GeV/c^{2};#",
4134 new TH1D(
"Ph_fromPions_Pt_reco_one_Outside",
4135 "Ph_fromPions_Pt_reco_one_Outside; P_{t} in GeV/c^{2};#",
4141 new TH1D(
"Ph_fromEtas_Pt_reco_one_Outside",
4142 "Ph_fromEtas_Pt_reco_one_Outside; P_{t} in GeV/c^{2};#",
4148 new TH1D(
"Ph_fromDalitz_Pt_reco_one_Outside",
4149 "Ph_fromDalitz_Pt_reco_one_Outside; P_{t} in GeV/c^{2};#",
4157 new TH1D(
"DP_InvMassReco_two_Outside",
4158 "DP_InvMassReco_two_Outside; invariant mass in GeV/c^{2};#",
4164 new TH1D(
"DP_OpeningAngleReco_two_Outside",
4165 "DP_OpeningAngleReco_two_Outside; angle [deg];#",
4171 new TH1D(
"DP_Pdg_two_Outside",
"DP_Pdg_two_Outside; Id;#", 1000, -10, 490);
4174 "DP_P_reco_two_Outside; P in GeV/c^{2};#",
4180 new TH1D(
"DP_Pt_reco_two_Outside",
4181 "DP_Pt_reco_two_Outside; P_{t} in GeV/c^{2};#",
4187 new TH1D(
"Ph_fromTarget_Pt_reco_two_Outside",
4188 "Ph_fromTarget_Pt_reco_two_Outside; P_{t} in GeV/c^{2};#",
4194 new TH1D(
"Ph_fromPions_Pt_reco_two_Outside",
4195 "Ph_fromPions_Pt_reco_two_Outside; P_{t} in GeV/c^{2};#",
4201 new TH1D(
"Ph_fromEtas_Pt_reco_two_Outside",
4202 "Ph_fromEtas_Pt_reco_two_Outside; P_{t} in GeV/c^{2};#",
4208 new TH1D(
"Ph_fromXi_Pt_reco_two_Outside",
4209 "Ph_fromXi_Pt_reco_two_Outside; P_{t} in GeV/c^{2};#",
4215 new TH1D(
"Ph_fromDalitz_Pt_reco_two_Outside",
4216 "Ph_fromDalitz_Pt_reco_two_Outside; P_{t} in GeV/c^{2};#",
4222 new TH1D(
"Ph_fromOther_Pt_reco_two_Outside",
4223 "Ph_fromOther_Pt_reco_two_Outside; P_{t} in GeV/c^{2};#",
4229 new TH1D(
"Ph_twoFromTarget_Pt_reco_two_Outside",
4230 "Ph_twoFromTarget_Pt_reco_two_Outside; P_{t} in GeV/c^{2};#",
4236 new TH1D(
"Ph_fromCombinatorial_Pt_reco_two_Outside",
4237 "Ph_fromCombinatorial_Pt_reco_two_Outside; P_{t} in GeV/c^{2};#",
4243 new TH1D(
"Ph_fromConversion_Pt_reco_two_Outside",
4244 "Ph_fromConversion_Pt_reco_two_Outside; P_{t} in GeV/c^{2};#",
4250 new TH1D(
"twoFromTarget_PDG_two_Outside",
4251 "twoFromTarget_PDG_two_Outside; Id;#",
4257 new TH1D(
"fromCombinatorial_PDG_two_Outside",
4258 "fromCombinatorial_PDG_two_Outside; Id;#",
4264 new TH1D(
"CombinatorialMotherPdg_two_Outside",
4265 "CombinatorialMotherPdg_two_Outside; Id;#",
4271 new TH1D(
"CombinatorialGrMotherPdg_two_Outside",
4272 "CombinatorialGrMotherPdg_two_Outside; Id;#",
4278 new TH1D(
"Electrons_two_Outside",
4279 "Electrons_two_Outside; P_{t} in GeV/c^{2};#",
4285 "Pions_two_Outside",
"Pions_two_Outside; P_{t} in GeV/c^{2};#", 30, 0, 3);
4288 new TH1D(
"PionElectron_two_Outside",
4289 "PionElectron_two_Outside; P_{t} in GeV/c^{2};#",
4295 new TH1D(
"elsePionOrElectron_two_Outside",
4296 "elsePionOrElectron_two_Outside; P_{t} in GeV/c^{2};#",
4302 new TH1D(
"DalitzAndConversion_Pt_two_Outside",
4303 "DalitzAndConversion_Pt_two_Outside; P_{t} in GeV/c^{2};#",
4309 new TH1D(
"DoubleConversion_Pt_two_Outside",
4310 "DoubleConversion_Pt_two_Outside; P_{t} in GeV/c^{2};#",
4316 new TH1D(
"fromFireball_P_two_Outside",
4317 "fromFireball_P_two_Outside; P_{t} in GeV/c^{2};#",
4323 new TH1D(
"twoFromTarget_P_two_Outside",
4324 "twoFromTarget_P_two_Outside; P_{t} in GeV/c^{2};#",
4330 new TH1D(
"fromCombinatorial_electron_P_two_Outside",
4331 "fromCombinatorial_electron_P_two_Outside; P_{t} in GeV/c^{2};#",
4337 "fromCombinatorial_NOTelectron_P_two_Outside",
4338 "fromCombinatorial_NOTelectron_P_two_Outside; P_{t} in GeV/c^{2};#",
4347 new TH1D(
"DP_InvMassReco_onetwo_Outside",
4348 "DP_InvMassReco_onetwo_Outside; invariant mass in GeV/c^{2};#",
4354 new TH1D(
"DP_OpeningAngleReco_onetwo_Outside",
4355 "DP_OpeningAngleReco_onetwo_Outside; angle [deg];#",
4361 "DP_Pdg_onetwo_Outside",
"DP_Pdg_onetwo_Outside; Id;#", 1000, -10, 490);
4364 new TH1D(
"DP_P_reco_onetwo_Outside",
4365 "DP_P_reco_onetwo_Outside; P in GeV/c^{2};#",
4371 new TH1D(
"DP_Pt_reco_onetwo_Outside",
4372 "DP_Pt_reco_onetwo_Outside; P_{t} in GeV/c^{2};#",
4378 new TH1D(
"Ph_fromTarget_Pt_reco_onetwo_Outside",
4379 "Ph_fromTarget_Pt_reco_onetwo_Outside; P_{t} in GeV/c^{2};#",
4385 new TH1D(
"Ph_fromPions_Pt_reco_onetwo_Outside",
4386 "Ph_fromPions_Pt_reco_onetwo_Outside; P_{t} in GeV/c^{2};#",
4392 new TH1D(
"Ph_fromEtas_Pt_reco_onetwo_Outside",
4393 "Ph_fromEtas_Pt_reco_onetwo_Outside; P_{t} in GeV/c^{2};#",
4399 new TH1D(
"Ph_fromXi_Pt_reco_onetwo_Outside",
4400 "Ph_fromXi_Pt_reco_onetwo_Outside; P_{t} in GeV/c^{2};#",
4406 new TH1D(
"Ph_fromDalitz_Pt_reco_onetwo_Outside",
4407 "Ph_fromDalitz_Pt_reco_onetwo_Outside; P_{t} in GeV/c^{2};#",
4413 new TH1D(
"Ph_fromOther_Pt_reco_onetwo_Outside",
4414 "Ph_fromOther_Pt_reco_onetwo_Outside; P_{t} in GeV/c^{2};#",
4420 new TH1D(
"Ph_twoFromTarget_Pt_reco_onetwo_Outside",
4421 "Ph_twoFromTarget_Pt_reco_onetwo_Outside; P_{t} in GeV/c^{2};#",
4428 "Ph_fromCombinatorial_Pt_reco_onetwo_Outside",
4429 "Ph_fromCombinatorial_Pt_reco_onetwo_Outside; P_{t} in GeV/c^{2};#",
4436 new TH1D(
"Ph_fromConversion_Pt_reco_onetwo_Outside",
4437 "Ph_fromConversion_Pt_reco_onetwo_Outside; P_{t} in GeV/c^{2};#",
4444 new TH1D(
"twoFromTarget_PDG_onetwo_Outside",
4445 "twoFromTarget_PDG_onetwo_Outside; Id;#",
4451 new TH1D(
"fromCombinatorial_PDG_onetwo_Outside",
4452 "fromCombinatorial_PDG_onetwo_Outside; Id;#",
4458 new TH1D(
"CombinatorialMotherPdg_onetwo_Outside",
4459 "CombinatorialMotherPdg_onetwo_Outside; Id;#",
4465 new TH1D(
"CombinatorialGrMotherPdg_onetwo_Outside",
4466 "CombinatorialGrMotherPdg_onetwo_Outside; Id;#",
4473 new TH1D(
"Electrons_onetwo_Outside",
4474 "Electrons_onetwo_Outside; P_{t} in GeV/c^{2};#",
4480 "Pions_onetwo_Outside; P_{t} in GeV/c^{2};#",
4486 new TH1D(
"PionElectron_onetwo_Outside",
4487 "PionElectron_onetwo_Outside; P_{t} in GeV/c^{2};#",
4493 new TH1D(
"elsePionOrElectron_onetwo_Outside",
4494 "elsePionOrElectron_onetwo_Outside; P_{t} in GeV/c^{2};#",
4500 new TH1D(
"DalitzAndConversion_Pt_onetwo_Outside",
4501 "DalitzAndConversion_Pt_onetwo_Outside; P_{t} in GeV/c^{2};#",
4507 new TH1D(
"DoubleConversion_Pt_onetwo_Outside",
4508 "DoubleConversion_Pt_onetwo_Outside; P_{t} in GeV/c^{2};#",
4514 new TH1D(
"fromFireball_P_onetwo_Outside",
4515 "fromFireball_P_onetwo_Outside; P_{t} in GeV/c^{2};#",
4521 new TH1D(
"twoFromTarget_P_onetwo_Outside",
4522 "twoFromTarget_P_onetwo_Outside; P_{t} in GeV/c^{2};#",
4528 "fromCombinatorial_electron_P_onetwo_Outside",
4529 "fromCombinatorial_electron_P_onetwo_Outside; P_{t} in GeV/c^{2};#",
4536 "fromCombinatorial_NOTelectron_P_onetwo_Outside",
4537 "fromCombinatorial_NOTelectron_P_onetwo_Outside; P_{t} in GeV/c^{2};#",
4547 new TH1D(
"DP_InvMassReco_all_Both",
4548 "DP_InvMassReco_all_Both; invariant mass in GeV/c^{2};#",
4554 new TH1D(
"DP_OpeningAngleReco_all_Both",
4555 "DP_OpeningAngleReco_all_Both; angle [deg];#",
4561 new TH1D(
"DP_Pdg_all_Both",
"DP_Pdg_all_Both; Id;#", 1000, -10, 490);
4564 "DP_P_reco_all_Both",
"DP_P_reco_all_Both; P in GeV/c^{2};#", 200, 0, 10);
4567 "DP_Pt_reco_all_Both; P_{t} in GeV/c^{2};#",
4573 new TH1D(
"Ph_fromTarget_Pt_reco_all_Both",
4574 "Ph_fromTarget_Pt_reco_all_Both; P_{t} in GeV/c^{2};#",
4580 new TH1D(
"Ph_fromPions_Pt_reco_all_Both",
4581 "Ph_fromPions_Pt_reco_all_Both; P_{t} in GeV/c^{2};#",
4587 new TH1D(
"Ph_fromEtas_Pt_reco_all_Both",
4588 "Ph_fromEtas_Pt_reco_all_Both; P_{t} in GeV/c^{2};#",
4594 new TH1D(
"Ph_fromDalitz_Pt_reco_all_Both",
4595 "Ph_fromDalitz_Pt_reco_all_Both; P_{t} in GeV/c^{2};#",
4601 new TH2D(
"Ph_pt_vs_rap_est_all_Both",
4602 "Ph_pt_vs_rap_est_all_Both; rapidity y; p_{t} in GeV/c ",
4613 new TH1D(
"DP_InvMassReco_zero_Both",
4614 "DP_InvMassReco_zero_Both; invariant mass in GeV/c^{2};#",
4620 new TH1D(
"DP_OpeningAngleReco_zero_Both",
4621 "DP_OpeningAngleReco_zero_Both; angle [deg];#",
4627 new TH1D(
"DP_Pdg_zero_Both",
"DP_Pdg_zero_Both; Id;#", 1000, -10, 490);
4630 "DP_P_reco_zero_Both",
"DP_P_reco_zero_Both; P in GeV/c^{2};#", 200, 0, 10);
4633 "DP_Pt_reco_zero_Both; P_{t} in GeV/c^{2};#",
4639 new TH1D(
"Ph_fromTarget_Pt_reco_zero_Both",
4640 "Ph_fromTarget_Pt_reco_zero_Both; P_{t} in GeV/c^{2};#",
4646 new TH1D(
"Ph_fromPions_Pt_reco_zero_Both",
4647 "Ph_fromPions_Pt_reco_zero_Both; P_{t} in GeV/c^{2};#",
4653 new TH1D(
"Ph_fromEtas_Pt_reco_zero_Both",
4654 "Ph_fromEtas_Pt_reco_zero_Both; P_{t} in GeV/c^{2};#",
4660 new TH1D(
"Ph_fromDalitz_Pt_reco_zero_Both",
4661 "Ph_fromDalitz_Pt_reco_zero_Both; P_{t} in GeV/c^{2};#",
4667 new TH2D(
"Ph_pt_vs_rap_est_zero_Both",
4668 "Ph_pt_vs_rap_est_zero_Both; rapidity y; p_{t} in GeV/c ",
4679 new TH1D(
"DP_InvMassReco_one_Both",
4680 "DP_InvMassReco_one_Both; invariant mass in GeV/c^{2};#",
4686 new TH1D(
"DP_OpeningAngleReco_one_Both",
4687 "DP_OpeningAngleReco_one_Both; angle [deg];#",
4693 new TH1D(
"DP_Pdg_one_Both",
"DP_Pdg_one_Both; Id;#", 1000, -10, 490);
4696 "DP_P_reco_one_Both",
"DP_P_reco_one_Both; P in GeV/c^{2};#", 200, 0, 10);
4699 "DP_Pt_reco_one_Both; P_{t} in GeV/c^{2};#",
4705 new TH1D(
"Ph_fromTarget_Pt_reco_one_Both",
4706 "Ph_fromTarget_Pt_reco_one_Both; P_{t} in GeV/c^{2};#",
4712 new TH1D(
"Ph_fromPions_Pt_reco_one_Both",
4713 "Ph_fromPions_Pt_reco_one_Both; P_{t} in GeV/c^{2};#",
4719 new TH1D(
"Ph_fromEtas_Pt_reco_one_Both",
4720 "Ph_fromEtas_Pt_reco_one_Both; P_{t} in GeV/c^{2};#",
4726 new TH1D(
"Ph_fromDalitz_Pt_reco_one_Both",
4727 "Ph_fromDalitz_Pt_reco_one_Both; P_{t} in GeV/c^{2};#",
4733 new TH2D(
"Ph_pt_vs_rap_est_one_Both",
4734 "Ph_pt_vs_rap_est_one_Both; rapidity y; p_{t} in GeV/c ",
4745 new TH1D(
"DP_InvMassReco_two_Both",
4746 "DP_InvMassReco_two_Both; invariant mass in GeV/c^{2};#",
4752 new TH1D(
"DP_OpeningAngleReco_two_Both",
4753 "DP_OpeningAngleReco_two_Both; angle [deg];#",
4759 new TH1D(
"DP_Pdg_two_Both",
"DP_Pdg_two_Both; Id;#", 1000, -10, 490);
4762 "DP_P_reco_two_Both",
"DP_P_reco_two_Both; P in GeV/c^{2};#", 200, 0, 10);
4765 "DP_Pt_reco_two_Both; P_{t} in GeV/c^{2};#",
4771 new TH1D(
"Ph_fromTarget_Pt_reco_two_Both",
4772 "Ph_fromTarget_Pt_reco_two_Both; P_{t} in GeV/c^{2};#",
4778 new TH1D(
"Ph_fromPions_Pt_reco_two_Both",
4779 "Ph_fromPions_Pt_reco_two_Both; P_{t} in GeV/c^{2};#",
4785 new TH1D(
"Ph_fromEtas_Pt_reco_two_Both",
4786 "Ph_fromEtas_Pt_reco_two_Both; P_{t} in GeV/c^{2};#",
4792 new TH1D(
"Ph_fromDalitz_Pt_reco_two_Both",
4793 "Ph_fromDalitz_Pt_reco_two_Both; P_{t} in GeV/c^{2};#",
4799 new TH1D(
"Ph_twoFromTarget_Pt_reco_two_Both",
4800 "Ph_twoFromTarget_Pt_reco_two_Both; P_{t} in GeV/c^{2};#",
4806 new TH1D(
"Ph_fromCombinatorial_Pt_reco_two_Both",
4807 "Ph_fromCombinatorial_Pt_reco_two_Both; P_{t} in GeV/c^{2};#",
4813 new TH1D(
"Ph_fromConversion_Pt_reco_two_Both",
4814 "Ph_fromConversion_Pt_reco_two_Both; P_{t} in GeV/c^{2};#",
4820 new TH2D(
"Ph_pt_vs_rap_est_two_Both",
4821 "Ph_pt_vs_rap_est_two_Both; rapidity y; p_{t} in GeV/c ",
4832 new TH1D(
"DP_InvMassReco_onetwo_Both",
4833 "DP_InvMassReco_onetwo_Both; invariant mass in GeV/c^{2};#",
4839 new TH1D(
"DP_OpeningAngleReco_onetwo_Both",
4840 "DP_OpeningAngleReco_onetwo_Both; angle [deg];#",
4846 new TH1D(
"DP_Pdg_onetwo_Both",
"DP_Pdg_onetwo_Both; Id;#", 1000, -10, 490);
4849 "DP_P_reco_onetwo_Both; P in GeV/c^{2};#",
4855 new TH1D(
"DP_Pt_reco_onetwo_Both",
4856 "DP_Pt_reco_onetwo_Both; P_{t} in GeV/c^{2};#",
4862 new TH1D(
"Ph_fromTarget_Pt_reco_onetwo_Both",
4863 "Ph_fromTarget_Pt_reco_onetwo_Both; P_{t} in GeV/c^{2};#",
4869 new TH1D(
"Ph_fromPions_Pt_reco_onetwo_Both",
4870 "Ph_fromPions_Pt_reco_onetwo_Both; P_{t} in GeV/c^{2};#",
4876 new TH1D(
"Ph_fromEtas_Pt_reco_onetwo_Both",
4877 "Ph_fromEtas_Pt_reco_onetwo_Both; P_{t} in GeV/c^{2};#",
4883 new TH1D(
"Ph_fromDalitz_Pt_reco_onetwo_Both",
4884 "Ph_fromDalitz_Pt_reco_onetwo_Both; P_{t} in GeV/c^{2};#",
4890 new TH1D(
"Ph_twoFromTarget_Pt_reco_onetwo_Both",
4891 "Ph_twoFromTarget_Pt_reco_onetwo_Both; P_{t} in GeV/c^{2};#",
4897 new TH1D(
"Ph_fromCombinatorial_Pt_reco_onetwo_Both",
4898 "Ph_fromCombinatorial_Pt_reco_onetwo_Both; P_{t} in GeV/c^{2};#",
4904 new TH1D(
"Ph_fromConversion_Pt_reco_onetwo_Both",
4905 "Ph_fromConversion_Pt_reco_onetwo_Both; P_{t} in GeV/c^{2};#",
4911 new TH2D(
"Ph_pt_vs_rap_est_onetwo_Both",
4912 "Ph_pt_vs_rap_est_onetwo_Both; rapidity y; p_{t} in GeV/c ",
4923 "DP_EMT_Pt_all_Target; P_{t} in GeV/c^{2};#",
4929 new TH1D(
"DP_EMT_Pt_zero_Target",
4930 "DP_EMT_Pt_zero_Target; P_{t} in GeV/c^{2};#",
4936 "DP_EMT_Pt_one_Target; P_{t} in GeV/c^{2};#",
4942 "DP_EMT_Pt_two_Target; P_{t} in GeV/c^{2};#",
4948 new TH1D(
"DP_EMT_Pt_onetwo_Target",
4949 "DP_EMT_Pt_onetwo_Target; P_{t} in GeV/c^{2};#",
4957 new TH1D(
"DP_EMT_Pt_all_Outside",
4958 "DP_EMT_Pt_all_Outside; P_{t} in GeV/c^{2};#",
4964 new TH1D(
"DP_EMT_Pt_zero_Outside",
4965 "DP_EMT_Pt_zero_Outside; P_{t} in GeV/c^{2};#",
4971 new TH1D(
"DP_EMT_Pt_one_Outside",
4972 "DP_EMT_Pt_one_Outside; P_{t} in GeV/c^{2};#",
4978 new TH1D(
"DP_EMT_Pt_two_Outside",
4979 "DP_EMT_Pt_two_Outside; P_{t} in GeV/c^{2};#",
4985 new TH1D(
"DP_EMT_Pt_onetwo_Outside",
4986 "DP_EMT_Pt_onetwo_Outside; P_{t} in GeV/c^{2};#",
4994 "DP_EMT_Pt_all_Both",
"DP_EMT_Pt_all_Both; P_{t} in GeV/c^{2};#", 30, 0, 3);
4997 "DP_EMT_Pt_zero_Both; P_{t} in GeV/c^{2};#",
5003 "DP_EMT_Pt_one_Both",
"DP_EMT_Pt_one_Both; P_{t} in GeV/c^{2};#", 30, 0, 3);
5006 "DP_EMT_Pt_two_Both",
"DP_EMT_Pt_two_Both; P_{t} in GeV/c^{2};#", 30, 0, 3);
5009 new TH1D(
"DP_EMT_Pt_onetwo_Both",
5010 "DP_EMT_Pt_onetwo_Both; P_{t} in GeV/c^{2};#",
5019 new TH1D(
"DP_LK_EMT_neg_Pt_all_Target",
5020 "DP_LK_EMT_neg_Pt_all_Target; P_{t} in GeV/c^{2};#",
5026 new TH1D(
"DP_LK_EMT_neg_Pt_zero_Target",
5027 "DP_LK_EMT_neg_Pt_zero_Target; P_{t} in GeV/c^{2};#",
5033 new TH1D(
"DP_LK_EMT_neg_Pt_one_Target",
5034 "DP_LK_EMT_neg_Pt_one_Target; P_{t} in GeV/c^{2};#",
5040 new TH1D(
"DP_LK_EMT_neg_Pt_two_Target",
5041 "DP_LK_EMT_neg_Pt_two_Target; P_{t} in GeV/c^{2};#",
5047 new TH1D(
"DP_LK_EMT_neg_Pt_onetwo_Target",
5048 "DP_LK_EMT_neg_Pt_onetwo_Target; P_{t} in GeV/c^{2};#",
5055 new TH1D(
"DP_LK_EMT_pos_Pt_all_Target",
5056 "DP_LK_EMT_pos_Pt_all_Target; P_{t} in GeV/c^{2};#",
5062 new TH1D(
"DP_LK_EMT_pos_Pt_zero_Target",
5063 "DP_LK_EMT_pos_Pt_zero_Target; P_{t} in GeV/c^{2};#",
5069 new TH1D(
"DP_LK_EMT_pos_Pt_one_Target",
5070 "DP_LK_EMT_pos_Pt_one_Target; P_{t} in GeV/c^{2};#",
5076 new TH1D(
"DP_LK_EMT_pos_Pt_two_Target",
5077 "DP_LK_EMT_pos_Pt_two_Target; P_{t} in GeV/c^{2};#",
5083 new TH1D(
"DP_LK_EMT_pos_Pt_onetwo_Target",
5084 "DP_LK_EMT_pos_Pt_onetwo_Target; P_{t} in GeV/c^{2};#",
5091 new TH1D(
"DP_LK_EMT_Pt_all_Target",
5092 "DP_LK_EMT_Pt_all_Target; P_{t} in GeV/c^{2};#",
5098 new TH1D(
"DP_LK_EMT_Pt_zero_Target",
5099 "DP_LK_EMT_Pt_zero_Target; P_{t} in GeV/c^{2};#",
5105 new TH1D(
"DP_LK_EMT_Pt_one_Target",
5106 "DP_LK_EMT_Pt_one_Target; P_{t} in GeV/c^{2};#",
5112 new TH1D(
"DP_LK_EMT_Pt_two_Target",
5113 "DP_LK_EMT_Pt_two_Target; P_{t} in GeV/c^{2};#",
5119 new TH1D(
"DP_LK_EMT_Pt_onetwo_Target",
5120 "DP_LK_EMT_Pt_onetwo_Target; P_{t} in GeV/c^{2};#",
5129 new TH1D(
"DP_LK_EMT_neg_Pt_all_Outside",
5130 "DP_LK_EMT_neg_Pt_all_Outside; P_{t} in GeV/c^{2};#",
5136 new TH1D(
"DP_LK_EMT_neg_Pt_zero_Outside",
5137 "DP_LK_EMT_neg_Pt_zero_Outside; P_{t} in GeV/c^{2};#",
5143 new TH1D(
"DP_LK_EMT_neg_Pt_one_Outside",
5144 "DP_LK_EMT_neg_Pt_one_Outside; P_{t} in GeV/c^{2};#",
5150 new TH1D(
"DP_LK_EMT_neg_Pt_two_Outside",
5151 "DP_LK_EMT_neg_Pt_two_Outside; P_{t} in GeV/c^{2};#",
5157 new TH1D(
"DP_LK_EMT_neg_Pt_onetwo_Outside",
5158 "DP_LK_EMT_neg_Pt_onetwo_Outside; P_{t} in GeV/c^{2};#",
5165 new TH1D(
"DP_LK_EMT_pos_Pt_all_Outside",
5166 "DP_LK_EMT_pos_Pt_all_Outside; P_{t} in GeV/c^{2};#",
5172 new TH1D(
"DP_LK_EMT_pos_Pt_zero_Outside",
5173 "DP_LK_EMT_pos_Pt_zero_Outside; P_{t} in GeV/c^{2};#",
5179 new TH1D(
"DP_LK_EMT_pos_Pt_one_Outside",
5180 "DP_LK_EMT_pos_Pt_one_Outside; P_{t} in GeV/c^{2};#",
5186 new TH1D(
"DP_LK_EMT_pos_Pt_two_Outside",
5187 "DP_LK_EMT_pos_Pt_two_Outside; P_{t} in GeV/c^{2};#",
5193 new TH1D(
"DP_LK_EMT_pos_Pt_onetwo_Outside",
5194 "DP_LK_EMT_pos_Pt_onetwo_Outside; P_{t} in GeV/c^{2};#",
5201 new TH1D(
"DP_LK_EMT_Pt_all_Outside",
5202 "DP_LK_EMT_Pt_all_Outside; P_{t} in GeV/c^{2};#",
5208 new TH1D(
"DP_LK_EMT_Pt_zero_Outside",
5209 "DP_LK_EMT_Pt_zero_Outside; P_{t} in GeV/c^{2};#",
5215 new TH1D(
"DP_LK_EMT_Pt_one_Outside",
5216 "DP_LK_EMT_Pt_one_Outside; P_{t} in GeV/c^{2};#",
5222 new TH1D(
"DP_LK_EMT_Pt_two_Outside",
5223 "DP_LK_EMT_Pt_two_Outside; P_{t} in GeV/c^{2};#",
5229 new TH1D(
"DP_LK_EMT_Pt_onetwo_Outside",
5230 "DP_LK_EMT_Pt_onetwo_Outside; P_{t} in GeV/c^{2};#",