Skip to content

Commit 29d5b02

Browse files
authored
[PWGLF] added histos for reweighting (#13041)
1 parent 3dd04aa commit 29d5b02

File tree

1 file changed

+21
-0
lines changed

1 file changed

+21
-0
lines changed

PWGLF/Tasks/Nuspex/antinucleiInJets.cxx

Lines changed: 21 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -338,6 +338,13 @@ struct AntinucleiInJets {
338338
registryMC.add("helium3_rec_tpc_ue", "helium3_rec_tpc_ue", HistType::kTH1F, {{nbins, 3 * min, 3 * max, "#it{p}_{T} (GeV/#it{c})"}});
339339
registryMC.add("antihelium3_rec_tpc_jet", "antihelium3_rec_tpc_jet", HistType::kTH1F, {{nbins, 3 * min, 3 * max, "#it{p}_{T} (GeV/#it{c})"}});
340340
registryMC.add("antihelium3_rec_tpc_ue", "antihelium3_rec_tpc_ue", HistType::kTH1F, {{nbins, 3 * min, 3 * max, "#it{p}_{T} (GeV/#it{c})"}});
341+
342+
// Generated spectra needed for reweighting
343+
registryMC.add("protonBar", "protonBar", HistType::kTH1F, {{5000, 0, 5, "#it{p}_{T} (GeV/#it{c})"}});
344+
registryMC.add("lambdaBar", "lambdaBar", HistType::kTH1F, {{5000, 0, 5, "#it{p}_{T} (GeV/#it{c})"}});
345+
registryMC.add("xiBar", "xiBar", HistType::kTH1F, {{5000, 0, 5, "#it{p}_{T} (GeV/#it{c})"}});
346+
registryMC.add("omegaBar", "omegaBar", HistType::kTH1F, {{5000, 0, 5, "#it{p}_{T} (GeV/#it{c})"}});
347+
registryMC.add("sigmaBar", "sigmaBar", HistType::kTH1F, {{5000, 0, 5, "#it{p}_{T} (GeV/#it{c})"}});
341348
}
342349

343350
// Systematic uncertainties (Data)
@@ -1261,6 +1268,7 @@ struct AntinucleiInJets {
12611268
case PDG_t::kProtonBar:
12621269
registryMC.fill(HIST("antip_gen_jet"), particle.pt());
12631270
registryMC.fill(HIST("antip_gen_ue"), particle.pt());
1271+
registryMC.fill(HIST("protonBar"), particle.pt());
12641272
break;
12651273
case o2::constants::physics::Pdg::kDeuteron:
12661274
registryMC.fill(HIST("deuteron_gen_jet"), particle.pt());
@@ -1278,6 +1286,19 @@ struct AntinucleiInJets {
12781286
registryMC.fill(HIST("antihelium3_gen_jet"), particle.pt());
12791287
registryMC.fill(HIST("antihelium3_gen_ue"), particle.pt());
12801288
break;
1289+
// Histograms for re-weighting
1290+
case PDG_t::kLambda0Bar:
1291+
registryMC.fill(HIST("lambdaBar"), particle.pt());
1292+
break;
1293+
case PDG_t::kXiPlusBar:
1294+
registryMC.fill(HIST("xiBar"), particle.pt());
1295+
break;
1296+
case PDG_t::kOmegaPlusBar:
1297+
registryMC.fill(HIST("omegaBar"), particle.pt());
1298+
break;
1299+
case PDG_t::kSigmaBarMinus:
1300+
registryMC.fill(HIST("sigmaBar"), particle.pt());
1301+
break;
12811302
}
12821303
}
12831304

0 commit comments

Comments
 (0)