From ffdca6aaa956592b03efaa7ca0758d5b2399fc83 Mon Sep 17 00:00:00 2001
From: loizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Mon, 21 Oct 2013 08:53:46 +0000
Subject: [PATCH] code debug flag (Ruediger)

---
 PWGJE/EMCALJetTasks/AliEmcalJetTask.cxx | 9 +++++++--
 PWGJE/EMCALJetTasks/AliEmcalJetTask.h   | 2 ++
 2 files changed, 9 insertions(+), 2 deletions(-)

diff --git a/PWGJE/EMCALJetTasks/AliEmcalJetTask.cxx b/PWGJE/EMCALJetTasks/AliEmcalJetTask.cxx
index 4bb637582f6..80bb965825c 100644
--- a/PWGJE/EMCALJetTasks/AliEmcalJetTask.cxx
+++ b/PWGJE/EMCALJetTasks/AliEmcalJetTask.cxx
@@ -63,6 +63,7 @@ AliEmcalJetTask::AliEmcalJetTask() :
   fIsMcPart(0),
   fIsEmcPart(0),
   fLegacyMode(kFALSE),
+  fCodeDebug(kFALSE),
   fJets(0),
   fEvent(0),
   fTracks(0),
@@ -103,6 +104,7 @@ AliEmcalJetTask::AliEmcalJetTask(const char *name) :
   fIsMcPart(0),
   fIsEmcPart(0),
   fLegacyMode(false),
+  fCodeDebug(kFALSE),
   fJets(0),
   fEvent(0),
   fTracks(0),
@@ -183,7 +185,7 @@ void AliEmcalJetTask::FindJets()
   fjw.SetGhostArea(fGhostArea);
   fjw.SetR(fRadius);
   fjw.SetAlgorithm(jalgo);
-  if(fRecombScheme>0)
+  if(fRecombScheme>=0)
     fjw.SetRecombScheme(static_cast<fastjet::RecombinationScheme>(fRecombScheme)); 
   fjw.SetMaxRap(fEtaMax);
   fjw.Clear();
@@ -247,7 +249,10 @@ void AliEmcalJetTask::FindJets()
 
       // offset of 100 for consistency with cluster ids
       AliDebug(2,Form("Track %d accepted (label = %d, pt = %f)", iTracks, t->GetLabel(), t->Pt()));
-      fjw.AddInputVector(t->Px(), t->Py(), t->Pz(), t->P(), iTracks + 100);  
+      if(fCodeDebug)
+        fjw.AddInputVector(t->Px(), t->Py(), t->Pz(), t->E(), iTracks + 100);
+      else
+        fjw.AddInputVector(t->Px(), t->Py(), t->Pz(), t->P(), iTracks + 100);  
     }
   }
 
diff --git a/PWGJE/EMCALJetTasks/AliEmcalJetTask.h b/PWGJE/EMCALJetTasks/AliEmcalJetTask.h
index 9fc29c175ce..a71dceec12b 100644
--- a/PWGJE/EMCALJetTasks/AliEmcalJetTask.h
+++ b/PWGJE/EMCALJetTasks/AliEmcalJetTask.h
@@ -77,6 +77,7 @@ class AliEmcalJetTask : public AliAnalysisTaskSE {
     }
   }
   void                   SetLegacyMode(Bool_t mode)       { fLegacyMode ^= mode; }
+  void                   SetCodeDebug(Bool_t val)         { fCodeDebug = val; }
 
   UInt_t                 GetJetType()                     { return fJetType; }
   Bool_t                 GetLegacyMode()                  { return fLegacyMode; }
@@ -115,6 +116,7 @@ class AliEmcalJetTask : public AliAnalysisTaskSE {
   Bool_t                 fIsMcPart;               //!=true if MC particles are given as input
   Bool_t                 fIsEmcPart;              //!=true if emcal particles are given as input (for clusters)
   Bool_t                 fLegacyMode;             //! if true, enable FJ 2.x behavior
+  Bool_t                 fCodeDebug;              // use nontested code changes 
   TClonesArray          *fJets;                   //!jet collection
   AliVEvent             *fEvent;                  //!current event
   TClonesArray          *fTracks;                 //!tracks collection
-- 
2.43.0