From 9e1514bbc0b4d71d523e0d78289a350528b787e9 Mon Sep 17 00:00:00 2001 From: Marta Razza Date: Fri, 10 Jul 2026 18:49:00 +0200 Subject: [PATCH] PWGHF/D2H: PID selections changed in taskDeuteronFromLb --- PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx | 18 ++++++++---------- 1 file changed, 8 insertions(+), 10 deletions(-) diff --git a/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx b/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx index db0af2da598..2c050c1e9c2 100644 --- a/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx +++ b/PWGHF/D2H/Tasks/taskDeuteronFromLb.cxx @@ -66,10 +66,9 @@ struct HfTaskDeuteronFromLb { Configurable cfgITScls{"cfgITScls", 2, "Minimum ITS clusters"}; Configurable cfgMaxPt{"cfgMaxPt", 5.0f, "Maximum pT cut"}; Configurable cfgMinPt{"cfgMinPt", 0.5f, "Minimum pT cut"}; - Configurable cfgTPCNsigma{"cfgTPCNsigma", 4.0f, "TPC n sigma for deuteron PID"}; - Configurable cfgTofNsigmaMin{"cfgTofNsigmaMin", 3.0f, "TOF n sigma min for deuteron PID"}; - Configurable cfgTofNsigmaMax{"cfgTofNsigmaMax", 4.0f, "TOF n sigma max for deuteron PID"}; - Configurable ptThresholdPid{"ptThresholdPid", 0.5f, "pT threshold to switch between 4 and 3 sigmas for TOF PID"}; + Configurable cfgTPCNsigma{"cfgTPCNsigma", 3.0f, "TPC n sigma for deuteron PID"}; + Configurable cfgTofNsigma{"cfgTofNsigma", 3.0f, "TOF n sigma for deuteron PID"}; + Configurable ptThresholdPid{"ptThresholdPid", 0.5f, "pT threshold to switch between TPC and TPC+TOF PID"}; Configurable cfgDCAmin{"cfgDCAmin", 0.05f, "Minimum DCA for deuteron PID"}; Configurable cfgDCAmax{"cfgDCAmax", 1000.0f, "Maximum DCA for deuteron PID"}; Configurable rapidityCut{"rapidityCut", 0.5f, "Rapidity cut"}; @@ -260,11 +259,10 @@ struct HfTaskDeuteronFromLb { } const bool isTPCDe = std::abs(track.tpcNSigmaDe()) < cfgTPCNsigma; - const bool isTOFDe_min = std::abs(track.tofNSigmaDe()) > cfgTofNsigmaMin; - const bool isTOFDe_max = std::abs(track.tofNSigmaDe()) < cfgTofNsigmaMax; + const bool isTOFDe = std::abs(track.tofNSigmaDe()) < cfgTofNsigma; if (track.pt() < ptThresholdPid) { - if (isTPCDe && isTOFDe_max) { + if (isTPCDe) { QAHistos.fill(HIST("Data/ptAntiDeuteron"), track.pt()); QAHistos.fill(HIST("Data/etaAntideuteron"), track.eta()); QAHistos.fill(HIST("Data/hDCAxyVsPt"), track.pt(), dca[0]); @@ -273,7 +271,7 @@ struct HfTaskDeuteronFromLb { QAHistos.fill(HIST("Data/hnSigmaTOFVsPt"), track.pt(), track.tofNSigmaDe()); } } else { - if (isTPCDe && isTOFDe_min && isTOFDe_max) { + if (isTPCDe && isTOFDe) { QAHistos.fill(HIST("Data/ptAntiDeuteron"), track.pt()); QAHistos.fill(HIST("Data/etaAntideuteron"), track.eta()); QAHistos.fill(HIST("Data/hDCAxyVsPt"), track.pt(), dca[0]); @@ -285,7 +283,7 @@ struct HfTaskDeuteronFromLb { } } } - PROCESS_SWITCH(HfTaskDeuteronFromLb, processData, "processData", false); + PROCESS_SWITCH(HfTaskDeuteronFromLb, processData, "processData", true); void processMC(MCCollisionCandidates::iterator const& collision, MCTrackCandidates const& tracks, o2::aod::McParticles const&) { @@ -423,7 +421,7 @@ struct HfTaskDeuteronFromLb { } } - PROCESS_SWITCH(HfTaskDeuteronFromLb, processGen, "processGen", true); + PROCESS_SWITCH(HfTaskDeuteronFromLb, processGen, "processGen", false); }; WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)