From 2d311e8bf884ed02cd12f5bb7da9f58cc72e64e2 Mon Sep 17 00:00:00 2001 From: tutripat <73981392+tutripat@users.noreply.github.com> Date: Thu, 21 Aug 2025 12:30:43 +0100 Subject: [PATCH 1/4] Update to include DCA cut --- .../GlobalEventProperties/dndeta-mft-pp.cxx | 17 +++++++++++++++++ 1 file changed, 17 insertions(+) diff --git a/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx b/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx index c90a3d89474..ce3bdbd6d90 100644 --- a/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx +++ b/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx @@ -91,6 +91,7 @@ struct PseudorapidityDensityMFT { "max allowed Z difference for reconstructed collisions (cm)"}; Configurable usePhiCut{"usePhiCut", true, "use azimuthal angle cut"}; + Configurable useDCAxyCut{"useDCAxyCut", false, "use DCAxy cut"}; Configurable cfgPhiCut{"cfgPhiCut", 0.1f, "Cut on azimuthal angle of MFT tracks"}; Configurable cfgPhiCut1{"cfgPhiCut1", 0.0f, @@ -108,6 +109,7 @@ struct PseudorapidityDensityMFT { Configurable cfgnEta2{"cfgnEta2", -1.0f, "Cut on eta1"}; Configurable cfgChi2NDFMax{"cfgChi2NDFMax", 2000.0f, "Max allowed chi2/NDF for MFT tracks"}; + Configurable maxDCAxy{"maxDCAxy", 2.0f, "Cut on dcaXY"}; HistogramRegistry registry{ "registry", @@ -588,6 +590,11 @@ struct PseudorapidityDensityMFT { if ((phi <= 0.02) || ((phi >= 3.10) && (phi <= 3.23)) || (phi >= 6.21)) continue; } + float dcaxy_cut = retrack.bestDCAXY(); + if (useDCAxyCut) { + if (dcaxy_cut > maxDCAxy) + continue; + } if ((cfgnEta1 < track.eta()) && (track.eta() < cfgnEta2) && track.nClusters() >= cfgnCluster && retrack.ambDegree() > 0 && chi2ndf < cfgChi2NDFMax && (phi > cfgPhiCut1 && phi < cfgPhiCut2)) { registry.fill(HIST("Tracks/2Danalysis/EtaZvtx"), track.eta(), z); } @@ -628,6 +635,11 @@ struct PseudorapidityDensityMFT { if ((phi <= 0.02) || ((phi >= 3.10) && (phi <= 3.23)) || (phi >= 6.21)) continue; } + float dcaxy_cut = retrack.bestDCAXY(); + if (useDCAxyCut) { + if (dcaxy_cut > maxDCAxy) + continue; + } if ((cfgnEta1 < track.eta()) && (track.eta() < cfgnEta2) && track.nClusters() >= cfgnCluster && retrack.ambDegree() > 0 && chi2ndf < cfgChi2NDFMax && (phi > cfgPhiCut1 && phi < cfgPhiCut2)) { registry.fill(HIST("Tracks/Control/Chi2NDF"), chi2ndf); registry.fill(HIST("Tracks/2Danalysis/EtaZvtx_sel8"), track.eta(), z); @@ -646,12 +658,17 @@ struct PseudorapidityDensityMFT { float ndf = std::max(2.0f * track.nClusters() - 5.0f, 1.0f); float chi2ndf = track.chi2() / ndf; float phi = track.phi(); + float dcaxy_cut = retrack.bestDCAXY(); o2::math_utils::bringTo02Pi(phi); if ((cfgnEta1 < track.eta()) && (track.eta() < cfgnEta2) && track.nClusters() >= cfgnCluster && chi2ndf < cfgChi2NDFMax && (phi > cfgPhiCut1 && phi < cfgPhiCut2)) { if (usePhiCut) { if ((phi <= 0.02) || ((phi >= 3.10) && (phi <= 3.23)) || (phi >= 6.21)) continue; } + if (useDCAxyCut) { + if (dcaxy_cut > maxDCAxy) + continue; + } registry.fill(HIST("TracksEtaZvtx"), track.eta(), z); if (midtracks.size() > 0 && retrack.ambDegree() > 0) { registry.fill(HIST("Tracks/EtaZvtx_gt0"), track.eta(), z); From 72eb665841ff03aac810193952bc660a9c5697b8 Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Thu, 21 Aug 2025 11:32:07 +0000 Subject: [PATCH 2/4] Please consider the following formatting changes --- PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx b/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx index ce3bdbd6d90..2dbd095a563 100644 --- a/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx +++ b/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx @@ -638,7 +638,7 @@ struct PseudorapidityDensityMFT { float dcaxy_cut = retrack.bestDCAXY(); if (useDCAxyCut) { if (dcaxy_cut > maxDCAxy) - continue; + continue; } if ((cfgnEta1 < track.eta()) && (track.eta() < cfgnEta2) && track.nClusters() >= cfgnCluster && retrack.ambDegree() > 0 && chi2ndf < cfgChi2NDFMax && (phi > cfgPhiCut1 && phi < cfgPhiCut2)) { registry.fill(HIST("Tracks/Control/Chi2NDF"), chi2ndf); @@ -667,7 +667,7 @@ struct PseudorapidityDensityMFT { } if (useDCAxyCut) { if (dcaxy_cut > maxDCAxy) - continue; + continue; } registry.fill(HIST("TracksEtaZvtx"), track.eta(), z); if (midtracks.size() > 0 && retrack.ambDegree() > 0) { From 659cfd6af6840c09a8e11fff5480daff639d24ff Mon Sep 17 00:00:00 2001 From: tutripat <73981392+tutripat@users.noreply.github.com> Date: Thu, 21 Aug 2025 18:09:34 +0100 Subject: [PATCH 3/4] Update dndeta-mft-pp.cxx --- PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx | 1 + 1 file changed, 1 insertion(+) diff --git a/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx b/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx index 2dbd095a563..0dec6b99ebc 100644 --- a/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx +++ b/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx @@ -41,6 +41,7 @@ #include #include #include +#include using namespace o2; using namespace o2::framework; From ffdf8a1cf9214e25c81c7fe455092f4b87a4e747 Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Thu, 21 Aug 2025 17:10:54 +0000 Subject: [PATCH 4/4] Please consider the following formatting changes --- PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx b/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx index 0dec6b99ebc..a37b4f8c368 100644 --- a/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx +++ b/PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx @@ -35,13 +35,13 @@ #include "TFile.h" +#include #include #include #include #include #include #include -#include using namespace o2; using namespace o2::framework;