[PWGHF] Use RecoDecay for kaon origin classification

Follow RecoDecay::getParticleOrigin with its default arguments and classify nonzero origins as HF, as suggested by Fabrizio. Remove the custom PDG and ancestry helpers without changing the remaining origin routing.

Co-authored-by: Fabrizio Grosa <fabrizio.grosa@cern.ch>

Author: FDUEnrich

PWGCF/FemtoDream/Core/femtoDreamUtils.h

@@ -105,44 +105,6 @@ inline float getMass(int pdgCode)
   return mass;
 }
 
-/// Identify conventional charm/beauty mesons and baryons from their PDG quark digits.
-inline bool isHeavyFlavorHadron(int pdgCode)
-{
-  constexpr int DigitBase = 10;
-  constexpr int HadronLimit = 10000000;
-  const auto code = std::abs(pdgCode);
-  if (code >= HadronLimit || code % DigitBase == 0) {
-    return false;
-  }
-  const auto quark1 = (code / RecoDecay::PdgDivisorBaryon) % DigitBase;
-  const auto quark2 = (code / RecoDecay::PdgDivisorMeson) % DigitBase;
-  const auto quark3 = (code / DigitBase) % DigitBase;
-  // Quarks, diquarks and generator-specific entries are not hadron decays.
-  if (quark1 > kTop || quark2 < kDown || quark2 > kTop || quark3 < kDown || quark3 > kTop) {
-    return false;
-  }
-  return quark1 == kCharm || quark1 == kBottom || quark2 == kCharm || quark2 == kBottom || quark3 == kCharm || quark3 == kBottom;
-}
-
-template <typename T>
-inline bool hasHeavyFlavorHadronAncestor(T const& particles, typename T::iterator const& particle)
-{
-  auto ancestor = particle;
-  std::vector<int64_t> motherIds;
-  while (ancestor.has_mothers()) {
-    motherIds.clear();
-    RecoDecay::getOrderedMotherIds(ancestor, motherIds);
-    if (motherIds.empty()) {
-      break;
-    }
-    ancestor = particles.rawIteratorAt(motherIds.front() - particles.offset());
-    if (isHeavyFlavorHadron(ancestor.pdgCode())) {
-      return true;
-    }
-  }
-  return false;
-}
-
 /// Classify true kaons from the correct collision for the DCA templates.
 template <typename T>
 inline int getKaonOrigin(T const& particles, typename T::iterator const& particle)
@@ -153,7 +115,7 @@ inline int getKaonOrigin(T const& particles, typename T::iterator const& particl
     return kMaterial;
   }
   // HF decay products can be physical primaries. Check their ancestry first.
-  if (hasHeavyFlavorHadronAncestor(particles, particle)) {
+  if (RecoDecay::getParticleOrigin(particles, particle) != RecoDecay::OriginType::None) {
     return kSecondaryFromHF;
   }
   if (particle.isPhysicalPrimary()) {