Add daughter particle handling for XicPlus and Xic0

Author: singhra1994

PWGHF/HFC/TableProducer/correlatorXicHadrons.cxx

@@ -1571,6 +1571,20 @@ struct HfCorrelatorXicHadrons {
       listDaughters.clear();
       const std::size_t nDaughtersExpected = IsXicPlus ? XicDecayDaughtersCount::XicPlusDaughtersCount : XicDecayDaughtersCount::Xic0DaughtersCount;
 
+      if (IsXicPlus) {
+        // Final state: p, pi-, pi-, pi+, pi+
+        std::array<int, XicDecayDaughtersCount::XicPlusDaughtersCount> const arrDaughXicPlusPDG = {kProton, kPiMinus, kPiMinus, kPiPlus, kPiPlus};
+        
+        // Pass -1 to automatically traverse the entire cascade down to the final state
+        RecoDecay::getDaughters(particle, &listDaughters, arrDaughXicPlusPDG);
+      } else {
+        // Final state: p, pi-, pi-, pi+
+        std::array<int, XicDecayDaughtersCount::Xic0DaughtersCount> const arrDaughXic0PDG = {kProton, kPiMinus, kPiMinus, kPiPlus};
+        
+        // Pass -1 to automatically traverse the entire cascade down to the final state
+        RecoDecay::getDaughters(particle, &listDaughters, arrDaughXic0PDG);
+      }
+
       int counterDaughters = 0;
       std::vector<int> prongsId(nDaughtersExpected);
       if (listDaughters.size() == nDaughtersExpected) {