From 90a608e09eab7dea32d1f824d2eedcde6a789a37 Mon Sep 17 00:00:00 2001 From: tutripat <73981392+tutripat@users.noreply.github.com> Date: Thu, 4 Jun 2026 13:42:44 +0200 Subject: [PATCH 1/3] Minor configurable changes --- .../PseudorapidityDensityMFT.cxx | 24 +++++++++++-------- 1 file changed, 14 insertions(+), 10 deletions(-) diff --git a/PWGLF/Tasks/GlobalEventProperties/PseudorapidityDensityMFT.cxx b/PWGLF/Tasks/GlobalEventProperties/PseudorapidityDensityMFT.cxx index c765142dbac..831e2747676 100644 --- a/PWGLF/Tasks/GlobalEventProperties/PseudorapidityDensityMFT.cxx +++ b/PWGLF/Tasks/GlobalEventProperties/PseudorapidityDensityMFT.cxx @@ -345,7 +345,8 @@ struct PseudorapidityDensityMFT { Configurable estimatorEta{"estimatorEta", 1.0, "eta range for INEL>0 sample definition"}; - + Configurable usePerCollisionSampleGt0Cut{"usePerCollisionSampleGt0Cut", true, "Require at least one central-barrel track in the per-collision sample"}; + Configurable useMidtracksAndPerCollisionSampleGt0Cut{"useMidtracksAndPerCollisionSampleGt0Cut", true, "Require both INEL>0 midtracks and at least one central-barrel track in the per-collision sample"}; Configurable useEvSel{"useEvSel", true, "use event selection"}; Configurable useNoSameBunchPileup{"useNoSameBunchPileup", true, "reject collisions in case of pileup with another collision in the same foundBC"}; Configurable useGoodItsLayersAll{"useGoodItsLayersAll", true, "all ITS layers are in a good state"}; @@ -1562,20 +1563,22 @@ struct PseudorapidityDensityMFT { registry.fill(HIST("Tracks/2Danalysis/EventsNtrkZvtx_all"), nTrk, z); fillDataCut(DataCutBin::VzWindow); // registry.fill(HIST("EventSelection"), static_cast(EventSelectionBin::Vz)); + const bool hasPerCollisionSample = perCollisionSample.size() > 0; + const bool hasMidtracks = midtracks.size() > 0; if (midtracks.size() > 0) { // registry.fill(HIST("EventSelection"), static_cast(EventSelectionBin::Sel8VzInelGt0)); registry.fill(HIST("EventsNtrkZvtx_gt0"), nTrk, z); registry.fill(HIST("Tracks/2Danalysis/EventsNtrkZvtx_sel8_inelgt0"), nTrk, z); eventsInel.insert(collision.globalIndex()); } - if (perCollisionSample.size() > 0) { - // registry.fill(HIST("EventSelection"), static_cast(EventSelectionBin::PerCollisionSampleGt0)); - fillDataCut(DataCutBin::PerCollisionSampleGt0); - } - if (midtracks.size() > 0 && perCollisionSample.size() > 0) { - // registry.fill(HIST("EventSelection"), static_cast(EventSelectionBin::MidtracksAndPerCollisionSampleGt0)); - fillDataCut(DataCutBin::InelGt0); - } + if (usePerCollisionSampleGt0Cut && !hasPerCollisionSample) { + return;} + fillDataCut(DataCutBin::PerCollisionSampleGt0); + + if (useMidtracksAndPerCollisionSampleGt0Cut && !(hasMidtracks && hasPerCollisionSample)) { + return;} + // registry.fill(HIST("EventSelection"), static_cast(EventSelectionBin::MidtracksAndPerCollisionSampleGt0)); + fillDataCut(DataCutBin::InelGt0); const auto passEventSelection = [&](auto const& collision) { struct EvSelStep { @@ -1626,6 +1629,7 @@ struct PseudorapidityDensityMFT { const float ndf = getTrackNdf(track); const float chi2ndf = track.chi2() / ndf; float phi = track.phi(); + float ptCut = track.pt(); o2::math_utils::bringTo02Pi(phi); const float dcaXyCut = retrack.bestDCAXY(); @@ -1641,7 +1645,7 @@ struct PseudorapidityDensityMFT { ((phi <= PhiVetoLow) || ((phi >= PhiVetoPiMin) && (phi <= PhiVetoPiMax)) || (phi >= PhiVetoHigh))) || - (useDCAxyCut && dcaXyCut > maxDCAxy); + (useDCAxyCut && dcaXyCut > maxDCAxy)|| (usePtCut && ptCut > cfgnPt); if constexpr (std::is_same_v>) { const float dcaZCut = retrack.bestDCAZ(); From 19f7d46334d27e4553dde404edb30f2e12e23c61 Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Thu, 4 Jun 2026 11:44:31 +0000 Subject: [PATCH 2/3] Please consider the following formatting changes --- .../GlobalEventProperties/PseudorapidityDensityMFT.cxx | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) diff --git a/PWGLF/Tasks/GlobalEventProperties/PseudorapidityDensityMFT.cxx b/PWGLF/Tasks/GlobalEventProperties/PseudorapidityDensityMFT.cxx index 831e2747676..be85e87e460 100644 --- a/PWGLF/Tasks/GlobalEventProperties/PseudorapidityDensityMFT.cxx +++ b/PWGLF/Tasks/GlobalEventProperties/PseudorapidityDensityMFT.cxx @@ -1572,11 +1572,13 @@ struct PseudorapidityDensityMFT { eventsInel.insert(collision.globalIndex()); } if (usePerCollisionSampleGt0Cut && !hasPerCollisionSample) { - return;} + return; + } fillDataCut(DataCutBin::PerCollisionSampleGt0); if (useMidtracksAndPerCollisionSampleGt0Cut && !(hasMidtracks && hasPerCollisionSample)) { - return;} + return; + } // registry.fill(HIST("EventSelection"), static_cast(EventSelectionBin::MidtracksAndPerCollisionSampleGt0)); fillDataCut(DataCutBin::InelGt0); @@ -1645,7 +1647,7 @@ struct PseudorapidityDensityMFT { ((phi <= PhiVetoLow) || ((phi >= PhiVetoPiMin) && (phi <= PhiVetoPiMax)) || (phi >= PhiVetoHigh))) || - (useDCAxyCut && dcaXyCut > maxDCAxy)|| (usePtCut && ptCut > cfgnPt); + (useDCAxyCut && dcaXyCut > maxDCAxy) || (usePtCut && ptCut > cfgnPt); if constexpr (std::is_same_v>) { const float dcaZCut = retrack.bestDCAZ(); From 32272dd2299dce48d59142b9a962c15702f79870 Mon Sep 17 00:00:00 2001 From: tutripat <73981392+tutripat@users.noreply.github.com> Date: Thu, 4 Jun 2026 15:27:41 +0200 Subject: [PATCH 3/3] Update PseudorapidityDensityMFT.cxx --- PWGLF/Tasks/GlobalEventProperties/PseudorapidityDensityMFT.cxx | 2 -- 1 file changed, 2 deletions(-) diff --git a/PWGLF/Tasks/GlobalEventProperties/PseudorapidityDensityMFT.cxx b/PWGLF/Tasks/GlobalEventProperties/PseudorapidityDensityMFT.cxx index be85e87e460..ed76d333570 100644 --- a/PWGLF/Tasks/GlobalEventProperties/PseudorapidityDensityMFT.cxx +++ b/PWGLF/Tasks/GlobalEventProperties/PseudorapidityDensityMFT.cxx @@ -83,7 +83,6 @@ AxisSpec dcaXyAxis = {6000, -30, 30}; AxisSpec dcaZAxis = {6000, -30, 30}; AxisSpec dcaXAxis = {6000, -30, 30}; AxisSpec dcaYAxis = {6000, -30, 30}; // previous AxisSpec dcaYAxis = {2000, -10, 10}; - // AxisSpec dcaXyAxis = {600, -0.15f, 0.15f}; // AxisSpec dcaZAxis = {600, -0.15f, 0.15f}; // AxisSpec dcaXAxis = {600, -0.15f, 0.15f}; @@ -91,7 +90,6 @@ AxisSpec dcaYAxis = {6000, -30, 30}; // previous AxisSpec dcaYAxis = {2000, -10 // bin width 0.0005 cm: range [-30, 30] cm => 60/0.0005 = 120000 bins // Keep bin width = 0.0005 cm (5 um): range [-1, 1] cm => 2.0/0.0005 = 4000 bins // AxisSpec axisBinsDCA = {600, -0.15f, 0.15f, "#it{dca}_{xy} (cm)"}; - AxisSpec centAxis = {{0, 10, 20, 30, 40, 50, 60, 70, 80, 100}}; // Vertex position axes (cm)