Skip to content

Commit 2494835

Browse files
authored
[DPG] Update eventSelectionQa.cxx - add dEdx vs occupancy QA histo (#10819)
1 parent c81e279 commit 2494835

File tree

1 file changed

+15
-0
lines changed

1 file changed

+15
-0
lines changed

DPG/Tasks/AOTEvent/eventSelectionQa.cxx

Lines changed: 15 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -356,6 +356,8 @@ struct EventSelectionQaTask {
356356
histos.add("occupancyQA/hITSTracks_ev1_vs_ev2_2coll_in_ROF", ";nITStracks event #1;nITStracks event #2", kTH2D, {{200, 0., 6000}, {200, 0., 6000}});
357357
histos.add("occupancyQA/hITSTracks_ev1_vs_ev2_2coll_in_ROF_UPC", ";nITStracks event #1;nITStracks event #2", kTH2D, {{41, -0.5, 40.5}, {41, -0.5, 40.5}});
358358
histos.add("occupancyQA/hITSTracks_ev1_vs_ev2_2coll_in_ROF_nonUPC", ";nITStracks event #1;nITStracks event #2", kTH2D, {{200, 0., 6000}, {200, 0., 6000}});
359+
360+
histos.add("occupancyQA/dEdx_vs_centr_vs_occup_narrow_p_win", "dE/dx", kTH3F, {{20, 0, 4000, "n PV tracks"}, {60, 0, 15000, "occupancy"}, {800, 0.0, 800.0, "dE/dx (a. u.)"}});
359361
}
360362
}
361363

@@ -1152,6 +1154,19 @@ struct EventSelectionQaTask {
11521154
histos.fill(HIST("occupancyQA/hOccupancyByFT0CvsByTracks"), occupancyByTracks, occupancyByFT0C);
11531155
histos.fill(HIST("occupancyQA/hNumTracksPV_vs_V0A_vs_occupancy"), multV0A, nPV, occupancyByTracks);
11541156
histos.fill(HIST("occupancyQA/hNumTracksPVTPC_vs_V0A_vs_occupancy"), multV0A, nContributorsAfterEtaTPCCuts, occupancyByTracks);
1157+
1158+
// dE/dx QA for a narrow pT bin
1159+
for (const auto& track : tracksGrouped) {
1160+
if (!track.isPVContributor())
1161+
continue;
1162+
if (std::fabs(track.eta()) < 0.8 && track.pt() > 0.2 && track.itsNCls() >= 5) {
1163+
float signedP = track.sign() * track.tpcInnerParam();
1164+
if (std::fabs(signedP) > 0.38 && std::fabs(signedP) < 0.4 && track.tpcNClsFound() > 70 && track.tpcNClsCrossedRows() > 80 && track.itsChi2NCl() < 36 && track.tpcChi2NCl() < 4) {
1165+
float dEdx = track.tpcSignal();
1166+
histos.fill(HIST("occupancyQA/dEdx_vs_centr_vs_occup_narrow_p_win"), nPV, occupancyByTracks, dEdx);
1167+
}
1168+
}
1169+
}
11551170
}
11561171
}
11571172

0 commit comments

Comments
 (0)