Skip to content

[PWGLF] Correction on histogram filling #15544

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
abmodak merged 4 commits into from
Mar 30, 2026
+5 −6
Merged

[PWGLF] Correction on histogram filling #15544

Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
a07790d
Correction for histograms filling
lucasjsilva Mar 27, 2026
91b669b
Correction on histogram filling
lucasjsilva Mar 27, 2026
fda4e1f
Correction on histogram filling
lucasjsilva Mar 27, 2026
6276ecc
Correction on histogram filling
lucasjsilva Mar 30, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
11 changes: 5 additions & 6 deletions PWGLF/Tasks/GlobalEventProperties/studyPnch.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -256,7 +256,6 @@ struct StudyPnch {
if (!isGenTrackSelected(track)) {
continue;
}
// Verify that the track belongs to the given MC collision
if (track.mcCollisionId() != McCol.globalIndex()) {
continue;
}
Expand All @@ -276,7 +275,6 @@ struct StudyPnch {
if (!isTrackSelected(track)) {
continue;
}
// Verify that the track belongs to the given MC collision
if (track.has_mcParticle()) {
auto particle = track.mcParticle();
if (particle.mcCollisionId() != McCol.mcCollisionId()) {
Expand All @@ -287,7 +285,7 @@ struct StudyPnch {
}
histos.fill(HIST("PhiVsEtaHist"), track.phi(), track.eta());
}
// Once all the frequencies have been counted, a loop can be made to fill the histogram
// Loop to fill the histogram without cloned tracks
for (const auto& [globalIndex, frequency] : recoFrequencies) {
histos.fill(HIST("hCountNTracks"), frequency);
// Fill histogram with not cloned tracks
Expand Down Expand Up @@ -355,18 +353,19 @@ struct StudyPnch {
if (!isEventSelected(RecCol)) {
continue;
}
// Evaluation of reconstructed collisions with more than 1 contributor
if (RecCol.globalIndex() != mcCollision.bestCollisionIndex()) {
continue;
}
auto recTracksPart = RecTracks.sliceBy(perCollision, RecCol.globalIndex());
auto multrec = countNTracksMcCol(recTracksPart, RecCol);
if (multrec > 0) {
histos.fill(HIST("hMultiplicityMCrec"), multrec);
histos.fill(HIST("hMultiplicityMCrecMod"), multrec);
}
auto multgen = countGenTracks(GenParticles, mcCollision);
if (multgen > 0 && multrec > 0) {
histos.fill(HIST("hMultiplicityMCgen"), multgen);
histos.fill(HIST("hResponseMatrix"), multrec, multgen);
histos.fill(HIST("hMultiplicityMCgenMod"), multgen);
histos.fill(HIST("hResponseMatrixMod"), multrec, multgen);
}
}
}
Expand Down
Loading