[PWGLF] Add an additional systematic uncertainty table, and also consider another method for removing shared tracks. (#18238)

Author: 15071832337

PWGLF/DataModel/LFSpincorrelationTables.h

@@ -53,6 +53,12 @@ DECLARE_SOA_COLUMN(ProtonPhi, protonPhi, float);                   //! Proton Ph
 DECLARE_SOA_COLUMN(ProtonIndex, protonIndex, int);                 //! Proton index
 DECLARE_SOA_COLUMN(PionIndex, pionIndex, int);                     //! Pion index
 DECLARE_SOA_COLUMN(DcaV0ToPV, dcaV0ToPV, float);                   //! DCA of V0 to primary vertex
+DECLARE_SOA_COLUMN(TpcPosNClsCrossedRows, tpcPosNClsCrossedRows, float);
+DECLARE_SOA_COLUMN(TpcNegNClsCrossedRows, tpcNegNClsCrossedRows, float);
+DECLARE_SOA_COLUMN(Ctau, ctau, float);
+DECLARE_SOA_COLUMN(TpcNsigmaPos, tpcNsigmaPos, float);
+DECLARE_SOA_COLUMN(TpcNsigmaNeg, tpcNsigmaNeg, float);
+DECLARE_SOA_COLUMN(MK0Short, mK0Short, float);
 
 } // namespace lambdapair
 DECLARE_SOA_TABLE(LambdaPairs, "AOD", "LAMBDAPAIR",
@@ -74,7 +80,13 @@ DECLARE_SOA_TABLE(LambdaPairs, "AOD", "LAMBDAPAIR",
                   lambdapair::ProtonPhi,
                   lambdapair::ProtonIndex,
                   lambdapair::PionIndex,
-                  lambdapair::DcaV0ToPV);
+                  lambdapair::DcaV0ToPV,
+                  lambdapair::TpcPosNClsCrossedRows,
+                  lambdapair::TpcNegNClsCrossedRows,
+                  lambdapair::Ctau,
+                  lambdapair::TpcNsigmaPos,
+                  lambdapair::TpcNsigmaNeg,
+                  lambdapair::MK0Short);
 
 using LambdaPair = LambdaPairs::iterator;
 
@@ -109,6 +121,8 @@ DECLARE_SOA_COLUMN(ProtonPhimc, protonPhimc, float);                   //! Proto
 DECLARE_SOA_COLUMN(ProtonIndexmc, protonIndexmc, int);                 //! Proton index in montecarlo
 DECLARE_SOA_COLUMN(PionIndexmc, pionIndexmc, int);                     //! Pion index in montecarlo
 DECLARE_SOA_COLUMN(DcaV0ToPVmc, dcaV0ToPVmc, float);                   //! DCA of V0 to primary vertex
+
+//! TPC n sigma
 } // namespace lambdapairmc
 DECLARE_SOA_TABLE(LambdaPairmcs, "AOD", "LAMBDAPAIRMC",
                   o2::soa::Index<>,

PWGLF/TableProducer/Strangeness/lambdaspincorrelation.cxx

@@ -299,6 +299,13 @@ struct lambdaspincorrelation {
     std::vector<int> negativeIndex = {};
     std::vector<float> dcaBetweenDaughter = {};
     std::vector<float> dcaV0ToPV = {};
+    std::vector<float> tpcPosNClsCrossedRows = {};
+    std::vector<float> tpcNegNClsCrossedRows = {};
+    std::vector<float> ctau = {};
+    std::vector<float> tpcNsigmaPos = {};
+    std::vector<float> tpcNsigmaNeg = {};
+    std::vector<float> mK0Short = {};
+
     int numbV0 = 0;
     // LOGF(info, "event collisions: (%d)", collision.index());
     auto centrality = collision.centFT0C();
@@ -350,13 +357,26 @@ struct lambdaspincorrelation {
           if (aLambdaTag) {
             histos.fill(HIST("hV0Info"), 3.5);
           }
+
+          float ctauLambda = v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * (o2::constants::physics::MassLambda);
+
           // LOGF(info, "v0 index2: (%d)", v0.index());
           auto postrack1 = v0.template posTrack_as<AllTrackCandidates>();
           auto negtrack1 = v0.template negTrack_as<AllTrackCandidates>();
           positiveIndex.push_back(postrack1.globalIndex());
           negativeIndex.push_back(negtrack1.globalIndex());
           v0Cospa.push_back(v0.v0cosPA());
           dcaV0ToPV.push_back(std::abs(v0.dcav0topv()));
+
+          tpcPosNClsCrossedRows.push_back(postrack1.tpcNClsCrossedRows());
+          tpcNegNClsCrossedRows.push_back(negtrack1.tpcNClsCrossedRows());
+          ctau.push_back(ctauLambda);
+
+          tpcNsigmaPos.push_back(std::abs(postrack1.tpcNSigmaPr()));
+          tpcNsigmaNeg.push_back(std::abs(negtrack1.tpcNSigmaPi()));
+
+          mK0Short.push_back(v0.mK0Short());
+
           v0Radius.push_back(v0.v0radius());
           dcaPositive.push_back(std::abs(v0.dcapostopv()));
           dcaNegative.push_back(std::abs(v0.dcanegtopv()));
@@ -393,7 +413,7 @@ struct lambdaspincorrelation {
           lambdaDummy = lambdaMother.at(i5);
           protonDummy = protonDaughter.at(i5);
           pionDummy = pionDaughter.at(i5);
-          lambdaPair(indexEvent, v0Status.at(i5), doubleStatus.at(i5), v0Cospa.at(i5), v0Radius.at(i5), dcaPositive.at(i5), dcaNegative.at(i5), dcaBetweenDaughter.at(i5), lambdaDummy.Pt(), lambdaDummy.Eta(), lambdaDummy.Phi(), lambdaDummy.M(), protonDummy.Pt(), protonDummy.Eta(), protonDummy.Phi(), positiveIndex.at(i5), negativeIndex.at(i5), dcaV0ToPV.at(i5));
+          lambdaPair(indexEvent, v0Status.at(i5), doubleStatus.at(i5), v0Cospa.at(i5), v0Radius.at(i5), dcaPositive.at(i5), dcaNegative.at(i5), dcaBetweenDaughter.at(i5), lambdaDummy.Pt(), lambdaDummy.Eta(), lambdaDummy.Phi(), lambdaDummy.M(), protonDummy.Pt(), protonDummy.Eta(), protonDummy.Phi(), positiveIndex.at(i5), negativeIndex.at(i5), dcaV0ToPV.at(i5), tpcPosNClsCrossedRows.at(i5), tpcNegNClsCrossedRows.at(i5), ctau.at(i5), tpcNsigmaPos.at(i5), tpcNsigmaNeg.at(i5), mK0Short.at(i5));
         }
       }
     }