From 8464240aad94ead2a45a7bc9709be4efe92460c8 Mon Sep 17 00:00:00 2001 From: Francesca Ercolessi Date: Tue, 4 Nov 2025 11:45:57 +0100 Subject: [PATCH 1/2] add histo for MC percentile --- PWGLF/Tasks/Strangeness/strangenessInJets.cxx | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/PWGLF/Tasks/Strangeness/strangenessInJets.cxx b/PWGLF/Tasks/Strangeness/strangenessInJets.cxx index 3bf011a7220..d07b1fa4372 100644 --- a/PWGLF/Tasks/Strangeness/strangenessInJets.cxx +++ b/PWGLF/Tasks/Strangeness/strangenessInJets.cxx @@ -231,6 +231,7 @@ struct StrangenessInJets { // Event counter registryMC.add("number_of_events_mc_gen", "number of gen events in mc", HistType::kTH1D, {{10, 0, 10, "Event Cuts"}}); + registryMC.add("number_of_events_vsmultiplicity", "number of events vs multiplicity", HistType::kTH1D, {{101, 0, 101, "Multiplicity percentile"}}); // Histograms for analysis switch (particleOfInterest) { @@ -268,6 +269,7 @@ struct StrangenessInJets { // Event counter registryMC.add("number_of_events_mc_rec", "number of rec events in mc", HistType::kTH1D, {{10, 0, 10, "Event Cuts"}}); + registryMC.add("number_of_events_vsmultiplicity", "number of events vs multiplicity", HistType::kTH1D, {{101, 0, 101, "Multiplicity percentile"}}); // Histograms for analysis switch (particleOfInterest) { @@ -1249,6 +1251,7 @@ struct StrangenessInJets { if (jetMinusBkg.pt() < minJetPt) continue; registryMC.fill(HIST("number_of_events_mc_gen"), 4.5); + registryMC.fill(HIST("number_of_events_vsmultiplicity"), genMultiplicity); // Set up two perpendicular cone axes for underlying event estimation TVector3 jetAxis(jet.px(), jet.py(), jet.pz()); @@ -1534,6 +1537,7 @@ struct StrangenessInJets { // Fill event counter for events with at least one selected jet registryMC.fill(HIST("number_of_events_mc_rec"), 4.5); + registryMC.fill(HIST("number_of_events_vsmultiplicity"), multiplicity); // Loop over selected jets for (int i = 0; i < static_cast(selectedJet.size()); i++) { From 3c8f3aca18bbe6da9f366eee848939b6b5c8512a Mon Sep 17 00:00:00 2001 From: Francesca Ercolessi Date: Tue, 4 Nov 2025 11:48:43 +0100 Subject: [PATCH 2/2] minor fix --- PWGLF/Tasks/Strangeness/strangenessInJets.cxx | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/PWGLF/Tasks/Strangeness/strangenessInJets.cxx b/PWGLF/Tasks/Strangeness/strangenessInJets.cxx index d07b1fa4372..f55f608a4f3 100644 --- a/PWGLF/Tasks/Strangeness/strangenessInJets.cxx +++ b/PWGLF/Tasks/Strangeness/strangenessInJets.cxx @@ -231,7 +231,7 @@ struct StrangenessInJets { // Event counter registryMC.add("number_of_events_mc_gen", "number of gen events in mc", HistType::kTH1D, {{10, 0, 10, "Event Cuts"}}); - registryMC.add("number_of_events_vsmultiplicity", "number of events vs multiplicity", HistType::kTH1D, {{101, 0, 101, "Multiplicity percentile"}}); + registryMC.add("number_of_events_vsmultiplicity_gen", "number of events vs multiplicity", HistType::kTH1D, {{101, 0, 101, "Multiplicity percentile"}}); // Histograms for analysis switch (particleOfInterest) { @@ -269,7 +269,7 @@ struct StrangenessInJets { // Event counter registryMC.add("number_of_events_mc_rec", "number of rec events in mc", HistType::kTH1D, {{10, 0, 10, "Event Cuts"}}); - registryMC.add("number_of_events_vsmultiplicity", "number of events vs multiplicity", HistType::kTH1D, {{101, 0, 101, "Multiplicity percentile"}}); + registryMC.add("number_of_events_vsmultiplicity_rec", "number of events vs multiplicity", HistType::kTH1D, {{101, 0, 101, "Multiplicity percentile"}}); // Histograms for analysis switch (particleOfInterest) { @@ -1251,7 +1251,7 @@ struct StrangenessInJets { if (jetMinusBkg.pt() < minJetPt) continue; registryMC.fill(HIST("number_of_events_mc_gen"), 4.5); - registryMC.fill(HIST("number_of_events_vsmultiplicity"), genMultiplicity); + registryMC.fill(HIST("number_of_events_vsmultiplicity_gen"), genMultiplicity); // Set up two perpendicular cone axes for underlying event estimation TVector3 jetAxis(jet.px(), jet.py(), jet.pz()); @@ -1537,7 +1537,7 @@ struct StrangenessInJets { // Fill event counter for events with at least one selected jet registryMC.fill(HIST("number_of_events_mc_rec"), 4.5); - registryMC.fill(HIST("number_of_events_vsmultiplicity"), multiplicity); + registryMC.fill(HIST("number_of_events_vsmultiplicity_rec"), multiplicity); // Loop over selected jets for (int i = 0; i < static_cast(selectedJet.size()); i++) {