Skip to content

Commit 66f048b

Browse files
altsybeealtsybee
authored
tune histograms, fixing issues
1 parent 5f48d15 commit 66f048b

File tree

1 file changed

+103
-42
lines changed

1 file changed

+103
-42
lines changed

‎DPG/Tasks/AOTEvent/lightIonsEvSelQa.cxx‎

Lines changed: 103 additions & 42 deletions
Original file line numberDiff line numberDiff line change
@@ -133,21 +133,19 @@ struct LightIonsEvSelQa {
133133
const AxisSpec axisMultFT0C{nBinsMultFwd, 0., static_cast<float>(nMaxMultFwd * 0.15), "mult FT0C"};
134134
const AxisSpec axisMultT0M{nBinsMultFwd * 2, 0., static_cast<float>(nMaxMultFwd * 0.4), "mult FT0M"};
135135

136-
const AxisSpec axisMultT0MlargeBins{nBinsMultFwd, 0., static_cast<float>(nMaxMultFwd * 0.75), "mult FT0M"};
137-
histos.add("multT0M_vs_multSumFT0", ";multT0M;multT0M_uncorrected", kTH2F, {axisMultT0MlargeBins, axisMultT0MlargeBins});
138-
139136
const AxisSpec axisVtxZ{800, -20., 20., ""};
140-
const AxisSpec axisBcDiff{600, -300., 300., "bc difference"};
137+
const AxisSpec axisBcDiff{601, -300.5, 300.5, "bc difference"};
141138

142139
const AxisSpec axisNcontrib{601, -0.5, 600.5, "n contributors"};
143140
const AxisSpec axisColTimeRes{1500, 0., 1500., "collision time resolution (ns)"};
144141

145142
AxisSpec axisVertexChi2{100, 0, 500, "Chi2 of vertex fit"};
146143
AxisSpec axisVertexChi2perContrib{100, 0, 10, "Chi2 of vertex fit"};
147144

148-
const AxisSpec axisTimeDiff{100, -10., 10., ""};
149-
const AxisSpec axisTimeSum{100, -10., 10., ""};
150-
const AxisSpec axisZNampl{200, 0, 3000, ""};
145+
const AxisSpec axisTimeZN{800, -20., 20., ""};
146+
const AxisSpec axisTimeDiff{150, -10., 10., ""};
147+
const AxisSpec axisTimeSum{150, -10., 10., ""};
148+
const AxisSpec axisZNampl{200, 0, 10000, ""};
151149

152150
histos.add("noSpecSelections/hBcColNoSel8", "", kTH1F, {axisBCs});
153151
histos.add("noSpecSelections/hBcOrigNoSel8", "", kTH1F, {axisBCs});
@@ -171,7 +169,10 @@ struct LightIonsEvSelQa {
171169
histos.add("noSpecSelections/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
172170
histos.add("noSpecSelections/hVertexChi2VsNcontrib", "", kTH2F, {axisNcontrib, axisVertexChi2perContrib});
173171
histos.add("noSpecSelections/hNPVvsNch", "", kTH2F, {axisNcontrib, axisNcontrib});
174-
histos.add("noSpecSelections/hTimeZAC", "ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});
172+
histos.add("noSpecSelections/hTimeZN_AC_sum_vs_diff", ";ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});
173+
histos.add("noSpecSelections/hTimeZN_A_vs_C", ";ZNA time (ns); ZNC time (ns)", kTH2F, {axisTimeZN, axisTimeZN});
174+
histos.add("noSpecSelections/hTimeZNA", ";ZNA time (ns)", kTH1F, {axisTimeZN});
175+
histos.add("noSpecSelections/hTimeZNC", ";ZNC time (ns)", kTH1F, {axisTimeZN});
175176

176177
histos.add("noPU/hBcColNoSel8", "", kTH1F, {axisBCs});
177178
histos.add("noPU/hBcOrigNoSel8", "", kTH1F, {axisBCs});
@@ -194,8 +195,11 @@ struct LightIonsEvSelQa {
194195
histos.add("noPU/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
195196
histos.add("noPU/hVertexChi2VsNcontrib", "", kTH2F, {axisNcontrib, axisVertexChi2perContrib});
196197
histos.add("noPU/hNPVvsNch", "", kTH2F, {axisNcontrib, axisNcontrib});
197-
histos.add("noPU/hTimeZAC", "ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});
198-
histos.add("noPU/hAmplZAC", "ZNC vs ZNA", kTH2F, {axisZNampl, axisZNampl});
198+
histos.add("noPU/hTimeZN_AC_sum_vs_diff", ";ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});
199+
histos.add("noPU/hTimeZN_A_vs_C", ";ZNA time (ns); ZNC time (ns)", kTH2F, {axisTimeZN, axisTimeZN});
200+
histos.add("noPU/hTimeZNA", ";ZNA time (ns)", kTH1F, {axisTimeZN});
201+
histos.add("noPU/hTimeZNC", ";ZNC time (ns)", kTH1F, {axisTimeZN});
202+
histos.add("noPU/hAmplZNAC", "ZNC vs ZNA", kTH2F, {axisZNampl, axisZNampl});
199203

200204
histos.add("noPU_pvTOFmatched/hBcColNoSel8", "", kTH1F, {axisBCs});
201205
histos.add("noPU_pvTOFmatched/hBcTVX", "", kTH1F, {axisBCs});
@@ -403,7 +407,10 @@ struct LightIonsEvSelQa {
403407
histos.add("noPU_lowMultCut/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
404408
histos.add("noPU_lowMultCut/hVertexChi2VsNcontrib", "", kTH2F, {axisNcontrib, axisVertexChi2perContrib});
405409
histos.add("noPU_lowMultCut/hNPVvsNch", "", kTH2F, {axisNcontrib, axisNcontrib});
406-
histos.add("noPU_lowMultCut/hTimeZAC", "ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});
410+
histos.add("noPU_lowMultCut/hTimeZN_AC_sum_vs_diff", ";ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});
411+
histos.add("noPU_lowMultCut/hTimeZN_A_vs_C", ";ZNA time (ns); ZNC time (ns)", kTH2F, {axisTimeZN, axisTimeZN});
412+
histos.add("noPU_lowMultCut/hTimeZNA", ";ZNA time (ns)", kTH1F, {axisTimeZN});
413+
histos.add("noPU_lowMultCut/hTimeZNC", ";ZNC time (ns)", kTH1F, {axisTimeZN});
407414

408415
histos.add("highMultCloudCut/hBcColNoSel8", "", kTH1F, {axisBCs});
409416
histos.add("highMultCloudCut/hBcOrigNoSel8", "", kTH1F, {axisBCs});
@@ -430,7 +437,10 @@ struct LightIonsEvSelQa {
430437
histos.add("noPU_highMultCloudCut/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
431438
histos.add("noPU_highMultCloudCut/hVertexChi2VsNcontrib", "", kTH2F, {axisNcontrib, axisVertexChi2perContrib});
432439
histos.add("noPU_highMultCloudCut/hNPVvsNch", "", kTH2F, {axisNcontrib, axisNcontrib});
433-
histos.add("noPU_highMultCloudCut/hTimeZAC", "ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});
440+
histos.add("noPU_highMultCloudCut/hTimeZN_AC_sum_vs_diff", ";ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});
441+
histos.add("noPU_highMultCloudCut/hTimeZN_A_vs_C", ";ZNA time (ns); ZNC time (ns)", kTH2F, {axisTimeZN, axisTimeZN});
442+
histos.add("noPU_highMultCloudCut/hTimeZNA", ";ZNA time (ns)", kTH1F, {axisTimeZN});
443+
histos.add("noPU_highMultCloudCut/hTimeZNC", ";ZNC time (ns)", kTH1F, {axisTimeZN});
434444

435445
histos.add("noPU_badVzDiff/hBcColNoSel8", "", kTH1F, {axisBCs});
436446
histos.add("noPU_badVzDiff/hBcTVX", "", kTH1F, {axisBCs});
@@ -570,7 +580,10 @@ struct LightIonsEvSelQa {
570580
histos.add("noPU_cutByVzDiff_pvTOF/hColTimeResVsNcontrib", "", kTH2F, {axisNcontrib, axisColTimeRes});
571581
histos.add("noPU_cutByVzDiff_pvTOF/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
572582
histos.add("noPU_cutByVzDiff_pvTOF/hNPVvsNch", "", kTH2F, {axisNcontrib, axisNcontrib});
573-
histos.add("noPU_cutByVzDiff_pvTOF/hTimeZAC", "ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});
583+
histos.add("noPU_cutByVzDiff_pvTOF/hTimeZN_AC_sum_vs_diff", ";ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});
584+
histos.add("noPU_cutByVzDiff_pvTOF/hTimeZN_A_vs_C", ";ZNA time (ns); ZNC time (ns)", kTH2F, {axisTimeZN, axisTimeZN});
585+
histos.add("noPU_cutByVzDiff_pvTOF/hTimeZNA", ";ZNA time (ns)", kTH1F, {axisTimeZN});
586+
histos.add("noPU_cutByVzDiff_pvTOF/hTimeZNC", ";ZNC time (ns)", kTH1F, {axisTimeZN});
574587

575588
histos.add("noPU_cutByVzDiff_noFT0activityNearby/hBcColNoSel8", "", kTH1F, {axisBCs});
576589
histos.add("noPU_cutByVzDiff_noFT0activityNearby/hBcTVX", "", kTH1F, {axisBCs});
@@ -606,18 +619,21 @@ struct LightIonsEvSelQa {
606619
histos.add("noPU_CutOnZNACtime/hColTimeResVsNcontrib", "", kTH2F, {axisNcontrib, axisColTimeRes});
607620
histos.add("noPU_CutOnZNACtime/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
608621
histos.add("noPU_CutOnZNACtime/hNPVvsNch", "", kTH2F, {axisNcontrib, axisNcontrib});
609-
histos.add("noPU_CutOnZNACtime/hTimeZAC", "ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});
610-
histos.add("noPU_CutOnZNACtime/hAmplZAC", "ZNC vs ZNA", kTH2F, {axisZNampl, axisZNampl});
622+
histos.add("noPU_CutOnZNACtime/hTimeZN_AC_sum_vs_diff", ";ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});
623+
histos.add("noPU_CutOnZNACtime/hAmplZNAC", "ZNC vs ZNA", kTH2F, {axisZNampl, axisZNampl});
624+
histos.add("noPU_CutOnZNACtime/hTimeZN_A_vs_C", ";ZNA time (ns); ZNC time (ns)", kTH2F, {axisTimeZN, axisTimeZN});
625+
histos.add("noPU_CutOnZNACtime/hTimeZNA", ";ZNA time (ns)", kTH1F, {axisTimeZN});
626+
histos.add("noPU_CutOnZNACtime/hTimeZNC", ";ZNC time (ns)", kTH1F, {axisTimeZN});
611627

612628
histos.add("noPU_AntiCutOnZNACtime/hBcFV0", "", kTH1F, {axisBCs});
613629
histos.add("noPU_AntiCutOnZNACtime/nTracksPV_vs_V0A", "", kTH2F, {axisMultV0A, axisNtracks});
614630
histos.add("noPU_AntiCutOnZNACtime/hVtxFT0MinusVtxColVsMultT0M", "", kTH2F, {axisVtxZ, axisMultT0M});
615-
histos.add("noPU_AntiCutOnZNACtime/hAmplZAC", "ZNC vs ZNA", kTH2F, {axisZNampl, axisZNampl});
631+
histos.add("noPU_AntiCutOnZNACtime/hAmplZNAC", "ZNC vs ZNA", kTH2F, {axisZNampl, axisZNampl});
616632

617633
histos.add("noPU_AntiCutOnZNAampl/hBcFV0", "", kTH1F, {axisBCs});
618634
histos.add("noPU_AntiCutOnZNAampl/hVtxFT0MinusVtxColVsMultT0M", "", kTH2F, {axisVtxZ, axisMultT0M});
619635
histos.add("noPU_AntiCutOnZNAampl/nTracksPV_vs_V0A", "", kTH2F, {axisMultV0A, axisNtracks});
620-
histos.add("noPU_AntiCutOnZNAampl/hAmplZAC", "ZNC vs ZNA", kTH2F, {axisZNampl, axisZNampl});
636+
histos.add("noPU_AntiCutOnZNAampl/hAmplZNAC", "ZNC vs ZNA", kTH2F, {axisZNampl, axisZNampl});
621637

622638
//
623639
histos.add("hNcontribColFromData", "", kTH1F, {axisNcontrib});
@@ -629,8 +645,12 @@ struct LightIonsEvSelQa {
629645
histos.add("MC/hMCdataBcDiffVsMult", "", kTH2F, {axisNcontrib, axisBcDiff});
630646
histos.add("MC/hMCdataFoundBcDiffVsMult", "", kTH2F, {axisNcontrib, axisBcDiff});
631647
histos.add("MCsel8/hMCdataVzDiff", "", kTH2F, {axisNcontrib, axisVtxZdiff});
648+
histos.add("MCsel8/hMCdataVzDiff_vertTRDmatched", "", kTH2F, {axisNcontrib, axisVtxZdiff});
649+
histos.add("MCsel8/hMCdataVzDiff_vertTOFmatched", "", kTH2F, {axisNcontrib, axisVtxZdiff});
632650
histos.add("MCsel8/hMCdataBcDiffVsMult", "", kTH2F, {axisNcontrib, axisBcDiff});
633651
histos.add("MCsel8/hMCdataFoundBcDiffVsMult", "", kTH2F, {axisNcontrib, axisBcDiff});
652+
histos.add("MCsel8/hMCdataFoundBcDiffVsMult_vertTRDmatched", "", kTH2F, {axisNcontrib, axisBcDiff});
653+
histos.add("MCsel8/hMCdataFoundBcDiffVsMult_vertTOFmatched", "", kTH2F, {axisNcontrib, axisBcDiff});
634654
}
635655

636656
Preslice<FullTracksIU> perCollision = aod::track::collisionId;
@@ -947,20 +967,13 @@ struct LightIonsEvSelQa {
947967
float multT0A = hasFT0 ? col.multFT0A() : 0;
948968
float multT0C = hasFT0 ? col.multFT0C() : 0;
949969
float multT0M = multT0A + multT0C;
950-
if (hasFT0) {
951-
float multT0A_uncorr = foundBC.ft0().sumAmpA();
952-
float multT0C_uncorr = foundBC.ft0().sumAmpC();
953-
float multT0M_uncorr = multT0A_uncorr + multT0C_uncorr;
954-
955-
histos.fill(HIST("multT0M_vs_multSumFT0"), multT0M, multT0M_uncorr);
956-
}
957970

958971
float timeZNA = foundBC.has_zdc() ? foundBC.zdc().timeZNA() : -999.f;
959972
float timeZNC = foundBC.has_zdc() ? foundBC.zdc().timeZNC() : -999.f;
960973
float znSum = timeZNA + timeZNC;
961974
float znDiff = timeZNA - timeZNC;
962975
// bool goodZNACtime = fabs(znSum) < 2 && fabs(znDiff) < 2;
963-
bool goodZNACtime = (timeZNA > -6 && timeZNA < 4) && (timeZNC > -6 && timeZNC < 4);
976+
bool goodZNACtime = (timeZNA > -5 && timeZNA < 2) && (timeZNC > -5 && timeZNC < 2);
964977

965978
float multZNA = foundBC.has_zdc() ? foundBC.zdc().energyCommonZNA() : -999;
966979
float multZNC = foundBC.has_zdc() ? foundBC.zdc().energyCommonZNC() : -999;
@@ -1030,14 +1043,20 @@ struct LightIonsEvSelQa {
10301043
histos.fill(HIST("noSpecSelections/hBcColNoSel8"), localBC);
10311044
histos.fill(HIST("noSpecSelections/hBcOrigNoSel8"), bcOriginal);
10321045
histos.fill(HIST("noSpecSelections/Vz"), vZ);
1033-
histos.fill(HIST("noSpecSelections/hTimeZAC"), znDiff, znSum);
1046+
histos.fill(HIST("noSpecSelections/hTimeZN_AC_sum_vs_diff"), znDiff, znSum);
1047+
histos.fill(HIST("noSpecSelections/hTimeZN_A_vs_C"), timeZNA, timeZNC);
1048+
histos.fill(HIST("noSpecSelections/hTimeZNA"), timeZNA);
1049+
histos.fill(HIST("noSpecSelections/hTimeZNC"), timeZNC);
10341050

10351051
if (noPU) {
10361052
histos.fill(HIST("noPU/hBcColNoSel8"), localBC);
10371053
histos.fill(HIST("noPU/hBcOrigNoSel8"), bcOriginal);
10381054
histos.fill(HIST("noPU/Vz"), vZ);
1039-
histos.fill(HIST("noPU/hTimeZAC"), znDiff, znSum);
1040-
histos.fill(HIST("noPU/hAmplZAC"), multZNA, multZNC);
1055+
histos.fill(HIST("noPU/hTimeZN_AC_sum_vs_diff"), znDiff, znSum);
1056+
histos.fill(HIST("noPU/hTimeZN_A_vs_C"), timeZNA, timeZNC);
1057+
histos.fill(HIST("noPU/hAmplZNAC"), multZNA, multZNC);
1058+
histos.fill(HIST("noPU/hTimeZNA"), timeZNA);
1059+
histos.fill(HIST("noPU/hTimeZNC"), timeZNC);
10411060
}
10421061
if (noPU && pvTOFmatched) {
10431062
histos.fill(HIST("noPU_pvTOFmatched/hBcColNoSel8"), localBC);
@@ -1093,7 +1112,10 @@ struct LightIonsEvSelQa {
10931112
if (noPU && underLine) {
10941113
histos.fill(HIST("noPU_lowMultCut/hBcColNoSel8"), localBC);
10951114
histos.fill(HIST("noPU_lowMultCut/hBcOrigNoSel8"), bcOriginal);
1096-
histos.fill(HIST("noPU_lowMultCut/hTimeZAC"), znDiff, znSum);
1115+
histos.fill(HIST("noPU_lowMultCut/hTimeZN_AC_sum_vs_diff"), znDiff, znSum);
1116+
histos.fill(HIST("noPU_lowMultCut/hTimeZN_A_vs_C"), timeZNA, timeZNC);
1117+
histos.fill(HIST("noPU_lowMultCut/hTimeZNA"), timeZNA);
1118+
histos.fill(HIST("noPU_lowMultCut/hTimeZNC"), timeZNC);
10971119
}
10981120
if (grassOnTheRight) {
10991121
histos.fill(HIST("highMultCloudCut/hBcColNoSel8"), localBC);
@@ -1102,11 +1124,17 @@ struct LightIonsEvSelQa {
11021124
if (noPU && grassOnTheRight) {
11031125
histos.fill(HIST("noPU_highMultCloudCut/hBcColNoSel8"), localBC);
11041126
histos.fill(HIST("noPU_highMultCloudCut/hBcOrigNoSel8"), bcOriginal);
1105-
histos.fill(HIST("noPU_highMultCloudCut/hTimeZAC"), znDiff, znSum);
1127+
histos.fill(HIST("noPU_highMultCloudCut/hTimeZN_AC_sum_vs_diff"), znDiff, znSum);
1128+
histos.fill(HIST("noPU_highMultCloudCut/hTimeZN_A_vs_C"), timeZNA, timeZNC);
1129+
histos.fill(HIST("noPU_highMultCloudCut/hTimeZNA"), timeZNA);
1130+
histos.fill(HIST("noPU_highMultCloudCut/hTimeZNC"), timeZNC);
11061131
}
11071132
if (noPU && !badVzDiff && pvTOFmatched) { // noPileup_cutByVzDiff_pvTOF_noFT0act
11081133
histos.fill(HIST("noPU_cutByVzDiff_pvTOF/hBcColNoSel8"), localBC);
1109-
histos.fill(HIST("noPU_cutByVzDiff_pvTOF/hTimeZAC"), znDiff, znSum);
1134+
histos.fill(HIST("noPU_cutByVzDiff_pvTOF/hTimeZN_AC_sum_vs_diff"), znDiff, znSum);
1135+
histos.fill(HIST("noPU_cutByVzDiff_pvTOF/hTimeZN_A_vs_C"), timeZNA, timeZNC);
1136+
histos.fill(HIST("noPU_cutByVzDiff_pvTOF/hTimeZNA"), timeZNA);
1137+
histos.fill(HIST("noPU_cutByVzDiff_pvTOF/hTimeZNC"), timeZNC);
11101138
}
11111139

11121140
// only here cut on sel8:
@@ -1298,7 +1326,10 @@ struct LightIonsEvSelQa {
12981326
histos.fill(HIST("noPU_CutOnZNACtime/hColBcDiffVsNcontrib"), nContributors, bcToClosestTVXdiff);
12991327
histos.fill(HIST("noPU_CutOnZNACtime/hColTimeResVsNcontrib"), nContributors, timeRes);
13001328
histos.fill(HIST("noPU_CutOnZNACtime/hNPVvsNch"), nPVtracks, nGlobalTracksAll);
1301-
histos.fill(HIST("noPU_CutOnZNACtime/hTimeZAC"), znDiff, znSum);
1329+
histos.fill(HIST("noPU_CutOnZNACtime/hTimeZN_AC_sum_vs_diff"), znDiff, znSum);
1330+
histos.fill(HIST("noPU_CutOnZNACtime/hTimeZN_A_vs_C"), timeZNA, timeZNC);
1331+
histos.fill(HIST("noPU_CutOnZNACtime/hTimeZNA"), timeZNA);
1332+
histos.fill(HIST("noPU_CutOnZNACtime/hTimeZNC"), timeZNC);
13021333
}
13031334

13041335
if (foundBC.has_ft0()) {
@@ -1499,7 +1530,7 @@ struct LightIonsEvSelQa {
14991530
histos.fill(HIST("noPU_CutOnZNACtime/hVtxFT0MinusVtxColVsMultT0M"), diffVz, multT0A + multT0C);
15001531
histos.fill(HIST("noPU_CutOnZNACtime/nTracksPV_vs_T0A"), multT0A, nPVtracks);
15011532
histos.fill(HIST("noPU_CutOnZNACtime/nTracksPV_vs_T0C"), multT0C, nPVtracks);
1502-
histos.fill(HIST("noPU_CutOnZNACtime/hAmplZAC"), multZNA, multZNC);
1533+
histos.fill(HIST("noPU_CutOnZNACtime/hAmplZNAC"), multZNA, multZNC);
15031534
}
15041535
}
15051536

@@ -1627,13 +1658,13 @@ struct LightIonsEvSelQa {
16271658
histos.fill(HIST("noPU_AntiCutOnZNACtime/hBcFV0"), localBC);
16281659
histos.fill(HIST("noPU_AntiCutOnZNACtime/nTracksPV_vs_V0A"), multV0A, nPVtracks);
16291660
histos.fill(HIST("noPU_AntiCutOnZNACtime/hVtxFT0MinusVtxColVsMultT0M"), diffVz, multT0A + multT0C);
1630-
histos.fill(HIST("noPU_AntiCutOnZNACtime/hAmplZAC"), multZNA, multZNC);
1661+
histos.fill(HIST("noPU_AntiCutOnZNACtime/hAmplZNAC"), multZNA, multZNC);
16311662
}
16321663
if (noPU && !cutZNACampl) {
16331664
histos.fill(HIST("noPU_AntiCutOnZNAampl/hBcFV0"), localBC);
16341665
histos.fill(HIST("noPU_AntiCutOnZNAampl/nTracksPV_vs_V0A"), multV0A, nPVtracks);
16351666
histos.fill(HIST("noPU_AntiCutOnZNAampl/hVtxFT0MinusVtxColVsMultT0M"), diffVz, multT0A + multT0C);
1636-
histos.fill(HIST("noPU_AntiCutOnZNAampl/hAmplZAC"), multZNA, multZNC);
1667+
histos.fill(HIST("noPU_AntiCutOnZNAampl/hAmplZNAC"), multZNA, multZNC);
16371668
}
16381669
}
16391670
if (foundBC.has_zdc()) {
@@ -1778,24 +1809,54 @@ struct LightIonsEvSelQa {
17781809
aod::McCollisions const&)
17791810
{
17801811
for (const auto& col : collisions) {
1812+
if (fabs(col.posZ()) > 10)
1813+
continue;
17811814
bool isSel8 = col.sel8();
17821815
if (col.has_mcCollision()) {
17831816
const auto mcCollision = col.mcCollision();
17841817
LOGP(debug, "col.posZ() = {}, mcCollision.posZ() = {}", col.posZ(), mcCollision.posZ());
1785-
histos.fill(HIST("MC/hMCdataVzDiff"), col.numContrib(), col.posZ() - mcCollision.posZ());
1818+
1819+
float diffVz = col.posZ() - mcCollision.posZ();
1820+
histos.fill(HIST("MC/hMCdataVzDiff"), col.numContrib(), diffVz);
17861821

17871822
auto bc = col.bc_as<BCsInfo>();
17881823
auto mcBc = mcCollision.bc_as<BCsInfo>();
17891824
auto foundBC = col.foundBC_as<BCsInfo>();
17901825

1791-
LOGP(debug, "col.bc().globalBC() = {}, mcCollision.bc().globalBC() = {}", col.bc().globalBC(), mcCollision.bc().globalBC());
1792-
histos.fill(HIST("MC/hMCdataBcDiffVsMult"), col.numContrib(), bc.globalBC() - mcBc.globalBC());
1793-
histos.fill(HIST("MC/hMCdataFoundBcDiffVsMult"), col.numContrib(), foundBC.globalBC() - mcBc.globalBC());
1826+
uint64_t globalBC = bc.globalBC();
1827+
uint64_t foundGlobalBC = foundBC.globalBC();
1828+
uint64_t mcGlobalBC = mcBc.globalBC();
1829+
1830+
int bcDiff = static_cast<int>(globalBC - mcGlobalBC);
1831+
int foundBcDiff = static_cast<int>(foundGlobalBC - mcGlobalBC);
1832+
1833+
// restrict bc diff range - to see values in the diff histograms
1834+
if (bcDiff > 300)
1835+
bcDiff = 300;
1836+
if (bcDiff < -300)
1837+
bcDiff = -300;
1838+
if (foundBcDiff > 300)
1839+
foundBcDiff = 300;
1840+
if (foundBcDiff < -300)
1841+
foundBcDiff = -300;
1842+
1843+
LOGP(debug, "globalBC = {}, mcGlobalBC = {}", globalBC, mcGlobalBC);
1844+
histos.fill(HIST("MC/hMCdataBcDiffVsMult"), col.numContrib(), bcDiff);
1845+
histos.fill(HIST("MC/hMCdataFoundBcDiffVsMult"), col.numContrib(), foundBcDiff);
17941846

17951847
if (isSel8) {
1796-
histos.fill(HIST("MCsel8/hMCdataVzDiff"), col.numContrib(), col.posZ() - mcCollision.posZ());
1797-
histos.fill(HIST("MCsel8/hMCdataBcDiffVsMult"), col.numContrib(), bc.globalBC() - mcBc.globalBC());
1798-
histos.fill(HIST("MCsel8/hMCdataFoundBcDiffVsMult"), col.numContrib(), foundBC.globalBC() - mcBc.globalBC());
1848+
histos.fill(HIST("MCsel8/hMCdataVzDiff"), col.numContrib(), diffVz);
1849+
histos.fill(HIST("MCsel8/hMCdataBcDiffVsMult"), col.numContrib(), bcDiff);
1850+
histos.fill(HIST("MCsel8/hMCdataFoundBcDiffVsMult"), col.numContrib(), foundBcDiff);
1851+
1852+
if (col.selection_bit(kIsVertexTRDmatched)) {
1853+
histos.fill(HIST("MCsel8/hMCdataVzDiff_vertTRDmatched"), col.numContrib(), diffVz);
1854+
histos.fill(HIST("MCsel8/hMCdataFoundBcDiffVsMult_vertTRDmatched"), col.numContrib(), foundBcDiff);
1855+
}
1856+
if (col.selection_bit(kIsVertexTOFmatched)) {
1857+
histos.fill(HIST("MCsel8/hMCdataVzDiff_vertTOFmatched"), col.numContrib(), diffVz);
1858+
histos.fill(HIST("MCsel8/hMCdataFoundBcDiffVsMult_vertTOFmatched"), col.numContrib(), foundBcDiff);
1859+
}
17991860
}
18001861
}
18011862
}

0 commit comments

Comments
 (0)