From c6bbbe280ba184c9c79b102f7569bcd75e4eb43e Mon Sep 17 00:00:00 2001
From: Peter Stratmann <peter.stratmann@cern.ch>
Date: Thu, 18 Dec 2025 12:45:52 +0100
Subject: [PATCH] Remove Pythia from hadron-photon correlation task

---
 PWGJE/Tasks/CMakeLists.txt              | 19 +++++--------------
 PWGJE/Tasks/hadronPhotonCorrelation.cxx | 24 ++++++++++++------------
 2 files changed, 17 insertions(+), 26 deletions(-)

diff --git a/PWGJE/Tasks/CMakeLists.txt b/PWGJE/Tasks/CMakeLists.txt
index 403f2c3590a..d341989737c 100644
--- a/PWGJE/Tasks/CMakeLists.txt
+++ b/PWGJE/Tasks/CMakeLists.txt
@@ -58,6 +58,10 @@ o2physics_add_dpl_workflow(photon-charged-trigger-producer
                     SOURCES photonChargedTriggerProducer.cxx
                     PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore
                     COMPONENT_NAME Analysis)
+o2physics_add_dpl_workflow(hadron-photon-correlation
+										SOURCES hadronPhotonCorrelation.cxx
+										PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore
+										COMPONENT_NAME Analysis)
 
 if(FastJet_FOUND)
     o2physics_add_dpl_workflow(jet-background-analysis
@@ -354,17 +358,4 @@ if(FastJet_FOUND)
                         SOURCES substructureDebug.cxx
                         PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::PWGJECore O2Physics::AnalysisCore
                         COMPONENT_NAME Analysis)
-endif()
-
-if(pythia_FOUND)
-  if(TARGET ROOT::EGPythia8)
-    set_target_properties(ROOT::EGPythia8 PROPERTIES
-        INTERFACE_LINK_LIBRARIES "ROOT::Core;ROOT::EG;ROOT::Graf;ROOT::Physics;${pythia_LIBRARY_SHARED}"
-    )
-
-    o2physics_add_dpl_workflow(hadron-photon-correlation
-                        SOURCES hadronPhotonCorrelation.cxx
-                        PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore ROOT::EGPythia8
-                        COMPONENT_NAME Analysis)
-  endif()
-endif()
+endif()
\ No newline at end of file
diff --git a/PWGJE/Tasks/hadronPhotonCorrelation.cxx b/PWGJE/Tasks/hadronPhotonCorrelation.cxx
index f0bfd076d11..f84d69d621e 100644
--- a/PWGJE/Tasks/hadronPhotonCorrelation.cxx
+++ b/PWGJE/Tasks/hadronPhotonCorrelation.cxx
@@ -40,10 +40,10 @@
 #include "Framework/runDataProcessing.h"
 
 #include <TClonesArray.h>
-#include <TLorentzVector.h> // o2-linter: disable= root/lorentz-vector (TLorentzVector is needed for TPythia8Decayer)
+// #include <TLorentzVector.h> // o2-linter: disable= root/lorentz-vector (TLorentzVector is needed for TPythia8Decayer)
 #include <TPDGCode.h>
-#include <TParticle.h>
-#include <TPythia8Decayer.h>
+// #include <TParticle.h>
+// #include <TPythia8Decayer.h>
 
 #include <cmath>
 #include <map>
@@ -112,9 +112,9 @@ struct HadronPhotonCorrelation {
   float ptMinAssoc;
   float ptMaxAssoc;
 
-  TPythia8Decayer* decayer = new TPythia8Decayer;
-  TLorentzVector* motherLV = new TLorentzVector(); // o2-linter: disable= root/lorentz-vector (TLorentzVector is needed for TPythia8Decayer)
-  TClonesArray* decayParticles = new TClonesArray("TParticle", 10);
+  // TPythia8Decayer* decayer = new TPythia8Decayer;
+  // TLorentzVector* motherLV = new TLorentzVector(); // o2-linter: disable= root/lorentz-vector (TLorentzVector is needed for TPythia8Decayer)
+  // TClonesArray* decayParticles = new TClonesArray("TParticle", 10);
 
   HistogramRegistry registry{"histogram registry"};
 
@@ -391,7 +391,7 @@ struct HadronPhotonCorrelation {
   }
 
   // Initialize Pythia8 decay particle
-  bool initDecayParticle(const TParticle* particle)
+  /*bool initDecayParticle(const TParticle* particle)
   {
 
     if (particle->GetMother(0) != 0) {
@@ -411,7 +411,7 @@ struct HadronPhotonCorrelation {
     }
 
     return true;
-  }
+  }*/
 
   // Initialize trigger tracks
   template <typename T>
@@ -991,7 +991,7 @@ struct HadronPhotonCorrelation {
       }
 
       // Use PYTHIA to simulate decay
-      decayer->Init();
+      /*decayer->Init();
       for (int pid = 1; pid <= nHadrons; pid++) {
         TParticlePDG* pdgParticle = nullptr;
 
@@ -1025,7 +1025,7 @@ struct HadronPhotonCorrelation {
             registry.fill(HIST("generated/charged/hCocktailPhotonCorrelGen"), track_trig.pt(), track_assoc.pt(), daughter->Pt(), daughter->Eta() - track_trig.eta(), dphi, pidCodes[pdgParticle->GetName()], collision.weight());
           }
         }
-      }
+      }*/
     }
     registry.fill(HIST("generated/charged/hPionMultGen"), nPions, collision.weight());
   }
@@ -1147,7 +1147,7 @@ struct HadronPhotonCorrelation {
       }
 
       // Use PYTHIA to simulate decay
-      decayer->Init();
+      /*decayer->Init();
       motherLV->SetPtEtaPhiM(track_assoc.pt(), track_assoc.eta(), track_assoc.phi(), pdgParticle->Mass());
       decayer->Decay(track_assoc.pdgCode(), motherLV);
       decayer->ImportParticles(decayParticles);
@@ -1167,7 +1167,7 @@ struct HadronPhotonCorrelation {
           float dphi = RecoDecay::constrainAngle(daughter->Phi() - track_trig.phi(), -PIHalf);
           registry.fill(HIST("generated/neutral/hCocktailPhotonCorrelGen"), track_trig.pt(), track_assoc.pt(), daughter->Pt(), daughter->Eta() - track_trig.eta(), dphi, pidCodes[pdgParticle->GetName()], collision.weight());
         }
-      }
+      }*/
     }
     registry.fill(HIST("generated/neutral/hNeutralMultGen"), nNeutrals, collision.weight());
   }