Skip to content

[PWGLF] Inclusion of process function for Monte Carlo process avoiding over count of generated tracks #14722

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

Open
lucasjsilva wants to merge 7 commits into
base: master
Choose a base branch
from
Open

[PWGLF] Inclusion of process function for Monte Carlo process avoiding over count of generated tracks #14722

Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
88a67a0
Added treatment for outliers in MC tracks
lucasjsilva Jan 15, 2026
e8078a4
Merge branch 'AliceO2Group:master' into master
lucasjsilva Jan 28, 2026
ea902fd
Treatment for outliers in generated MC tracks
lucasjsilva Jan 28, 2026
32a4ca5
Process Monte Carlo with correct count of generated tracks
lucasjsilva Jan 30, 2026
8c1debc
Please consider the following formatting changes
alibuild Jan 30, 2026
04e1f91
Merge pull request #2 from alibuild/alibot-cleanup-14722
lucasjsilva Jan 30, 2026
4e3f129
Merge branch 'master' into processMCOutlier
lucasjsilva Jan 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
29 changes: 27 additions & 2 deletions PWGLF/Tasks/GlobalEventProperties/studyPnch.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -125,7 +125,7 @@ struct StudyPnch {
x->SetBinLabel(6, "INEL > 0");
x->SetBinLabel(7, "|vz| < 10");

if (doprocessData || doprocessCorrelation || doprocessMonteCarlo) {
if (doprocessData || doprocessCorrelation || doprocessMonteCarlo || doprocessTreatedMonteCarlo) {
histos.add("PhiVsEtaHist", "PhiVsEtaHist", kTH2F, {axisPhi, axisEta}, false);
}
if (doprocessData) {
Expand All @@ -143,6 +143,11 @@ struct StudyPnch {
histos.add("hMultiplicityMCgen", "hMultiplicityMCgen", kTH1F, {axisMult}, true);
histos.add("hResponseMatrix", "hResponseMatrix", kTH2F, {axisMult, axisMult}, true);
}
if (doprocessTreatedMonteCarlo) {
histos.add("hMultiplicityTreatMCrec", "hMultiplicityTreatMCrec", kTH1F, {axisMult}, true);
histos.add("hMultiplicityTreatMCgen", "hMultiplicityTreatMCgen", kTH1F, {axisMult}, true);
histos.add("hResponseMatrixTreat", "hResponseMatrixTreat", kTH2F, {axisMult, axisMult}, true);
}
if (doprocessEvtLossSigLossMC) {
histos.add("MCEventHist", "MCEventHist", kTH1F, {axisEvent}, false);
auto hstat = histos.get<TH1>(HIST("MCEventHist"));
Expand Down Expand Up @@ -264,7 +269,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();
Expand Down Expand Up @@ -322,6 +326,26 @@ struct StudyPnch {
}
}

void processTreatedMonteCarlo(ColMCTrueTable::iterator const& mcCollision, ColMCRecTable const& RecCols, TrackMCTrueTable const& GenParticles, FilTrackMCRecTable const& RecTracks)
{
// Count generated tracks at each iterator
auto multgen = countGenTracks(GenParticles, mcCollision);
histos.fill(HIST("hMultiplicityTreatMCgen"), multgen);
for (const auto& RecCol : RecCols) {
if (!isEventSelected(RecCol)) {
continue;
}
// Verify that the reconstructed collision corresponds to the given MC collision
if (RecCol.mcCollisionId() != mcCollision.globalIndex()) {
continue;
}
auto recTracksPart = RecTracks.sliceBy(perCollision, RecCol.globalIndex());
auto multrec = countNTracksMcCol(recTracksPart, RecCol);
histos.fill(HIST("hMultiplicityTreatMCrec"), multrec);
histos.fill(HIST("hResponseMatrixTreat"), multrec, multgen);
}
}

void processEvtLossSigLossMC(soa::Join<ColMCTrueTable, aod::MultMCExtras>::iterator const& mcCollision, ColMCRecTable const& RecCols, TrackMCTrueTable const& GenParticles)
{
if (isApplyInelgt0 && !mcCollision.isInelGt0()) {
Expand Down Expand Up @@ -359,6 +383,7 @@ struct StudyPnch {
PROCESS_SWITCH(StudyPnch, processData, "process data CentFT0C", false);
PROCESS_SWITCH(StudyPnch, processCorrelation, "do correlation study in data", false);
PROCESS_SWITCH(StudyPnch, processMonteCarlo, "process MC CentFT0C", false);
PROCESS_SWITCH(StudyPnch, processTreatedMonteCarlo, "process Treated MC CentFT0C", false);
PROCESS_SWITCH(StudyPnch, processEvtLossSigLossMC, "process Signal Loss, Event Loss", false);
};

Expand Down
Loading