[PWGEM] Update PCMPCM ML analysis (#16705)

Author: amaringarcia

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);
             }

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)