[PWGCF,PWGHF] added table for sys studies (#11102)

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

Author: singhra1994

PWGCF/DataModel/FemtoDerived.h

@@ -323,6 +323,13 @@ DECLARE_SOA_TABLE(FDResultsHF, "AOD", "FDRESULTSHF", //! table to store results
                   fdhf::FlagMc,
                   fdhf::OriginMcRec);
 
+DECLARE_SOA_TABLE(FDResultsHFTrkInfo, "AOD", "FDRESULTSHFTRKINFO", //! table to store results for HF femtoscopy
+                  femtodreamparticle::TPCNClsFound,
+                  track::TPCNClsFindable,
+                  femtodreamparticle::TPCNClsCrossedRows,
+                  femtodreamparticle::TPCNSigmaPr,
+                  femtodreamparticle::TOFNSigmaPr);
+
 DECLARE_SOA_TABLE(FDHfCandMC, "AOD", "FDHFCANDMC", //! Table for reconstructed MC charm hadron candidates
                   o2::soa::Index<>,
                   fdhf::FlagMc,

PWGHF/HFC/TableProducer/femtoDreamProducer.cxx

@@ -434,7 +434,8 @@ struct HfFemtoDreamProducer {
     bool isSelectedMlLcToPKPi = true;
     bool isSelectedMlLcToPiKP = true;
     for (const auto& candidate : candidates) {
-
+      outputMlPKPi = {-1.0f, -1.0f, -1.0f};
+      outputMlPiKP = {-1.0f, -1.0f, -1.0f};
       auto trackPos1 = candidate.template prong0_as<TrackType>(); // positive daughter (negative for the antiparticles)
       auto trackNeg = candidate.template prong1_as<TrackType>();  // negative daughter (positive for the antiparticles)
       auto trackPos2 = candidate.template prong2_as<TrackType>(); // positive daughter (negative for the antiparticles)

PWGHF/HFC/Tasks/taskCharmHadronsFemtoDream.cxx

@@ -196,6 +196,7 @@ struct HfTaskCharmHadronsFemtoDream {
   SliceCache cache;
   Preslice<aod::FDParticles> perCol = aod::femtodreamparticle::fdCollisionId;
   Produces<o2::aod::FDResultsHF> fillFemtoResult;
+  Produces<o2::aod::FDResultsHFTrkInfo> fillFemtoResultTrkInfo;
 
   void init(InitContext& /*context*/)
   {
@@ -242,6 +243,17 @@ struct HfTaskCharmHadronsFemtoDream {
     registryMixQa.fill(HIST("MixingQA/hSECollisionPool"), col.posZ(), col.multNtr());
   }
 
+  template <typename Part1>
+  void fillTableTrkInfo(Part1 const& p1)
+  {
+    fillFemtoResultTrkInfo(
+      p1.tpcNClsFound(),
+      p1.tpcNClsFindable(),
+      p1.tpcNClsCrossedRows(),
+      p1.tpcNSigmaPr(),
+      p1.tofNSigmaPr());
+  }
+
   /// This function processes the same event and takes care of all the histogramming
   template <bool isMc, typename PartitionType, typename CandType, typename TableTracks, typename Collision>
   void doSameEvent(PartitionType& sliceTrk1, CandType& sliceCharmHad, TableTracks const& parts, Collision const& col)
@@ -331,6 +343,7 @@ struct HfTaskCharmHadronsFemtoDream {
         processType,
         charmHadMc,
         originType);
+      fillTableTrkInfo(p1);
 
       sameEventCont.setPair<isMc, true>(p1, p2, col.multNtr(), col.multV0M(), use4D, extendedPlots, smearingByOrigin);
     }
@@ -427,6 +440,7 @@ struct HfTaskCharmHadronsFemtoDream {
           processType,
           charmHadMc,
           originType);
+        fillTableTrkInfo(p1);
 
         // if constexpr (!isMc) mixedEventCont.setPair<isMc, true>(p1, p2, collision1.multNtr(), collision1.multV0M(), use4D, extendedPlots, smearingByOrigin);
         mixedEventCont.setPair<isMc, true>(p1, p2, collision1.multNtr(), collision1.multV0M(), use4D, extendedPlots, smearingByOrigin);