From 5cad9e80626b731c99535ec27f43a6c8aff49772 Mon Sep 17 00:00:00 2001
From: feisenhu <eisenhut@ikf.uni-frankfurt.de>
Date: Mon, 13 May 2024 13:24:19 +0200
Subject: [PATCH] [PWGEM] MCtemplates: add a boolean to run MCGenPairs

---
 PWGDQ/Core/HistogramsLibrary.cxx     | 2 +-
 PWGEM/Dilepton/Tasks/MCtemplates.cxx | 7 ++++---
 2 files changed, 5 insertions(+), 4 deletions(-)

diff --git a/PWGDQ/Core/HistogramsLibrary.cxx b/PWGDQ/Core/HistogramsLibrary.cxx
index 2a378fc80f4..80127198f11 100644
--- a/PWGDQ/Core/HistogramsLibrary.cxx
+++ b/PWGDQ/Core/HistogramsLibrary.cxx
@@ -972,7 +972,7 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
         int nbins_ptee = sizeof(ptee_bins) / sizeof(*ptee_bins) - 1;
 
         // binning for phiv:
-        // every 0.2 GeV/c from 0 to 10 GeV/c
+        // steps of size pi/100
         double phiv_bins[101];
         for (int i = 0; i <= 100; i++)
           phiv_bins[i] = TMath::Pi() / 100. * i;
diff --git a/PWGEM/Dilepton/Tasks/MCtemplates.cxx b/PWGEM/Dilepton/Tasks/MCtemplates.cxx
index b6e06ddd446..0fde7ec3508 100644
--- a/PWGEM/Dilepton/Tasks/MCtemplates.cxx
+++ b/PWGEM/Dilepton/Tasks/MCtemplates.cxx
@@ -414,6 +414,7 @@ struct AnalysisSameEventPairing {
   Configurable<std::string> fConfigMCRecSignals{"cfgBarrelMCRecSignals", "", "Comma separated list of MC signals (reconstructed)"};
   Configurable<std::string> fConfigMCGenSignals{"cfgBarrelMCGenSignals", "", "Comma separated list of MC signals (generated)"};
   Configurable<std::string> fConfigAddSEPHistogram{"cfgAddSEPHistogram", "", "Comma separated list of histograms"};
+  Configurable<bool> fConfigRunMCGenPair{"cfgRunMCGenPair", false, "Do pairing of true MC particles"};
   Configurable<bool> fPropToPCA{"cfgPropToPCA", false, "Propagate tracks to secondary vertex"};
   Configurable<bool> fConfigDoSecVtxProp{"cfgDoSecVtxProp", false, "Propagate tracks to secondary vertex"};
   // TODO: here we specify signals, however signal decisions are precomputed and stored in mcReducedFlags
@@ -627,7 +628,7 @@ struct AnalysisSameEventPairing {
     runPairing<VarManager::kDecayToEE, gkTrackFillMap>(tracks, tracks);
     auto groupedMCTracks = tracksMC.sliceBy(perReducedMcEvent, event.reducedMCevent().globalIndex());
     groupedMCTracks.bindInternalIndicesTo(&tracksMC);
-    runMCGenPair(groupedMCTracks);
+    if(fConfigRunMCGenPair) runMCGenPair(groupedMCTracks);
   }
 
   void processDecayToEESkimmedWithCov(soa::Filtered<MyEventsVtxCovSelected>::iterator const& event,
@@ -642,7 +643,7 @@ struct AnalysisSameEventPairing {
     runPairing<VarManager::kDecayToEE, gkTrackFillMapWithCov>(tracks, tracks);
     auto groupedMCTracks = tracksMC.sliceBy(perReducedMcEvent, event.reducedMCevent().globalIndex());
     groupedMCTracks.bindInternalIndicesTo(&tracksMC);
-    runMCGenPair(groupedMCTracks);
+    if(fConfigRunMCGenPair) runMCGenPair(groupedMCTracks);
   }
 
   void processDecayToEEAOD(soa::Filtered<MyEventsSelectedAOD>::iterator const& event,
@@ -657,7 +658,7 @@ struct AnalysisSameEventPairing {
     runPairing<VarManager::kDecayToEE, gkTrackFillMapAOD>(tracks, tracks);
     auto groupedMCTracks = tracksMC.sliceBy(perMcCollision, event.mcCollision().globalIndex());
     groupedMCTracks.bindInternalIndicesTo(&tracksMC);
-    runMCGenPair(groupedMCTracks);
+    if(fConfigRunMCGenPair) runMCGenPair(groupedMCTracks);
   }
 
   void processDummy(MyEvents&)