From b2b3c06cf01ba51a7ffe11557a33c41d49a75587 Mon Sep 17 00:00:00 2001 From: lorber98 <95907752+lorber98@users.noreply.github.com> Date: Wed, 24 Dec 2025 16:53:27 +0100 Subject: [PATCH 1/3] Added configurable to select FT0C or FT0M --- .../Strangeness/strangenessInJetsIons.cxx | 41 ++++++++++++++++--- 1 file changed, 35 insertions(+), 6 deletions(-) diff --git a/PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx b/PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx index 7d23099b6ff..e53d60c9ecb 100644 --- a/PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx +++ b/PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx @@ -76,11 +76,15 @@ using namespace o2::constants::math; using std::array; // Define convenient aliases for joined AOD tables -using SelCollisions = soa::Join; +using SelCollisions = soa::Join; using SimCollisions = soa::Join; + aod::CentFT0Cs, + aod::CentFT0Ms>; using DaughterTracks = soa::Join; @@ -123,6 +127,7 @@ struct StrangenessInJetsIons { Configurable deltaEtaEdge{"deltaEtaEdge", 0.05, "eta gap from detector edge"}; Configurable cfgSkimmedProcessing{"cfgSkimmedProcessing", false, "Enable processing of skimmed data"}; Configurable triggerName{"triggerName", "fOmega", "Software trigger name"}; + Configurable centrEstimator{"centrEstimator", "FT0C", "Select centrality estimator. Options: FT0C, FT0M"}; // Event selection Configurable requireNoSameBunchPileup{"requireNoSameBunchPileup", true, "Require kNoSameBunchPileup selection"}; @@ -204,7 +209,17 @@ struct StrangenessInJetsIons { // Define binning and axis specifications for multiplicity, eta, pT, PID, and invariant mass histograms std::vector multBinning = {0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100}; - AxisSpec multAxis = {multBinning, "FT0C percentile"}; + std::string multAxTitle; + std::string centrEstimatorStr = static_cast(centrEstimator); + if (centrEstimatorStr == "FT0C") { + multAxTitle = "FT0C percentile"; + } else if (centrEstimatorStr == "FT0M") { + multAxTitle = "FT0M percentile"; + } else { + LOG(fatal) << "Centrality estimator " << centrEstimatorStr << " not available. Exit." << endl; + } + AxisSpec multAxis = {multBinning, multAxTitle}; + const AxisSpec ptAxis{100, 0.0, 10.0, "#it{p}_{T} (GeV/#it{c})"}; const AxisSpec invMassK0sAxis{200, 0.44, 0.56, "m_{#pi#pi} (GeV/#it{c}^{2})"}; const AxisSpec invMassLambdaAxis{200, 1.09, 1.14, "m_{p#pi} (GeV/#it{c}^{2})"}; @@ -1140,7 +1155,12 @@ struct StrangenessInJetsIons { registryData.fill(HIST("number_of_events_data"), 7.5); // Event multiplicity - const float multiplicity = collision.centFT0C(); + float multiplicity; + if (static_cast(centrEstimator) == "FT0C") { + multiplicity = collision.centFT0C(); + } else { + multiplicity = collision.centFT0M(); + } // Fill event multiplicity registryData.fill(HIST("number_of_events_vsmultiplicity"), multiplicity); @@ -1351,7 +1371,11 @@ struct StrangenessInJetsIons { if (!selectRecoEvent(recoColl)) continue; - multiplicity = recoColl.centFT0C(); + if (static_cast(centrEstimator) == "FT0C") { + multiplicity = recoColl.centFT0C(); + } else { + multiplicity = recoColl.centFT0M(); + } // LOGF(info, " MC GEN index: %d", collision.globalIndex()); // LOGF(info, " MC RECO index: %d", recoColl.mcCollisionId()); // LOGF(info, " multiplicity: %f", multiplicity); @@ -1668,7 +1692,12 @@ struct StrangenessInJetsIons { registryMC.fill(HIST("number_of_events_mc_rec"), 4.5); // Event multiplicity - const float multiplicity = collision.centFT0C(); + float multiplicity; + if (static_cast(centrEstimator) == "FT0C") { + multiplicity = collision.centFT0C(); + } else { + multiplicity = collision.centFT0M(); + } // Number of V0 and cascades per collision auto v0sPerColl = fullV0s.sliceBy(perCollisionV0, collision.globalIndex()); From 2951c3f4f55197c256332e3ce5c67b9babd5ea14 Mon Sep 17 00:00:00 2001 From: lorber98 <95907752+lorber98@users.noreply.github.com> Date: Wed, 24 Dec 2025 22:09:03 +0100 Subject: [PATCH 2/3] Added configurable for centrality axis --- PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx b/PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx index e53d60c9ecb..4134d023193 100644 --- a/PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx +++ b/PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx @@ -122,7 +122,7 @@ struct StrangenessInJetsIons { Configurable> particleOfInterest{"particleOfInterest", {0, 0, 0, 0, 0}, "Particles to study: [K0 and Lambda, Xi and Omega, Pion, Kaon, Proton]"}; Configurable minJetPt{"minJetPt", 10.0, "Minimum reconstructed pt of the jet (GeV/c)"}; - Configurable rJet{"rJet", 0.3, "Jet resolution parameter (R)"}; + Configurable rJet{"rJet", 0.4, "Jet resolution parameter (R)"}; Configurable zVtx{"zVtx", 10.0, "Maximum z-vertex position"}; Configurable deltaEtaEdge{"deltaEtaEdge", 0.05, "eta gap from detector edge"}; Configurable cfgSkimmedProcessing{"cfgSkimmedProcessing", false, "Enable processing of skimmed data"}; @@ -183,6 +183,9 @@ struct StrangenessInJetsIons { Configurable deltaMassOmega{"deltaMassOmega", 0.02f, "Mass window for Omega rejection"}; Configurable deltaMassLambda{"deltaMassLambda", 0.02f, "Mass window for Lambda inclusion"}; + // Axes + ConfigurableAxis multBinning{"multBinning", {VARIABLE_WIDTH, 0, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100}, "Binning of the centrality axis"}; + struct : ConfigurableGroup { ConfigurableAxis longLivedBinsNsigma{"longLivedBinsNsigma", {200, -10.f, 10.f}, "Binning of nSigma axis"}; ConfigurableAxis longLivedBinsPt{"longLivedBinsPt", {VARIABLE_WIDTH, -5.0, -4.8, -4.6, -4.4, -4.2, -4.0, -3.8, -3.6, -3.4, -3.2, -3.0, -2.8, -2.6, -2.4, -2.2, -2.0, -1.9, -1.8, -1.7, -1.6, -1.5, -1.4, -1.3, -1.2, -1.1, -1.0, -0.95, -0.9, -0.85, -0.8, -0.75, -0.7, -0.65, -0.6, -0.55, -0.5, -0.45, -0.4, -0.35, -0.3, -0.25, -0.2, -0.18, -0.16, -0.14, -0.12, -0.1, 0.0, 0.1, 0.12, 0.14, 0.16, 0.18, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.2, 4.4, 4.6, 4.8, 5.0}, "Binning of the pT axis"}; @@ -208,7 +211,6 @@ struct StrangenessInJetsIons { } // Define binning and axis specifications for multiplicity, eta, pT, PID, and invariant mass histograms - std::vector multBinning = {0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100}; std::string multAxTitle; std::string centrEstimatorStr = static_cast(centrEstimator); if (centrEstimatorStr == "FT0C") { From e605e61db41e5e192f02eee1872f670bba3448b4 Mon Sep 17 00:00:00 2001 From: lorber98 <95907752+lorber98@users.noreply.github.com> Date: Wed, 14 Jan 2026 12:19:49 +0100 Subject: [PATCH 3/3] Converted configurable for centrality estimator from string to int --- PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx | 15 +++++++-------- 1 file changed, 7 insertions(+), 8 deletions(-) diff --git a/PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx b/PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx index 4134d023193..6981b638514 100644 --- a/PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx +++ b/PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx @@ -127,7 +127,7 @@ struct StrangenessInJetsIons { Configurable deltaEtaEdge{"deltaEtaEdge", 0.05, "eta gap from detector edge"}; Configurable cfgSkimmedProcessing{"cfgSkimmedProcessing", false, "Enable processing of skimmed data"}; Configurable triggerName{"triggerName", "fOmega", "Software trigger name"}; - Configurable centrEstimator{"centrEstimator", "FT0C", "Select centrality estimator. Options: FT0C, FT0M"}; + Configurable centrEstimator{"centrEstimator", 0, "Select centrality estimator. Options: 0 = FT0C, 1 = FT0M"}; // Event selection Configurable requireNoSameBunchPileup{"requireNoSameBunchPileup", true, "Require kNoSameBunchPileup selection"}; @@ -212,13 +212,12 @@ struct StrangenessInJetsIons { // Define binning and axis specifications for multiplicity, eta, pT, PID, and invariant mass histograms std::string multAxTitle; - std::string centrEstimatorStr = static_cast(centrEstimator); - if (centrEstimatorStr == "FT0C") { + if (centrEstimator == 0) { multAxTitle = "FT0C percentile"; - } else if (centrEstimatorStr == "FT0M") { + } else if (centrEstimator == 1) { multAxTitle = "FT0M percentile"; } else { - LOG(fatal) << "Centrality estimator " << centrEstimatorStr << " not available. Exit." << endl; + LOG(fatal) << "Centrality estimator " << centrEstimator << " not available. Exit." << endl; } AxisSpec multAxis = {multBinning, multAxTitle}; @@ -1158,7 +1157,7 @@ struct StrangenessInJetsIons { // Event multiplicity float multiplicity; - if (static_cast(centrEstimator) == "FT0C") { + if (centrEstimator == 0) { multiplicity = collision.centFT0C(); } else { multiplicity = collision.centFT0M(); @@ -1373,7 +1372,7 @@ struct StrangenessInJetsIons { if (!selectRecoEvent(recoColl)) continue; - if (static_cast(centrEstimator) == "FT0C") { + if (centrEstimator == 0) { multiplicity = recoColl.centFT0C(); } else { multiplicity = recoColl.centFT0M(); @@ -1695,7 +1694,7 @@ struct StrangenessInJetsIons { // Event multiplicity float multiplicity; - if (static_cast(centrEstimator) == "FT0C") { + if (centrEstimator == 0) { multiplicity = collision.centFT0C(); } else { multiplicity = collision.centFT0M();