From ecd476733e979a9c8af6d558683db9389c8165bb Mon Sep 17 00:00:00 2001
From: mcosenti <mcosenti@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Fri, 14 Sep 2012 12:41:05 +0000
Subject: [PATCH] including setter to enable debugging printouts

---
 .../AliAnalysisTaskEMCALIsoPhoton.cxx         | 29 ++++++++++++++++---
 .../AliAnalysisTaskEMCALIsoPhoton.h           |  2 ++
 2 files changed, 27 insertions(+), 4 deletions(-)

diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx
index 15f3641dd7a..be0f947ff9a 100644
--- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx
+++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx
@@ -21,11 +21,8 @@
 #include "AliESDHeader.h"
 #include "AliESDInputHandler.h"
 #include "AliESDUtils.h"
-#include "AliESDpid.h"
 #include "AliESDtrack.h"
 #include "AliESDtrackCuts.h"
-#include "AliESDv0.h"
-#include "AliKFParticle.h"
 #include "AliMCEvent.h"
 #include "AliMCEventHandler.h"
 #include "AliStack.h"
@@ -48,6 +45,7 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton() :
   fTrigBit("kEMC7"),
   fIsTrain(0),
   fIsMc(0),
+  fDebug(0),
   fExoticCut(0.97),
   fIsoConeR(0.4),
   fNDimensions(7),
@@ -83,6 +81,7 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton(const char *name) :
   fTrigBit("kEMC7"),
   fIsTrain(0),
   fIsMc(0),
+  fDebug(0),
   fExoticCut(0.97),
   fIsoConeR(0.4),
   fNDimensions(7),
@@ -183,6 +182,8 @@ void AliAnalysisTaskEMCALIsoPhoton::UserExec(Option_t *)
   }
   if(fIsMc)
     isSelected = kTRUE;
+  if(fDebug)
+    printf("isSelected = %d, fIsMC=%d\n", isSelected, fIsMc);
   if(!isSelected )
         return; 
 
@@ -193,6 +194,8 @@ void AliAnalysisTaskEMCALIsoPhoton::UserExec(Option_t *)
   }
   
   fEvtSel->Fill(0);
+  if(fDebug)
+    printf("fESD is ok\n");
   
   AliESDVertex *pv = (AliESDVertex*)fESD->GetPrimaryVertex();
   if(!pv) 
@@ -200,6 +203,8 @@ void AliAnalysisTaskEMCALIsoPhoton::UserExec(Option_t *)
   fPVtxZ->Fill(pv->GetZ());
   if(TMath::Abs(pv->GetZ())>15)
     return;
+  if(fDebug)
+    printf("passed vertex cut\n");
 
   fEvtSel->Fill(1);
 
@@ -234,14 +239,26 @@ void AliAnalysisTaskEMCALIsoPhoton::UserExec(Option_t *)
   
     
   FillClusHists(); 
+  if(fDebug)
+    printf("passed calling of FillClusHists\n");
 
   fCaloClusters->Clear();
   fSelPrimTracks->Clear();
 
   fMCEvent = MCEvent();
-  if(fMCEvent)
+  if(fMCEvent){
+    if(fDebug)
+      cout<<"MCevent exists\n";
     fStack = (AliStack*)fMCEvent->Stack();
+  }
+  else{
+    if(fDebug)
+      cout<<"ERROR: NO MC EVENT!!!!!!\n";
+  }
   FillMcHists();
+  if(fDebug)
+    printf("passed calling of FillMcHists\n");
+
 
   PostData(1, fOutputList);
 }      
@@ -256,6 +273,8 @@ void AliAnalysisTaskEMCALIsoPhoton::FillClusHists()
   const Int_t nclus = fCaloClusters->GetEntries();
   if(nclus==0)
     return;
+  if(fDebug)
+    printf("Inside FillClusHists and there are %d clusters\n",nclus);
   Double_t maxE = 0;
   Int_t nclus10 = 0;
   for(Int_t ic=0;ic<nclus;ic++){
@@ -477,6 +496,8 @@ void AliAnalysisTaskEMCALIsoPhoton ::FillMcHists()
   if(!fStack)
     return;
   Int_t nTracks = fStack->GetNtrack();
+  if(fDebug)
+    printf("Inside FillMcHists and there are %d mcparts\n",nTracks);
   for(Int_t iTrack=0;iTrack<nTracks;iTrack++){
     TParticle *mcp = static_cast<TParticle*>(fStack->Particle(iTrack));  
     if(!mcp)
diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h
index 6b28a8f2425..bdedcc449b9 100644
--- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h
+++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h
@@ -44,6 +44,7 @@ class AliAnalysisTaskEMCALIsoPhoton : public AliAnalysisTaskSE {
   void                   SetPrimTrackCuts(AliESDtrackCuts *c)   { fPrTrCuts           = c;       }
   void                   SetTrainMode(Bool_t t)                 { fIsTrain            = t;       }
   void                   SetMcMode(Bool_t mc)                   { fIsMc               = mc;      }
+  void                   SetDebugOn(Bool_t d)                   { fDebug              = d;       }
   
  protected:
   TRefArray             *fCaloClusters;          //!pointer to EMCal clusters
@@ -57,6 +58,7 @@ class AliAnalysisTaskEMCALIsoPhoton : public AliAnalysisTaskSE {
   TString                fTrigBit;               // string to the trigger bit name
   Bool_t                 fIsTrain;               // variable to set train mode
   Bool_t                 fIsMc;                  // variable to set mc mode
+  Bool_t                 fDebug;                 // variable to set on/off debugging printouts
   Double_t               fExoticCut;             // variable to set the cut on exotic clusters
   Double_t               fIsoConeR;              // variable to set the isolation cone radius
   Int_t                  fNDimensions;           // variable to set the number of dimensions of n-sparse
-- 
2.39.3