[PWGLF] V0 analysis vs multiplicity: remove unbound index (#9630)

Author: ercolessi

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,

PWGLF/TableProducer/Strangeness/v0qaanalysis.cxx

@@ -237,7 +237,7 @@ struct LfV0qaanalysis {
           TMath::Abs(v0.negTrack_as<DauTracks>().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(),

PWGLF/Tasks/Strangeness/v0postprocessing.cxx

@@ -44,6 +44,7 @@ struct v0postprocessing {
   Configurable<float> ntpcsigma{"ntpcsigma", 5, "N sigma TPC"};
   Configurable<float> ntpcsigmaMC{"ntpcsigmaMC", 100, "N sigma TPC for MC"};
   Configurable<float> etadau{"etadau", 0.8, "Eta Daughters"};
+  Configurable<float> minITShits{"minITShits", 2, "min ITS hits"};
   Configurable<bool> isMC{"isMC", 1, "isMC"};
   Configurable<bool> evSel{"evSel", 1, "evSel"};
   Configurable<bool> 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());