[PWGHF] Refactor XiC0 and OmegaC0 QA workflows

PWGHF/Core/DecayChannels.h

@@ -27,7 +27,6 @@ namespace o2::hf_decay
 
 // TODO
 // - HF cascades (Λc+ → p K0short)
-// - HF cascades to LF cascades (Ωc0/Ξc0 → Ξ+ π−, Ξc+ → Ξ+ π− π+)
 // - Σc
 
 using HfDecayChannel = int8_t;
@@ -296,6 +295,21 @@ enum DecayChannelMain : HfDecayChannel {
 };
 } // namespace hf_cand_sigmac
 
+// @brief XiC0 and OmegaC0 candidates: main channels
+namespace hf_cand_casc_bach_qa
+{
+enum DecayChannelMain : HfDecayChannel {
+  // XiC0
+  Xic0ToXiPi = 1,   // Xi- π+ → Λ π- π+ → p π- π- π+
+  Xic0ToOmegaK = 2, // Ω- K+ → Λ K- K+ → p π- K- K+
+  // OmegaC0
+  Omegac0ToOmegaPi = 3,         // Ω- π+ → Λ π- π+ → p π- π- π+
+  Omegac0ToOmegaK = 4,          // Ω- K+ → Λ K- K+ → p π- K- K+
+  Omegac0ToXiPi = 5,            // Xi- π+ → Λ π- π+ → p π- π- π+
+  NChannelsMain = Omegac0ToXiPi // last channel
+};
+} // namespace hf_cand_casc_bach_qa
+
 } // namespace o2::hf_decay
 
 #endif // PWGHF_CORE_DECAYCHANNELS_H_

PWGHF/Core/HfHelper.h

@@ -491,12 +491,84 @@ struct HfHelper {
     return candidate.m(std::array{o2::constants::physics::MassXiMinus, o2::constants::physics::MassPiPlus});
   }
 
+  template <typename T>
+  static auto invMassOmegaczeroToXiPi(const T& candidate)
+  {
+    return candidate.m(std::array{o2::constants::physics::MassXiMinus, o2::constants::physics::MassPiPlus});
+  }
+
   template <typename T>
   static auto invMassOmegaczeroToOmegaPi(const T& candidate)
   {
     return candidate.m(std::array{o2::constants::physics::MassOmegaMinus, o2::constants::physics::MassPiPlus});
   }
 
+  template <typename T>
+  static auto invMassOmegaczeroToOmegaK(const T& candidate)
+  {
+    return candidate.m(std::array{o2::constants::physics::MassOmegaMinus, o2::constants::physics::MassKPlus});
+  }
+
+  template <typename T>
+  static auto yXiczero(const T& candidate)
+  {
+    return candidate.y(o2::constants::physics::MassXiC0);
+  }
+
+  template <typename T>
+  static auto eXiczero(const T& candidate)
+  {
+    return candidate.e(o2::constants::physics::MassXiC0);
+  }
+
+  template <typename T>
+  static auto ctXiczero(const T& candidate)
+  {
+    return candidate.ct(o2::constants::physics::MassXiC0);
+  }
+
+  template <typename T>
+  static auto cosThetaStarXiczeroToXiPi(const T& candidate)
+  {
+    return candidate.cosThetaStar(std::array{o2::constants::physics::MassXiMinus, o2::constants::physics::MassPiPlus}, o2::constants::physics::MassXiC0, 1);
+  }
+
+  template <typename T>
+  static auto yOmegaczero(const T& candidate)
+  {
+    return candidate.y(o2::constants::physics::MassOmegaC0);
+  }
+
+  template <typename T>
+  static auto eOmegaczero(const T& candidate)
+  {
+    return candidate.e(o2::constants::physics::MassOmegaC0);
+  }
+
+  template <typename T>
+  static auto ctOmegaczero(const T& candidate)
+  {
+    return candidate.ct(o2::constants::physics::MassOmegaC0);
+  }
+
+  template <typename T>
+  static auto cosThetaStarOmegaczeroToXiPi(const T& candidate)
+  {
+    return candidate.cosThetaStar(std::array{o2::constants::physics::MassXiMinus, o2::constants::physics::MassPiPlus}, o2::constants::physics::MassOmegaC0, 1);
+  }
+
+  template <typename T>
+  static auto cosThetaStarOmegaczeroToOmegaPi(const T& candidate)
+  {
+    return candidate.cosThetaStar(std::array{o2::constants::physics::MassOmegaMinus, o2::constants::physics::MassPiPlus}, o2::constants::physics::MassOmegaC0, 1);
+  }
+
+  template <typename T>
+  static auto cosThetaStarOmegaczeroToOmegaK(const T& candidate)
+  {
+    return candidate.cosThetaStar(std::array{o2::constants::physics::MassOmegaMinus, o2::constants::physics::MassKPlus}, o2::constants::physics::MassOmegaC0, 1);
+  }
+
   // hf_cand_casc_lf_3prong
 
   template <typename T>