[PWGLF] Update on Xi1530 and Reso-init (#15668)

Co-authored-by: ALICE Action Bot <alibuild@cern.ch>

Author: mjkim525

PWGLF/DataModel/LFResonanceTables.h

@@ -231,6 +231,9 @@ DECLARE_SOA_COLUMN(DaughterTOFNSigmaNegPr10, daughterTOFNSigmaNegPr10, int8_t);
 DECLARE_SOA_COLUMN(DaughterTOFNSigmaBachPi10, daughterTOFNSigmaBachPi10, int8_t); //! TOF PID x10 of the bachelor daughter as Pion
 DECLARE_SOA_COLUMN(DaughterTOFNSigmaBachKa10, daughterTOFNSigmaBachKa10, int8_t); //! TOF PID x10 of the bachelor daughter as Kaon
 DECLARE_SOA_COLUMN(DaughterTOFNSigmaBachPr10, daughterTOFNSigmaBachPr10, int8_t); //! TOF PID x10 of the bachelor daughter as Proton
+DECLARE_SOA_COLUMN(NCrossedRowsPos, nCrossedRowsPos, uint8_t);                    //! Number of TPC crossed rows of the positive daughter
+DECLARE_SOA_COLUMN(NCrossedRowsNeg, nCrossedRowsNeg, uint8_t);                    //! Number of TPC crossed rows of the negative daughter
+DECLARE_SOA_COLUMN(NCrossedRowsBach, nCrossedRowsBach, uint8_t);                  //! Number of TPC crossed rows of the bachelor daughter
 // For MC
 DECLARE_SOA_INDEX_COLUMN(McParticle, mcParticle); //! Index of the corresponding MC particle
 DECLARE_SOA_COLUMN(IsPhysicalPrimary, isPhysicalPrimary, bool);
@@ -643,6 +646,8 @@ DECLARE_SOA_TABLE(ResoV0s, "AOD", "RESOV0",
                   v0data::DCAPosToPV,
                   v0data::DCANegToPV,
                   v0data::DCAV0ToPV,
+                  resodaughter::NCrossedRowsPos,
+                  resodaughter::NCrossedRowsNeg,
                   resodaughter::MLambda,
                   resodaughter::MAntiLambda,
                   resodaughter::MK0Short,
@@ -710,6 +715,9 @@ DECLARE_SOA_TABLE(ResoCascades, "AOD", "RESOCASCADE",
                   cascdata::DCAXYCascToPV,
                   cascdata::DCAZCascToPV,
                   cascdata::Sign,
+                  resodaughter::NCrossedRowsPos,
+                  resodaughter::NCrossedRowsNeg,
+                  resodaughter::NCrossedRowsBach,
                   resodaughter::MLambda,
                   resodaughter::MXi,
                   resodaughter::TransRadius,

PWGLF/TableProducer/Resonances/resonanceInitializer.cxx

@@ -796,6 +796,8 @@ struct ResonanceInitializer {
                v0.dcapostopv(),
                v0.dcanegtopv(),
                v0.dcav0topv(),
+               static_cast<uint8_t>(v0.template posTrack_as<TrackType>().tpcNClsCrossedRows()),
+               static_cast<uint8_t>(v0.template negTrack_as<TrackType>().tpcNClsCrossedRows()),
                v0.mLambda(),
                v0.mAntiLambda(),
                v0.mK0Short(),
@@ -858,6 +860,9 @@ struct ResonanceInitializer {
                     casc.dcaXYCascToPV(),
                     casc.dcaZCascToPV(),
                     casc.sign(),
+                    static_cast<uint8_t>(casc.template posTrack_as<TrackType>().tpcNClsCrossedRows()),
+                    static_cast<uint8_t>(casc.template negTrack_as<TrackType>().tpcNClsCrossedRows()),
+                    static_cast<uint8_t>(casc.template bachelor_as<TrackType>().tpcNClsCrossedRows()),
                     casc.mLambda(),
                     casc.mXi(),
                     casc.v0radius(), casc.cascradius(), casc.x(), casc.y(), casc.z());
@@ -1736,6 +1741,9 @@ struct ResonanceInitializer {
     resoSpheroCollisions(computeSpherocity(tracks, trackSphMin, trackSphDef));
     resoEvtPlCollisions(0, 0, 0, 0);
     fillMCCollision<false>(collision, mcParticles);
+    // Loop over all MC particles
+    auto mcParts = selectedMCParticles->sliceBy(perMcCollision, collision.mcCollision().globalIndex());
+    fillMCParticles(mcParts, mcParticles);
 
     // Loop over tracks
     if (FilterForDerivedTables.cfgBypassNoPairV0s && (V0s.size() < 1)) {
@@ -1746,10 +1754,6 @@ struct ResonanceInitializer {
       fillMicroTracks<true>(collision, tracks);
     }
     fillV0s<true>(collision, V0s, tracks);
-
-    // Loop over all MC particles
-    auto mcParts = selectedMCParticles->sliceBy(perMcCollision, collision.mcCollision().globalIndex());
-    fillMCParticles(mcParts, mcParticles);
   }
   PROCESS_SWITCH(ResonanceInitializer, processTrackV0MC, "Process for MC", false);
 
@@ -1825,6 +1829,9 @@ struct ResonanceInitializer {
       mult = mcCollision.multMCNParticlesEta10();
 
     fillMCCollision<false>(collision, mcParticles, impactpar, mult);
+    // Loop over all MC particles
+    auto mcParts = selectedMCParticles->sliceBy(perMcCollision, mcId);
+    fillMCParticles(mcParts, mcParticles);
 
     // Loop over tracks
     if (FilterForDerivedTables.cfgBypassNoPairV0s && (V0s.size() < 1)) {
@@ -1839,16 +1846,11 @@ struct ResonanceInitializer {
     }
     fillV0s<true>(collision, V0s, tracks);
     fillCascades<true>(collision, Cascades, tracks);
-
-    // Loop over all MC particles
-    auto mcParts = selectedMCParticles->sliceBy(perMcCollision, mcId);
-    fillMCParticles(mcParts, mcParticles);
   }
   PROCESS_SWITCH(ResonanceInitializer, processTrackV0CascMC, "Process for MC", false);
 
   //  Following the discussions at the PAG meeting (https://indico.cern.ch/event/1583408/)
   //  we have introduced an auxiliary task that, when the resonanceInitializer.cxx is used,
-  // Only consider N_rec / N_gen i.e. not consider level of N_gen at least once
   void processMCgen(soa::Join<aod::McCollisions, aod::McCentFT0Ms, aod::MultMCExtras>::iterator const& mcCollision,
                     aod::McParticles const& mcParticles,
                     const soa::SmallGroups<o2::soa::Join<ResoEvents001, aod::McCollisionLabels>>& collisions,

PWGLF/TableProducer/Resonances/resonanceModuleInitializer.cxx

@@ -1034,6 +1034,8 @@ struct ResonanceDaughterInitializer {
                v0.dcapostopv(),
                v0.dcanegtopv(),
                v0.dcav0topv(),
+               0,
+               0,
                v0.mLambda(),
                v0.mAntiLambda(),
                v0.mK0Short(),
@@ -1206,6 +1208,9 @@ struct ResonanceDaughterInitializer {
                     casc.dcaXYCascToPV(),
                     casc.dcaZCascToPV(),
                     casc.sign(),
+                    0,
+                    0,
+                    0,
                     casc.mLambda(),
                     casc.mXi(),
                     casc.v0radius(), casc.cascradius(), casc.x(), casc.y(), casc.z());