[PWGHF] Add derived data creator for Lc in pK0s. Farewell!

PWGHF/DataModel/CandidateReconstructionTables.h

@@ -495,6 +495,8 @@ DECLARE_SOA_DYNAMIC_COLUMN(PtV0Neg, ptV0Neg, //! pt of the negative V0 daughter
                            [](float px, float py) { return RecoDecay::pt(px, py); });
 DECLARE_SOA_DYNAMIC_COLUMN(CtV0, ctV0, //! c*t of the V0
                            [](float xVtxP, float yVtxP, float zVtxP, float xVtxS, float yVtxS, float zVtxS, float px, float py, float pz, double m) -> float { return RecoDecay::ct(std::array{px, py, pz}, RecoDecay::distance(std::array{xVtxP, yVtxP, zVtxP}, std::array{xVtxS, yVtxS, zVtxS}), m); });
+DECLARE_SOA_DYNAMIC_COLUMN(DecayLengthV0, decayLengthV0, //!
+                           [](float xVtxP, float yVtxP, float zVtxP, float xVtxS, float yVtxS, float zVtxS) -> float { return RecoDecay::distance(std::array{xVtxP, yVtxP, zVtxP}, std::array{xVtxS, yVtxS, zVtxS}); });
 DECLARE_SOA_COLUMN(FlagMcMatchRec, flagMcMatchRec, int8_t); //! reconstruction level
 DECLARE_SOA_COLUMN(FlagMcMatchGen, flagMcMatchGen, int8_t); //! generator level
 DECLARE_SOA_COLUMN(OriginMcRec, originMcRec, int8_t);       //! particle origin, reconstruction level
@@ -552,7 +554,8 @@ DECLARE_SOA_TABLE(HfCandCascBase, "AOD", "HFCANDCASCBASE", //!
                   v0data::legacy::MAntiLambda<v0data::PxPos, v0data::PyPos, v0data::PzPos, v0data::PxNeg, v0data::PyNeg, v0data::PzNeg>,
                   v0data::legacy::MK0Short<v0data::PxPos, v0data::PyPos, v0data::PzPos, v0data::PxNeg, v0data::PyNeg, v0data::PzNeg>,
                   v0data::MGamma<v0data::PxPos, v0data::PyPos, v0data::PzPos, v0data::PxNeg, v0data::PyNeg, v0data::PzNeg>,
-                  hf_cand_casc::CtV0<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand::ZSecondaryVertex, hf_cand_casc::V0X, hf_cand_casc::V0Y, hf_cand_casc::V0Z, hf_cand::PxProng1, hf_cand::PyProng1, hf_cand::PzProng1>);
+                  hf_cand_casc::CtV0<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand::ZSecondaryVertex, hf_cand_casc::V0X, hf_cand_casc::V0Y, hf_cand_casc::V0Z, hf_cand::PxProng1, hf_cand::PyProng1, hf_cand::PzProng1>,
+                  hf_cand_casc::DecayLengthV0<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand::ZSecondaryVertex, hf_cand_casc::V0X, hf_cand_casc::V0Y, hf_cand_casc::V0Z>);
 //                  ,
 //                  v0data::MLambda<v0data::PxPos, v0data::PyPos, v0data::PzPos, v0data::PxNeg, v0data::PyNeg, v0data::PzNeg>,
 //                  v0data::MAntiLambda<v0data::PxPos, v0data::PyPos, v0data::PzPos, v0data::PxNeg, v0data::PyNeg, v0data::PzNeg>,

PWGHF/DataModel/DerivedTables.h

@@ -307,9 +307,10 @@ DECLARE_SOA_COLUMN(FlagMcDecayChanGen, flagMcDecayChanGen, int8_t); //! resonant
 DECLARE_TABLES_2P(D0, "D0", d0, 2);
 DECLARE_TABLES_3P(Lc, "LC", lc, 3);
 DECLARE_TABLES_3P(Dplus, "DP", dplus, 4);
-DECLARE_TABLES_3P(Ds, "DS", ds, 9);
 DECLARE_TABLES_3P(Bplus, "BP", bplus, 5);
 DECLARE_TABLES_3P(Dstar, "DST", dstar, 6);
+DECLARE_TABLES_3P(Ds, "DS", ds, 9);
+DECLARE_TABLES_3P(LcToK0sP, "LCC", lc_to_k0s_p, 10);
 // Workaround for the existing B0 macro in termios.h
 #pragma push_macro("B0")
 #undef B0
@@ -355,6 +356,20 @@ DECLARE_SOA_COLUMN(PtProngXi, ptProngXi, float);
 DECLARE_SOA_COLUMN(PtProngPi0, ptProngPi0, float);                                     //! transverse momentum of the first pion prong
 DECLARE_SOA_COLUMN(PtProngPi1, ptProngPi1, float);                                     //! transverse momentum of the second pion prong
 DECLARE_SOA_COLUMN(RSecondaryVertex, rSecondaryVertex, float);                         //! distance of the secondary vertex from the z axis
+// Lc± → K0s p±
+DECLARE_SOA_COLUMN(V0CosPA, v0cosPA, float);               //! cosine of the V0 pointing angle
+DECLARE_SOA_COLUMN(CtV0, ctV0, float);                     //! V0 proper lifetime times c
+DECLARE_SOA_COLUMN(DecayLengthV0, decayLengthV0, float);   //! V0 decay length
+DECLARE_SOA_COLUMN(PtV0Pos, ptV0Pos, float);               //! transverse momentum of V0 pos prong
+DECLARE_SOA_COLUMN(PtV0Neg, ptV0Neg, float);               //! transverse momentum of V0 neg prong
+DECLARE_SOA_COLUMN(DCAPosToPV, dcapostopv, float);         //! V0 pos prong impact param wrt prim vtx in xy
+DECLARE_SOA_COLUMN(DCANegToPV, dcanegtopv, float);         //! V0 neg prong impact param wrt prim vtx in xy
+DECLARE_SOA_COLUMN(DCAV0Daughters, dcaV0daughters, float); //! V0 daughters dca
+DECLARE_SOA_COLUMN(V0Radius, v0radius, float);             //! transverse distance of the V0 decay vertex from the beam line
+DECLARE_SOA_COLUMN(MLambda, mLambda, float);
+DECLARE_SOA_COLUMN(MAntiLambda, mAntiLambda, float);
+DECLARE_SOA_COLUMN(MK0Short, mK0Short, float);
+DECLARE_SOA_COLUMN(MGamma, mGamma, float);
 // D*± → D0(bar) π±
 DECLARE_SOA_COLUMN(SignProng1, signProng1, int8_t);
 // TOF
@@ -801,6 +816,72 @@ DECLARE_SOA_TABLE_STAGED(HfLcMcs, "HFLCMC", //! Table with MC candidate info
                          hf_cand_mc::IsCandidateSwapped,
                          o2::soa::Marker<MarkerLc>);
 
+// ----------------
+// Lc to K0sP
+// ----------------
+
+DECLARE_SOA_TABLE_STAGED(HfLcToK0sPPars, "HFLCTOK0SPPAR", //! Table with candidate properties used for selection
+                         hf_cand::Chi2PCA,
+                         hf_cand_par::Cpa,
+                         hf_cand_par::CpaXY,
+                         hf_cand_par::DecayLength,
+                         hf_cand_par::DecayLengthXY,
+                         hf_cand_par::PtProng0,
+                         hf_cand_par::PtProng1,
+                         hf_cand::ImpactParameter0,
+                         hf_cand::ImpactParameter1,
+                         hf_cand_par::V0Radius,
+                         hf_cand_par::V0CosPA,
+                         hf_cand_par::MLambda,
+                         hf_cand_par::MAntiLambda,
+                         hf_cand_par::MK0Short,
+                         hf_cand_par::MGamma,
+                         hf_cand_par::DCAV0Daughters,
+                         hf_cand_par::PtV0Pos,
+                         hf_cand_par::PtV0Neg,
+                         hf_cand_par::DecayLengthV0,
+                         hf_cand_par::DCAPosToPV,
+                         hf_cand_par::DCANegToPV,
+                         hf_cand_par::NSigTpcPr0,
+                         hf_cand_par::NSigTofPr0,
+                         hf_cand_par::CtV0,
+                         o2::soa::Marker<MarkerLcToK0sP>);
+
+DECLARE_SOA_TABLE_STAGED(HfLcToK0sPParEs, "HFLCTOK0SPPARE", //! Table with additional candidate properties used for selection
+                         hf_cand::XSecondaryVertex,
+                         hf_cand::YSecondaryVertex,
+                         hf_cand::ZSecondaryVertex,
+                         hf_cand::ErrorDecayLength,
+                         hf_cand::ErrorDecayLengthXY,
+                         hf_cand_par::RSecondaryVertex,
+                         hf_cand_par::PProng0,
+                         hf_cand_par::PProng1,
+                         hf_cand::PxProng0,
+                         hf_cand::PyProng0,
+                         hf_cand::PzProng0,
+                         hf_cand::PxProng1,
+                         hf_cand::PyProng1,
+                         hf_cand::PzProng1,
+                         v0data::PxPos,
+                         v0data::PyPos,
+                         v0data::PzPos,
+                         v0data::PxNeg,
+                         v0data::PyNeg,
+                         v0data::PzNeg,
+                         hf_cand::ErrorImpactParameter0,
+                         hf_cand::ErrorImpactParameter1,
+                         hf_cand_par::Ct,
+                         o2::soa::Marker<MarkerLcToK0sP>);
+
+DECLARE_SOA_TABLE_STAGED(HfLcToK0sPMls, "HFLCTOK0SPML", //! Table with candidate selection ML scores
+                         hf_cand_mc::MlScores,
+                         o2::soa::Marker<MarkerLcToK0sP>);
+
+DECLARE_SOA_TABLE_STAGED(HfLcToK0sPMcs, "HFLCTOK0SPMC", //! Table with MC candidate info
+                         hf_cand_mc::FlagMcMatchRec,
+                         hf_cand_mc::OriginMcRec,
+                         o2::soa::Marker<MarkerLcToK0sP>);
+
 // ----------------
 // D+
 // ----------------

PWGHF/TableProducer/CMakeLists.txt

@@ -330,6 +330,11 @@ o2physics_add_dpl_workflow(derived-data-creator-lc-to-p-k-pi
                     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
                     COMPONENT_NAME Analysis)
 
+o2physics_add_dpl_workflow(derived-data-creator-lc-to-k0s-p
+                    SOURCES derivedDataCreatorLcToK0sP.cxx
+                    PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
+                    COMPONENT_NAME Analysis)
+
 o2physics_add_dpl_workflow(derived-data-creator-xic-to-xi-pi-pi
                     SOURCES derivedDataCreatorXicToXiPiPi.cxx
                     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore