From 47162b1fa61ce53d846893f868b06fb4f47730b5 Mon Sep 17 00:00:00 2001 From: chiarapinto Date: Fri, 19 Sep 2025 11:36:56 +0200 Subject: [PATCH] [PWGLF] added histos for reweighting --- PWGLF/Tasks/Nuspex/antinucleiInJets.cxx | 21 +++++++++++++++++++++ 1 file changed, 21 insertions(+) diff --git a/PWGLF/Tasks/Nuspex/antinucleiInJets.cxx b/PWGLF/Tasks/Nuspex/antinucleiInJets.cxx index e2b2459bc2b..20d48d727d0 100644 --- a/PWGLF/Tasks/Nuspex/antinucleiInJets.cxx +++ b/PWGLF/Tasks/Nuspex/antinucleiInJets.cxx @@ -338,6 +338,13 @@ struct AntinucleiInJets { registryMC.add("helium3_rec_tpc_ue", "helium3_rec_tpc_ue", HistType::kTH1F, {{nbins, 3 * min, 3 * max, "#it{p}_{T} (GeV/#it{c})"}}); registryMC.add("antihelium3_rec_tpc_jet", "antihelium3_rec_tpc_jet", HistType::kTH1F, {{nbins, 3 * min, 3 * max, "#it{p}_{T} (GeV/#it{c})"}}); registryMC.add("antihelium3_rec_tpc_ue", "antihelium3_rec_tpc_ue", HistType::kTH1F, {{nbins, 3 * min, 3 * max, "#it{p}_{T} (GeV/#it{c})"}}); + + // Generated spectra needed for reweighting + registryMC.add("protonBar", "protonBar", HistType::kTH1F, {{5000, 0, 5, "#it{p}_{T} (GeV/#it{c})"}}); + registryMC.add("lambdaBar", "lambdaBar", HistType::kTH1F, {{5000, 0, 5, "#it{p}_{T} (GeV/#it{c})"}}); + registryMC.add("xiBar", "xiBar", HistType::kTH1F, {{5000, 0, 5, "#it{p}_{T} (GeV/#it{c})"}}); + registryMC.add("omegaBar", "omegaBar", HistType::kTH1F, {{5000, 0, 5, "#it{p}_{T} (GeV/#it{c})"}}); + registryMC.add("sigmaBar", "sigmaBar", HistType::kTH1F, {{5000, 0, 5, "#it{p}_{T} (GeV/#it{c})"}}); } // Systematic uncertainties (Data) @@ -1261,6 +1268,7 @@ struct AntinucleiInJets { case PDG_t::kProtonBar: registryMC.fill(HIST("antip_gen_jet"), particle.pt()); registryMC.fill(HIST("antip_gen_ue"), particle.pt()); + registryMC.fill(HIST("protonBar"), particle.pt()); break; case o2::constants::physics::Pdg::kDeuteron: registryMC.fill(HIST("deuteron_gen_jet"), particle.pt()); @@ -1278,6 +1286,19 @@ struct AntinucleiInJets { registryMC.fill(HIST("antihelium3_gen_jet"), particle.pt()); registryMC.fill(HIST("antihelium3_gen_ue"), particle.pt()); break; + // Histograms for re-weighting + case PDG_t::kLambda0Bar: + registryMC.fill(HIST("lambdaBar"), particle.pt()); + break; + case PDG_t::kXiPlusBar: + registryMC.fill(HIST("xiBar"), particle.pt()); + break; + case PDG_t::kOmegaPlusBar: + registryMC.fill(HIST("omegaBar"), particle.pt()); + break; + case PDG_t::kSigmaBarMinus: + registryMC.fill(HIST("sigmaBar"), particle.pt()); + break; } }