BmnRoot 26.10.0-18-g4e3414f9
Software framework of the BM@N experiment
Loading...
Searching...
No Matches
BmnLambdaQa.cxx
Go to the documentation of this file.
1/**
2 * \file BmnLambdaQa.cxx
3 * \brief FairTask for MC simulated lambda reconstruction performance calculation.
4 * \author Andrey Lebedev <andrey.lebedev@gsi.de> - original author for CBM experiment
5 * \author Sergey Merts <sergey.merts@gmail.com> - modifications for BMN experiment
6 * \author Alexander Lytaev <sas-lyt@ya.ru> - modifications for BMN experiment
7 * \date 2018 July
8 */
9
10#include "BmnLambdaQa.h"
11
12#include "BmnGemTrack.h"
13#include "BmnGlobalTrack.h"
14#include "BmnHistManager.h"
15#include "BmnLambdaQaReport.h"
16#include "BmnParticlePair.h"
17#include "BmnSiliconPoint.h"
18#include "BmnSimulationReport.h"
19#include "CbmMCTrack.h"
20#include "FairRootManager.h"
21
22#include <TChain.h>
23#include <TClonesArray.h>
24#include <TString.h>
25#include <boost/assign/list_of.hpp>
26#include <cmath>
27#include <fstream>
28#include <sstream>
29
30using boost::assign::list_of;
31using std::binary_search;
32using std::cout;
33using std::map;
34using std::vector;
35using namespace TMath;
36
37Int_t BmnLambdaQa::fCurrentEvent = 0;
38
39BmnLambdaQa::BmnLambdaQa()
40 : fHM(NULL)
41 , fOutputDir("./")
42 , fNOfBins(50)
43 , fLambdaPRangeMin(0.)
44 , fLambdaPRangeMax(4.)
45 , fLambdaEtaRangeMin(0.)
46 , fLambdaEtaRangeMax(5.)
47 , fLambdaThetaRangeMin(0.)
48 , fLambdaThetaRangeMax(45.)
49 , fLambdaThetaRRangeMin(0.)
50 , fLambdaThetaRRangeMax(2.5)
51 , fProtonsPRangeMin(0.)
52 , fProtonsPRangeMax(4.)
53 , fProtonsEtaRangeMin(0.)
54 , fProtonsEtaRangeMax(5.)
55 , fProtonsThetaRangeMin(0.)
56 , fProtonsThetaRangeMax(45.)
57 , fProtonsThetaRRangeMin(0.)
58 , fProtonsThetaRRangeMax(2.5)
59 , fMesonsPRangeMin(0.)
60 , fMesonsPRangeMax(2.)
61 , fMesonsEtaRangeMin(0.)
62 , fMesonsEtaRangeMax(4.)
63 , fMesonsThetaRangeMin(0.)
64 , fMesonsThetaRangeMax(45.)
65 , fMesonsThetaRRangeMin(0.)
66 , fMesonsThetaRRangeMax(2.5)
67 , fInvMassMin(1.07)
68 , fInvMassMax(1.22)
69 , fUseMCFile(kFALSE)
70 , fUseRecoFile(kFALSE)
71 , fMCTracks(NULL)
72{
73 // Initialize cut arrays ...
74 Double_t val = 0.;
75 for (Int_t i = 0; i < 2; i++) {
76 val = (i == 0) ? -DBL_MAX : DBL_MAX;
77 for (Int_t j = 0; j < 2; j++) {
78 val = (j == 0) ? -DBL_MAX : DBL_MAX;
79 fMom[i][j] = val;
80 fEta[i][j] = val;
81 fDCA[i][j] = val;
82 fDCA12[j] = val;
83 fPath[j] = val;
84 }
85 }
86}
87
88BmnLambdaQa::~BmnLambdaQa()
89{
90 if (fHM)
91 delete fHM;
92}
93
94InitStatus BmnLambdaQa::Init()
95{
96
97 fHM = new BmnHistManager();
98
99 ReadDataBranches();
100 fNLambdas = 0;
101 fNReconstructable = 0;
102
103 fNOfParticlePairs = 0;
104 fNOfParticlePairsWithMatchedLambda = 0;
105
106 fNOfParticlePairsMC = 0;
107 fNOfLambdasParticlePairsMC = 0;
108 fNOfParticlePairsMCAll = 0;
109
110 CreateH1("numberOfLambdas_LambdaQa", "", "", 1, 0, 1.);
111 CreateH1("numberOfReconstructableLambdas_LambdaQa", "", "", 1, 0, 1.);
112 CreateNumberOfReconstructableLambdaHistograms();
113 CreateNumberOfLambdaHistograms();
114 CreateLambdaRecEfficiencyHistograms();
115 CreateTwoDimensionalRecEfficiencyHistograms();
116 CreateNumberOLambdaDecayProtonsHistograms();
117 CreateNumberOfLambdaDecayMesonsHistograms();
118 CreateRecProtonsRecEfficiencyHistograms();
119 CreateRecMesonsRecEfficiencyHistograms();
120 CreateTwoDimensionalRecProtonsRecEfficiencyHistograms();
121 CreateTwoDimensionalRecMesonsRecEfficiencyHistograms();
122 CreateH1("numberOfMCReconstructedParticlePairsQA", "", "", 1, 0, 1.);
123 CreateH1("numberOfMCReconstructedLambdasQA", "", "", 1, 0, 1.);
124 CreateH1("numberOfMCReconstructedAllParticlePairsQA", "", "", 1, 0, 1.);
125 CreateTwoDimensionalReconstructedParticlePairsFromMCDataHistograms();
126 CreateReconstructedParticlePairsFromMCDataHistograms();
127 CreateTwoDimensionalReconstructedParticlePairsFromMCDataWOCutsHistograms();
128 CreateReconstructedParticlePairsFromMCDataWOCutsHistograms();
129 CreateH1("numberOfReconstructedParticlePairWithMatchedLambdaQA", "", "", 1, 0, 1.);
130 CreateH1("numberOfNotReconstructedLambdasQA", "", "", 1, 0, 1.);
131 CreateNumberOfReconstrcutedLambdaHistograms();
132 CreateReconstructionEfficiencyHistograms();
133 CreateTwoDimensionalRealRecEfficiencyHistograms();
134 CreateNumberOfReconstrcutedRecProtonHistograms();
135 CreateTwoDimensionalNumberOfReconstructedProtonsHistograms();
136 CreateRecProtonsReconstructionEfficiencyHistograms();
137 CreateNumberOfReconstrcutedRecMesonsHistograms();
138 CreateTwoDimensionalNumberOfReconstructedMesonsHistograms();
139 CreateRecMesonsReconstructionEfficiencyHistograms();
140 CreateNumberOfNotReconstructedLambdaHistograms();
141 CreateH1("numberOfReconstructedParticlePairsQA", "", "", 1, 0, 1.);
142 CreateTwoDimensionalReconstructedParticlePairsHistograms();
143 CreateReconstructedParticlePairsHistograms();
144 CreateTwoDimensionalReconstructedParticlePairsWOCutsHistograms();
145 CreateReconstructedParticlePairsWOCutsHistograms();
146
147 return kSUCCESS;
148}
149
150void BmnLambdaQa::Exec(Option_t* opt)
151{
152 if (fCurrentEvent % 100 == 0)
153 cout << "Event# " << fCurrentEvent << endl;
154 fCurrentEvent++;
155
156 vector<Int_t> indicesOfLambdaProtons;
157 vector<Int_t> indicesOfLambdaMesons;
158
159 /* array: 0-th index is mother lambda's Id, 1-st is proton's Id, 2-nd is negative pi-meson's Id*/
160 vector<array<int, 3>> indicesOfLamda;
161
162 for (Int_t iTrack = 0; iTrack < fMCTracks->GetEntriesFast(); iTrack++) {
163 CbmMCTrack* track = (CbmMCTrack*)fMCTracks->At(iTrack);
164
165 if (NULL == track) {
166 cout << "Monte-Carlo simulation branch error" << endl;
167 return;
168 }
169
170 Double_t P = ((CbmMCTrack*)fMCTracks->At(iTrack))->GetP();
171 Double_t P_x = ((CbmMCTrack*)fMCTracks->At(iTrack))->GetPx();
172 Double_t P_y = ((CbmMCTrack*)fMCTracks->At(iTrack))->GetPy();
173 Double_t P_z = ((CbmMCTrack*)fMCTracks->At(iTrack))->GetPz();
174 Double_t pRapidity = 0.5 * Log((P + P_z) / (P - P_z));
175 Double_t rapidity = ((CbmMCTrack*)fMCTracks->At(iTrack))->GetRapidity();
176
177 if (3122 == track->GetPdgCode()) {
178 fNLambdas++;
179
180 fHM->H1("numberOfLambdas_LambdaQa")->Fill(1);
181
182 fHM->H1("simNLambda_P_sim")->Fill(P);
183 fHM->H1("simNLambda_eta_sim")->Fill(pRapidity);
184 fHM->H1("simNLambda_theta_sim")->Fill(TMath::ACos(P_z / P) * TMath::RadToDeg());
185 fHM->H1("simNLambda_theta_r_sim")->Fill(rapidity);
186 fHM->H1("simNLambda_P_x_sim")->Fill(P_x);
187 fHM->H1("simNLambda_P_y_sim")->Fill(P_y);
188 fHM->H1("simNLambda_P_z_sim")->Fill(P_z);
189 fHM->H2("simNLambda_eta_P_sim")->Fill(P, pRapidity);
190 }
191
192 if (2212 == track->GetPdgCode()) {
193 int motherId = track->GetMotherId();
194
195 if (motherId < 0)
196 continue;
197
198 if (3122 == ((CbmMCTrack*)fMCTracks->At(motherId))->GetPdgCode()) {
199 indicesOfLambdaProtons.push_back(iTrack);
200
201 fHM->H1("NProtons_P_sim")->Fill(P);
202 fHM->H1("NProtons_eta_sim")->Fill(pRapidity);
203 fHM->H1("NProtons_theta_sim")->Fill(TMath::ACos(P_z / P) * TMath::RadToDeg());
204 fHM->H1("NProtons_theta_r_sim")->Fill(rapidity);
205 fHM->H2("simNRecProtons_eta_P_sim")->Fill(pRapidity, P);
206
207 Double_t nSiliconHitsOnTrackP = 0;
208
209 for (Int_t hit = 0; hit < fSiliconPoints->GetEntriesFast(); hit++) {
210 if (((BmnSiliconPoint*)fSiliconPoints->At(hit))->GetTrackID() == iTrack)
211 nSiliconHitsOnTrackP++;
212 }
213
214 if (track->GetNPoints(kGEM) + nSiliconHitsOnTrackP >= 4) {
215 fHM->H1("NRecProtons_P_sim")->Fill(P);
216 fHM->H1("NRecProtons_eta_sim")->Fill(pRapidity);
217 fHM->H1("NRecProtons_theta_sim")->Fill(TMath::ACos(P_z / P) * TMath::RadToDeg());
218 fHM->H1("NRecProtons_theta_r_sim")->Fill(rapidity);
219 fHM->H2("simNReconstructableRecProtons_eta_P_sim")->Fill(pRapidity, P);
220 }
221 }
222 }
223
224 if (-211 == track->GetPdgCode()) {
225 int motherIdMeson = track->GetMotherId();
226
227 if (motherIdMeson < 0) {
228 continue;
229 }
230
231 if (3122 == ((CbmMCTrack*)fMCTracks->At(motherIdMeson))->GetPdgCode()) {
232 indicesOfLambdaMesons.push_back(iTrack);
233
234 fHM->H1("NMesons_P_sim")->Fill(P);
235 fHM->H1("NMesons_eta_sim")->Fill(pRapidity);
236 fHM->H1("NMesons_theta_sim")->Fill(TMath::ACos(P_z / P) * TMath::RadToDeg());
237 fHM->H1("NMesons_theta_r_sim")->Fill(rapidity);
238 fHM->H2("simNRecMesons_eta_P_sim")->Fill(pRapidity, P);
239
240 Double_t nSiliconHitsOnTrackM = 0;
241
242 for (Int_t hit = 0; hit < fSiliconPoints->GetEntriesFast(); hit++) {
243 if (((BmnSiliconPoint*)fSiliconPoints->At(hit))->GetTrackID() == iTrack)
244 nSiliconHitsOnTrackM++;
245 }
246
247 if (track->GetNPoints(kGEM) + nSiliconHitsOnTrackM >= 4) {
248 fHM->H1("NRecMesons_P_sim")->Fill(P);
249 fHM->H1("NRecMesons_eta_sim")->Fill(pRapidity);
250 fHM->H1("NRecMesons_theta_sim")->Fill(TMath::ACos(P_z / P) * TMath::RadToDeg());
251 fHM->H1("NRecMesons_theta_r_sim")->Fill(rapidity);
252 fHM->H2("simNReconstructableRecMesons_eta_P_sim")->Fill(pRapidity, P);
253 }
254 }
255 }
256 }
257
258 for (size_t cnt1 = 0; cnt1 < indicesOfLambdaProtons.size(); cnt1++) {
259 for (size_t cnt2 = 0; cnt2 < indicesOfLambdaMesons.size(); cnt2++) {
260 if (((CbmMCTrack*)fMCTracks->At(indicesOfLambdaProtons[cnt1]))->GetMotherId()
261 == ((CbmMCTrack*)fMCTracks->At(indicesOfLambdaMesons[cnt2]))->GetMotherId())
262 {
263 array<int, 3> tmp = {((CbmMCTrack*)fMCTracks->At(indicesOfLambdaProtons[cnt1]))->GetMotherId(),
264 indicesOfLambdaProtons[cnt1], indicesOfLambdaMesons[cnt2]};
265 indicesOfLamda.push_back(tmp);
266 }
267 }
268 }
269
270 for (size_t cnt = 0; cnt < indicesOfLamda.size(); cnt++) {
271 Int_t nSiliconHitsOnTrackP = 0;
272 Int_t nSiliconHitsOnTrackM = 0;
273 Int_t nSSDHitsOnTrackP = 0;
274 Int_t nSSDHitsOnTrackM = 0;
275 Int_t nGEMHitsOnTrackP = 0;
276 Int_t nGEMHitsOnTrackM = 0;
277
278 for (Int_t hit = 0; hit < fSiliconPoints->GetEntriesFast(); hit++) {
279 Int_t trId = ((BmnSiliconPoint*)fSiliconPoints->At(hit))->GetTrackID();
280 if (trId == indicesOfLamda[cnt][1])
281 nSiliconHitsOnTrackP++;
282 else if (trId == indicesOfLamda[cnt][2])
283 nSiliconHitsOnTrackM++;
284 }
285
286 nGEMHitsOnTrackP = ((CbmMCTrack*)fMCTracks->At(indicesOfLamda[cnt][1]))->GetNPoints(kGEM);
287 nGEMHitsOnTrackM = ((CbmMCTrack*)fMCTracks->At(indicesOfLamda[cnt][2]))->GetNPoints(kGEM);
288
289 if (nGEMHitsOnTrackP + nSiliconHitsOnTrackP + nSSDHitsOnTrackP >= 4
290 && nGEMHitsOnTrackM + nSiliconHitsOnTrackM + nSSDHitsOnTrackM >= 4)
291 {
292 fNReconstructable++;
293
294 Double_t P = ((CbmMCTrack*)fMCTracks->At(indicesOfLamda[cnt][0]))->GetP();
295 Double_t P_x = ((CbmMCTrack*)fMCTracks->At(indicesOfLamda[cnt][0]))->GetPx();
296 Double_t P_y = ((CbmMCTrack*)fMCTracks->At(indicesOfLamda[cnt][0]))->GetPy();
297 Double_t P_z = ((CbmMCTrack*)fMCTracks->At(indicesOfLamda[cnt][0]))->GetPz();
298 Double_t pRapidity = 0.5 * Log((P + P_z) / (P - P_z));
299 Double_t rapidity = ((CbmMCTrack*)fMCTracks->At(indicesOfLamda[cnt][0]))->GetRapidity();
300
301 fHM->H1("numberOfReconstructableLambdas_LambdaQa")->Fill(1);
302
303 fHM->H1("simNReconstructableLambda_P_sim")->Fill(P);
304 fHM->H1("simNReconstructableLambda_eta_sim")->Fill(pRapidity);
305 fHM->H1("simNReconstructableLambda_theta_r_sim")->Fill(rapidity);
306 fHM->H1("simNReconstructableLambda_theta_sim")->Fill(TMath::ACos(P_z / P) * TMath::RadToDeg());
307 fHM->H1("simNReconstructableLambda_P_x_sim")->Fill(P_x);
308 fHM->H1("simNReconstructableLambda_P_y_sim")->Fill(P_y);
309 fHM->H1("simNReconstructableLambda_P_z_sim")->Fill(P_z);
310 fHM->H2("simNReconstructableLambda_eta_P_sim")->Fill(P, pRapidity);
311 }
312 }
313
314 if (fParticlePair_MC)
315 for (Int_t iPair = 0; iPair < fParticlePair_MC->GetEntriesFast(); iPair++) {
316 BmnParticlePair* pair = (BmnParticlePair*)fParticlePair_MC->UncheckedAt(iPair);
317 Double_t momPart1 = pair->GetMomPart1();
318 Double_t momPart2 = pair->GetMomPart2();
319 Double_t etaPart1 = pair->GetEtaPart1();
320 Double_t etaPart2 = pair->GetEtaPart2();
321 Double_t dca1 = pair->GetDCA1();
322 Double_t dca2 = pair->GetDCA2();
323 Double_t dca12 = pair->GetDCA12();
324 Double_t path = pair->GetPath();
325
326 // Apply cuts ...
327 if (momPart1 < fMom[0][0] || momPart1 > fMom[0][1] || momPart2 < fMom[1][0] || momPart2 > fMom[1][1])
328 continue;
329
330 if (etaPart1 < fEta[0][0] || etaPart1 > fEta[0][1] || etaPart2 < fEta[1][0] || etaPart2 > fEta[1][1])
331 continue;
332
333 // Geom. cuts applied ...
334 if (dca1 < fDCA[0][0] || dca1 > fDCA[0][1] || dca2 < fDCA[1][0] || dca2 > fDCA[1][1])
335 continue;
336
337 if (dca12 < fDCA12[0] || dca12 > fDCA12[1])
338 continue;
339
340 if (path < fPath[0] || path > fPath[1])
341 continue;
342
343 fNOfParticlePairsMC++;
344 fHM->H1("numberOfMCReconstructedAllParticlePairsQA")->Fill(1);
345 fHM->H1("NPairsRecoFromMCInvMass")->Fill(pair->GetInvMass());
346
347 fHM->H2("NPairsRecoFromMCInvMassMomProton")->Fill(momPart1, pair->GetInvMass());
348 fHM->H2("NPairsRecoFromMCInvMassMomMeson")->Fill(momPart2, pair->GetInvMass());
349 fHM->H2("NPairsRecoFromMCInvMassEtaProton")->Fill(etaPart1, pair->GetInvMass());
350 fHM->H2("NPairsRecoFromMCInvMassEtaMeson")->Fill(etaPart2, pair->GetInvMass());
351 fHM->H2("NPairsRecoFromMCInvMassDCA1")->Fill(dca1, pair->GetInvMass());
352 fHM->H2("NPairsRecoFromMCInvMassDCA2")->Fill(dca2, pair->GetInvMass());
353 fHM->H2("NPairsRecoFromMCInvMassDCA12X")->Fill(dca12, pair->GetInvMass());
354 fHM->H2("NPairsRecoFromMCInvMassPath")->Fill(path, pair->GetInvMass());
355
356 CbmMCTrack* posMCTrack = (CbmMCTrack*)fMCTracks->At(pair->GetMCTrackIdPart1());
357 CbmMCTrack* negMCTrack = (CbmMCTrack*)fMCTracks->At(pair->GetMCTrackIdPart2());
358 if (!negMCTrack || !posMCTrack)
359 continue;
360 if ((posMCTrack->GetMotherId() != negMCTrack->GetMotherId()) || (posMCTrack->GetMotherId() < 0)) {
361 continue;
362 }
363 CbmMCTrack* lambda = (CbmMCTrack*)fMCTracks->At(posMCTrack->GetMotherId());
364 if (lambda->GetPdgCode() != 3122) {
365 continue;
366 }
367 fNOfLambdasParticlePairsMC++;
368 fHM->H1("numberOfMCReconstructedLambdasQA")->Fill(1);
369 }
370
371 if (fParticlePair_MC)
372 for (Int_t iPair = 0; iPair < fParticlePair_MC->GetEntriesFast(); iPair++) {
373 fNOfParticlePairsMCAll++;
374 fHM->H1("numberOfMCReconstructedParticlePairsQA")->Fill(1);
375
376 BmnParticlePair* pair = (BmnParticlePair*)fParticlePair_MC->At(iPair);
377
378 fHM->H1("NPairsRecoFromMCWOCutsInvMass")->Fill(pair->GetInvMass());
379
380 Double_t momPart1 = pair->GetMomPart1();
381 Double_t momPart2 = pair->GetMomPart2();
382 Double_t etaPart1 = pair->GetEtaPart1();
383 Double_t etaPart2 = pair->GetEtaPart2();
384 Double_t dca1 = pair->GetDCA1();
385 Double_t dca2 = pair->GetDCA2();
386 Double_t dca12x = pair->GetDCA12();
387 Double_t dca12y = pair->GetDCA12();
388 Double_t path = pair->GetPath();
389
390 fHM->H2("NPairsRecoFromMCWOCutsInvMassMomProton")->Fill(momPart1, pair->GetInvMass());
391 fHM->H2("NPairsRecoFromMCWOCutsInvMassMomMeson")->Fill(momPart2, pair->GetInvMass());
392 fHM->H2("NPairsRecoFromMCWOCutsInvMassEtaProton")->Fill(etaPart1, pair->GetInvMass());
393 fHM->H2("NPairsRecoFromMCWOCutsInvMassEtaMeson")->Fill(etaPart2, pair->GetInvMass());
394
395 fHM->H2("NPairsRecoFromMCWOCutsInvMassDCA1")->Fill(dca1, pair->GetInvMass());
396 fHM->H2("NPairsRecoFromMCWOCutsInvMassDCA2")->Fill(dca2, pair->GetInvMass());
397 fHM->H2("NPairsRecoFromMCWOCutsInvMassDCA12X")->Fill(dca12x, pair->GetInvMass());
398 fHM->H2("NPairsRecoFromMCWOCutsInvMassDCA12Y")->Fill(dca12y, pair->GetInvMass());
399 fHM->H2("NPairsRecoFromMCWOCutsInvMassDCA12XsubY")->Fill(abs(dca12x - dca12y), pair->GetInvMass());
400 fHM->H2("NPairsRecoFromMCWOCutsInvMassPath")->Fill(path, pair->GetInvMass());
401
402 CbmMCTrack* posMCTrack = (CbmMCTrack*)fMCTracks->At(pair->GetMCTrackIdPart1());
403 CbmMCTrack* negMCTrack = (CbmMCTrack*)fMCTracks->At(pair->GetMCTrackIdPart2());
404 if (!negMCTrack || !posMCTrack)
405 continue;
406 if ((posMCTrack->GetMotherId() != negMCTrack->GetMotherId()) || (posMCTrack->GetMotherId() < 0)) {
407 continue;
408 }
409 CbmMCTrack* lambda = (CbmMCTrack*)fMCTracks->At(posMCTrack->GetMotherId());
410 if (lambda->GetPdgCode() != 3122) {
411 continue;
412 }
413 }
414
415 /// RECO
416 if (fParticlePair_RECO)
417 for (Int_t iPair = 0; iPair < fParticlePair_RECO->GetEntriesFast(); iPair++) {
418 BmnParticlePair* pair = (BmnParticlePair*)fParticlePair_RECO->UncheckedAt(iPair);
419
420 CbmMCTrack* posPartTr = (CbmMCTrack*)fMCTracks->At(pair->GetMCTrackIdPart1());
421 if (!posPartTr) {
422 cout << "proton track obtaining error" << endl;
423 return;
424 }
425
426 CbmMCTrack* negPartTr = (CbmMCTrack*)fMCTracks->At(pair->GetMCTrackIdPart2());
427 if (!negPartTr) {
428 cout << "meson track obtaining error" << endl;
429 return;
430 }
431
432 if (-1 == posPartTr->GetMotherId() || -1 == negPartTr->GetMotherId()
433 || negPartTr->GetMotherId() != posPartTr->GetMotherId())
434 continue;
435
436 CbmMCTrack* primaryPartTr = (CbmMCTrack*)fMCTracks->At(posPartTr->GetMotherId());
437 if (!primaryPartTr) {
438 cout << "lambda track obtaining error" << endl;
439 return;
440 }
441
442 if (primaryPartTr->GetPdgCode() != 3122)
443 continue;
444
445 fNOfParticlePairsWithMatchedLambda++;
446 fHM->H1("numberOfReconstructedParticlePairWithMatchedLambdaQA")->Fill(1);
447
448 // Double_t momLambdaXMC = primaryPartTr->GetPx();
449 // Double_t momLambdaYMC = primaryPartTr->GetPy();
450 Double_t momLambdaZMC = primaryPartTr->GetPz();
451 Double_t momLambdaMC = primaryPartTr->GetP();
452 Double_t thetaMC = TMath::ACos(momLambdaZMC / momLambdaMC) * TMath::RadToDeg();
453 Double_t pRapidityLMC = 0.5 * Log((momLambdaMC + momLambdaZMC) / (momLambdaMC - momLambdaZMC));
454 Double_t thetaR = primaryPartTr->GetRapidity();
455
456 fHM->H1("recNLambda_P_rec")->Fill(momLambdaMC);
457 fHM->H1("recNLambda_theta_rec")->Fill(thetaMC);
458 fHM->H1("recNLambda_eta_rec")->Fill(pRapidityLMC);
459 fHM->H1("recNLambda_theta_r_rec")->Fill(thetaR);
460 fHM->H2("recNLambda_eta_P_sim")->Fill(pRapidityLMC, momLambdaMC);
461
462 Double_t momProtMC = posPartTr->GetP();
463 Double_t momProtZMC = posPartTr->GetPz();
464 Double_t pRapidityPMC = 0.5 * Log((momProtMC + momProtZMC) / (momProtMC - momProtZMC));
465 Double_t thetaRProt = posPartTr->GetRapidity();
466
467 fHM->H1("recNRecProtons_P_rec")->Fill(momProtMC);
468 fHM->H1("recNRecProtons_theta_rec")->Fill(TMath::ACos(momProtZMC / momProtMC) * TMath::RadToDeg());
469 fHM->H1("recNRecProtons_eta_rec")->Fill(pRapidityPMC);
470 fHM->H1("recNRecProtons_theta_r_rec")->Fill(thetaRProt);
471 fHM->H2("recNRecProtons_eta_P_sim")->Fill(pair->GetEtaPart1(), momProtMC);
472
473 Double_t momMesMC = negPartTr->GetP();
474 Double_t momMesZMC = negPartTr->GetPz();
475 Double_t pRapidityMMC = 0.5 * Log((momMesMC + momMesZMC) / (momMesMC - momMesZMC));
476 Double_t thetaRMes = negPartTr->GetRapidity();
477
478 fHM->H1("recNRecMesons_P_rec")->Fill(momMesMC);
479 fHM->H1("recNRecMesons_theta_rec")->Fill(TMath::ACos(momMesZMC / momMesMC) * TMath::RadToDeg());
480 fHM->H1("recNRecMesons_eta_rec")->Fill(pRapidityMMC);
481 fHM->H1("recNRecMesons_theta_r_rec")->Fill(thetaRMes);
482 fHM->H2("recNRecMesons_eta_P_sim")->Fill(pair->GetEtaPart1(), momMesMC);
483 }
484
485 if (fParticlePair_RECO)
486 for (Int_t iPair = 0; iPair < fParticlePair_RECO->GetEntriesFast(); iPair++) {
487 BmnParticlePair* pair = (BmnParticlePair*)fParticlePair_RECO->UncheckedAt(iPair);
488 Double_t momPart1 = pair->GetMomPart1();
489 Double_t momPart2 = pair->GetMomPart2();
490 Double_t etaPart1 = pair->GetEtaPart1();
491 Double_t etaPart2 = pair->GetEtaPart2();
492 Double_t dca1 = pair->GetDCA1();
493 Double_t dca2 = pair->GetDCA2();
494 Double_t dca12x = pair->GetDCA12();
495 Double_t dca12y = pair->GetDCA12();
496 Double_t path = pair->GetPath();
497
498 // Apply cuts ...
499 if (momPart1 < fMom[0][0] || momPart1 > fMom[0][1] || momPart2 < fMom[1][0] || momPart2 > fMom[1][1])
500 continue;
501
502 if (etaPart1 < fEta[0][0] || etaPart1 > fEta[0][1] || etaPart2 < fEta[1][0] || etaPart2 > fEta[1][1])
503 continue;
504
505 // Geom. cuts applied ...
506 if (dca1 < fDCA[0][0] || dca1 > fDCA[0][1] || dca2 < fDCA[1][0] || dca2 > fDCA[1][1])
507 continue;
508
509 if (dca12x < fDCA12[0] || dca12x > fDCA12[1])
510 continue;
511
512 if (path < fPath[0] || path > fPath[1])
513 continue;
514
515 fNOfParticlePairs++;
516 fHM->H1("numberOfReconstructedParticlePairsQA")->Fill(1);
517 fHM->H1("NPairsRecoInvMass")->Fill(pair->GetInvMass());
518
519 fHM->H2("NPairsRecoInvMassMomProton")->Fill(momPart1, pair->GetInvMass());
520 fHM->H2("NPairsRecoInvMassMomMeson")->Fill(momPart2, pair->GetInvMass());
521 fHM->H2("NPairsRecoInvMassEtaProton")->Fill(etaPart1, pair->GetInvMass());
522 fHM->H2("NPairsRecoInvMassEtaMeson")->Fill(etaPart2, pair->GetInvMass());
523 fHM->H2("NPairsRecoInvMassDCA1")->Fill(dca1, pair->GetInvMass());
524 fHM->H2("NPairsRecoInvMassDCA2")->Fill(dca2, pair->GetInvMass());
525 fHM->H2("NPairsRecoInvMassDCA12X")->Fill(dca12x, pair->GetInvMass());
526 fHM->H2("NPairsRecoInvMassDCA12Y")->Fill(dca12y, pair->GetInvMass());
527 fHM->H2("NPairsRecoInvMassDCA12XsubY")->Fill(abs(dca12x - dca12y), pair->GetInvMass());
528 fHM->H2("NPairsRecoInvMassPath")->Fill(path, pair->GetInvMass());
529 }
530
531 if (fParticlePair_RECO)
532 for (Int_t iPair = 0; iPair < fParticlePair_RECO->GetEntriesFast(); iPair++) {
533 BmnParticlePair* pair = (BmnParticlePair*)fParticlePair_RECO->At(iPair);
534
535 fHM->H1("NPairsRecoInvMassWOCuts")->Fill(pair->GetInvMass());
536
537 Double_t momPart1 = pair->GetMomPart1();
538 Double_t momPart2 = pair->GetMomPart2();
539 Double_t etaPart1 = pair->GetEtaPart1();
540 Double_t etaPart2 = pair->GetEtaPart2();
541 Double_t dca1 = pair->GetDCA1();
542 Double_t dca2 = pair->GetDCA2();
543 Double_t dca12x = pair->GetDCA12();
544 Double_t dca12y = pair->GetDCA12();
545 Double_t path = pair->GetPath();
546
547 fHM->H2("NPairsRecoInvMassMomProtonWOCuts")->Fill(momPart1, pair->GetInvMass());
548 fHM->H2("NPairsRecoInvMassMomMesonWOCuts")->Fill(momPart2, pair->GetInvMass());
549 fHM->H2("NPairsRecoInvMassEtaProtonWOCuts")->Fill(etaPart1, pair->GetInvMass());
550 fHM->H2("NPairsRecoInvMassEtaMesonWOCuts")->Fill(etaPart2, pair->GetInvMass());
551
552 fHM->H2("NPairsRecoInvMassDCA1WOCuts")->Fill(dca1, pair->GetInvMass());
553 fHM->H2("NPairsRecoInvMassDCA2WOCuts")->Fill(dca2, pair->GetInvMass());
554 fHM->H2("NPairsRecoInvMassDCA12XWOCuts")->Fill(dca12x, pair->GetInvMass());
555 fHM->H2("NPairsRecoInvMassDCA12YWOCuts")->Fill(dca12y, pair->GetInvMass());
556 fHM->H2("NPairsRecoInvMassDCA12XsubYWOCuts")->Fill(abs(dca12x - dca12y), pair->GetInvMass());
557 fHM->H2("NPairsRecoInvMassPathWOCuts")->Fill(path, pair->GetInvMass());
558 }
559}
560
561void BmnLambdaQa::Finish()
562{
563 fHM->WriteToFile();
564
565 vector<TClonesArray*> isPairBranches;
566 isPairBranches.push_back(fParticlePair_MC);
567 isPairBranches.push_back(fParticlePair_RECO);
568
569 BmnParticlePairsInfo* pairInfo = new BmnParticlePairsInfo();
570 pairInfo->setMomPart1Min(fMom[0][0]);
571 pairInfo->setMomPart1Max(fMom[0][1]);
572 pairInfo->setMomPart2Min(fMom[1][0]);
573 pairInfo->setMomPart2Max(fMom[1][1]);
574
575 pairInfo->setEtaPart1Min(fEta[0][0]);
576 pairInfo->setEtaPart1Max(fEta[0][1]);
577 pairInfo->setEtaPart2Min(fEta[1][0]);
578 pairInfo->setEtaPart2Max(fEta[1][1]);
579
580 pairInfo->setDCAPart1Min(fDCA[0][0]);
581 pairInfo->setDCAPart1Max(fDCA[0][1]);
582 pairInfo->setDCAPart2Min(fDCA[1][0]);
583 pairInfo->setDCAPart2Max(fDCA[1][1]);
584
585 pairInfo->setDCA12Min(fDCA12[0]);
586 pairInfo->setDCA12Max(fDCA12[1]);
587 pairInfo->setPathMin(fPath[0]);
588 pairInfo->setPathMax(fPath[1]);
589 BmnSimulationReport* report = new BmnLambdaQaReport(fUseMCFile, fUseRecoFile, pairInfo, isPairBranches);
590 report->Create(fHM, fOutputDir);
591 delete report;
592 delete pairInfo;
593}
594
595void BmnLambdaQa::ReadDataBranches()
596{
597 FairRootManager* ioman = FairRootManager::Instance();
598
599 if (!ioman)
600 Fatal("Init", "BmnRootManager is not instantiated");
601
602 fMCTracks = (TClonesArray*)ioman->GetObject("MCTrack");
603 fGlobalTracks = (TClonesArray*)ioman->GetObject("BmnGlobalTrack");
604 fSiliconPoints = (TClonesArray*)ioman->GetObject("SiliconPoint");
605
606 fParticlePair_MC = (TClonesArray*)ioman->GetObject("ParticlePair_MC");
607 fParticlePair_RECO = (TClonesArray*)ioman->GetObject("ParticlePair_RECO");
608
609 if (fMCTracks)
610 fUseMCFile = kTRUE;
611
612 if (fGlobalTracks)
613 fUseRecoFile = kTRUE;
614}
615
616void BmnLambdaQa::CreateH1(const string& name,
617 const string& xTitle,
618 const string& yTitle,
619 Int_t nofBins,
620 Double_t minBin,
621 Double_t maxBin)
622{
623 TH1F* h = new TH1F(name.c_str(), string(";" + xTitle + ";" + yTitle).c_str(), nofBins, minBin, maxBin);
624 fHM->Add(name, h);
625}
626
627void BmnLambdaQa::CreateH2(const string& name,
628 const string& xTitle,
629 const string& yTitle,
630 const string& zTitle,
631 Int_t nofBinsX,
632 Double_t minBinX,
633 Double_t maxBinX,
634 Int_t nofBinsY,
635 Double_t minBinY,
636 Double_t maxBinY)
637{
638 TH2F* h = new TH2F(name.c_str(), (";" + xTitle + ";" + yTitle + ";" + zTitle).c_str(), nofBinsX, minBinX, maxBinX,
639 nofBinsY, minBinY, maxBinY);
640 fHM->Add(name, h);
641}
642
643void BmnLambdaQa::CreateNumberOfReconstructableLambdaHistograms()
644{
645 CreateH1("simNReconstructableLambda_P_sim", "P_{sim}, GeV/c", "number of reconstructable #Lambda^{0}", fNOfBins,
646 fLambdaPRangeMin, fLambdaPRangeMax);
647 CreateH1("simNReconstructableLambda_eta_sim", "#eta_{sim}", "number of reconstructable #Lambda^{0}", fNOfBins,
648 fLambdaEtaRangeMin, fLambdaEtaRangeMax);
649 CreateH1("simNReconstructableLambda_theta_sim", "#theta_{sim}, degrees", "number of reconstructable #Lambda^{0}",
650 fNOfBins, fLambdaThetaRangeMin, fLambdaThetaRangeMax);
651 CreateH1("simNReconstructableLambda_theta_r_sim", "Y^{r}_{sim}", "number of reconstructable #Lambda^{0}", fNOfBins,
652 fLambdaThetaRRangeMin, fLambdaThetaRRangeMax);
653 CreateH1("simNReconstructableLambda_P_x_sim", "P^{x}_{sim}, GeV/c", "number of reconstructable #Lambda^{0}",
654 fNOfBins, -2, 2);
655 CreateH1("simNReconstructableLambda_P_y_sim", "P^{y}_{sim}, GeV/c", "number of reconstructable #Lambda^{0}",
656 fNOfBins, -2, 2);
657 CreateH1("simNReconstructableLambda_P_z_sim", "P^{z}_{sim}, GeV/c", "number of reconstructable #Lambda^{0}",
658 fNOfBins, -1, 4);
659}
660
661void BmnLambdaQa::CreateNumberOfLambdaHistograms()
662{
663 CreateH1("simNLambda_P_sim", "P_{sim}, GeV/c", "number of #Lambda^{0}", fNOfBins, fLambdaPRangeMin,
664 fLambdaPRangeMax);
665 CreateH1("simNLambda_eta_sim", "#eta_{sim}", "number of #Lambda^{0}", fNOfBins, fLambdaEtaRangeMin,
666 fLambdaEtaRangeMax);
667 CreateH1("simNLambda_theta_sim", "#theta_{sim}, degrees", "number of #Lambda^{0}", fNOfBins, fLambdaThetaRangeMin,
668 fLambdaThetaRangeMax);
669 CreateH1("simNLambda_theta_r_sim", "Y^{r}_{sim}", "number of #Lambda^{0}", fNOfBins, fLambdaThetaRRangeMin,
670 fLambdaThetaRRangeMax);
671 CreateH1("simNLambda_P_x_sim", "P^{x}_{sim}, GeV/c", "number of #Lambda^{0}", fNOfBins, -2, 2);
672 CreateH1("simNLambda_P_y_sim", "P^{y}_{sim}, GeV/c", "number of #Lambda^{0}", fNOfBins, -2, 2);
673 CreateH1("simNLambda_P_z_sim", "P^{z}_{sim}, GeV/c", "number of #Lambda^{0}", fNOfBins, -1, 4);
674}
675
676void BmnLambdaQa::CreateLambdaRecEfficiencyHistograms()
677{
678 CreateH1("simLambdaRecEff_P_sim", "P_{sim}, GeV/c", "geometrical efficiency, %", fNOfBins, fLambdaPRangeMin,
679 fLambdaPRangeMax);
680 CreateH1("simLambdaRecEff_eta_sim", "#eta_{sim}", "geometrical efficiency, % efficiency, %", fNOfBins,
681 fLambdaEtaRangeMin, fLambdaEtaRangeMax);
682 CreateH1("simLambdaRecEff_theta_sim", "#theta_{sim}, degrees", "geometrical efficiency, %", fNOfBins,
683 fLambdaThetaRangeMin, fLambdaThetaRangeMax);
684 CreateH1("simLambdaRecEff_theta_r_sim", "Y^{r}_{sim}", "geometrical efficiency, %", fNOfBins, fLambdaThetaRRangeMin,
685 fLambdaThetaRRangeMax);
686 CreateH1("simLambdaRecEff_P_x_sim", "P^{x}_{sim}, GeV/c", "geometrical efficiency, %", fNOfBins, -2, 2);
687 CreateH1("simLambdaRecEff_P_y_sim", "P^{y}_{sim}, GeV/c", "geometrical efficiency, %", fNOfBins, -2, 2);
688 CreateH1("simLambdaRecEff_P_z_sim", "P^{z}_{sim}, GeV/c", "geometrical efficiency, %", fNOfBins, -1, 4);
689}
690
691void BmnLambdaQa::CreateTwoDimensionalRecEfficiencyHistograms()
692{
693 CreateH2("simNReconstructableLambda_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "number of #Lambda^{0}", fNOfBins,
694 fLambdaPRangeMin, fLambdaPRangeMax, fNOfBins, fLambdaEtaRangeMin, fLambdaEtaRangeMax);
695 CreateH2("simNLambda_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "number of #Lambda^{0}", fNOfBins,
696 fLambdaPRangeMin, fLambdaPRangeMax, fNOfBins, fLambdaEtaRangeMin, fLambdaEtaRangeMax);
697 CreateH2("simEffRecLambda_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "geometrical efficiency, %", fNOfBins,
698 fLambdaPRangeMin, fLambdaPRangeMax, fNOfBins, fLambdaEtaRangeMin, fLambdaEtaRangeMax);
699}
700
701void BmnLambdaQa::CreateNumberOLambdaDecayProtonsHistograms()
702{
703 CreateH1("NProtons_P_sim", "P_{sim}, GeV/c", "number of lambda decay p", fNOfBins, fProtonsPRangeMin,
704 fProtonsPRangeMax);
705 CreateH1("NProtons_eta_sim", "#eta_{sim}", "number of lambda decay p", fNOfBins, fProtonsEtaRangeMin,
706 fProtonsEtaRangeMax);
707 CreateH1("NProtons_theta_sim", "#theta_{sim}, degrees", "number of lambda decay p", fNOfBins, fProtonsThetaRangeMin,
708 fProtonsThetaRangeMax);
709 CreateH1("NProtons_theta_r_sim", "Y^{r}_{sim}", "number of lambda decay p", fNOfBins, fProtonsThetaRRangeMin,
710 fProtonsThetaRRangeMax);
711 CreateH1("NRecProtons_P_sim", "P_{sim}, GeV/c", "number of reconstructable lambda decay p", fNOfBins,
712 fProtonsPRangeMin, fProtonsPRangeMax);
713 CreateH1("NRecProtons_eta_sim", "#eta_{sim}", "number of reconstructable lambda decay p", fNOfBins,
714 fProtonsEtaRangeMin, fProtonsEtaRangeMax);
715 CreateH1("NRecProtons_theta_sim", "#theta_{sim}, degrees", "number of reconstructable lambda decay p", fNOfBins,
716 fProtonsThetaRangeMin, fProtonsThetaRangeMax);
717 CreateH1("NRecProtons_theta_r_sim", "Y^{r}_{sim}, degrees", "number of reconstructable lambda decay p", fNOfBins,
718 fProtonsThetaRRangeMin, fProtonsThetaRRangeMax);
719}
720
721void BmnLambdaQa::CreateRecProtonsRecEfficiencyHistograms()
722{
723 CreateH1("simRecProtonsRecEff_P_sim", "P_{sim}, GeV/c", "geometrical efficiency, %", fNOfBins, fProtonsPRangeMin,
724 fProtonsPRangeMax);
725 CreateH1("simRecProtonsRecEff_eta_sim", "#eta_{sim}", "geometrical efficiency, %", fNOfBins, fProtonsEtaRangeMin,
726 fProtonsEtaRangeMax);
727 CreateH1("simRecProtonsRecEff_theta_sim", "#theta_{sim}, degrees", "geometrical efficiency, %", fNOfBins,
728 fProtonsThetaRangeMin, fProtonsThetaRangeMax);
729 CreateH1("simRecProtonsRecEff_theta_r_sim", "Y^{r}_{sim}, degrees", "geometrical efficiency, %", fNOfBins,
730 fProtonsThetaRRangeMin, fProtonsThetaRRangeMax);
731}
732
733void BmnLambdaQa::CreateTwoDimensionalRecProtonsRecEfficiencyHistograms()
734{
735 CreateH2("simNRecProtons_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "number of lambda decay p", fNOfBins,
736 fProtonsEtaRangeMin, fProtonsEtaRangeMax, fNOfBins, fProtonsPRangeMin, fProtonsPRangeMax);
737 CreateH2("simNReconstructableRecProtons_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "number of lambda decay p",
738 fNOfBins, fProtonsEtaRangeMin, fProtonsEtaRangeMax, fNOfBins, fProtonsPRangeMin, fProtonsPRangeMax);
739 CreateH2("simEffRecRecProtons_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "reconstruction efficiency, %", fNOfBins,
740 fProtonsEtaRangeMin, fProtonsEtaRangeMax, fNOfBins, fProtonsPRangeMin, fProtonsPRangeMax);
741}
742
743void BmnLambdaQa::CreateNumberOfLambdaDecayMesonsHistograms()
744{
745 CreateH1("NMesons_P_sim", "P_{sim}, GeV/c", "number of lambda decay #pi^{-}", fNOfBins, fMesonsPRangeMin,
746 fMesonsPRangeMax);
747 CreateH1("NMesons_eta_sim", "#eta_{sim}", "number of lambda decay #pi^{-}", fNOfBins, fMesonsEtaRangeMin,
748 fMesonsEtaRangeMax);
749 CreateH1("NMesons_theta_sim", "#theta_{sim}, degrees", "number of lambda decay #pi^{-}", fNOfBins,
750 fMesonsThetaRangeMin, fMesonsThetaRangeMax);
751 CreateH1("NMesons_theta_r_sim", "Y^{r}_{sim}, degrees", "number of lambda decay #pi^{-}", fNOfBins,
752 fMesonsThetaRRangeMin, fMesonsThetaRRangeMax);
753 CreateH1("NRecMesons_P_sim", "P_{sim}, GeV/c", "number of reconstructable lambda decay #pi^{-}", fNOfBins,
754 fMesonsPRangeMin, fMesonsPRangeMax);
755 CreateH1("NRecMesons_eta_sim", "#eta_{sim}", "number of reconstructable lambda decay #pi^{-}", fNOfBins,
756 fMesonsEtaRangeMin, fMesonsEtaRangeMax);
757 CreateH1("NRecMesons_theta_sim", "#theta_{sim}, degrees", "number of reconstructable lambda decay #pi^{-}",
758 fNOfBins, fMesonsThetaRangeMin, fMesonsThetaRangeMax);
759 CreateH1("NRecMesons_theta_r_sim", "Y^{r}_{sim}, degrees", "number of reconstructable lambda decay #pi^{-}",
760 fNOfBins, fMesonsThetaRRangeMin, fMesonsThetaRRangeMax);
761}
762
763void BmnLambdaQa::CreateRecMesonsRecEfficiencyHistograms()
764{
765 CreateH1("simRecMesonsRecEff_P_sim", "P_{sim}, GeV/c", "geometrical efficiency, %", fNOfBins, fMesonsPRangeMin,
766 fMesonsPRangeMax);
767 CreateH1("simRecMesonsRecEff_eta_sim", "#eta_{sim}", "geometrical efficiency, %", fNOfBins, fMesonsEtaRangeMin,
768 fMesonsEtaRangeMax);
769 CreateH1("simRecMesonsRecEff_theta_sim", "#theta_{sim}, degrees", "geometrical efficiency, %", fNOfBins,
770 fMesonsThetaRangeMin, fMesonsThetaRangeMax);
771 CreateH1("simRecMesonsRecEff_theta_r_sim", "Y^{r}_{sim}, degrees", "geometrical efficiency, %", fNOfBins,
772 fMesonsThetaRRangeMin, fMesonsThetaRRangeMax);
773}
774
775void BmnLambdaQa::CreateTwoDimensionalRecMesonsRecEfficiencyHistograms()
776{
777 CreateH2("simNRecMesons_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "number of lambda decay #pi^{-}", fNOfBins,
778 fMesonsEtaRangeMin, fMesonsEtaRangeMax, fNOfBins, fMesonsPRangeMin, fMesonsPRangeMax);
779 CreateH2("simNReconstructableRecMesons_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "number of lambda decay #pi^{-}",
780 fNOfBins, fMesonsEtaRangeMin, fMesonsEtaRangeMax, fNOfBins, fMesonsPRangeMin, fMesonsPRangeMax);
781 CreateH2("simEffRecRecMesons_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "geometrical efficiency, %", fNOfBins,
782 fMesonsEtaRangeMin, fMesonsEtaRangeMax, fNOfBins, fMesonsPRangeMin, fMesonsPRangeMax);
783}
784
785void BmnLambdaQa::CreateNumberOfReconstrcutedLambdaHistograms()
786{
787 CreateH1("recNLambda_P_rec", "P_{sim}, GeV/c", "number of #Lambda^{0} reconstructed", fNOfBins, fLambdaPRangeMin,
788 fLambdaPRangeMax);
789 CreateH1("recNLambda_theta_rec", "#theta_{sim}, degrees", "number of #Lambda^{0} reconstructed", fNOfBins,
790 fLambdaThetaRangeMin, fLambdaThetaRangeMax);
791 CreateH1("recNLambda_eta_rec", "#eta_{sim}", "number of #Lambda^{0} reconstructed", fNOfBins, fLambdaEtaRangeMin,
792 fLambdaEtaRangeMax);
793 CreateH1("recNLambda_theta_r_rec", "Y^{r}_{sim}", "number of #Lambda^{0} reconstructed", fNOfBins,
794 fLambdaThetaRRangeMin, fLambdaThetaRRangeMax);
795}
796
797void BmnLambdaQa::CreateReconstructionEfficiencyHistograms()
798{
799 CreateH1("recEffLambda_P_rec", "P_{sim}, GeV/c", "reconstruction efficiency, %", fNOfBins, fLambdaPRangeMin,
800 fLambdaPRangeMax);
801 CreateH1("recEffLambda_theta_rec", "#theta_{sim}, degrees", "reconstruction efficiency, %", fNOfBins,
802 fLambdaThetaRangeMin, fLambdaThetaRangeMax);
803 CreateH1("recEffLambda_eta_rec", "#eta_{sim}", "reconstruction efficiency, %", fNOfBins, fLambdaEtaRangeMin,
804 fLambdaEtaRangeMax);
805 CreateH1("recEffLambda_theta_r_rec", "Y^{r}_{sim}", "reconstruction efficiency, %", fNOfBins, fLambdaThetaRRangeMin,
806 fLambdaThetaRRangeMax);
807}
808
809void BmnLambdaQa::CreateTwoDimensionalRealRecEfficiencyHistograms()
810{
811 CreateH2("recNLambda_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "number of reconstructed #Lambda^{0}", fNOfBins,
812 fLambdaEtaRangeMin, fLambdaEtaRangeMax, fNOfBins, fLambdaPRangeMin, fLambdaPRangeMax);
813 CreateH2("recEffRecLambda_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "efficiency, %", fNOfBins, fLambdaEtaRangeMin,
814 fLambdaEtaRangeMax, fNOfBins, fLambdaPRangeMin, fLambdaPRangeMax);
815}
816
817void BmnLambdaQa::CreateNumberOfReconstrcutedRecProtonHistograms()
818{
819 CreateH1("recNRecProtons_P_rec", "P_{sim}, GeV/c", "number of lambda decay p reconstructed", fNOfBins,
820 fProtonsPRangeMin, fProtonsPRangeMax);
821 CreateH1("recNRecProtons_theta_rec", "#theta_{sim}, degrees", "number of lambda decay p reconstructed", fNOfBins,
822 fProtonsThetaRangeMin, fProtonsThetaRangeMax);
823 CreateH1("recNRecProtons_eta_rec", "#eta_{sim}", "number of lambda decay p reconstructed", fNOfBins,
824 fProtonsEtaRangeMin, fProtonsEtaRangeMax);
825 CreateH1("recNRecProtons_theta_r_rec", "Y^{r}_{sim}", "number of lambda decay p reconstructed", fNOfBins,
826 fProtonsThetaRRangeMin, fProtonsThetaRRangeMax);
827}
828
829void BmnLambdaQa::CreateRecProtonsReconstructionEfficiencyHistograms()
830{
831 CreateH1("recEffRecProtons_P_rec", "P_{sim}, GeV/c", "reconstruction efficiency, %", fNOfBins, fProtonsPRangeMin,
832 fProtonsPRangeMax);
833 CreateH1("recEffRecProtons_theta_rec", "#theta_{sim}, degrees", "reconstruction efficiency, %", fNOfBins,
834 fProtonsThetaRangeMin, fProtonsThetaRangeMax);
835 CreateH1("recEffRecProtons_eta_rec", "#eta_{sim}", "reconstruction efficiency, %", fNOfBins, fProtonsEtaRangeMin,
836 fProtonsEtaRangeMax);
837 CreateH1("recEffRecProtons_theta_r_rec", "Y^{r}_{sim}", "reconstruction efficiency, %", fNOfBins,
838 fProtonsThetaRRangeMin, fProtonsThetaRRangeMax);
839}
840
841void BmnLambdaQa::CreateTwoDimensionalNumberOfReconstructedProtonsHistograms()
842{
843 CreateH2("recNRecProtons_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "number of lambda decay p", fNOfBins,
844 fProtonsEtaRangeMin, fProtonsEtaRangeMax, fNOfBins, fProtonsPRangeMin, fProtonsPRangeMax);
845 CreateH2("recEffRecProtons_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "efficiency, %", fNOfBins,
846 fProtonsEtaRangeMin, fProtonsEtaRangeMax, fNOfBins, fProtonsPRangeMin, fProtonsPRangeMax);
847}
848
849void BmnLambdaQa::CreateNumberOfReconstrcutedRecMesonsHistograms()
850{
851 CreateH1("recNRecMesons_P_rec", "P_{sim}, GeV/c", "number of lambda decay p reconstructed", fNOfBins,
852 fMesonsPRangeMin, fMesonsPRangeMax);
853 CreateH1("recNRecMesons_theta_rec", "#theta_{sim}, degrees", "number of lambda decay p reconstructed", fNOfBins,
854 fMesonsThetaRangeMin, fMesonsThetaRangeMax);
855 CreateH1("recNRecMesons_eta_rec", "#eta_{sim}", "number of lambda decay p reconstructed", fNOfBins,
856 fMesonsEtaRangeMin, fMesonsEtaRangeMax);
857 CreateH1("recNRecMesons_theta_r_rec", "Y^{r}_{sim}", "number of lambda decay p reconstructed", fNOfBins,
858 fMesonsThetaRRangeMin, fMesonsThetaRRangeMax);
859}
860
861void BmnLambdaQa::CreateRecMesonsReconstructionEfficiencyHistograms()
862{
863 CreateH1("recEffRecMesons_P_rec", "P_{sim}, GeV/c", "reconstruction efficiency, %", fNOfBins, fMesonsPRangeMin,
864 fMesonsPRangeMax);
865 CreateH1("recEffRecMesons_theta_rec", "#theta_{sim}, degrees", "reconstruction efficiency, %", fNOfBins,
866 fMesonsThetaRangeMin, fMesonsThetaRangeMax);
867 CreateH1("recEffRecMesons_eta_rec", "#eta_{sim}", "reconstruction efficiency, %", fNOfBins, fMesonsEtaRangeMin,
868 fMesonsEtaRangeMax);
869 CreateH1("recEffRecMesons_theta_r_rec", "Y^{r}_{sim}", "reconstruction efficiency, %", fNOfBins,
870 fMesonsThetaRRangeMin, fMesonsThetaRRangeMax);
871}
872
873void BmnLambdaQa::CreateTwoDimensionalNumberOfReconstructedMesonsHistograms()
874{
875 CreateH2("recNRecMesons_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "number of lambda decay p", fNOfBins,
876 fMesonsEtaRangeMin, fMesonsEtaRangeMax, fNOfBins, fMesonsPRangeMin, fMesonsPRangeMax);
877 CreateH2("recEffRecMesons_eta_P_sim", "P_{sim}, GeV/c", "#eta_{sim}", "efficiency, %", fNOfBins, fMesonsEtaRangeMin,
878 fMesonsEtaRangeMax, fNOfBins, fMesonsPRangeMin, fMesonsPRangeMax);
879}
880
881void BmnLambdaQa::CreateReconstructedParticlePairsFromMCDataHistograms()
882{
883 CreateH1("NPairsRecoFromMCInvMass", "m_{inv}, GeV/c^2", "number of particle pairs", fNOfBins, 1.07, 1.22);
884}
885
886void BmnLambdaQa::CreateTwoDimensionalReconstructedParticlePairsFromMCDataHistograms()
887{
888 CreateH2("NPairsRecoFromMCInvMassMomProton", "P^{p}_{sim}, GeV/c", "m_{inv}, GeV/c^{2}", "number of particle pairs",
889 fNOfBins, fProtonsPRangeMin, fProtonsPRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
890 CreateH2("NPairsRecoFromMCInvMassEtaProton", "#eta^{p}_{sim}", "m_{inv}, GeV/c^{2}", "number of particle pairs",
891 fNOfBins, fProtonsEtaRangeMin, fProtonsEtaRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
892 CreateH2("NPairsRecoFromMCInvMassDCA1", "DCA^{p}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs", fNOfBins,
893 fDCA1RangeMin, fDCA1RangeMax, fNOfBins, fInvMassMin, fInvMassMax);
894 CreateH2("NPairsRecoFromMCInvMassMomMeson", "P^{#pi^{-}}_{sim}, GeV/c", "m_{inv}, GeV/c^{2}",
895 "number of particle pairs", fNOfBins, fMesonsPRangeMin, fMesonsPRangeMax, fNOfBins, fInvMassMin,
896 fInvMassMax);
897 CreateH2("NPairsRecoFromMCInvMassEtaMeson", "#eta^{#pi^{-}}_{sim}", "m_{inv}, GeV/c^{2}",
898 "number of particle pairs", fNOfBins, fMesonsEtaRangeMin, fMesonsEtaRangeMax, fNOfBins, fInvMassMin,
899 fInvMassMax);
900 CreateH2("NPairsRecoFromMCInvMassDCA2", "DCA^{#pi^{-}}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs",
901 fNOfBins, fDCA2RangeMin, fDCA2RangeMax, fNOfBins, fInvMassMin, fInvMassMax);
902 CreateH2("NPairsRecoFromMCInvMassDCA12X", "DCA12_{X}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs",
903 fNOfBins, fDCA12XRangeMin, fDCA12XRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
904 CreateH2("NPairsRecoFromMCInvMassDCA12Y", "DCA12_{Y}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs",
905 fNOfBins, fDCA12YRangeMin, fDCA12YRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
906 CreateH2("NPairsRecoFromMCInvMassDCA12XsubY", "|DCA12_{X} - DCA12_{Y}|, cm", "m_{inv}, GeV/c^{2}",
907 "number of particle pairs", fNOfBins, fDCA12XRangeMin, fDCA12XRangeMax, fNOfBins, fInvMassMin,
908 fInvMassMax);
909 CreateH2("NPairsRecoFromMCInvMassPath", "path, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs", fNOfBins,
910 fPathRangeMin, fPathRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
911}
912
913void BmnLambdaQa::CreateReconstructedParticlePairsFromMCDataWOCutsHistograms()
914{
915 CreateH1("NPairsRecoFromMCWOCutsInvMass", "m_{inv}, GeV/c^2", "number of particle pairs", fNOfBins, 1.07, 1.22);
916}
917
918void BmnLambdaQa::CreateTwoDimensionalReconstructedParticlePairsFromMCDataWOCutsHistograms()
919{
920 CreateH2("NPairsRecoFromMCWOCutsInvMassMomProton", "P^{p}_{sim}, GeV/c", "m_{inv}, GeV/c^{2}",
921 "number of particle pairs", fNOfBins, fProtonsPRangeMin, fProtonsPRangeMax, fNOfBins, fInvMassMin,
922 fInvMassMax);
923 CreateH2("NPairsRecoFromMCWOCutsInvMassEtaProton", "#eta^{p}_{sim}", "m_{inv}, GeV/c^{2}",
924 "number of particle pairs", fNOfBins, fProtonsEtaRangeMin, fProtonsEtaRangeMax, fNOfBins, fInvMassMin,
925 fInvMassMax);
926 CreateH2("NPairsRecoFromMCWOCutsInvMassDCA1", "DCA^{p}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs",
927 fNOfBins, fDCA1RangeMin, fDCA1RangeMax, fNOfBins, fInvMassMin, fInvMassMax);
928 CreateH2("NPairsRecoFromMCWOCutsInvMassMomMeson", "P^{#pi^{-}}_{sim}, GeV/c", "m_{inv}, GeV/c^{2}",
929 "number of particle pairs", fNOfBins, fMesonsPRangeMin, fMesonsPRangeMax, fNOfBins, fInvMassMin,
930 fInvMassMax);
931 CreateH2("NPairsRecoFromMCWOCutsInvMassEtaMeson", "#eta^{#pi^{-}}_{sim}", "m_{inv}, GeV/c^{2}",
932 "number of particle pairs", fNOfBins, fMesonsEtaRangeMin, fMesonsEtaRangeMax, fNOfBins, fInvMassMin,
933 fInvMassMax);
934 CreateH2("NPairsRecoFromMCWOCutsInvMassDCA2", "DCA^{#pi^{-}}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs",
935 fNOfBins, fDCA2RangeMin, fDCA2RangeMax, fNOfBins, fInvMassMin, fInvMassMax);
936 CreateH2("NPairsRecoFromMCWOCutsInvMassDCA12X", "DCA12_{X}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs",
937 fNOfBins, fDCA12XRangeMin, fDCA12XRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
938 CreateH2("NPairsRecoFromMCWOCutsInvMassDCA12Y", "DCA12_{Y}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs",
939 fNOfBins, fDCA12YRangeMin, fDCA12YRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
940 CreateH2("NPairsRecoFromMCWOCutsInvMassDCA12XsubY", "|DCA12_{X} - DCA12_{Y}|, cm", "m_{inv}, GeV/c^{2}",
941 "number of particle pairs", fNOfBins, fDCA12XRangeMin, fDCA12XRangeMax, fNOfBins, fInvMassMin,
942 fInvMassMax);
943 CreateH2("NPairsRecoFromMCWOCutsInvMassPath", "path, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs",
944 fNOfBins, fPathRangeMin, fPathRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
945}
946
947void BmnLambdaQa::CreateNumberOfNotReconstructedLambdaHistograms()
948{
949 CreateH1("recNNotRecLambda_P_rec", "P_{sim}, GeV/c", "number of #Lambda^{0} not reconstructed", fNOfBins,
950 fLambdaPRangeMin, fLambdaPRangeMax);
951 CreateH1("recNNotRecLambda_theta_rec", "#theta_{sim}, degrees", "number of #Lambda^{0} not reconstructed", fNOfBins,
952 fLambdaThetaRangeMin, fLambdaThetaRangeMax);
953 CreateH1("recNNotRecLambda_eta_rec", "#eta_{sim}", "number of #Lambda^{0} not reconstructed", fNOfBins,
954 fLambdaEtaRangeMin, fLambdaEtaRangeMax);
955 CreateH1("recNNotRecLambda_theta_r_rec", "Y^{r}_{sim}", "number of #Lambda^{0} not reconstructed", fNOfBins,
956 fLambdaThetaRRangeMin, fLambdaThetaRRangeMax);
957}
958
959void BmnLambdaQa::CreateReconstructedParticlePairsHistograms()
960{
961 CreateH1("NPairsRecoInvMass", "m_{inv}, GeV/c^2", "number of particle pairs", fNOfBins, fInvMassMin, fInvMassMax);
962}
963
964void BmnLambdaQa::CreateTwoDimensionalReconstructedParticlePairsHistograms()
965{
966 CreateH2("NPairsRecoInvMassMomProton", "P^{p}_{reco}, GeV/c", "m_{inv}, GeV/c^{2}", "number of particle pairs",
967 fNOfBins, fProtonsPRangeMin, fProtonsPRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
968 CreateH2("NPairsRecoInvMassEtaProton", "#eta^{p}_{reco}", "m_{inv}, GeV/c^{2}", "number of particle pairs",
969 fNOfBins, fProtonsEtaRangeMin, fProtonsEtaRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
970 CreateH2("NPairsRecoInvMassDCA1", "DCA^{p}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs", fNOfBins,
971 fDCA1RangeMin, fDCA1RangeMax, fNOfBins, fInvMassMin, fInvMassMax);
972 CreateH2("NPairsRecoInvMassMomMeson", "P^{#pi^{-}}_{reco}, GeV/c", "m_{inv}, GeV/c^{2}", "number of particle pairs",
973 fNOfBins, fMesonsPRangeMin, fMesonsPRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
974 CreateH2("NPairsRecoInvMassEtaMeson", "#eta^{#pi^{-}}_{reco}", "m_{inv}, GeV/c^{2}", "number of particle pairs",
975 fNOfBins, fMesonsEtaRangeMin, fMesonsEtaRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
976 CreateH2("NPairsRecoInvMassDCA2", "DCA^{#pi^{-}}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs", fNOfBins,
977 fDCA2RangeMin, fDCA2RangeMax, fNOfBins, fInvMassMin, fInvMassMax);
978 CreateH2("NPairsRecoInvMassDCA12X", "DCA12_{X}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs", fNOfBins,
979 fDCA12XRangeMin, fDCA12XRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
980 CreateH2("NPairsRecoInvMassDCA12Y", "DCA12_{Y}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs", fNOfBins,
981 fDCA12YRangeMin, fDCA12YRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
982 CreateH2("NPairsRecoInvMassDCA12XsubY", "|DCA12_{X} - DCA12_{Y}|, cm", "m_{inv}, GeV/c^{2}",
983 "number of particle pairs", fNOfBins, fDCA12XRangeMin, fDCA12XRangeMax, fNOfBins, fInvMassMin,
984 fInvMassMax);
985 CreateH2("NPairsRecoInvMassPath", "path, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs", fNOfBins,
986 fPathRangeMin, fPathRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
987}
988
989void BmnLambdaQa::CreateReconstructedParticlePairsWOCutsHistograms()
990{
991 CreateH1("NPairsRecoInvMassWOCuts", "m_{inv}, GeV/c^2", "number of particle pairs", fNOfBins, fInvMassMin,
992 fInvMassMax);
993}
994
995void BmnLambdaQa::CreateTwoDimensionalReconstructedParticlePairsWOCutsHistograms()
996{
997 CreateH2("NPairsRecoInvMassMomProtonWOCuts", "P^{p}_{reco}, GeV/c", "m_{inv}, GeV/c^{2}",
998 "number of particle pairs", fNOfBins, fProtonsPRangeMin, fProtonsPRangeMax, fNOfBins, fInvMassMin,
999 fInvMassMax);
1000 CreateH2("NPairsRecoInvMassEtaProtonWOCuts", "#eta^{p}_{reco}", "m_{inv}, GeV/c^{2}", "number of particle pairs",
1001 fNOfBins, fProtonsEtaRangeMin, fProtonsEtaRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
1002 CreateH2("NPairsRecoInvMassDCA1WOCuts", "DCA^{p}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs", fNOfBins,
1003 fDCA1RangeMin, fDCA1RangeMax, fNOfBins, fInvMassMin, fInvMassMax);
1004 CreateH2("NPairsRecoInvMassMomMesonWOCuts", "P^{#pi^{-}}_{reco}, GeV/c", "m_{inv}, GeV/c^{2}",
1005 "number of particle pairs", fNOfBins, fMesonsPRangeMin, fMesonsPRangeMax, fNOfBins, fInvMassMin,
1006 fInvMassMax);
1007 CreateH2("NPairsRecoInvMassEtaMesonWOCuts", "#eta^{#pi^{-}}_{reco}", "m_{inv}, GeV/c^{2}",
1008 "number of particle pairs", fNOfBins, fMesonsEtaRangeMin, fMesonsEtaRangeMax, fNOfBins, fInvMassMin,
1009 fInvMassMax);
1010 CreateH2("NPairsRecoInvMassDCA2WOCuts", "DCA^{#pi^{-}}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs",
1011 fNOfBins, fDCA2RangeMin, fDCA2RangeMax, fNOfBins, fInvMassMin, fInvMassMax);
1012 CreateH2("NPairsRecoInvMassDCA12XWOCuts", "DCA12_{X}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs",
1013 fNOfBins, fDCA12XRangeMin, fDCA12XRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
1014 CreateH2("NPairsRecoInvMassDCA12YWOCuts", "DCA12_{Y}, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs",
1015 fNOfBins, fDCA12YRangeMin, fDCA12YRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
1016 CreateH2("NPairsRecoInvMassDCA12XsubYWOCuts", "|DCA12_{X} - DCA12_{Y}|, cm", "m_{inv}, GeV/c^{2}",
1017 "number of particle pairs", fNOfBins, fDCA12XRangeMin, fDCA12XRangeMax, fNOfBins, fInvMassMin,
1018 fInvMassMax);
1019 CreateH2("NPairsRecoInvMassPathWOCuts", "path, cm", "m_{inv}, GeV/c^{2}", "number of particle pairs", fNOfBins,
1020 fPathRangeMin, fPathRangeMax, fNOfBins, fInvMassMin, fInvMassMax);
1021}
i
int i
Definition P4_F32vec4.h:22
BmnLambdaQaReport.h
visualising simulated lambda reconstruction performance calculation.
BmnLambdaQa.h
FairTask for MC simulated lambda reconstruction performance calculation.
Int_t
int Int_t
kFALSE
constexpr Bool_t kFALSE
Double_t
double Double_t
kTRUE
constexpr Bool_t kTRUE
Option_t
const char Option_t
Fatal
void Fatal(const char *location, const char *msgfmt,...)
kGEM
@ kGEM
Definition BmnDetectorList.h:15
BmnHistManager
Histogram manager.
Definition BmnHistManager.h:42
BmnLambdaQa::CreateTwoDimensionalRecProtonsRecEfficiencyHistograms
void CreateTwoDimensionalRecProtonsRecEfficiencyHistograms()
Definition BmnLambdaQa.cxx:733
BmnLambdaQa::~BmnLambdaQa
virtual ~BmnLambdaQa()
Destructor.
Definition BmnLambdaQa.cxx:88
BmnLambdaQa::CreateRecMesonsReconstructionEfficiencyHistograms
void CreateRecMesonsReconstructionEfficiencyHistograms()
Definition BmnLambdaQa.cxx:861
BmnLambdaQa::fMesonsThetaRangeMax
Int_t fMesonsThetaRangeMax
Definition BmnLambdaQa.h:185
BmnLambdaQa::CreateReconstructedParticlePairsHistograms
void CreateReconstructedParticlePairsHistograms()
Definition BmnLambdaQa.cxx:959
BmnLambdaQa::CreateReconstructedParticlePairsFromMCDataHistograms
void CreateReconstructedParticlePairsFromMCDataHistograms()
Definition BmnLambdaQa.cxx:881
BmnLambdaQa::fNOfParticlePairsMC
Int_t fNOfParticlePairsMC
Definition BmnLambdaQa.h:160
BmnLambdaQa::fMesonsThetaRRangeMax
Int_t fMesonsThetaRRangeMax
Definition BmnLambdaQa.h:187
BmnLambdaQa::CreateNumberOfLambdaDecayMesonsHistograms
void CreateNumberOfLambdaDecayMesonsHistograms()
Definition BmnLambdaQa.cxx:743
BmnLambdaQa::CreateTwoDimensionalReconstructedParticlePairsHistograms
void CreateTwoDimensionalReconstructedParticlePairsHistograms()
Definition BmnLambdaQa.cxx:964
BmnLambdaQa::CreateH2
void CreateH2(const string &name, const string &xTitle, const string &yTitle, const string &zTitle, Int_t nofBinsX, Double_t minBinX, Double_t maxBinX, Int_t nofBinsY, Double_t minBinY, Double_t maxBinY)
Definition BmnLambdaQa.cxx:627
BmnLambdaQa::CreateTwoDimensionalRecEfficiencyHistograms
void CreateTwoDimensionalRecEfficiencyHistograms()
Definition BmnLambdaQa.cxx:691
BmnLambdaQa::fGlobalTracks
TClonesArray * fGlobalTracks
Definition BmnLambdaQa.h:205
BmnLambdaQa::fDCA
Double_t fDCA[2][2]
Definition BmnLambdaQa.h:221
BmnLambdaQa::fNOfParticlePairsWithMatchedLambda
Int_t fNOfParticlePairsWithMatchedLambda
Definition BmnLambdaQa.h:159
BmnLambdaQa::fNOfParticlePairs
Int_t fNOfParticlePairs
Definition BmnLambdaQa.h:158
BmnLambdaQa::fMesonsThetaRRangeMin
Int_t fMesonsThetaRRangeMin
Definition BmnLambdaQa.h:186
BmnLambdaQa::fEta
Double_t fEta[2][2]
Definition BmnLambdaQa.h:217
BmnLambdaQa::fDCA12YRangeMin
Int_t fDCA12YRangeMin
Definition BmnLambdaQa.h:196
BmnLambdaQa::fDCA12XRangeMax
Int_t fDCA12XRangeMax
Definition BmnLambdaQa.h:195
BmnLambdaQa::fMesonsThetaRangeMin
Int_t fMesonsThetaRangeMin
Definition BmnLambdaQa.h:184
BmnLambdaQa::fLambdaThetaRangeMax
Int_t fLambdaThetaRangeMax
Definition BmnLambdaQa.h:169
BmnLambdaQa::fUseRecoFile
Bool_t fUseRecoFile
Definition BmnLambdaQa.h:202
BmnLambdaQa::Finish
virtual void Finish()
Derived from FairTask.
Definition BmnLambdaQa.cxx:561
BmnLambdaQa::fMesonsEtaRangeMin
Int_t fMesonsEtaRangeMin
Definition BmnLambdaQa.h:182
BmnLambdaQa::fLambdaThetaRRangeMin
Int_t fLambdaThetaRRangeMin
Definition BmnLambdaQa.h:170
BmnLambdaQa::Init
virtual InitStatus Init()
Derived from FairTask.
Definition BmnLambdaQa.cxx:94
BmnLambdaQa::CreateRecProtonsReconstructionEfficiencyHistograms
void CreateRecProtonsReconstructionEfficiencyHistograms()
Definition BmnLambdaQa.cxx:829
BmnLambdaQa::CreateNumberOfNotReconstructedLambdaHistograms
void CreateNumberOfNotReconstructedLambdaHistograms()
Definition BmnLambdaQa.cxx:947
BmnLambdaQa::fDCA1RangeMax
Int_t fDCA1RangeMax
Definition BmnLambdaQa.h:191
BmnLambdaQa::CreateTwoDimensionalNumberOfReconstructedMesonsHistograms
void CreateTwoDimensionalNumberOfReconstructedMesonsHistograms()
Definition BmnLambdaQa.cxx:873
BmnLambdaQa::fNOfParticlePairsMCAll
Int_t fNOfParticlePairsMCAll
Definition BmnLambdaQa.h:161
BmnLambdaQa::CreateH1
void CreateH1(const string &name, const string &xTitle, const string &yTitle, Int_t nofBins, Double_t minBin, Double_t maxBin)
Definition BmnLambdaQa.cxx:616
BmnLambdaQa::fMCTracks
TClonesArray * fMCTracks
Definition BmnLambdaQa.h:204
BmnLambdaQa::fNReconstructable
Int_t fNReconstructable
Definition BmnLambdaQa.h:157
BmnLambdaQa::fProtonsPRangeMax
Int_t fProtonsPRangeMax
Definition BmnLambdaQa.h:173
BmnLambdaQa::CreateNumberOfReconstrcutedLambdaHistograms
void CreateNumberOfReconstrcutedLambdaHistograms()
Definition BmnLambdaQa.cxx:785
BmnLambdaQa::CreateTwoDimensionalRealRecEfficiencyHistograms
void CreateTwoDimensionalRealRecEfficiencyHistograms()
Definition BmnLambdaQa.cxx:809
BmnLambdaQa::fDCA2RangeMin
Int_t fDCA2RangeMin
Definition BmnLambdaQa.h:192
BmnLambdaQa::fUseMCFile
Bool_t fUseMCFile
Definition BmnLambdaQa.h:201
BmnLambdaQa::CreateNumberOfLambdaHistograms
void CreateNumberOfLambdaHistograms()
Definition BmnLambdaQa.cxx:661
BmnLambdaQa::fDCA1RangeMin
Int_t fDCA1RangeMin
Definition BmnLambdaQa.h:190
BmnLambdaQa::fProtonsThetaRangeMax
Int_t fProtonsThetaRangeMax
Definition BmnLambdaQa.h:177
BmnLambdaQa::CreateReconstructedParticlePairsFromMCDataWOCutsHistograms
void CreateReconstructedParticlePairsFromMCDataWOCutsHistograms()
Definition BmnLambdaQa.cxx:913
BmnLambdaQa::fLambdaThetaRRangeMax
Int_t fLambdaThetaRRangeMax
Definition BmnLambdaQa.h:171
BmnLambdaQa::fLambdaEtaRangeMax
Int_t fLambdaEtaRangeMax
Definition BmnLambdaQa.h:167
BmnLambdaQa::fHM
BmnHistManager * fHM
Definition BmnLambdaQa.h:103
BmnLambdaQa::CreateTwoDimensionalReconstructedParticlePairsFromMCDataHistograms
void CreateTwoDimensionalReconstructedParticlePairsFromMCDataHistograms()
Definition BmnLambdaQa.cxx:886
BmnLambdaQa::fPathRangeMax
Int_t fPathRangeMax
Definition BmnLambdaQa.h:199
BmnLambdaQa::fNLambdas
Int_t fNLambdas
Definition BmnLambdaQa.h:156
BmnLambdaQa::fProtonsPRangeMin
Int_t fProtonsPRangeMin
Definition BmnLambdaQa.h:172
BmnLambdaQa::fPath
Double_t fPath[2]
Definition BmnLambdaQa.h:220
BmnLambdaQa::fMesonsPRangeMin
Int_t fMesonsPRangeMin
Definition BmnLambdaQa.h:180
BmnLambdaQa::fNOfLambdasParticlePairsMC
Int_t fNOfLambdasParticlePairsMC
Definition BmnLambdaQa.h:162
BmnLambdaQa::fLambdaThetaRangeMin
Int_t fLambdaThetaRangeMin
Definition BmnLambdaQa.h:168
BmnLambdaQa::fProtonsEtaRangeMax
Int_t fProtonsEtaRangeMax
Definition BmnLambdaQa.h:175
BmnLambdaQa::fParticlePair_RECO
TClonesArray * fParticlePair_RECO
Definition BmnLambdaQa.h:211
BmnLambdaQa::fDCA12XRangeMin
Int_t fDCA12XRangeMin
Definition BmnLambdaQa.h:194
BmnLambdaQa::CreateReconstructionEfficiencyHistograms
void CreateReconstructionEfficiencyHistograms()
Definition BmnLambdaQa.cxx:797
BmnLambdaQa::fDCA2RangeMax
Int_t fDCA2RangeMax
Definition BmnLambdaQa.h:193
BmnLambdaQa::CreateReconstructedParticlePairsWOCutsHistograms
void CreateReconstructedParticlePairsWOCutsHistograms()
Definition BmnLambdaQa.cxx:989
BmnLambdaQa::fPathRangeMin
Int_t fPathRangeMin
Definition BmnLambdaQa.h:198
BmnLambdaQa::fParticlePair_MC
TClonesArray * fParticlePair_MC
Definition BmnLambdaQa.h:210
BmnLambdaQa::ReadDataBranches
void ReadDataBranches()
Definition BmnLambdaQa.cxx:595
BmnLambdaQa::CreateTwoDimensionalRecMesonsRecEfficiencyHistograms
void CreateTwoDimensionalRecMesonsRecEfficiencyHistograms()
Definition BmnLambdaQa.cxx:775
BmnLambdaQa::fProtonsEtaRangeMin
Int_t fProtonsEtaRangeMin
Definition BmnLambdaQa.h:174
BmnLambdaQa::fMesonsEtaRangeMax
Int_t fMesonsEtaRangeMax
Definition BmnLambdaQa.h:183
BmnLambdaQa::CreateRecProtonsRecEfficiencyHistograms
void CreateRecProtonsRecEfficiencyHistograms()
Definition BmnLambdaQa.cxx:721
BmnLambdaQa::fMom
Double_t fMom[2][2]
Definition BmnLambdaQa.h:216
BmnLambdaQa::fProtonsThetaRRangeMin
Int_t fProtonsThetaRRangeMin
Definition BmnLambdaQa.h:178
BmnLambdaQa::CreateTwoDimensionalReconstructedParticlePairsWOCutsHistograms
void CreateTwoDimensionalReconstructedParticlePairsWOCutsHistograms()
Definition BmnLambdaQa.cxx:995
BmnLambdaQa::fDCA12
Double_t fDCA12[2]
Definition BmnLambdaQa.h:222
BmnLambdaQa::CreateNumberOLambdaDecayProtonsHistograms
void CreateNumberOLambdaDecayProtonsHistograms()
Definition BmnLambdaQa.cxx:701
BmnLambdaQa::fInvMassMin
Int_t fInvMassMin
Definition BmnLambdaQa.h:188
BmnLambdaQa::fOutputDir
string fOutputDir
Definition BmnLambdaQa.h:104
BmnLambdaQa::fInvMassMax
Int_t fInvMassMax
Definition BmnLambdaQa.h:189
BmnLambdaQa::CreateLambdaRecEfficiencyHistograms
void CreateLambdaRecEfficiencyHistograms()
Definition BmnLambdaQa.cxx:676
BmnLambdaQa::CreateTwoDimensionalReconstructedParticlePairsFromMCDataWOCutsHistograms
void CreateTwoDimensionalReconstructedParticlePairsFromMCDataWOCutsHistograms()
Definition BmnLambdaQa.cxx:918
BmnLambdaQa::fSiliconPoints
TClonesArray * fSiliconPoints
Definition BmnLambdaQa.h:206
BmnLambdaQa::CreateNumberOfReconstrcutedRecProtonHistograms
void CreateNumberOfReconstrcutedRecProtonHistograms()
Definition BmnLambdaQa.cxx:817
BmnLambdaQa::fLambdaPRangeMax
Int_t fLambdaPRangeMax
Definition BmnLambdaQa.h:165
BmnLambdaQa::fLambdaEtaRangeMin
Int_t fLambdaEtaRangeMin
Definition BmnLambdaQa.h:166
BmnLambdaQa::fCurrentEvent
static Int_t fCurrentEvent
Definition BmnLambdaQa.h:213
BmnLambdaQa::fProtonsThetaRangeMin
Int_t fProtonsThetaRangeMin
Definition BmnLambdaQa.h:176
BmnLambdaQa::Exec
virtual void Exec(Option_t *opt)
Derived from FairTask.
Definition BmnLambdaQa.cxx:150
BmnLambdaQa::fDCA12YRangeMax
Int_t fDCA12YRangeMax
Definition BmnLambdaQa.h:197
BmnLambdaQa::fMesonsPRangeMax
Int_t fMesonsPRangeMax
Definition BmnLambdaQa.h:181
BmnLambdaQa::CreateNumberOfReconstrcutedRecMesonsHistograms
void CreateNumberOfReconstrcutedRecMesonsHistograms()
Definition BmnLambdaQa.cxx:849
BmnLambdaQa::fLambdaPRangeMin
Int_t fLambdaPRangeMin
Definition BmnLambdaQa.h:164
BmnLambdaQa::CreateNumberOfReconstructableLambdaHistograms
void CreateNumberOfReconstructableLambdaHistograms()
Definition BmnLambdaQa.cxx:643
BmnLambdaQa::fProtonsThetaRRangeMax
Int_t fProtonsThetaRRangeMax
Definition BmnLambdaQa.h:179
BmnLambdaQa::CreateRecMesonsRecEfficiencyHistograms
void CreateRecMesonsRecEfficiencyHistograms()
Definition BmnLambdaQa.cxx:763
BmnLambdaQa::CreateTwoDimensionalNumberOfReconstructedProtonsHistograms
void CreateTwoDimensionalNumberOfReconstructedProtonsHistograms()
Definition BmnLambdaQa.cxx:841
BmnParticlePair::GetMomPart1
Double_t GetMomPart1()
Definition BmnParticlePair.h:168
BmnParticlePair::GetInvMass
Double_t GetInvMass()
Definition BmnParticlePair.h:154
BmnParticlePair::GetEtaPart2
Double_t GetEtaPart2()
Definition BmnParticlePair.h:174
BmnParticlePair::GetEtaPart1
Double_t GetEtaPart1()
Definition BmnParticlePair.h:172
BmnParticlePair::GetMomPart2
Double_t GetMomPart2()
Definition BmnParticlePair.h:170
BmnParticlePairsInfo::setDCAPart2Min
void setDCAPart2Min(Double_t val)
Definition BmnParticlePairsInfo.cxx:88
BmnParticlePairsInfo::setEtaPart1Max
void setEtaPart1Max(Double_t val)
Definition BmnParticlePairsInfo.cxx:63
BmnParticlePairsInfo::setEtaPart2Max
void setEtaPart2Max(Double_t val)
Definition BmnParticlePairsInfo.cxx:73
BmnParticlePairsInfo::setMomPart2Min
void setMomPart2Min(Double_t val)
Definition BmnParticlePairsInfo.cxx:48
BmnParticlePairsInfo::setDCA12Min
void setDCA12Min(Double_t val)
Definition BmnParticlePairsInfo.cxx:98
BmnParticlePairsInfo::setMomPart1Min
void setMomPart1Min(Double_t val)
Definition BmnParticlePairsInfo.cxx:38
BmnParticlePairsInfo::setPathMin
void setPathMin(Double_t val)
Definition BmnParticlePairsInfo.cxx:108
BmnParticlePairsInfo::setDCAPart1Max
void setDCAPart1Max(Double_t val)
Definition BmnParticlePairsInfo.cxx:83
BmnParticlePairsInfo::setDCAPart1Min
void setDCAPart1Min(Double_t val)
Definition BmnParticlePairsInfo.cxx:78
BmnParticlePairsInfo::setDCA12Max
void setDCA12Max(Double_t val)
Definition BmnParticlePairsInfo.cxx:103
BmnParticlePairsInfo::setPathMax
void setPathMax(Double_t val)
Definition BmnParticlePairsInfo.cxx:113
BmnParticlePairsInfo::setMomPart1Max
void setMomPart1Max(Double_t val)
Definition BmnParticlePairsInfo.cxx:43
BmnParticlePairsInfo::setEtaPart1Min
void setEtaPart1Min(Double_t val)
Definition BmnParticlePairsInfo.cxx:58
BmnParticlePairsInfo::setMomPart2Max
void setMomPart2Max(Double_t val)
Definition BmnParticlePairsInfo.cxx:53
BmnParticlePairsInfo::setEtaPart2Min
void setEtaPart2Min(Double_t val)
Definition BmnParticlePairsInfo.cxx:68
BmnParticlePairsInfo::setDCAPart2Max
void setDCAPart2Max(Double_t val)
Definition BmnParticlePairsInfo.cxx:93
BmnSimulationReport
Base class for simulation reports.
Definition BmnSimulationReport.h:34
BmnSimulationReport::Create
void Create(BmnHistManager *histManager, const string &outputDir)
Main function which creates report data.
Definition BmnSimulationReport.cxx:35
CbmMCTrack
CbmMCTrack.h.
Definition CbmMCTrack.h:26
CbmMCTrack::GetNPoints
Long64_t GetNPoints(DetectorId detId) const
Accessors to the number of MCPoints in the detectors.
Definition CbmMCTrack.cxx:151
CbmMCTrack::GetPz
Double_t GetPz() const
Definition CbmMCTrack.h:60
CbmMCTrack::GetRapidity
Double_t GetRapidity() const
Definition CbmMCTrack.cxx:141
CbmMCTrack::GetMotherId
Int_t GetMotherId() const
Definition CbmMCTrack.h:57
CbmMCTrack::GetP
Double_t GetP() const
Definition CbmMCTrack.h:68
CbmMCTrack::GetPdgCode
Int_t GetPdgCode() const
Accessors.
Definition CbmMCTrack.h:56
abs
RVec< PromoteType< T > > abs(const RVec< T > &v)
TMath::ACos
Double_t ACos(Double_t)
TMath::Log
Double_t Log(Double_t x)
TMath::RadToDeg
constexpr Double_t RadToDeg()