Skip to content

[PWGLF] Add histograms for closure test #15287

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
alcaliva merged 4 commits into from
Mar 7, 2026
+85 −0
Merged

[PWGLF] Add histograms for closure test #15287

Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
1d03dcf
[PWGLF] Add histograms for closure test
alcaliva Mar 6, 2026
128af53
fix var
alcaliva Mar 6, 2026
21d1681
Please consider the following formatting changes
alibuild Mar 6, 2026
eac851d
Merge pull request #99 from alibuild/alibot-cleanup-15287
alcaliva Mar 6, 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
85 changes: 85 additions & 0 deletions PWGLF/Tasks/Nuspex/antinucleiInJets.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -137,6 +137,10 @@ struct JetMatching {

struct AntinucleiInJets {

// Random engine
std::mt19937 rng;
std::uniform_int_distribution<int> generateRandomNr{0, 1};

// Histogram registries for data, MC, quality control, multiplicity and correlations
HistogramRegistry registryData{"registryData", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};
HistogramRegistry registryMC{"registryMC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};
Expand Down Expand Up @@ -276,6 +280,9 @@ struct AntinucleiInJets {
auto timeSeed = std::chrono::high_resolution_clock::now().time_since_epoch().count();
mRand.SetSeed(timeSeed);

// Set seed of random engine
rng.seed(12345);

// Load reweighting histograms from CCDB if antinuclei efficiency processing is enabled
if (doprocessAntinucleiEfficiency || doprocessJetsMCgen || doprocessJetsMCrec) {
ccdb->setURL(urlToCcdb.value);
Expand Down Expand Up @@ -389,6 +396,12 @@ struct AntinucleiInJets {
registryMC.add("antiproton_gen_ue", "antiproton_gen_ue", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_gen_full", "antiproton_gen_full", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});

// Generated spectra of antiprotons for closure test
registryMC.add("antiproton_gen_jet_data", "antiproton_gen_jet_data", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_gen_ue_data", "antiproton_gen_ue_data", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_gen_jet_mc", "antiproton_gen_jet_mc", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_gen_ue_mc", "antiproton_gen_ue_mc", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});

// Normalization histogram
registryMC.add("antiproton_deltay_deltaphi_jet", "antiproton_deltay_deltaphi_jet", HistType::kTH2F, {{2000, -1.0, 1.0, "#Delta#it{y}"}, {2000, 0.0, 2.0, "#Delta#phi"}});
registryMC.add("antiproton_deltay_deltaphi_ue", "antiproton_deltay_deltaphi_ue", HistType::kTH2F, {{2000, -1.0, 1.0, "#Delta#it{y}"}, {2000, 0.0, 2.0, "#Delta#phi"}});
Expand Down Expand Up @@ -416,6 +429,16 @@ struct AntinucleiInJets {
registryMC.add("antiproton_rec_tpc_full", "antiproton_rec_tpc_full", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_rec_tof_full", "antiproton_rec_tof_full", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});

// Reconstructed spectra of antiprotons for closure test
registryMC.add("antiproton_rec_tpc_jet_data", "antiproton_rec_tpc_jet_data", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_rec_tof_jet_data", "antiproton_rec_tof_jet_data", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_rec_tpc_ue_data", "antiproton_rec_tpc_ue_data", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_rec_tof_ue_data", "antiproton_rec_tof_ue_data", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_rec_tpc_jet_mc", "antiproton_rec_tpc_jet_mc", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_rec_tof_jet_mc", "antiproton_rec_tof_jet_mc", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_rec_tpc_ue_mc", "antiproton_rec_tpc_ue_mc", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_rec_tof_ue_mc", "antiproton_rec_tof_ue_mc", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});

// Fraction of primary antiprotons
registryMC.add("antiproton_prim_jet", "antiproton_prim_jet", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_incl_jet", "antiproton_incl_jet", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
Expand Down Expand Up @@ -2044,6 +2067,10 @@ struct AntinucleiInJets {
// Get particles in this MC collision
const auto mcParticlesThisMcColl = mcParticles.sliceBy(mcParticlesPerMcCollision, collision.globalIndex());

// Generate random number to decide if simulated event has to be considered as data or MC in the closure test
int sample = generateRandomNr(rng);
bool isPseudoData = (sample == 0);

// Loop over MC particles
for (const auto& particle : mcParticlesThisMcColl) {

Expand Down Expand Up @@ -2134,6 +2161,13 @@ struct AntinucleiInJets {
// Fill histogram for generated antiprotons
registryMC.fill(HIST("antiproton_gen_jet"), particle.pt(), weightJet);

// Fill histograms for generated antiprotons for closure test
if (isPseudoData) {
registryMC.fill(HIST("antiproton_gen_jet_data"), particle.pt());
} else {
registryMC.fill(HIST("antiproton_gen_jet_mc"), particle.pt());
}

// Fill 2d (pt,eta) distribution of antiprotons
registryMC.fill(HIST("antiproton_eta_pt_jet"), particle.pt(), particle.eta(), weightJet);
}
Expand Down Expand Up @@ -2182,12 +2216,24 @@ struct AntinucleiInJets {
// Fill histogram for antiprotons in the UE
registryMC.fill(HIST("antiproton_gen_ue"), protonVec.Pt(), weightUe);

// Fill histograms for generated antiprotons for closure test
if (sample == 0) {
registryMC.fill(HIST("antiproton_gen_ue_data"), protonVec.Pt());
} else {
registryMC.fill(HIST("antiproton_gen_ue_mc"), protonVec.Pt());
}

// Fill 2d (pt,eta) distribution of antiprotons
registryMC.fill(HIST("antiproton_eta_pt_ue"), protonVec.Pt(), protonVec.Eta(), weightUe);
}
}
if (isAtLeastOneJetSelected) {
registryMC.fill(HIST("genEvents"), 3.5);
if (isPseudoData) {
registryMC.fill(HIST("genEvents"), 4.5);
} else {
registryMC.fill(HIST("genEvents"), 5.5);
}
}

// Shrink large vectors
Expand Down Expand Up @@ -2269,6 +2315,10 @@ struct AntinucleiInJets {
// Get tracks in this MC collision
const auto mcTracksThisMcColl = mcTracks.sliceBy(mcTracksPerMcCollision, collision.globalIndex());

// Generate random number to decide if simulated event has to be considered as data or MC in the closure test
int sample = generateRandomNr(rng);
bool isPseudoData = (sample == 0);

// Loop over reconstructed tracks
int id(-1);
for (auto const& track : mcTracksThisMcColl) {
Expand Down Expand Up @@ -2433,8 +2483,23 @@ struct AntinucleiInJets {
// Fill histograms (TPC and TOF) only for selected candidates
if (passedItsPidProt && nsigmaTPCPr > minNsigmaTpc && nsigmaTPCPr < maxNsigmaTpc) {
registryMC.fill(HIST("antiproton_rec_tpc_jet"), pt, weightJet);

// Fill histograms for reconstructed antiprotons for closure test
if (isPseudoData) {
registryMC.fill(HIST("antiproton_rec_tpc_jet_data"), pt);
} else {
registryMC.fill(HIST("antiproton_rec_tpc_jet_mc"), pt);
}

if (track.hasTOF() && nsigmaTOFPr > minNsigmaTof && nsigmaTOFPr < maxNsigmaTof) {
registryMC.fill(HIST("antiproton_rec_tof_jet"), pt, weightJet);

// Fill histograms for reconstructed antiprotons for closure test
if (isPseudoData) {
registryMC.fill(HIST("antiproton_rec_tof_jet_data"), pt);
} else {
registryMC.fill(HIST("antiproton_rec_tof_jet_mc"), pt);
}
}
}
}
Expand Down Expand Up @@ -2515,14 +2580,34 @@ struct AntinucleiInJets {
// Fill histograms (TPC and TOF) only for selected candidates
if (passedItsPidProt && nsigmaTPCPr > minNsigmaTpc && nsigmaTPCPr < maxNsigmaTpc) {
registryMC.fill(HIST("antiproton_rec_tpc_ue"), pt, weightUe);

// Fill histograms for reconstructed antiprotons for closure test
if (isPseudoData) {
registryMC.fill(HIST("antiproton_rec_tpc_ue_data"), pt);
} else {
registryMC.fill(HIST("antiproton_rec_tpc_ue_mc"), pt);
}

if (track.hasTOF() && nsigmaTOFPr > minNsigmaTof && nsigmaTOFPr < maxNsigmaTof) {
registryMC.fill(HIST("antiproton_rec_tof_ue"), pt, weightUe);

// Fill histograms for reconstructed antiprotons for closure test
if (isPseudoData) {
registryMC.fill(HIST("antiproton_rec_tof_ue_data"), pt);
} else {
registryMC.fill(HIST("antiproton_rec_tof_ue_mc"), pt);
}
}
}
}
}
if (isAtLeastOneJetSelected) {
registryMC.fill(HIST("recEvents"), 9.5);
if (isPseudoData) {
registryMC.fill(HIST("recEvents"), 10.5);
} else {
registryMC.fill(HIST("recEvents"), 11.5);
}
}

// Shrink large vectors
Expand Down
Loading