From 916e1bea170200d8cf06c9967cced6b90ea7594d Mon Sep 17 00:00:00 2001 From: fmazzasc Date: Fri, 19 Jan 2024 18:32:23 +0100 Subject: [PATCH 1/2] Keep bachelor if compatible with one of the two cascades --- PWGLF/Tasks/nonPromptCascade.cxx | 17 +++++++++++------ 1 file changed, 11 insertions(+), 6 deletions(-) diff --git a/PWGLF/Tasks/nonPromptCascade.cxx b/PWGLF/Tasks/nonPromptCascade.cxx index cd7f79e8073..bf6251924cb 100644 --- a/PWGLF/Tasks/nonPromptCascade.cxx +++ b/PWGLF/Tasks/nonPromptCascade.cxx @@ -724,23 +724,28 @@ struct NonPromptCascadeTask { // QA PID float nSigmaTPC[nParticles]{bachelor.tpcNSigmaKa(), bachelor.tpcNSigmaPi(), protonTrack.tpcNSigmaPr(), pionTrack.tpcNSigmaPi()}; + bool isBachelorSurvived = false; if (isOmega) { if (bachelor.hasTPC()) { LOG(debug) << "TPCSignal bachelor " << bachelor.sign() << "/" << bachelor.tpcInnerParam() << "/" << bachelor.tpcSignal(); - if (nSigmaTPC[0] < cfgCutsPID->get(0u, 0u) || nSigmaTPC[0] > cfgCutsPID->get(0u, 1u)) { - continue; + if (nSigmaTPC[0] > cfgCutsPID->get(0u, 0u) && nSigmaTPC[0] < cfgCutsPID->get(0u, 1u)) { + registry.fill(HIST("h_PIDcutsOmega"), 3, massOmega); + isBachelorSurvived = true; } } - registry.fill(HIST("h_PIDcutsOmega"), 3, massOmega); } if (bachelor.hasTPC()) { LOG(debug) << "TPCSignal bachelor " << bachelor.sign() << "/" << bachelor.tpcInnerParam() << "/" << bachelor.tpcSignal(); - if (nSigmaTPC[1] < cfgCutsPID->get(1u, 0u) || nSigmaTPC[1] > cfgCutsPID->get(1u, 1u)) { - continue; + if (nSigmaTPC[1] > cfgCutsPID->get(1u, 0u) && nSigmaTPC[1] < cfgCutsPID->get(1u, 1u)) { + registry.fill(HIST("h_PIDcutsXi"), 3, massXi); + isBachelorSurvived = true; } } - registry.fill(HIST("h_PIDcutsXi"), 3, massXi); + + if (!isBachelorSurvived) { + continue; + } LOG(debug) << "TPCSignal protonTrack " << protonTrack.sign() << "/" << protonTrack.tpcInnerParam() << "/" << protonTrack.tpcSignal(); if (nSigmaTPC[2] < cfgCutsPID->get(2u, 0u) || nSigmaTPC[2] > cfgCutsPID->get(2u, 1u)) { From c39353285dae53aef972bcae206270aefbd958e9 Mon Sep 17 00:00:00 2001 From: fmazzasc Date: Fri, 19 Jan 2024 18:45:39 +0100 Subject: [PATCH 2/2] Fix omega histogram filling --- PWGLF/Tasks/nonPromptCascade.cxx | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/PWGLF/Tasks/nonPromptCascade.cxx b/PWGLF/Tasks/nonPromptCascade.cxx index bf6251924cb..799bb6c2173 100644 --- a/PWGLF/Tasks/nonPromptCascade.cxx +++ b/PWGLF/Tasks/nonPromptCascade.cxx @@ -752,8 +752,10 @@ struct NonPromptCascadeTask { continue; } + if (isOmega) { + registry.fill(HIST("h_PIDcutsOmega"), 4, massOmega); + } registry.fill(HIST("h_PIDcutsXi"), 4, massXi); - registry.fill(HIST("h_PIDcutsOmega"), 4, massOmega); LOG(debug) << "TPCSignal ntrack " << pionTrack.sign() << "/" << pionTrack.tpcInnerParam() << "/" << pionTrack.tpcSignal(); if (nSigmaTPC[3] < cfgCutsPID->get(3u, 0u) || nSigmaTPC[3] > cfgCutsPID->get(3u, 1u)) {