From 5e5fe182e8e9eb9933e0d14144f67b698df47f75 Mon Sep 17 00:00:00 2001 From: Francesca Ercolessi Date: Thu, 30 Jan 2025 17:06:15 +0100 Subject: [PATCH 1/2] remove unbounded index --- PWGLF/DataModel/v0qaanalysis.h | 5 ++--- PWGLF/TableProducer/Strangeness/v0qaanalysis.cxx | 4 ++-- 2 files changed, 4 insertions(+), 5 deletions(-) diff --git a/PWGLF/DataModel/v0qaanalysis.h b/PWGLF/DataModel/v0qaanalysis.h index dfbc50b8a32..efe82ed730f 100644 --- a/PWGLF/DataModel/v0qaanalysis.h +++ b/PWGLF/DataModel/v0qaanalysis.h @@ -22,7 +22,6 @@ namespace o2::aod namespace myv0candidates { -DECLARE_SOA_INDEX_COLUMN(Collision, collision); DECLARE_SOA_COLUMN(V0Pt, v0pt, float); DECLARE_SOA_COLUMN(RapLambda, raplambda, float); DECLARE_SOA_COLUMN(RapK0Short, rapk0short, float); @@ -66,8 +65,8 @@ DECLARE_SOA_COLUMN(QtArm, qtarm, float); } // namespace myv0candidates -DECLARE_SOA_TABLE(MyV0Candidates, "AOD", "MYV0CANDIDATES", o2::soa::Index<>, - myv0candidates::CollisionId, myv0candidates::V0Pt, myv0candidates::RapLambda, myv0candidates::RapK0Short, +DECLARE_SOA_TABLE(MyV0Candidates, "AOD", "MYV0CANDIDATES", + myv0candidates::V0Pt, myv0candidates::RapLambda, myv0candidates::RapK0Short, myv0candidates::MassLambda, myv0candidates::MassAntiLambda, myv0candidates::MassK0Short, myv0candidates::V0Radius, myv0candidates::V0CosPA, myv0candidates::V0DCAPosToPV, myv0candidates::V0DCANegToPV, myv0candidates::V0DCAV0Daughters, diff --git a/PWGLF/TableProducer/Strangeness/v0qaanalysis.cxx b/PWGLF/TableProducer/Strangeness/v0qaanalysis.cxx index 228850ece38..34b19005e4d 100644 --- a/PWGLF/TableProducer/Strangeness/v0qaanalysis.cxx +++ b/PWGLF/TableProducer/Strangeness/v0qaanalysis.cxx @@ -237,7 +237,7 @@ struct LfV0qaanalysis { TMath::Abs(v0.negTrack_as().eta()) < etadau) { // Fill table - myv0s(v0.globalIndex(), v0.pt(), v0.yLambda(), v0.yK0Short(), + myv0s(v0.pt(), v0.yLambda(), v0.yK0Short(), v0.mLambda(), v0.mAntiLambda(), v0.mK0Short(), v0.v0radius(), v0.v0cosPA(), v0.dcapostopv(), v0.dcanegtopv(), v0.dcaV0daughters(), @@ -363,7 +363,7 @@ struct LfV0qaanalysis { ) { // Fill table - myv0s(v0.globalIndex(), v0.pt(), v0.yLambda(), v0.yK0Short(), + myv0s(v0.pt(), v0.yLambda(), v0.yK0Short(), v0.mLambda(), v0.mAntiLambda(), v0.mK0Short(), v0.v0radius(), v0.v0cosPA(), v0.dcapostopv(), v0.dcanegtopv(), v0.dcaV0daughters(), From 7f7465fa20e8cc53b94629a6e41fb169060913bb Mon Sep 17 00:00:00 2001 From: Francesca Ercolessi Date: Thu, 30 Jan 2025 17:16:28 +0100 Subject: [PATCH 2/2] add MC dau check and ITS cl cut --- PWGLF/Tasks/Strangeness/v0postprocessing.cxx | 11 ++++++++--- 1 file changed, 8 insertions(+), 3 deletions(-) diff --git a/PWGLF/Tasks/Strangeness/v0postprocessing.cxx b/PWGLF/Tasks/Strangeness/v0postprocessing.cxx index d66f2814e1e..7d768b96ed3 100644 --- a/PWGLF/Tasks/Strangeness/v0postprocessing.cxx +++ b/PWGLF/Tasks/Strangeness/v0postprocessing.cxx @@ -44,6 +44,7 @@ struct v0postprocessing { Configurable ntpcsigma{"ntpcsigma", 5, "N sigma TPC"}; Configurable ntpcsigmaMC{"ntpcsigmaMC", 100, "N sigma TPC for MC"}; Configurable etadau{"etadau", 0.8, "Eta Daughters"}; + Configurable minITShits{"minITShits", 2, "min ITS hits"}; Configurable isMC{"isMC", 1, "isMC"}; Configurable evSel{"evSel", 1, "evSel"}; Configurable hasTOF2Leg{"hasTOF2Leg", 0, "hasTOF2Leg"}; @@ -131,6 +132,10 @@ struct v0postprocessing { continue; if (TMath::Abs(candidate.v0negeta()) > etadau) continue; + if (TMath::Abs(candidate.v0positshits()) < minITShits) + continue; + if (TMath::Abs(candidate.v0negitshits()) < minITShits) + continue; if (TMath::Abs(candidate.v0dcanegtopv()) < dcanegtopv) continue; if (TMath::Abs(candidate.v0dcapostopv()) < dcapostopv) @@ -244,7 +249,7 @@ struct v0postprocessing { TMath::Abs(candidate.masslambda() - o2::constants::physics::MassLambda0) > v0rejK0s && TMath::Abs(candidate.ntpcsigmanegpi()) <= ntpcsigmaMC && TMath::Abs(candidate.ntpcsigmapospi()) <= ntpcsigmaMC && - (candidate.pdgcode() == 310)) { + (candidate.pdgcode() == 310) && candidate.isdauk0short()) { registry.fill(HIST("hArmenterosPodolanski"), candidate.alpha(), candidate.qtarm()); @@ -275,7 +280,7 @@ struct v0postprocessing { if (TMath::Abs(candidate.ntpcsigmanegpi()) <= ntpcsigmaMC && TMath::Abs(candidate.ntpcsigmapospr()) <= ntpcsigmaMC && candidate.ctaulambda() < ctauLambda && TMath::Abs(candidate.masslambda() - o2::constants::physics::MassLambda0) < 0.075 && - candidate.pdgcode() == 3122) { + candidate.pdgcode() == 3122 && candidate.isdaulambda()) { registry.fill(HIST("hMassLambda_MC"), candidate.masslambda()); registry.fill(HIST("hMassVsPtLambda_MC"), candidate.v0pt(), candidate.multft0m(), candidate.masslambda()); @@ -295,7 +300,7 @@ struct v0postprocessing { if (TMath::Abs(candidate.ntpcsigmanegpr()) <= ntpcsigmaMC && TMath::Abs(candidate.ntpcsigmapospi()) <= ntpcsigmaMC && candidate.ctauantilambda() < ctauLambda && TMath::Abs(candidate.massantilambda() - o2::constants::physics::MassLambda0) < 0.075 && - candidate.pdgcode() == -3122) { + candidate.pdgcode() == -3122 && candidate.isdauantilambda()) { registry.fill(HIST("hMassAntiLambda_MC"), candidate.massantilambda()); registry.fill(HIST("hMassVsPtAntiLambda_MC"), candidate.v0pt(), candidate.multft0m(), candidate.massantilambda());