code debug flag (Ruediger)
authorloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 21 Oct 2013 08:53:46 +0000 (08:53 +0000)
committerloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 21 Oct 2013 08:53:46 +0000 (08:53 +0000)
PWGJE/EMCALJetTasks/AliEmcalJetTask.cxx
PWGJE/EMCALJetTasks/AliEmcalJetTask.h

index 4bb6375..80bb965 100644 (file)
@@ -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);  
     }
   }
 
index 9fc29c1..a71dcee 100644 (file)
@@ -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