From 1059b372c69927a16d2ad22032a685d0fe1adbfb Mon Sep 17 00:00:00 2001
From: Ana Marin <a.marin@gsi.de>
Date: Wed, 17 Jun 2026 11:18:42 +0200
Subject: [PATCH] Update PCMPCM ML analysis

---
 PWGEM/PhotonMeson/Core/Pi0EtaToGammaGamma.h   | 19 +++++++++++++++++++
 .../Tasks/Pi0EtaToGammaGammaMCPCMPCMML.cxx    |  2 +-
 2 files changed, 20 insertions(+), 1 deletion(-)

diff --git a/PWGEM/PhotonMeson/Core/Pi0EtaToGammaGamma.h b/PWGEM/PhotonMeson/Core/Pi0EtaToGammaGamma.h
index f022be669ef..2a84548f9bc 100644
--- a/PWGEM/PhotonMeson/Core/Pi0EtaToGammaGamma.h
+++ b/PWGEM/PhotonMeson/Core/Pi0EtaToGammaGamma.h
@@ -978,6 +978,25 @@ struct Pi0EtaToGammaGamma {
               if (std::fabs(v12.Rapidity()) > maxY) {
                 continue;
               }
+              // as photon has mass= 0 e = p
+              float alphaMeson = std::fabs(g1.p() - g2.p()) / (g1.p() + g2.p());
+              float alphaCut = 999.f;
+              switch (static_cast<AlphaMesonCutOption>(cfgAlphaMesonCut.value)) {
+                case AlphaMesonCutOption::Off:
+                  break;
+                case AlphaMesonCutOption::SpecificValue:
+                  alphaCut = cfgAlphaMeson;
+                  break;
+                case AlphaMesonCutOption::PTDependent: {
+                  alphaCut = cfgAlphaMesonA * std::tanh(cfgAlphaMesonB * v12.pt());
+                  break;
+                }
+                default:
+                  LOGF(error, "Invalid option for alpha meson cut. No alpha cut will be applied.");
+              }
+              if (alphaMeson > alphaCut) {
+                continue;
+              }
 
               fRegistry.fill(HIST("Pair/mix/hs"), v12.M(), v12.Pt(), weight);
             }
diff --git a/PWGEM/PhotonMeson/Tasks/Pi0EtaToGammaGammaMCPCMPCMML.cxx b/PWGEM/PhotonMeson/Tasks/Pi0EtaToGammaGammaMCPCMPCMML.cxx
index c5b06741e75..1f2f76d6499 100644
--- a/PWGEM/PhotonMeson/Tasks/Pi0EtaToGammaGammaMCPCMPCMML.cxx
+++ b/PWGEM/PhotonMeson/Tasks/Pi0EtaToGammaGammaMCPCMPCMML.cxx
@@ -26,7 +26,7 @@ using namespace o2::aod;
 using namespace o2::framework;
 using namespace o2::aod::pwgem::photonmeson::photonpair;
 
-using MyV0Photons = o2::soa::Filtered<o2::soa::Join<o2::aod::V0PhotonsKF, o2::aod::V0PhotonsPhiVPsi, o2::aod::V0KFEMEventIds, o2::aod::V0PhotonsKFPrefilterBitDerived>>;
+using MyV0Photons = o2::soa::Filtered<o2::soa::Join<o2::aod::V0PhotonsKF, o2::aod::V0PhotonsPhiVPsi, o2::aod::V0KFEMEventIds, o2::aod::V0PhotonsKFPrefilterBitDerived, o2::aod::V0PhotonOmegaMBWeights>>;
 using MyMCV0Legs = soa::Join<aod::V0Legs, aod::V0LegMCLabels>;
 
 WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)