replace printf's by AliDebug/AliInfo/AliWarning
authorgconesab <gustavo.conesa.balbastre@cern.ch>
Tue, 21 Oct 2014 17:40:49 +0000 (19:40 +0200)
committergconesab <gustavo.conesa.balbastre@cern.ch>
Tue, 21 Oct 2014 17:41:13 +0000 (19:41 +0200)
15 files changed:
PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrBaseClass.cxx
PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrBaseClass.h
PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.cxx
PWG/CaloTrackCorrBase/AliAnalysisTaskCaloTrackCorrelation.cxx
PWG/CaloTrackCorrBase/AliAnalysisTaskCaloTrackCorrelationM.cxx
PWG/CaloTrackCorrBase/AliAnalysisTaskCounter.cxx
PWG/CaloTrackCorrBase/AliCaloPID.cxx
PWG/CaloTrackCorrBase/AliCaloTrackAODReader.cxx
PWG/CaloTrackCorrBase/AliCaloTrackESDReader.cxx
PWG/CaloTrackCorrBase/AliCaloTrackMCReader.cxx
PWG/CaloTrackCorrBase/AliCaloTrackReader.cxx
PWG/CaloTrackCorrBase/AliCalorimeterUtils.cxx
PWG/CaloTrackCorrBase/AliIsolationCut.cxx
PWG/CaloTrackCorrBase/AliIsolationCut.h
PWG/CaloTrackCorrBase/AliMCAnalysisUtils.cxx

index ad383ba..0df5034 100755 (executable)
@@ -37,7 +37,6 @@
 #include "AliVCaloCells.h" 
 #include "AliAODEvent.h"
 #include "AliAODHandler.h"
-#include "AliAnalysisManager.h"
 #include "AliAODPWG4Particle.h"
 
 ClassImp(AliAnaCaloTrackCorrBaseClass)
@@ -151,7 +150,7 @@ void AliAnaCaloTrackCorrBaseClass::ConnectInputOutputAODBranches()
   //Recover ouput and input AOD pointers for each event in the maker
        
   //Delta AODs
-  if(fDebug > 3) printf("AliAnaCaloTrackCorrBaseClass::ConnectInputOutputAODBranches() - Connect Input with name: <%s>; Connect output with name <%s>\n",fInputAODName.Data(),fOutputAODName.Data());
+  AliDebug(3,Form("AliAnaCaloTrackCorrBaseClass::ConnectInputOutputAODBranches() - Connect Input with name: <%s>; Connect output with name <%s>\n",fInputAODName.Data(),fOutputAODName.Data()));
   
   //Get the AOD handler, if output AOD is created use it, if not get the branches from the input which should be deltaAODs
   AliAODHandler* aodHandler = 0x0;
@@ -165,46 +164,54 @@ void AliAnaCaloTrackCorrBaseClass::ConnectInputOutputAODBranches()
     fOutputAODBranch =  (TClonesArray *) (fReader->GetAODBranchList())->FindObject(fOutputAODName);
     fInputAODBranch  =  (TClonesArray *) (fReader->GetAODBranchList())->FindObject(fInputAODName);     
   }
-  else if (aodHandler->GetExtensions()) { 
-    
+  else if (aodHandler->GetExtensions())
+  {
     AliAODExtension *ext = (AliAODExtension*)aodHandler->GetExtensions()->FindObject(GetReader()->GetDeltaAODFileName()); 
-    if(ext){
+    if(ext)
+    {
       AliAODEvent *aodEvent = ext->GetAOD(); 
       if(fNewAOD)fOutputAODBranch = (TClonesArray*) aodEvent->FindListObject(fOutputAODName);
       fInputAODBranch = (TClonesArray*) aodEvent->FindListObject(fInputAODName);         
       if(!fOutputAODBranch && fNewAOD) fOutputAODBranch =  (TClonesArray *) fReader->GetOutputEvent()->FindListObject(fOutputAODName);
       if(!fInputAODBranch)  fInputAODBranch  =  (TClonesArray *) fReader->GetOutputEvent()->FindListObject(fInputAODName);
     }
-    else{//If no Delta AODs, kept in standard branch, to revise. 
-      if(fNewAOD && fReader->GetOutputEvent()) {
+    else
+    { // If no Delta AODs, kept in standard branch, to revise.
+      if(fNewAOD && fReader->GetOutputEvent())
+      {
         fOutputAODBranch =  (TClonesArray *) fReader->GetOutputEvent()->FindListObject(fOutputAODName);
         fInputAODBranch  =  (TClonesArray *) fReader->GetOutputEvent()->FindListObject(fInputAODName); 
       }
-      else {
+      else
+      {
         fInputAODBranch  =  (TClonesArray *) fReader->GetInputEvent()->FindListObject(fInputAODName);  
         if(!fInputAODBranch && fReader->GetOutputEvent() ) 
           fInputAODBranch  =  (TClonesArray *) fReader->GetOutputEvent()->FindListObject(fInputAODName);//Try the output event.
       }
     }
   }
-  else{ //If no Delta AODs, kept in standard branch
-    if(fNewAOD && fReader->GetOutputEvent()) {
+  else
+  { // If no Delta AODs, kept in standard branch
+    if(fNewAOD && fReader->GetOutputEvent())
+    {
       fOutputAODBranch =  (TClonesArray *) fReader->GetOutputEvent()->FindListObject(fOutputAODName);
       fInputAODBranch  =  (TClonesArray *) fReader->GetOutputEvent()->FindListObject(fInputAODName);   
     }
-    else{ 
+    else
+    {
       fInputAODBranch  =  (TClonesArray *) fReader->GetInputEvent()->FindListObject(fInputAODName);
       if(!fInputAODBranch && fReader->GetOutputEvent())  
         fInputAODBranch  =  (TClonesArray *) fReader->GetOutputEvent()->FindListObject(fInputAODName);//Try the output event.
     }
   }
   
-  if(GetDebug() > 1){
-    if(fNewAOD && !fOutputAODBranch) 
-      printf(" AliAnaCaloTrackCorrBaseClass::ConnectInputOutputAODBranches() - Output Branch <%s>, not found!\n",fOutputAODName.Data());
-    if(!fNewAOD && !fInputAODBranch) 
-      printf(" AliAnaCaloTrackCorrBaseClass::ConnectInputOutputAODBranches() - Input Branch  <%s>, not found!\n",fInputAODName.Data());
-  }
+//  if(GetDebug() > 1)
+//  {
+//    if(fNewAOD && !fOutputAODBranch) 
+//      AliInfo(Form("Output Branch <%s>, not found!\n",fOutputAODName.Data()));
+//    if(!fNewAOD && !fInputAODBranch) 
+//      AliInfo(Form("Input Branch  <%s>, not found!\n",fInputAODName.Data()));
+//  }
 }
 
 //_____________________________________________________________________________________
@@ -237,7 +244,8 @@ TClonesArray * AliAnaCaloTrackCorrBaseClass::GetAODBranch(const TString & aodNam
        //Recover ouput and input AOD pointers for each event in the maker
        
        //Delta AODs
-       if(fDebug > 3) printf("AliAnaCaloTrackCorrBaseClass::GetAODBranch() - Get Input Branch with name: <%s>; \n",aodName.Data());
+  
+       AliDebug(3,Form("AliAnaCaloTrackCorrBaseClass::GetAODBranch() - Get Input Branch with name: <%s>; \n",aodName.Data()));
        
   //Get the AOD handler, if output AOD is created use it, if not get the branches from the input which should be deltaAODs
   AliAODHandler* aodHandler = 0x0;
@@ -347,7 +355,8 @@ TString  AliAnaCaloTrackCorrBaseClass::GetBaseParametersList()
   parList+=onePar ;
   snprintf(onePar,buffersize,"fInputAODName  =%s Input AOD name \n",fInputAODName.Data()) ;
   parList+=onePar ;    
-  if(fNewAOD){
+  if(fNewAOD)
+  {
     snprintf(onePar,buffersize,"fOutputAODName  =%s Output AOD name \n",fOutputAODName.Data()) ;
     parList+=onePar ;  
     snprintf(onePar,buffersize,"fOutputAODClassName  =%s Output AOD class name \n",fOutputAODClassName.Data()) ;
@@ -367,8 +376,8 @@ TClonesArray * AliAnaCaloTrackCorrBaseClass::GetCreateOutputAODBranch()
 {
   //Create AOD branch filled in the analysis
   
-  printf("Create AOD branch of %s objects and with name < %s >\n",
-         fOutputAODClassName.Data(),fOutputAODName.Data()) ;
+  AliInfo(Form("Create AOD branch of %s objects and with name < %s >\n",
+          fOutputAODClassName.Data(),fOutputAODName.Data())) ;
   
   TClonesArray * aodBranch = new TClonesArray(fOutputAODClassName, 0);
   aodBranch->SetName(fOutputAODName);
@@ -427,7 +436,7 @@ Int_t AliAnaCaloTrackCorrBaseClass::GetTrackMultiplicityBin() const
     if(trackMult >= fTrackMultBins[ibin] && trackMult < fTrackMultBins[ibin+1]) return ibin;
   }
   
-  printf("AliAnaCaloTrackCorrBaseClass::GetTrackMultiplicityBin() - Bin not found for track multiplicity %d\n",trackMult);
+  AliWarning(Form("Bin not found for track multiplicity %d",trackMult));
   
   return -1;
 }
@@ -455,7 +464,7 @@ Int_t AliAnaCaloTrackCorrBaseClass::GetEventCentralityBin() const
       if(curCentrBin==GetNCentrBin())
       {
         curCentrBin = GetNCentrBin()-1;
-        printf("AliAnaCaloTrackCorrBaseClass::GetEventCentralityBin() - Centrality = %d, put it in last bin \n",GetEventCentrality());
+        AliDebug(1,Form("Centrality = %d, put it in last bin \n",GetEventCentrality()));
       }
     }
     else
@@ -464,9 +473,8 @@ Int_t AliAnaCaloTrackCorrBaseClass::GetEventCentralityBin() const
       if(curCentrBin==GetNCentrBin()) curCentrBin = GetNCentrBin()-1;
     }
     
-    if(GetDebug() > 0 )
-      printf("AliAnaCaloTrackCorrBaseClass::GetEventCentralityBin() - %d, centrality %d, n bins %d, max bin from centrality %d\n",
-             curCentrBin, GetEventCentrality(), GetNCentrBin(), GetReader()->GetCentralityOpt());
+    AliDebug(1,Form("Current CentrBin %d, centrality %d, n bins %d, max bin from centrality %d",
+                    curCentrBin, GetEventCentrality(), GetNCentrBin(), GetReader()->GetCentralityOpt()));
   }
   
   return curCentrBin;
@@ -486,16 +494,15 @@ Int_t AliAnaCaloTrackCorrBaseClass::GetEventRPBin() const
     
     if(epAngle < 0 || epAngle >TMath::Pi())
     { 
-      printf("AliAnaCaloTrackCorrBaseClass::GetEventRPBin() - Wrong event plane angle : %f \n",epAngle);
+      AliWarning(Form("Wrong event plane angle : %f \n",epAngle));
       return -1;
     }
     
     curRPBin = TMath::Nint(epAngle*(GetNRPBin()-1)/TMath::Pi());
-    if(curRPBin >= GetNRPBin()) printf("RP Bin %d out of range %d\n",curRPBin,GetNRPBin());
+    if(curRPBin >= GetNRPBin()) printf("RP Bin %d out of range %d",curRPBin,GetNRPBin());
     
-    if(GetDebug() > 0 )
-      printf("AliAnaCaloTrackCorrBaseClass::GetEventRPBin() - %d, bin float %f, angle %f, n bins %d\n", 
-             curRPBin,epAngle*(GetNRPBin()-1)/TMath::Pi(),epAngle,GetNRPBin());
+    AliDebug(1,Form("Current RP bin %d, bin float %f, angle %f, n bins %d",
+                    curRPBin,epAngle*(GetNRPBin()-1)/TMath::Pi(),epAngle,GetNRPBin()));
   }  
   
   return curRPBin ;
@@ -513,9 +520,8 @@ Int_t AliAnaCaloTrackCorrBaseClass::GetEventVzBin() const
   
   Int_t curZvertBin = (Int_t)(0.5*GetNZvertBin()*(v[2]+GetZvertexCut())/GetZvertexCut());
   
-  if(GetDebug() > 0 )
-    printf("AliAnaCaloTrackCorrBaseClass::GetEventVzBin() - %d, vz %2.2f, n bins %d \n",
-           curZvertBin, v[2], GetNZvertBin()); 
+  AliDebug(1,Form("AliAnaCaloTrackCorrBaseClass::GetEventVzBin() - %d, vz %2.2f, n bins %d",
+                  curZvertBin, v[2], GetNZvertBin()));
   
   return curZvertBin;
 }
@@ -547,20 +553,38 @@ Int_t AliAnaCaloTrackCorrBaseClass::GetEventMixBin() const
   
   Int_t eventBin = GetEventMixBin(iCen, iVz, iRP);
   
-  if(GetDebug() > 0)
-    printf("AliAnaCaloTrackCorrBaseClass::GetEventMixBin() - Bins : cent %d, vz %d, RP %d, event %d/%d\n",
-           iCen,iVz, iRP, eventBin, GetNZvertBin()*GetNRPBin()*GetNCentrBin());
+  AliDebug(1,Form("Bins : cent %d, vz %d, RP %d, event %d/%d",
+                  iCen,iVz, iRP, eventBin, GetNZvertBin()*GetNRPBin()*GetNCentrBin()));
   
   return eventBin;
   
 }
 
+//____________________________________________
+void AliAnaCaloTrackCorrBaseClass::InitDebug()
+{
+  // Init once the debugging level, if requested
+  
+  // Activate debug level in analysis
+
+  printf("Debug levels: Ana %d, MC %d, Iso %d\n",fDebug,GetMCAnalysisUtils()->GetDebug(),GetIsolationCut()->GetDebug());
+
+  if( fDebug >= 0 )
+    (AliAnalysisManager::GetAnalysisManager())->AddClassDebug(this->ClassName(),fDebug);
+  
+  if( GetMCAnalysisUtils()->GetDebug() >= 0 )
+    (AliAnalysisManager::GetAnalysisManager())->AddClassDebug(GetMCAnalysisUtils()->ClassName(),GetMCAnalysisUtils()->GetDebug());
+  
+  if( GetIsolationCut()->GetDebug() >= 0 )
+    (AliAnalysisManager::GetAnalysisManager())->AddClassDebug(GetIsolationCut()   ->ClassName(),GetIsolationCut()->GetDebug());
+}
+
 //_________________________________________________
 void AliAnaCaloTrackCorrBaseClass::InitParameters()
 { 
   //Initialize the parameters of the analysis.
   fDataMC              = kFALSE;
-  fDebug               = -1;
+  fDebug               = 0;
   fCheckCaloPID        = kTRUE ;
   fCheckFidCut         = kFALSE ;
   fCheckRealCaloAcc    = kFALSE ;
index 47a0929..1880127 100755 (executable)
@@ -39,6 +39,8 @@ class AliEMCALGeometry;
 class AliPHOSGeoUtils;
 class AliCentrality;
 class AliEventplane;
+#include "AliAnalysisManager.h"
+#include "AliLog.h"
 
 //jets
 class AliAODJetEventBackground;
@@ -53,7 +55,8 @@ public:
   
   virtual TList *        GetCreateOutputObjects()               { return (new TList)          ; }
   
-  virtual void           Init() {;}
+  virtual void           Init()                                 { ; }
+  virtual void           InitDebug()      ;
   virtual void           InitParameters() ;
   
   virtual void           FillEventMixPool()                     { ; }
index 3c0d894..fc81bc3 100755 (executable)
@@ -38,6 +38,7 @@
 #include "AliAODEvent.h"
 #include "AliAnaCaloTrackCorrBaseClass.h"
 #include "AliAnaCaloTrackCorrMaker.h"
+#include "AliLog.h"
 
 ClassImp(AliAnaCaloTrackCorrMaker)
 
@@ -76,7 +77,7 @@ fhClusterTriggerBCEventBC(0),           fhClusterTriggerBCEventBCUnMatch(0),
 fhClusterTriggerBCExoticEventBC(0),     fhClusterTriggerBCExoticEventBCUnMatch(0)
 {
   //Default Ctor
-  if(fAnaDebug > 1 ) printf("*** Analysis Maker Constructor *** \n");
+  AliDebug(1,"*** Analysis Maker Constructor ***");
   
   for(Int_t i = 0; i < 3; i++)
   {
@@ -190,8 +191,7 @@ void    AliAnaCaloTrackCorrMaker::AddAnalysis(TObject* ana, Int_t n)
   }
   else
   {
-    printf("AliAnaCaloTrackCorrMaker::AddAnalysis() - AnalysisContainer not initialized\n");
-    abort();
+    AliFatal("AnalysisContainer not initialized");
   }
 }
 
@@ -369,7 +369,7 @@ void AliAnaCaloTrackCorrMaker::FillTriggerControlHistograms()
   
   Int_t eventBC = fReader->GetInputEvent()->GetBunchCrossNumber();
   if(eventBC%4 < 0 || eventBC%4 > 3 )
-    printf("AliAnaCaloTrackCorrMaker::ProcessEvent() - STRANGE: Trigger BC %d - Event BC %d, modulo4 %d \n",triggerBC,eventBC,eventBC%4);
+    AliWarning(Form("STRANGE: Trigger BC %d - Event BC %d, modulo4 %d",triggerBC,eventBC,eventBC%4));
   
   if(triggerMatch)
   {
@@ -778,7 +778,7 @@ TList *AliAnaCaloTrackCorrMaker::GetOutputContainer()
   
   if(!fAnalysisContainer || fAnalysisContainer->GetEntries()==0)
   {
-    printf("AliAnaCaloTrackCorrMaker::GetOutputContainer() - Analysis job list not initialized!!!\n");
+    AliWarning("Analysis job list not initialized!!!");
     return fOutputContainer;
   }
   
@@ -828,27 +828,42 @@ void AliAnaCaloTrackCorrMaker::Init()
   //Init container histograms and other common variables
   // Fill the output list of histograms during the CreateOutputObjects stage.
   
+  // Activate debug level in maker
+  if( fAnaDebug >= 0 )
+    (AliAnalysisManager::GetAnalysisManager())->AddClassDebug(this->ClassName(),fAnaDebug);
+  
   //Initialize reader
   GetReader()->Init();
   GetReader()->SetCaloUtils(GetCaloUtils()); // pass the calo utils pointer to the reader
-       
+
+  // Activate debug level in reader
+  if( fReader->GetDebug() >= 0 )
+    (AliAnalysisManager::GetAnalysisManager())->AddClassDebug(fReader->ClassName(), fReader->GetDebug());
+  
+  // Activate debug level in calo utils
+  if( fCaloUtils->GetDebug() >= 0 )
+    (AliAnalysisManager::GetAnalysisManager())->AddClassDebug(fCaloUtils->ClassName(), fCaloUtils->GetDebug());
   
   if(!fAnalysisContainer || fAnalysisContainer->GetEntries()==0)
   {
-    printf("AliAnaCaloTrackCorrMaker::GetOutputInit() - Analysis job list not initialized!!!\n");
+    AliWarning("Analysis job list not initialized");
     return;
   }
   
   for(Int_t iana = 0; iana <  fAnalysisContainer->GetEntries(); iana++)
   {
-    
     AliAnaCaloTrackCorrBaseClass * ana =  ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
     
     ana->SetReader(fReader);       // Set Reader for each analysis
     ana->SetCaloUtils(fCaloUtils); // Set CaloUtils for each analysis
     
     ana->Init();
+    ana->InitDebug();
     
+//    // Activate debug level in analysis
+//    if( ana->GetDebug() >= 0 )
+//      (AliAnalysisManager::GetAnalysisManager())->AddClassDebug(ana->ClassName(),ana->GetDebug());
+
   }//Loop on analysis defined
   
 }
@@ -901,20 +916,11 @@ void AliAnaCaloTrackCorrMaker::ProcessEvent(Int_t iEntry, const char * currentFi
   //Process analysis for this event
   
   if(fMakeHisto && !fOutputContainer)
-  {
-    printf("AliAnaCaloTrackCorrMaker::ProcessEvent() - Histograms not initialized\n");
-    abort();
-  }
+    AliFatal("Histograms not initialized");
   
-  if(fAnaDebug >= 0 )
-  {
-    printf("***  AliAnaCaloTrackCorrMaker::ProcessEvent() Event %d   ***  \n",iEntry);
-    if(fAnaDebug > 1 )
-    {
-      printf("AliAnaCaloTrackCorrMaker::ProcessEvent() - Current File Name : %s\n", currentFileName);
-      //printf("fAODBranchList %p, entries %d\n",fAODBranchList,fAODBranchList->GetEntries());
-    }
-  }
+  AliDebug(1,Form("***  AliAnaCaloTrackCorrMaker::ProcessEvent() Event %d   ***",iEntry));
+  AliDebug(2,Form("Current File Name : %s", currentFileName));
+  //printf("fAODBranchList %p, entries %d\n",fAODBranchList,fAODBranchList->GetEntries());
   
   //Each event needs an empty branch
   TList * aodList = fReader->GetAODBranchList();
@@ -956,7 +962,7 @@ void AliAnaCaloTrackCorrMaker::ProcessEvent(Int_t iEntry, const char * currentFi
   
   if(!ok)
   {
-    if(fAnaDebug >= 1 )printf("*** Skip event *** %d \n",iEntry);
+    AliDebug(1,Form("*** Skip event *** %d",iEntry));
     fReader->ResetLists();
     return ;
   }
@@ -973,7 +979,7 @@ void AliAnaCaloTrackCorrMaker::ProcessEvent(Int_t iEntry, const char * currentFi
   
   //Loop on analysis algorithms
   
-  if(fAnaDebug > 0 ) printf("*** Begin analysis *** \n");
+  AliDebug(1,"*** Begin analysis ***");
   
   Int_t nana = fAnalysisContainer->GetEntries() ;
   for(Int_t iana = 0; iana <  nana; iana++)
@@ -1003,7 +1009,7 @@ void AliAnaCaloTrackCorrMaker::ProcessEvent(Int_t iEntry, const char * currentFi
   // do not fill control histograms for a non requested triggered event
   if(!fReader->IsEventTriggerAtSEOn() && !isTrigger)
   {
-    if(fAnaDebug > 0 ) printf("AliAnaCaloTrackMaker::ProcessEvent() - *** End analysis, MB for mixing *** \n");
+    AliDebug(1,"*** End analysis, MB for mixing ***");
     return;
   }
   
@@ -1012,7 +1018,7 @@ void AliAnaCaloTrackCorrMaker::ProcessEvent(Int_t iEntry, const char * currentFi
   //printf(">>>>>>>>>> AFTER >>>>>>>>>>>\n");
   //gObjectTable->Print();
        
-  if(fAnaDebug > 0 ) printf("AliAnaCaloTrackMaker::ProcessEvent() - *** End analysis *** \n");
+  AliDebug(1,"*** End analysis ***");
   
 }
 
@@ -1024,7 +1030,7 @@ void AliAnaCaloTrackCorrMaker::Terminate(TList * outputList)
   
   if (!outputList)
   {
-    Error("Terminate", "No output list");
+    AliError("No output list");
     return;
   }
   
index c00a9d3..5dbe9b0 100755 (executable)
@@ -39,6 +39,7 @@
 #include "AliPDG.h"
 #include "AliAnalysisManager.h"
 #include "AliInputEventHandler.h"
+#include "AliLog.h"
 
 ClassImp(AliAnalysisTaskCaloTrackCorrelation)
 
@@ -88,12 +89,12 @@ AliAnalysisTaskCaloTrackCorrelation::~AliAnalysisTaskCaloTrackCorrelation()
 void AliAnalysisTaskCaloTrackCorrelation::UserCreateOutputObjects()
 {
   // Create the output container
-  if (DebugLevel() > 1) printf("AliAnalysisTaskCaloTrackCorrelation::UserCreateOutputObjects() - Begin\n");
+  AliDebug(1,"Begin");
   
   //Get list of aod arrays, add each aod array to analysis frame
   TList * list = fAna->FillAndGetAODBranchList(); //Loop the analysis and create the list of branches
   
-  if (DebugLevel() >= 1) printf("AliAnalysisTaskCaloTrackCorrelation::UserCreateOutputObjects() - n AOD branches %d\n",list->GetEntries());
+  AliDebug(1,Form("n AOD branches %d",list->GetEntries()));
   
   //Put the delta AODs in output file, std or delta
   if((fAna->GetReader())->WriteDeltaAODToFile())
@@ -111,11 +112,11 @@ void AliAnalysisTaskCaloTrackCorrelation::UserCreateOutputObjects()
   OpenFile(1);
   fOutputContainer = fAna->GetOutputContainer();
   
-  if (DebugLevel() >= 1) printf("AliAnalysisTaskCaloTrackCorrelation::UserCreateOutputObjects() - n histograms %d\n",fOutputContainer->GetEntries());
+  AliDebug(1,Form("n histograms %d",fOutputContainer->GetEntries()));
   
   fOutputContainer->SetOwner(kTRUE);
   
-  if (DebugLevel() > 1) printf("AliAnalysisTaskCaloTrackCorrelation::UserCreateOutputObjects() - End\n");
+  AliDebug(1,"End");
   
   PostData(1,fOutputContainer);
   
@@ -136,21 +137,24 @@ void AliAnalysisTaskCaloTrackCorrelation::Init()
 {
   // Initialization
   
-  if (DebugLevel() > 1) printf("AliAnalysisTaskCaloTrackCorrelation::Init() - Begin\n");
+  AliDebug(1,"Begin");
+  
+  if( fDebug >= 0 )
+    (AliAnalysisManager::GetAnalysisManager())->AddClassDebug(this->ClassName(),fDebug);
   
   // Call configuration file if specified
   
   if (fConfigName.Length())
   {
-    printf("AliAnalysisTaskCaloTrackCorrelation::Init() - ### Configuration file is %s.C ###\n", fConfigName.Data());
+    AliInfo(Form("### Configuration file is %s.C ###", fConfigName.Data()));
     gROOT->LoadMacro(fConfigName+".C");
     fAna = (AliAnaCaloTrackCorrMaker*) gInterpreter->ProcessLine("ConfigAnalysis()");
   }
   
   if(!fAna)
   {
-    printf("AliAnalysisTaskCaloTrackCorrelation::Init() - Analysis maker pointer not initialized, no analysis specified, STOP !\n");
-    abort();
+    AliFatal("Analysis maker pointer not initialized, no analysis specified, STOP!");
+    return; // coverity
   }
   
   // Add different generator particles to PDG Data Base
@@ -170,7 +174,7 @@ void AliAnalysisTaskCaloTrackCorrelation::Init()
   // Selected Trigger
   if(fAna->GetReader()->IsEventTriggerAtSEOn()) fAna->GetReader()->SetEventTriggerMask(GetCollisionCandidates());
   
-  if (DebugLevel() > 1) printf("AliAnalysisTaskCaloTrackCorrelation::Init() - End\n");
+  AliDebug(1,"End");
   
 }
 
@@ -180,14 +184,14 @@ void AliAnalysisTaskCaloTrackCorrelation::UserExec(Option_t */*option*/)
 {
   // Execute analysis for current event
   
-  if (DebugLevel() > 1) printf("AliAnalysisTaskCaloTrackCorrelation::UserExec() - Begin\n");
+  AliDebug(1,"Begin");
   
   //Get the type of data, check if type is correct
   Int_t  datatype = fAna->GetReader()->GetDataType();
   if(datatype != AliCaloTrackReader::kESD && datatype != AliCaloTrackReader::kAOD &&
      datatype != AliCaloTrackReader::kMC)
   {
-    printf("AliAnalysisTaskCaloTrackCorrelation::UserExec() - Wrong type of data\n");
+    AliError("Wrong type of data");
     return ;
   }
   
@@ -197,10 +201,12 @@ void AliAnalysisTaskCaloTrackCorrelation::UserExec(Option_t */*option*/)
   fAna->ProcessEvent((Int_t) Entry(), CurrentFileName());
   
   //printf("AliAnalysisTaskCaloTrackCorrelation::Current Event %d; Current File Name : %s\n",(Int_t) Entry(), CurrentFileName());
-  if (DebugLevel() > 1) printf("AliAnalysisTaskCaloTrackCorrelation::UserExec() - End\n");
+  AliDebug(1,"End");
   
   PostData(1, fOutputContainer);
   
+  AliDebug(1,"End");
+  
   //gObjectTable->Print();
   
 }
@@ -237,8 +243,7 @@ void AliAnalysisTaskCaloTrackCorrelation::FinishTaskOutput()
   TH2F *histBin0 = dynamic_cast<TH2F*>(inputH->GetStatistics("BIN0"));
   
   if(histStat)fOutputContainer->Add(histStat);
-  else if(DebugLevel()>0)
-    printf("AliAnalysisTaskCaloTrackCorrelation::FinishTaskOutput() - Stat histogram not available check, \n if ESDs, that AliPhysicsSelection was on, \n if AODs, if EventStat_temp.root exists \n");
+  else AliDebug(0,"Stat histogram not available check, \n if ESDs, that AliPhysicsSelection was on, \n if AODs, if EventStat_temp.root exists");
   
   if(histBin0)fOutputContainer->Add(histBin0);
   
index 9b3f0c7..2d6f56e 100755 (executable)
@@ -43,6 +43,7 @@
 #include "AliAODEvent.h"
 #include "AliESDEvent.h"
 #include "AliAnalysisDataSlot.h"
+#include "AliLog.h"
 
 ClassImp(AliAnalysisTaskCaloTrackCorrelationM)
 
@@ -93,12 +94,13 @@ AliAnalysisTaskCaloTrackCorrelationM::~AliAnalysisTaskCaloTrackCorrelationM()
 void AliAnalysisTaskCaloTrackCorrelationM::UserCreateOutputObjects()
 {
   // Create the output container
-  if (DebugLevel() > 1) printf("AliAnalysisTaskCaloTrackCorrelationM::UserCreateOutputObjects() - Begin\n");
+  AliDebug(1,"Begin");
   
   //Get list of aod arrays, add each aod array to analysis frame 
   TClonesArray *array = 0;
   TList * list = fAna->FillAndGetAODBranchList();
-  if (DebugLevel() >= 1) printf("AliAnalysisTaskCaloTrackCorrelationM::UserCreateOutputObjects() - n AOD branches %d\n",list->GetEntries());
+
+  AliDebug(1,Form("n AOD branches %d",list->GetEntries()));
   
   //Put the delta AODs in output file, std or delta
   if((fAna->GetReader())->WriteDeltaAODToFile())
@@ -116,11 +118,11 @@ void AliAnalysisTaskCaloTrackCorrelationM::UserCreateOutputObjects()
   OpenFile(1);
   fOutputContainer = fAna->GetOutputContainer();
   
-  if (DebugLevel() >= 1) printf("AliAnalysisTaskCaloTrackCorrelationM::UserCreateOutputObjects() - n histograms %d\n",fOutputContainer->GetEntries());
+  AliDebug(1,Form("n histograms %d",fOutputContainer->GetEntries()));
   
   fOutputContainer->SetOwner(kTRUE);
   
-  if (DebugLevel() > 1) printf("AliAnalysisTaskCaloTrackCorrelationM::UserCreateOutputObjects() - End\n");
+  AliDebug(1,"End");
   
   PostData(1,fOutputContainer);
        
@@ -141,7 +143,7 @@ void AliAnalysisTaskCaloTrackCorrelationM::Init()
 {
   // Initialization
   
-  if (DebugLevel() > 1) printf("AliAnalysisTaskCaloTrackCorrelationM::Init() - Begin\n");
+  AliDebug(1,"Begin");
   
   fInputEvent = new AliMixedEvent() ; 
   
@@ -149,15 +151,15 @@ void AliAnalysisTaskCaloTrackCorrelationM::Init()
   
   if (fConfigName.Length())
   {
-    printf("AliAnalysisTaskCaloTrackCorrelationM::Init() - ### Configuration file is %s.C ###\n", fConfigName.Data());
+    AliInfo(Form("### Configuration file is %s.C ###", fConfigName.Data()));
     gROOT->LoadMacro(fConfigName+".C");
     fAna = (AliAnaCaloTrackCorrMaker*) gInterpreter->ProcessLine("ConfigAnalysis()");
   }
   
   if(!fAna)
   {
-    printf("AliAnalysisTaskCaloTrackCorrelationM::Init() - Analysis maker pointer not initialized, no analysis specified, STOP !\n");
-    abort();
+    AliFatal("Analysis maker pointer not initialized, no analysis specified, STOP!");
+    return; // coverity
   }
   
   // Add different generator particles to PDG Data Base 
@@ -174,7 +176,7 @@ void AliAnalysisTaskCaloTrackCorrelationM::Init()
   if((fAna->GetReader())->GetDeltaAODFileName()!="")
          AliAnalysisManager::GetAnalysisManager()->RegisterExtraFile((fAna->GetReader())->GetDeltaAODFileName());
   
-  if (DebugLevel() > 1) printf("AliAnalysisTaskCaloTrackCorrelationM::Init() - End\n");
+  AliDebug(1,"End");
   
 }
 
@@ -184,13 +186,13 @@ void AliAnalysisTaskCaloTrackCorrelationM::UserExec(Option_t */*option*/)
 {
   // Execute analysis for current event
   
-  if (DebugLevel() > 1) printf("AliAnalysisTaskCaloTrackCorrelationM::UserExec() - Begin\n");
+  AliDebug(1,"Begin");
   
   //Get the type of data, check if type is correct
   Int_t  datatype = fAna->GetReader()->GetDataType();
   if(datatype != AliCaloTrackReader::kESD && datatype != AliCaloTrackReader::kAOD &&
      datatype != AliCaloTrackReader::kMC){
-    printf("AliAnalysisTaskCaloTrackCorrelationM::UserExec() - Wrong type of data\n");
+    AliError("Wrong type of data");
     return ;
   }
   
@@ -221,13 +223,14 @@ void AliAnalysisTaskCaloTrackCorrelationM::UserExec(Option_t */*option*/)
   fAna->ProcessEvent((Int_t) Entry(), CurrentFileName());
   
   //printf("AliAnalysisTaskCaloTrackCorrelationM::Current Event %d; Current File Name : %s\n",(Int_t) Entry(), CurrentFileName());
-  if (DebugLevel() > 1) printf("AliAnalysisTaskCaloTrackCorrelationM::UserExec() - End\n");
   
   PostData(1, fOutputContainer);
        
   AliAnalysisDataSlot *out0 = GetOutputSlot(0);
   if (out0 && out0->IsConnected()) PostData(0, fTreeA);  
   
+  AliDebug(1,"End");
+
   //gObjectTable->Print();
   
 }
index 44069fe..10a72fa 100644 (file)
@@ -482,8 +482,7 @@ void AliAnalysisTaskCounter::FinishTaskOutput()
   
   if(histStat)
     fOutputContainer->Add(histStat);
-  else if(DebugLevel() > 1)
-    printf("AliAnalysisTaskCounter::FinishTaskOutput() - Stat histogram not available check, \n if ESDs, that AliPhysicsSelection was on, \n if AODs, if EventStat_temp.root exists \n");
+  else AliInfo("Stat histogram not available check, \n if ESDs, that AliPhysicsSelection was on, \n if AODs, if EventStat_temp.root exists");
 
   if(histBin0)
     fOutputContainer->Add(histBin0); 
@@ -521,7 +520,7 @@ Bool_t AliAnalysisTaskCounter::Notify()
   
   if(!fh1Xsec||!fh1Trials)
   {
-    Printf("%s%d No Histogram fh1Xsec",(char*)__FILE__,__LINE__);
+    AliInfo(Form("%s%d No Histogram fh1Xsec",(char*)__FILE__,__LINE__));
     return kFALSE;
   }
   
@@ -538,9 +537,9 @@ Bool_t AliAnalysisTaskCounter::Notify()
   
   fh1Trials->Fill("#sum{ntrials}",trials);
   
-  printf("AliAnalysisTaskCounter::Notify() - xs %f, trial %f, avg trials %f\n",xsection,trials, fAvgTrials);
+  AliInfo(Form("xs %f, trial %f, avg trials %f\n",xsection,trials, fAvgTrials));
   
-  if(fDebug) Printf("Reading File %s",fInputHandler->GetTree()->GetCurrentFile()->GetName());
+  AliDebug(1,Form("Reading File %s",fInputHandler->GetTree()->GetCurrentFile()->GetName()));
   
   return kTRUE;
 }
index 70cd4a9..d6b44b9 100755 (executable)
@@ -53,6 +53,7 @@
 #include "AliAODPWG4Particle.h"
 #include "AliCalorimeterUtils.h"
 #include "AliVEvent.h"
+#include "AliLog.h"
 
 // ---- Detector ----
 #include "AliEMCALPIDUtils.h"
@@ -202,13 +203,16 @@ void AliCaloPID::InitParameters()
   fPHOSPhotonWeightFormulaExpression = "0.98*(x<40)+ 0.68*(x>=100)+(x>=40 && x<100)*(0.98+x*(6e-3)-x*x*(2e-04)+x*x*x*(1.1e-06))";
   fPHOSPi0WeightFormulaExpression    = "0.98*(x<65)+ 0.915*(x>=100)+(x>=65 && x-x*(1.95e-3)-x*x*(4.31e-05)+x*x*x*(3.61e-07))"   ;
   
-  if(fRecalculateBayesian){
-    if(fParticleFlux == kLow){
-      printf("AliCaloPID::Init() - SetLOWFluxParam\n");
+  if(fRecalculateBayesian)
+  {
+    if(fParticleFlux == kLow)
+    {
+      AliInfo("SetLOWFluxParam");
       fEMCALPIDUtils->SetLowFluxParam() ;
     }
-    else if (fParticleFlux == kHigh){
-      printf("AliCaloPID::Init() - SetHIGHFluxParam\n");
+    else if (fParticleFlux == kHigh)
+    {
+      AliInfo("SetHighFluxParam");
       fEMCALPIDUtils->SetHighFluxParam() ;
     }
   }
@@ -586,13 +590,13 @@ Int_t AliCaloPID::GetIdentifiedParticleType(AliVCluster * cluster)
   // Calculate PID SS from data, do not use bayesian weights
   // -------------------------------------------------------
   
-  if(fDebug > 0) printf("AliCaloPID::GetIdentifiedParticleType: EMCAL %d?, E %3.2f, l0 %3.2f, l1 %3.2f, disp %3.2f, tof %1.11f, distCPV %3.2f, distToBC %1.1f, NMax %d\n",
-                        cluster->IsEMCAL(),energy,lambda0,cluster->GetM20(),cluster->GetDispersion(),cluster->GetTOF(), 
-                        cluster->GetEmcCpvDistance(), cluster->GetDistanceToBadChannel(),cluster->GetNExMax());
+  AliDebug(1,Form("EMCAL %d?, E %3.2f, l0 %3.2f, l1 %3.2f, disp %3.2f, tof %1.11f, distCPV %3.2f, distToBC %1.1f, NMax %d",
+                  cluster->IsEMCAL(),energy,lambda0,cluster->GetM20(),cluster->GetDispersion(),cluster->GetTOF(),
+                  cluster->GetEmcCpvDistance(), cluster->GetDistanceToBadChannel(),cluster->GetNExMax()));
   
   if(cluster->IsEMCAL())
   {
-    if(fDebug > 0) printf("AliCaloPID::GetIdentifiedParticleType() - EMCAL SS %f <%f < %f?\n",fEMCALL0CutMin, lambda0, fEMCALL0CutMax);
+    AliDebug(1,Form("EMCAL SS %f <%f < %f?",fEMCALL0CutMin, lambda0, fEMCALL0CutMax));
     
     if(lambda0 < fEMCALL0CutMax && lambda0 > fEMCALL0CutMin) return kPhoton ;
     else                                                     return kNeutralUnknown ; 
@@ -669,7 +673,7 @@ Int_t AliCaloPID::GetIdentifiedParticleTypeFromBayesWeights(Bool_t isEMCAL, Doub
     else                                                     pdg = kNeutralUnknown ;
   }
   
-  if(fDebug > 0)printf("AliCaloPID::GetIdentifiedParticleType:Final Pdg: %d, cluster energy %2.2f \n", pdg,energy);
+  AliDebug(1,Form("Final Pdg: %d, cluster energy %2.2f", pdg,energy));
 
   return pdg ;
   
@@ -702,13 +706,12 @@ Int_t AliCaloPID::GetIdentifiedParticleTypeFromClusterSplitting(AliVCluster* clu
   //If too low number of cells, skip it
   if ( nc < fSplitMinNCells)  return kNeutralUnknown ; 
   
-  if(fDebug > 0) printf("\t pass nCells cut\n");
+  AliDebug(2,"\t pass nCells cut");
   
   // Get Number of local maxima
   nMax  = caloutils->GetNumberOfLocalMaxima(cluster, cells, absIdList, maxEList) ; 
   
-  if(fDebug > 0) printf("AliCaloPID::GetIdentifiedParticleTypeFromClusterSplitting() - Cluster : E %1.1f, M02 %1.2f, NLM %d, N Cells %d\n",
-                        eClus,m02,nMax,nc);
+  AliDebug(1,Form("Cluster : E %1.1f, M02 %1.2f, NLM %d, N Cells %d",eClus,m02,nMax,nc));
 
   //---------------------------------------------------------------------
   // Get the 2 max indeces and do inv mass
@@ -787,9 +790,9 @@ Int_t AliCaloPID::GetIdentifiedParticleTypeFromClusterSplitting(AliVCluster* clu
   
   if(absId2<0 || absId1<0) 
   {
-    if(fDebug > 0) printf("AliCaloPID::GetIdentifiedParticleTypeFromClusterSplitting() - Bad index for local maxima : N max %d, i1 %d, i2 %d, cluster E %2.2f, ncells %d, m02 %2.2f\n",
-                          nMax,absId1,absId2,eClus,nc,m02);
-    return kNeutralUnknown ; 
+    AliDebug(1,Form("Bad index for local maxima : N max %d, i1 %d, i2 %d, cluster E %2.2f, ncells %d, m02 %2.2f",
+                    nMax,absId1,absId2,eClus,nc,m02));
+    return kNeutralUnknown ;
   }
   
   //---------------------------------------------------------------------
@@ -828,7 +831,7 @@ Int_t AliCaloPID::GetIdentifiedParticleTypeFromClusterSplitting(AliVCluster* clu
   else          splitFracCut = fSplitEFracMin[2];
   if((e1+e2)/eClus < splitFracCut) return kNeutralUnknown ;
 
-  if(fDebug > 0) printf("\t pass Split E frac cut\n");
+  AliDebug(2,"\t pass Split E frac cut");
   
   // Consider sub-clusters with minimum energy
   Float_t minECut = fSubClusterEMin[2];
@@ -840,7 +843,7 @@ Int_t AliCaloPID::GetIdentifiedParticleTypeFromClusterSplitting(AliVCluster* clu
     return kNeutralUnknown ;
   }
 
-  if(fDebug > 0) printf("\t pass min sub-cluster E cut\n");
+  AliDebug(2,"\t pass min sub-cluster E cut");
   
   // Asymmetry of cluster
   Float_t asy =-10;
@@ -849,7 +852,7 @@ Int_t AliCaloPID::GetIdentifiedParticleTypeFromClusterSplitting(AliVCluster* clu
   if( !IsInPi0SplitAsymmetryRange(eClus,asy,nMax) ) return kNeutralUnknown ;
   
   
-  if (fDebug>0) printf("\t pass asymmetry cut\n");
+  AliDebug(2,"\t pass asymmetry cut");
   
   Bool_t pi0OK = kFALSE;
   Bool_t etaOK = kFALSE;
@@ -864,10 +867,10 @@ Int_t AliCaloPID::GetIdentifiedParticleTypeFromClusterSplitting(AliVCluster* clu
   if(nMax > 2) energy = e1+e2; // In case of NLM>2 use mass cut for NLM=2 but for the split sum not the cluster energy that is not the pi0 E.
   
   // Check the mass, and set an ID to the splitted cluster
-  if     ( conOK && mass < fMassPhoMax && mass > fMassPhoMin ) { if(fDebug > 0) printf("\t Split Conv \n"); return kPhoton ; }
-  else if( etaOK && mass < fMassEtaMax && mass > fMassEtaMin ) { if(fDebug > 0) printf("\t Split Eta \n");  return kEta    ; }
-  else if( pi0OK && IsInPi0SplitMassRange(energy,mass,nMax)  ) { if(fDebug > 0) printf("\t Split Pi0 \n");  return kPi0    ; }
-  else                                                                                                      return kNeutralUnknown ;
+  if     ( conOK && mass < fMassPhoMax && mass > fMassPhoMin ) { AliDebug(2,"\t Split Conv"); return kPhoton ; }
+  else if( etaOK && mass < fMassEtaMax && mass > fMassEtaMin ) { AliDebug(2,"\t Split Eta" ); return kEta    ; }
+  else if( pi0OK && IsInPi0SplitMassRange(energy,mass,nMax)  ) { AliDebug(2,"\t Split Pi0" ); return kPi0    ; }
+  else                                                                                        return kNeutralUnknown ;
   
 }
 
@@ -879,67 +882,67 @@ TString  AliCaloPID::GetPIDParametersList()
   TString parList ; //this will be list of parameters used for this analysis.
   const Int_t buffersize = 255;
   char onePar[buffersize] ;
-  snprintf(onePar,buffersize,"--- AliCaloPID ---\n") ;
+  snprintf(onePar,buffersize,"--- AliCaloPID ---") ;
   parList+=onePar ;    
   if(fUseBayesianWeights)
   {
-    snprintf(onePar,buffersize,"fEMCALPhotonWeight =%2.2f (EMCAL bayesian weight for photons)\n",fEMCALPhotonWeight) ;
+    snprintf(onePar,buffersize,"fEMCALPhotonWeight =%2.2f (EMCAL bayesian weight for photons)",fEMCALPhotonWeight) ;
     parList+=onePar ;
-    snprintf(onePar,buffersize,"fEMCALPi0Weight =%2.2f (EMCAL bayesian weight for pi0)\n",fEMCALPi0Weight) ;
+    snprintf(onePar,buffersize,"fEMCALPi0Weight =%2.2f (EMCAL bayesian weight for pi0)",fEMCALPi0Weight) ;
     parList+=onePar ;
-    snprintf(onePar,buffersize,"fEMCALElectronWeight =%2.2f(EMCAL bayesian weight for electrons)\n",fEMCALElectronWeight) ;
+    snprintf(onePar,buffersize,"fEMCALElectronWeight =%2.2f(EMCAL bayesian weight for electrons)",fEMCALElectronWeight) ;
     parList+=onePar ;
-    snprintf(onePar,buffersize,"fEMCALChargeWeight =%2.2f (EMCAL bayesian weight for charged hadrons)\n",fEMCALChargeWeight) ;
+    snprintf(onePar,buffersize,"fEMCALChargeWeight =%2.2f (EMCAL bayesian weight for charged hadrons)",fEMCALChargeWeight) ;
     parList+=onePar ;
-    snprintf(onePar,buffersize,"fEMCALNeutralWeight =%2.2f (EMCAL bayesian weight for neutral hadrons)\n",fEMCALNeutralWeight) ;
+    snprintf(onePar,buffersize,"fEMCALNeutralWeight =%2.2f (EMCAL bayesian weight for neutral hadrons)",fEMCALNeutralWeight) ;
     parList+=onePar ;
-    snprintf(onePar,buffersize,"fPHOSPhotonWeight =%2.2f (PHOS bayesian weight for photons)\n",fPHOSPhotonWeight) ;
+    snprintf(onePar,buffersize,"fPHOSPhotonWeight =%2.2f (PHOS bayesian weight for photons)",fPHOSPhotonWeight) ;
     parList+=onePar ;
-    snprintf(onePar,buffersize,"fPHOSPi0Weight =%2.2f (PHOS bayesian weight for pi0)\n",fPHOSPi0Weight) ;
+    snprintf(onePar,buffersize,"fPHOSPi0Weight =%2.2f (PHOS bayesian weight for pi0)",fPHOSPi0Weight) ;
     parList+=onePar ;
-    snprintf(onePar,buffersize,"fPHOSElectronWeight =%2.2f(PHOS bayesian weight for electrons)\n",fPHOSElectronWeight) ;
+    snprintf(onePar,buffersize,"fPHOSElectronWeight =%2.2f(PHOS bayesian weight for electrons)",fPHOSElectronWeight) ;
     parList+=onePar ;
-    snprintf(onePar,buffersize,"fPHOSChargeWeight =%2.2f (PHOS bayesian weight for charged hadrons)\n",fPHOSChargeWeight) ;
+    snprintf(onePar,buffersize,"fPHOSChargeWeight =%2.2f (PHOS bayesian weight for charged hadrons)",fPHOSChargeWeight) ;
     parList+=onePar ;
-    snprintf(onePar,buffersize,"fPHOSNeutralWeight =%2.2f (PHOS bayesian weight for neutral hadrons)\n",fPHOSNeutralWeight) ;
+    snprintf(onePar,buffersize,"fPHOSNeutralWeight =%2.2f (PHOS bayesian weight for neutral hadrons)",fPHOSNeutralWeight) ;
     parList+=onePar ;
     
     if(fPHOSWeightFormula)
     {
-      snprintf(onePar,buffersize,"PHOS Photon Weight Formula: %s\n",fPHOSPhotonWeightFormulaExpression.Data() ) ;
+      snprintf(onePar,buffersize,"PHOS Photon Weight Formula: %s",fPHOSPhotonWeightFormulaExpression.Data() ) ;
       parList+=onePar;
-      snprintf(onePar,buffersize,"PHOS Pi0    Weight Formula: %s\n",fPHOSPi0WeightFormulaExpression.Data()    ) ;
+      snprintf(onePar,buffersize,"PHOS Pi0    Weight Formula: %s",fPHOSPi0WeightFormulaExpression.Data()    ) ;
       parList+=onePar;   
     }
   }
   else
   {
-    snprintf(onePar,buffersize,"EMCAL: fEMCALL0CutMin =%2.2f, fEMCALL0CutMax =%2.2f  (Cut on Shower Shape) \n",fEMCALL0CutMin, fEMCALL0CutMax) ;
+    snprintf(onePar,buffersize,"EMCAL: fEMCALL0CutMin =%2.2f, fEMCALL0CutMax =%2.2f  (Cut on Shower Shape)",fEMCALL0CutMin, fEMCALL0CutMax) ;
     parList+=onePar ;
-    snprintf(onePar,buffersize,"EMCAL: fEMCALDEtaCut =%2.2f, fEMCALDPhiCut =%2.2f  (Cut on track matching) \n",fEMCALDEtaCut, fEMCALDPhiCut) ;
+    snprintf(onePar,buffersize,"EMCAL: fEMCALDEtaCut =%2.2f, fEMCALDPhiCut =%2.2f  (Cut on track matching)",fEMCALDEtaCut, fEMCALDPhiCut) ;
     parList+=onePar ;
-    snprintf(onePar,buffersize,"fTOFCut  =%e (Cut on TOF, used in PID evaluation) \n",fTOFCut) ;
+    snprintf(onePar,buffersize,"fTOFCut  =%e (Cut on TOF, used in PID evaluation)",fTOFCut) ;
     parList+=onePar ;  
-    snprintf(onePar,buffersize,"fPHOSRCut =%2.2f, fPHOSDispersionCut =%2.2f  (Cut on Shower Shape and CPV) \n",fPHOSRCut,fPHOSDispersionCut) ;
+    snprintf(onePar,buffersize,"fPHOSRCut =%2.2f, fPHOSDispersionCut =%2.2f  (Cut on Shower Shape and CPV)",fPHOSRCut,fPHOSDispersionCut) ;
     parList+=onePar ;
     
   }
   
   if(fUseSimpleM02Cut)
   {
-    snprintf(onePar,buffersize,"%2.2f< M02 < %2.2f \n",    fSplitM02MinCut, fSplitM02MaxCut) ;
+    snprintf(onePar,buffersize,"%2.2f< M02 < %2.2f",    fSplitM02MinCut, fSplitM02MaxCut) ;
     parList+=onePar ;
   }
-  snprintf(onePar,buffersize,"fMinNCells =%d \n",        fSplitMinNCells) ;
+  snprintf(onePar,buffersize,"fMinNCells =%d",        fSplitMinNCells) ;
   parList+=onePar ;
   if(fUseSimpleMassCut)
   {
-    snprintf(onePar,buffersize,"pi0 : %2.1f < m <%2.1f\n", fMassPi0Min,fMassPi0Max);
+    snprintf(onePar,buffersize,"pi0 : %2.1f < m <%2.1f", fMassPi0Min,fMassPi0Max);
     parList+=onePar ;
   }
-  snprintf(onePar,buffersize,"eta : %2.1f < m <%2.1f\n", fMassEtaMin,fMassEtaMax);
+  snprintf(onePar,buffersize,"eta : %2.1f < m <%2.1f", fMassEtaMin,fMassEtaMax);
   parList+=onePar ;
-  snprintf(onePar,buffersize,"conv: %2.1f < m <%2.1f\n", fMassPhoMin,fMassPhoMax);
+  snprintf(onePar,buffersize,"conv: %2.1f < m <%2.1f", fMassPhoMin,fMassPhoMax);
   parList+=onePar ;
   
   
@@ -1043,45 +1046,43 @@ void AliCaloPID::SetPIDBits(AliVCluster * cluster,
   
   //Set PID pdg
   ph->SetIdentifiedParticleType(GetIdentifiedParticleType(cluster));
-  
-  if(fDebug > 0)
-  { 
-    printf("AliCaloPID::SetPIDBits: TOF %e, Lambda0 %2.2f, Lambda1 %2.2f\n",tof , l0, l1);     
-    printf("AliCaloPID::SetPIDBits: pdg %d, bits: TOF %d, Dispersion %d, Charge %d\n",
-           ph->GetIdentifiedParticleType(), ph->GetTOFBit() , ph->GetDispBit() , ph->GetChargedBit()); 
-  }
+  AliDebug(1,Form("TOF %e, Lambda0 %2.2f, Lambda1 %2.2f",tof , l0, l1));
+  AliDebug(1,Form("pdg %d, bits: TOF %d, Dispersion %d, Charge %d",
+           ph->GetIdentifiedParticleType(), ph->GetTOFBit() , ph->GetDispBit() , ph->GetChargedBit()));
 }
 
 //_________________________________________________________
 Bool_t AliCaloPID::IsTrackMatched(AliVCluster* cluster,
-                                  AliCalorimeterUtils * cu, 
-                                  AliVEvent* event) const 
+                                  AliCalorimeterUtils * cu,
+                                  AliVEvent* event) const
 {
   //Check if there is any track attached to this cluster
   
   Int_t nMatches = cluster->GetNTracksMatched();
   AliVTrack * track = 0;
-
+  
   if(nMatches > 0)
   {
     //In case of ESDs, by default without match one entry with negative index, no match, reject.
     if(!strcmp("AliESDCaloCluster",Form("%s",cluster->ClassName())))
-    {    
+    {
       Int_t iESDtrack = cluster->GetTrackMatchedIndex();
       if(iESDtrack >= 0) track = dynamic_cast<AliVTrack*> (event->GetTrack(iESDtrack));
       else return kFALSE;
       
       if (!track)
       {
-        if(fDebug > 0) printf("AliCaloPID::IsTrackMatched() - Null matched track in ESD when index is OK!\n");
+        AliWarning("Null matched track in ESD when index is OK!");
         return kFALSE;
       }
-    }      
-    else { // AOD
+    }
+    else
+    { // AOD
       track = dynamic_cast<AliVTrack*> (cluster->GetTrackMatched(0));
       if (!track)
       {
-        if(fDebug > 0) printf("AliCaloPID::IsTrackMatched() - Null matched track in AOD!\n");
+        AliWarning("Null matched track in AOD!");
         return kFALSE;
       }
     }
@@ -1095,8 +1096,8 @@ Bool_t AliCaloPID::IsTrackMatched(AliVCluster* cluster,
       dR = 2000., dZ = 2000.;
       cu->GetEMCALRecoUtils()->GetMatchedResiduals(cluster->GetID(),dZ,dR);
     }
-        
-    if(cluster->IsPHOS()) 
+    
+    if(cluster->IsPHOS())
     {
       Int_t charge = track->Charge();
       Double_t mf  = event->GetMagneticField();
@@ -1106,25 +1107,25 @@ Bool_t AliCaloPID::IsTrackMatched(AliVCluster* cluster,
     }//PHOS
     else //EMCAL
     {
-    if(fDebug > 1) 
-        printf("AliCaloPID::IsTrackMatched - EMCAL dR %f < %f, dZ %f < %f \n",dR, fEMCALDPhiCut, dZ, fEMCALDEtaCut);
       
-      if(TMath::Abs(dR) < fEMCALDPhiCut && 
+      AliDebug(1,Form("EMCAL dR %f < %f, dZ %f < %f ",dR, fEMCALDPhiCut, dZ, fEMCALDEtaCut));
+      
+      if(TMath::Abs(dR) < fEMCALDPhiCut &&
          TMath::Abs(dZ) < fEMCALDEtaCut)   return kTRUE;
       else                                 return kFALSE;
       
-    }//EMCAL cluster 
+    }//EMCAL cluster
     
     
   } // more than 1 match, at least one track in array
   else return kFALSE;
-    
+  
 }
 
 //___________________________________________________________________________________________________
-Float_t AliCaloPID::TestPHOSDispersion(const Double_t pt, const Double_t l1, const Double_t l2) const 
+Float_t AliCaloPID::TestPHOSDispersion(const Double_t pt, const Double_t l1, const Double_t l2) const
 {
-  //Check if cluster photon-like. Uses photon cluster parameterization in real pp data 
+  //Check if cluster photon-like. Uses photon cluster parameterization in real pp data
   //Returns distance in sigmas. Recommended cut 2.5
   
   Double_t l2Mean  = 1.53126+9.50835e+06/(1.+1.08728e+07*pt+1.73420e+06*pt*pt) ;
@@ -1132,13 +1133,13 @@ Float_t AliCaloPID::TestPHOSDispersion(const Double_t pt, const Double_t l1, con
   Double_t l2Sigma = 6.48260e-02+7.60261e+10/(1.+1.53012e+11*pt+5.01265e+05*pt*pt)+9.00000e-03*pt;
   Double_t l1Sigma = 4.44719e-04+6.99839e-01/(1.+1.22497e+00*pt+6.78604e-07*pt*pt)+9.00000e-03*pt;
   Double_t c       =-0.35-0.550*TMath::Exp(-0.390730*pt) ;
-  Double_t r2      = 0.5*  (l1-l1Mean)*(l1-l1Mean)/l1Sigma/l1Sigma + 
-                     0.5*  (l2-l2Mean)*(l2-l2Mean)/l2Sigma/l2Sigma +
-                     0.5*c*(l1-l1Mean)*(l2-l2Mean)/l1Sigma/l2Sigma ;
+  Double_t r2      = 0.5*  (l1-l1Mean)*(l1-l1Mean)/l1Sigma/l1Sigma +
+  0.5*  (l2-l2Mean)*(l2-l2Mean)/l2Sigma/l2Sigma +
+  0.5*c*(l1-l1Mean)*(l2-l2Mean)/l1Sigma/l2Sigma ;
   
-  if(fDebug > 0) printf("AliCaloPID::TestPHOSDispersion() - PHOS SS R %f < %f?\n", TMath::Sqrt(r2), fPHOSDispersionCut);
+  AliDebug(1,Form("PHOS SS R %f < %f?", TMath::Sqrt(r2), fPHOSDispersionCut));
   
-  return TMath::Sqrt(r2) ; 
+  return TMath::Sqrt(r2) ;
   
 }
 
@@ -1186,9 +1187,9 @@ Float_t AliCaloPID::TestPHOSChargedVeto(const Double_t dx,  const Double_t dz, c
   Double_t rz = (dz-meanZ)/sz ;
   Double_t rx = (dx-meanX)/sx ;
   
-  if(fDebug > 0) 
-    printf("AliCaloPID::TestPHOSDispersion() - PHOS Matching R %f < %f\n",TMath::Sqrt(rx*rx+rz*rz), fPHOSRCut);
+  AliDebug(1,Form("PHOS Matching R %f < %f",TMath::Sqrt(rx*rx+rz*rz), fPHOSRCut));
   
   return TMath::Sqrt(rx*rx+rz*rz) ;
   
 }
+
index eea178e..e97583f 100755 (executable)
@@ -31,6 +31,7 @@
 #include "AliAnalysisManager.h"
 #include "AliMixedEvent.h"
 #include "AliAODEvent.h"
+#include "AliLog.h"
 
 ClassImp(AliCaloTrackAODReader)
 
@@ -127,10 +128,9 @@ Bool_t AliCaloTrackAODReader::SelectTrack(AliVTrack* track, Double_t pTrack[3])
   
   if(!aodtrack) return kFALSE;
   
-  
-  if(fDebug > 2 ) printf("AliCaloTrackAODReader::FillInputCTS():AOD track type: %d (primary %d), hybrid? %d \n",
-                         aodtrack->GetType(),AliAODTrack::kPrimary,
-                         aodtrack->IsHybridGlobalConstrainedGlobal());
+  AliDebug(2,Form("AOD track type: %d (primary %d), hybrid? %d",
+                  aodtrack->GetType(),AliAODTrack::kPrimary,
+                  aodtrack->IsHybridGlobalConstrainedGlobal()));
   
   // Hybrid?
   if (fSelectHybridTracks && fTrackFilterMaskComplementary == 0)
@@ -162,7 +162,7 @@ Bool_t AliCaloTrackAODReader::SelectTrack(AliVTrack* track, Double_t pTrack[3])
     Double_t frac = Double_t(aodtrack->GetTPCnclsS()) / Double_t(aodtrack->GetTPCncls());
     if (frac > fCutTPCSharedClustersFraction)
     {
-      if (fDebug > 2 )printf("\t Reject track, shared cluster fraction %f > %f\n",frac, fCutTPCSharedClustersFraction);
+      AliDebug(2,Form("\t Reject track, shared cluster fraction %f > %f",frac, fCutTPCSharedClustersFraction));
       return kFALSE ;
     }
   }
@@ -172,12 +172,12 @@ Bool_t AliCaloTrackAODReader::SelectTrack(AliVTrack* track, Double_t pTrack[3])
   {
     if ( aodtrack->GetType()!= AliAODTrack::kPrimary )
     {
-      if (fDebug > 2 ) printf("\t Remove not primary track\n");
+      AliDebug(2,"\t Remove not primary track");
       return kFALSE ;
     }
   }
   
-  if (fDebug > 2 ) printf("\t accepted track! \n");
+  AliDebug(2,"\t accepted track!");
   
   track->GetPxPyPz(pTrack) ;
   
@@ -212,8 +212,8 @@ void AliCaloTrackAODReader::SetInputOutputMCEvent(AliVEvent* input,
     }
     else
     {
-      printf("AliCaloTrackAODReader::SetInputOutputMCEvent() - MultiEventHandler is NULL");
-      abort();
+      AliFatal("MultiEventHandler is NULL");
+      return;
     }
   }
   if        (strcmp(input->GetName(),"AliESDEvent") == 0) 
@@ -250,7 +250,7 @@ void AliCaloTrackAODReader::SetInputOutputMCEvent(AliVEvent* input,
   }
   else
   { 
-    AliFatal(Form("AliCaloTrackAODReader::SetInputOutputMCEvent() - STOP : Wrong data format: %s\n",input->GetName()));
+    AliFatal(Form("STOP : Wrong data format: %s",input->GetName()));
   }
   
   SetMC(mc);
index e9a02f9..bc251d5 100755 (executable)
@@ -33,6 +33,7 @@
 #include "AliMixedEvent.h"
 #include "AliESDEvent.h"
 #include "AliESDtrackCuts.h"
+#include "AliLog.h"
 
 ClassImp(AliCaloTrackESDReader)
 
@@ -192,8 +193,8 @@ void AliCaloTrackESDReader::SetInputOutputMCEvent(AliVEvent* esd,
     }
     else
     {
-      printf("AliCaloTrackESDReader::SetInputOutputMCEvent() - MultiEventHandler is NULL");
-      abort();
+      AliFatal("MultiEventHandler is NULL");
+      return;
     }
   }
   if (strcmp(esd->GetName(),"AliESDEvent") == 0)
@@ -203,7 +204,7 @@ void AliCaloTrackESDReader::SetInputOutputMCEvent(AliVEvent* esd,
   
   if(!tesd)
   {
-    AliFatal(Form("AliCaloTrackESDReader::SetInputOutputMCEvent() - STOP ::Wrong reader, here only ESDs. Input name: %s != AliESDEvent \n",esd->GetName()));
+    AliFatal(Form("STOP ::Wrong reader, here only ESDs. Input name: %s != AliESDEvent",esd->GetName()));
   }
   
   SetInputEvent(esd);
index 4178e0e..f662353 100755 (executable)
@@ -42,6 +42,7 @@
 #include "AliAODEvent.h"
 #include "AliFiducialCut.h"
 #include "AliMCAnalysisUtils.h"
+#include "AliLog.h"
 
 ClassImp(AliCaloTrackMCReader)
 
@@ -206,9 +207,9 @@ void  AliCaloTrackMCReader::FillCalorimeters(Int_t & iParticle, TParticle* parti
   AliAODCaloCluster *calo = new AliAODCaloCluster(index,1,labels,fMomentum.E(), x, NULL, ttype, 0);
   
   SetCaloClusterPID(pdg,calo) ;
-  if(fDebug > 3 && fMomentum.Pt() > 0.2)
-    printf("AliCaloTrackMCReader::FillCalorimeters() - PHOS %d?, EMCAL %d? : Selected cluster pdg %d, E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
-           fFillPHOS, fFillEMCAL, pdg,fMomentum.E(), fMomentum.Pt(), fMomentum.Phi()*TMath::RadToDeg(),fMomentum.Eta());
+  AliDebug(3,Form("PHOS %d?, EMCAL %d? : Selected cluster pdg %d, E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f",
+                  fFillPHOS, fFillEMCAL, pdg,fMomentum.E(), fMomentum.Pt(), fMomentum.Phi()*TMath::RadToDeg(),fMomentum.Eta()));
   
   //reference the selected object to the list
   if(fFillPHOS)fPHOSClusters ->Add(calo);
@@ -273,9 +274,8 @@ Bool_t AliCaloTrackMCReader::FillInputEvent(Int_t iEntry,
           
           if(TMath::Abs(pdg) == 11 && GetStack()->Particle(particle->GetFirstMother())->GetPdgCode()==22) continue ;
           
-          if(fDebug > 3 && pt > 0.2)
-            printf("AliCaloTrackMCReader::FillInputEvent() - CTS : Selected tracks E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
-                   en,pt,phi*TMath::RadToDeg(),eta);
+          AliDebug(2,Form("CTS : Selected tracks E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f",
+                          en,pt,phi*TMath::RadToDeg(),eta));
           
           x[0] = particle->Vx(); x[1] = particle->Vy(); x[2] = particle->Vz();
           p[0] = particle->Px(); p[1] = particle->Py(); p[2] = particle->Pz();
@@ -304,6 +304,8 @@ Bool_t AliCaloTrackMCReader::FillInputEvent(Int_t iEntry,
         //Fill particle/calocluster arrays
         if(!fDecayPi0) 
         {
+          AliDebug(2,Form("Calo : Selected tracks E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f",
+                          en,pt,phi*TMath::RadToDeg(),eta));
           FillCalorimeters(iParticle, particle);
         }
         else 
@@ -329,7 +331,12 @@ Bool_t AliCaloTrackMCReader::FillInputEvent(Int_t iEntry,
               FillCalorimeters(iParticle,pPhoton2);
             }//pt cut
           }//pi0
-          else FillCalorimeters(iParticle,particle); //Add the rest
+          else
+          {
+            AliDebug(2,Form("Calo : Selected tracks E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f",
+                            en,pt,phi*TMath::RadToDeg(),eta));
+            FillCalorimeters(iParticle,particle); //Add the rest
+          }
         }
       }//neutral particles
     } //particle with correct status
index 4a9c0dc..d98fe3b 100755 (executable)
@@ -47,6 +47,7 @@
 #include "AliInputEventHandler.h"
 #include "AliAODMCParticle.h"
 #include "AliStack.h"
+#include "AliLog.h"
 
 // ---- Detectors ----
 #include "AliPHOSGeoUtils.h"
@@ -309,12 +310,12 @@ Bool_t AliCaloTrackReader::CheckEventTriggers()
     if(eventType!=7)
     return kFALSE; //Only physics event, do not use for simulated events!!!
     
-    if(fDebug > 0)
-      printf("AliCaloTrackReader::CheckEventTriggers() - FiredTriggerClass <%s>, selected class <%s>, compare name %d\n",
-             GetFiredTriggerClasses().Data(),fFiredTriggerClassName.Data(), GetFiredTriggerClasses().Contains(fFiredTriggerClassName));
+   AliDebug(1,Form("FiredTriggerClass <%s>, selected class <%s>, compare name %d",
+                   GetFiredTriggerClasses().Data(),fFiredTriggerClassName.Data(),
+                   GetFiredTriggerClasses().Contains(fFiredTriggerClassName)));
     
     if( !GetFiredTriggerClasses().Contains(fFiredTriggerClassName) ) return kFALSE;
-    else if(fDebug > 0) printf("AliCaloTrackReader::CheckEventTriggers() - Accepted triggered event\n");
+    else AliDebug(1,"Accepted triggered event");
   }
   else if(fAnaLED)
   {
@@ -332,11 +333,11 @@ Bool_t AliCaloTrackReader::CheckEventTriggers()
     //   kSystemSoftwareTriggerEvent   = 12, // SYSTEM_SOFTWARE_TRIGGER_EVENT
     //   kDetectorSoftwareTriggerEvent = 13  // DETECTOR_SOFTWARE_TRIGGER_EVENT
     
-         if(eventType!=7 && fDebug > 1 )printf("AliCaloTrackReader::CheckEventTriggers() - DO LED, Event Type <%d>, 8 Calibration \n",  eventType);
-         if(eventType!=8)return kFALSE;
+         //if(eventType!=7 && fDebug > 1 )printf("AliCaloTrackReader::CheckEventTriggers() - DO LED, Event Type <%d>, 8 Calibration \n",  eventType);
+         if(eventType!=8) return kFALSE;
   }
   
-  if(fDebug > 0) printf("AliCaloTrackReader::CheckEventTriggers() - Pass Trigger name rejection \n");
+  AliDebug(1,"Pass Trigger name rejection");
 
   //-----------------------------------------------------------------
   // In case of mixing analysis, select here the trigger of the event
@@ -360,7 +361,7 @@ Bool_t AliCaloTrackReader::CheckEventTriggers()
     if(!isTrigger && !isMB) return kFALSE;
     
     //printf("Selected triggered event : %s\n",GetFiredTriggerClasses().Data());
-    if(fDebug > 0) printf("AliCaloTrackReader::CheckEventTriggers() - Pass uninteresting triggered events rejection in case of mixing analysis \n");
+    AliDebug(0,"Pass uninteresting triggered events rejection in case of mixing analysis");
   }
 
   //-------------------------------------------------------------------------------------
@@ -374,7 +375,7 @@ Bool_t AliCaloTrackReader::CheckEventTriggers()
   
   if(!okA || !okR) return kFALSE;
   
-  if(fDebug > 0) printf("AliCaloTrackReader::CheckEventTriggers() - Pass event bit rejection \n");
+  AliDebug(1,"Pass event bit rejection");
   
   //----------------------------------------------------------------------
   // Do not count events that were likely triggered by an exotic cluster
@@ -407,17 +408,16 @@ Bool_t AliCaloTrackReader::CheckEventTriggers()
     // If requested, remove badly triggeed events, but only when the EMCal trigger bit is set
     if(fRemoveBadTriggerEvents)
     {
-     if(fDebug > 0)
-      printf("AliCaloTrackReader::CheckEventTriggers() - ACCEPT triggered event? \n exotic? %d - bad cell %d - bad Max cell %d - BC %d  - Matched %d\n",
-             fIsExoticEvent,fIsBadCellEvent, fIsBadMaxCellEvent, fTriggerClusterBC,fIsTriggerMatch);
+     AliDebug(1,Form("ACCEPT triggered event? \n exotic? %d - bad cell %d - bad Max cell %d - BC %d  - Matched %d\n",
+                     fIsExoticEvent,fIsBadCellEvent, fIsBadMaxCellEvent, fTriggerClusterBC,fIsTriggerMatch));
       if     (fIsExoticEvent)         return kFALSE;
       else if(fIsBadCellEvent)        return kFALSE;
       else if(fRemoveUnMatchedTriggers && !fIsTriggerMatch) return kFALSE ;
       else if(fTriggerClusterBC != 0) return kFALSE;
-      if(fDebug > 0) printf("\t *** YES for %s\n",GetFiredTriggerClasses().Data());
+      AliDebug(1,Form("\t *** YES for %s",GetFiredTriggerClasses().Data()));
     }
     
-    if(fDebug > 0) printf("AliCaloTrackReader::CheckEventTriggers() - Pass EMCal triggered event rejection \n");
+   AliDebug(1,"Pass EMCal triggered event rejection \n");
   }
   
   //-------------------------------------------------------------------------------------
@@ -425,7 +425,7 @@ Bool_t AliCaloTrackReader::CheckEventTriggers()
   //-------------------------------------------------------------------------------------
   if (GetFiredTriggerClasses().Contains("FAST")  && !GetFiredTriggerClasses().Contains("ALL") && !fAcceptFastCluster)
   {
-    if(fDebug > 0)  printf("AliCaloTrackReader::CheckEventTriggers() - Do not count events from fast cluster, trigger name %s\n",fFiredTriggerClassName.Data());
+    AliDebug(1,Form("Do not count events from fast cluster, trigger name %s\n",fFiredTriggerClassName.Data()));
     return kFALSE;
   }
   
@@ -441,7 +441,7 @@ Bool_t AliCaloTrackReader::CheckEventTriggers()
 
     if(reject) return kFALSE;
     
-    if(fDebug > 0) printf("AliCaloTrackReader::CheckEventTriggers() - Pass LED event rejection \n");
+    AliDebug(1,"Pass LED event rejection");
 
   }// Remove LED events
   
@@ -463,8 +463,7 @@ Bool_t AliCaloTrackReader::ComparePtHardAndJetPt()
     Int_t nTriggerJets =  pygeh->NTriggerJets();
     Float_t ptHard = pygeh->GetPtHard();
     
-    if(fDebug > 1)
-      printf("AliCaloTrackReader::ComparePtHardAndJetPt() - Njets: %d, pT Hard %f\n",nTriggerJets, ptHard);
+    AliDebug(1,Form("Njets: %d, pT Hard %f",nTriggerJets, ptHard));
     
     Float_t tmpjet[]={0,0,0,0};
     for(Int_t ijet = 0; ijet< nTriggerJets; ijet++)
@@ -472,14 +471,13 @@ Bool_t AliCaloTrackReader::ComparePtHardAndJetPt()
       pygeh->TriggerJet(ijet, tmpjet);
       jet = new TParticle(94, 21, -1, -1, -1, -1, tmpjet[0],tmpjet[1],tmpjet[2],tmpjet[3], 0,0,0,0);
       
-      if(fDebug > 1)
-        printf("AliCaloTrackReader::ComparePtHardAndJetPt() - jet %d; pycell jet pT %f\n",ijet, jet->Pt());
+      AliDebug(1,Form("jet %d; pycell jet pT %f",ijet, jet->Pt()));
       
       //Compare jet pT and pt Hard
       if(jet->Pt() > fPtHardAndJetPtFactor * ptHard)
       {
-        printf("AliCaloTrackReader::ComparePtHardAndJetPt() - Reject jet event with : pT Hard %2.2f, pycell jet pT %2.2f, rejection factor %1.1f\n",
-               ptHard, jet->Pt(), fPtHardAndJetPtFactor);
+        AliInfo(Form("Reject jet event with : pT Hard %2.2f, pycell jet pT %2.2f, rejection factor %1.1f\n",
+                     ptHard, jet->Pt(), fPtHardAndJetPtFactor));
         return kFALSE;
       }
     }
@@ -510,7 +508,7 @@ Bool_t AliCaloTrackReader::ComparePtHardAndClusterPt()
       
       if(ecluster > fPtHardAndClusterPtFactor*ptHard)
       {
-        printf("AliCaloTrackReader::ComparePtHardAndClusterPt() - Reject : ecluster %2.2f, calo %d, factor %2.2f, ptHard %f\n",ecluster,clus->GetType(),fPtHardAndClusterPtFactor,ptHard);
+        AliInfo(Form("Reject : ecluster %2.2f, calo %d, factor %2.2f, ptHard %f",ecluster,clus->GetType(),fPtHardAndClusterPtFactor,ptHard));
         
         return kFALSE;
       }
@@ -530,7 +528,7 @@ AliStack* AliCaloTrackReader::GetStack() const
     return fMC->Stack();
   else
   {
-    if(fDebug > 1) printf("AliCaloTrackReader::GetStack() - Stack is not available\n");
+    AliDebug(1,"Stack is not available");
     return 0x0 ;
   }
 }
@@ -545,7 +543,7 @@ AliHeader* AliCaloTrackReader::GetHeader() const
   }
   else
   {
-    printf("AliCaloTrackReader::Header is not available\n");
+    AliInfo("Header is not available");
     return 0x0 ;
   }
 }
@@ -687,7 +685,7 @@ TClonesArray* AliCaloTrackReader::GetAODMCParticles() const
 {
   //Return list of particles in AOD, implemented in AliCaloTrackAODReader.
   
-  printf("AliCaloTrackReader::GetAODMCParticles() - Input are not AODs\n");
+  AliInfo("Input are not AODs");
   
   return NULL ;
 }
@@ -697,7 +695,7 @@ AliAODMCHeader* AliCaloTrackReader::GetAODMCHeader() const
 {
   //Return MC header in AOD, implemented in AliCaloTrackAODReader.
   
-  printf("AliCaloTrackReader::GetAODMCHeader() - Input are not AODs\n");
+  AliInfo("Input are not AODs");
   
   return NULL ;
 }
@@ -754,13 +752,14 @@ Int_t AliCaloTrackReader::GetVertexBC(const AliVVertex * vtx)
 void AliCaloTrackReader::Init()
 {
   //Init reader. Method to be called in AliAnaCaloTrackCorrMaker
-    
+  
   if(fReadStack && fReadAODMCParticles)
   {
-    printf("AliCaloTrackReader::Init() - Cannot access stack and mcparticles at the same time, change them \n");
+    AliInfo("Cannot access stack and mcparticles at the same time, change them");
     fReadStack          = kFALSE;
     fReadAODMCParticles = kFALSE;
   }
+  
 }
 
 //_______________________________________
@@ -1067,7 +1066,7 @@ Bool_t AliCaloTrackReader::FillInputEvent(Int_t iEntry, const char * /*curFileNa
   //fCurrentFileName = TString(currentFileName);
   if(!fInputEvent)
   {
-    if(fDebug >= 0) printf("AliCaloTrackReader::FillInputEvent() - Input event not available, skip event analysis\n");
+    AliInfo("Input event not available, skip event analysis");
     return kFALSE;
   }
   
@@ -1080,7 +1079,7 @@ Bool_t AliCaloTrackReader::FillInputEvent(Int_t iEntry, const char * /*curFileNa
   Bool_t accept = CheckEventTriggers();
   if(!accept) return kFALSE;
   
-  if(fDebug > 0) printf("AliCaloTrackReader::FillInputEvent() - Pass Event trigger selection \n");
+  AliDebug(1,"Pass Event trigger selection");
   
   //---------------------------------------------------------------------------
   // In case of analysis of events with jets, skip those with jet pt > 5 pt hard
@@ -1090,13 +1089,13 @@ Bool_t AliCaloTrackReader::FillInputEvent(Int_t iEntry, const char * /*curFileNa
   if(fComparePtHardAndJetPt)
   {
     if(!ComparePtHardAndJetPt()) return kFALSE ;
-    if(fDebug > 0) printf("AliCaloTrackReader::FillInputEvent() - Pass Pt Hard - Jet rejection \n");
+    AliDebug(1,"Pass Pt Hard - Jet rejection");
   }
   
   if(fComparePtHardAndClusterPt)
   {
     if(!ComparePtHardAndClusterPt()) return kFALSE ;
-    if(fDebug > 0) printf("AliCaloTrackReader::FillInputEvent() - Pass Pt Hard - Cluster rejection \n");
+    AliDebug(1,"Pass Pt Hard - Cluster rejection");
   }
   
   //------------------------------------------------------
@@ -1115,7 +1114,8 @@ Bool_t AliCaloTrackReader::FillInputEvent(Int_t iEntry, const char * /*curFileNa
       
       if(timeStampFrac < fTimeStampEventFracMin || timeStampFrac > fTimeStampEventFracMax) return kFALSE;
     }
-    if(fDebug > 0) printf("AliCaloTrackReader::FillInputEvent() - Pass Time Stamp rejection \n");
+    
+    AliDebug(1,"Pass Time Stamp rejection");
   }
   
   //------------------------------------------------------
@@ -1135,7 +1135,7 @@ Bool_t AliCaloTrackReader::FillInputEvent(Int_t iEntry, const char * /*curFileNa
         TMath::Abs(fVertex[0][2] ) < 1.e-6    ) return kFALSE;
   }
   
-  if(fDebug > 0) printf("AliCaloTrackReader::FillInputEvent() - Pass primary vertex rejection \n");
+  AliDebug(1,"Pass primary vertex rejection");
   
   //printf("Reader : IsPileUp %d, Multi %d\n",IsPileUpFromSPD(),fInputEvent->IsPileupFromSPDInMultBins());
   
@@ -1147,7 +1147,7 @@ Bool_t AliCaloTrackReader::FillInputEvent(Int_t iEntry, const char * /*curFileNa
     //printf("pile-up %d, %d, %2.2f, %2.2f, %2.2f, %2.2f\n",bPileup, (Int_t) fPileUpParamSPD[0], fPileUpParamSPD[1], fPileUpParamSPD[2], fPileUpParamSPD[3], fPileUpParamSPD[4]);
     if(bPileup) return kFALSE;
     
-    if(fDebug > 0) printf("AliCaloTrackReader::FillInputEvent() - Pass Pile-Up event rejection \n");
+    AliDebug(1,"Pass Pile-Up event rejection");
   }
   
   if(fDoV0ANDEventSelection)
@@ -1160,10 +1160,11 @@ Bool_t AliCaloTrackReader::FillInputEvent(Int_t iEntry, const char * /*curFileNa
 
     if(!bV0AND)
     {
-      printf("AliCaloTrackReader::FillInputEvent() - Reject event by V0AND\n");
+      AliDebug(1,"Reject event by V0AND");
       return kFALSE;
     }
-    if(fDebug > 0) printf("AliCaloTrackReader::FillInputEvent() - Pass V0AND event rejection \n");
+    
+    AliDebug(1,"Pass V0AND event rejection");
   }
 
   //------------------------------------------------------
@@ -1177,7 +1178,7 @@ Bool_t AliCaloTrackReader::FillInputEvent(Int_t iEntry, const char * /*curFileNa
     Int_t cen = GetEventCentrality();
     if(cen > fCentralityBin[1] || cen < fCentralityBin[0]) return kFALSE; //reject events out of bin.
     
-    if(fDebug > 0) printf("AliCaloTrackReader::FillInputEvent() - Pass centrality rejection \n");
+    AliDebug(1,"Pass centrality rejection");
   }
   
   //---------------------------------------------------------------------------
@@ -1204,14 +1205,14 @@ Bool_t AliCaloTrackReader::FillInputEvent(Int_t iEntry, const char * /*curFileNa
     //Accept events with at least one track
     if(fTrackMult == 0 && fDoRejectNoTrackEvents) return kFALSE ;
     
-    if(fDebug > 0) printf("AliCaloTrackReader::FillInputEvent() - Pass rejection of null track events \n");
+    AliDebug(1,"Pass rejection of null track events");
   }
   
   if(fDoVertexBCEventSelection)
   {
     if(fVertexBC != 0 && fVertexBC != AliVTrack::kTOFBCNA) return kFALSE ;
     
-    if(fDebug > 0) printf("AliCaloTrackReader::FillInputEvent() - Pass rejection of events with vertex at BC!=0 \n");
+    AliDebug(1,"Pass rejection of events with vertex at BC!=0");
   }
   
   if(fFillEMCALCells)
@@ -1234,7 +1235,7 @@ Bool_t AliCaloTrackReader::FillInputEvent(Int_t iEntry, const char * /*curFileNa
   if(fFillInputBackgroundJetBranch)
     FillInputBackgroundJets();
 
-  if(fDebug > 0) printf("AliCaloTrackReader::FillInputEvent() - Event accepted for analysis \n");
+  AliDebug(1,"Event accepted for analysis");
 
   return kTRUE ;
 }
@@ -1251,7 +1252,7 @@ Int_t AliCaloTrackReader::GetEventCentrality() const
     else if(fCentralityOpt==20)  return GetCentrality()->GetCentralityClass5(fCentralityClass); // 20 bins max
     else
     {
-      printf("AliCaloTrackReader::GetEventCentrality() - Unknown centrality option %d, use 10, 20 or 100\n",fCentralityOpt);
+      AliInfo(Form("Unknown centrality option %d, use 10, 20 or 100\n",fCentralityOpt));
       return -1;
     }
   }
@@ -1270,14 +1271,14 @@ Double_t AliCaloTrackReader::GetEventPlaneAngle() const
     
     if(GetEventPlaneMethod()=="Q" && (ep < 0 || ep > TMath::Pi()))
     {
-      if(fDebug > 0 ) printf("AliCaloTrackReader::GetEventPlaneAngle() -  Bad EP for <Q> method : %f\n",ep);
+      AliDebug(1,Form("Bad EP for <Q> method : %f\n",ep));
       return -1000;
     }
     else if(GetEventPlaneMethod().Contains("V0")  )
     {
       if((ep > TMath::Pi()/2 || ep < -TMath::Pi()/2))
       {
-        if(fDebug > 0 ) printf("AliCaloTrackReader::GetEventPlaneAngle() -  Bad EP for <%s> method : %f\n",GetEventPlaneMethod().Data(), ep);
+        AliDebug(1,Form("Bad EP for <%s> method : %f\n",GetEventPlaneMethod().Data(), ep));
         return -1000;
       }
       
@@ -1285,18 +1286,18 @@ Double_t AliCaloTrackReader::GetEventPlaneAngle() const
       
     }
     
-    //printf("AliCaloTrackReader::GetEventPlaneAngle() = %f\n",ep);
-    if(fDebug > 0 )
-    {
-      if     (ep > TMath::Pi()) printf("AliCaloTrackReader::GetEventPlaneAngle() - Too large angle = %f\n",ep);
-      else if(ep < 0          ) printf("AliCaloTrackReader::GetEventPlaneAngle() - Negative angle = %f\n" ,ep);
-    }
+    AliDebug(3,Form("Event plane angle %f",ep));
+    
+//    if(fDebug > 0 )
+//    {
+//      if     (ep > TMath::Pi()) printf("AliCaloTrackReader::GetEventPlaneAngle() - Too large angle = %f\n",ep);
+//      else if(ep < 0          ) printf("AliCaloTrackReader::GetEventPlaneAngle() - Negative angle = %f\n" ,ep);
+//    }
     
     return ep;
   }
   else
   {
-    if(fDataType!=kMC && fDebug > 0) printf("AliCaloTrackReader::GetEventPlaneAngle() -  No EP pointer\n");
     return -1000;
   }
   
@@ -1323,7 +1324,6 @@ void AliCaloTrackReader::GetVertex(Double_t vertex[3], Int_t evtIndex) const
 //________________________________________
 void AliCaloTrackReader::FillVertexArray()
 {
-  
   //Fill data member with vertex
   //In case of Mixed event, multiple vertices
   
@@ -1357,7 +1357,7 @@ void AliCaloTrackReader::FillVertexArray()
       }
       else
       {
-        printf("AliCaloTrackReader::FillVertexArray() - NULL primary vertex\n");
+        AliWarning("NULL primary vertex");
         fVertex[0][0]=0.;   fVertex[0][1]=0.;   fVertex[0][2]=0.;
       }//Primary vertex pointer do not exist
       
@@ -1366,8 +1366,7 @@ void AliCaloTrackReader::FillVertexArray()
       fVertex[0][0]=0.;   fVertex[0][1]=0.;   fVertex[0][2]=0.;
     }
     
-    if(fDebug > 1)
-      printf("AliCaloTrackReader::FillVertexArray() - Single Event Vertex : %f,%f,%f\n",fVertex[0][0],fVertex[0][1],fVertex[0][2]);
+    AliDebug(1,Form("Single Event Vertex : %f,%f,%f\n",fVertex[0][0],fVertex[0][1],fVertex[0][2]));
     
   } else
   { // MultiEvent analysis
@@ -1378,9 +1377,7 @@ void AliCaloTrackReader::FillVertexArray()
       else
          AliWarning("No vertex found");
       
-      if(fDebug > 1)
-        printf("AliCaloTrackReader::FillVertexArray() - Multi Event %d Vertex : %f,%f,%f\n",
-               iev,fVertex[iev][0],fVertex[iev][1],fVertex[iev][2]);
+      AliDebug(1,Form("Multi Event %d Vertex : %f,%f,%f",iev,fVertex[iev][0],fVertex[iev][1],fVertex[iev][2]));
     }
   }
   
@@ -1391,7 +1388,7 @@ void AliCaloTrackReader::FillInputCTS()
 {
   //Return array with Central Tracking System (CTS) tracks
   
-  if(fDebug > 2 ) printf("AliCaloTrackReader::FillInputCTS()\n");
+  AliDebug(1,"Begin");
   
   Double_t pTrack[3] = {0,0,0};
   
@@ -1492,9 +1489,8 @@ void AliCaloTrackReader::FillInputCTS()
     
     if(fCheckFidCut && !fFiducialCut->IsInFiducialCut(fMomentum.Eta(),fMomentum.Phi(),kCTS)) continue;
     
-    if(fDebug > 2 && fMomentum.Pt() > 0.1)
-      printf("AliCaloTrackReader::FillInputCTS() - Selected tracks pt %3.2f, phi %3.2f, eta %3.2f\n",
-             fMomentum.Pt(),fMomentum.Phi()*TMath::RadToDeg(),fMomentum.Eta());
+    AliDebug(2,Form("Selected tracks pt %3.2f, phi %3.2f, eta %3.2f",
+                    fMomentum.Pt(),fMomentum.Phi()*TMath::RadToDeg(),fMomentum.Eta()));
     
     if (fMixedEvent)  track->SetID(itrack);
     
@@ -1509,8 +1505,7 @@ void AliCaloTrackReader::FillInputCTS()
     else      fVertexBC = AliVTrack::kTOFBCNA ;
   }
   
-  if(fDebug > 1)
-    printf("AliCaloTrackReader::FillInputCTS()   - aod entries %d, input tracks %d, pass status %d, multipliticy %d\n", fCTSTracks->GetEntriesFast(), nTracks, nstatus, fTrackMult);//fCTSTracksNormalInputEntries);
+  AliDebug(1,Form("AOD entries %d, input tracks %d, pass status %d, multipliticy %d", fCTSTracks->GetEntriesFast(), nTracks, nstatus, fTrackMult));//fCTSTracksNormalInputEntries);
   
 }
 
@@ -1529,9 +1524,9 @@ void AliCaloTrackReader::FillInputEMCALAlgorithm(AliVCluster * clus, Int_t iclus
   
   clus->GetMomentum(fMomentum, fVertex[vindex]);
   
-  if( (fDebug > 2 && fMomentum.E() > 0.1) || fDebug > 10 )
-    printf("AliCaloTrackReader::FillInputEMCALAlgorithm() - Input cluster E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
-           fMomentum.E(),fMomentum.Pt(),fMomentum.Phi()*TMath::RadToDeg(),fMomentum.Eta());
+  //if( (fDebug > 2 && fMomentum.E() > 0.1) || fDebug > 10 )
+  AliDebug(2,Form("Input cluster E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f",
+                  fMomentum.E(),fMomentum.Pt(),fMomentum.Phi()*TMath::RadToDeg(),fMomentum.Eta()));
 
   if(fRecalculateClusters)
   {
@@ -1585,9 +1580,9 @@ void AliCaloTrackReader::FillInputEMCALAlgorithm(AliVCluster * clus, Int_t iclus
   
   if(!goodCluster)
   {
-    if( (fDebug > 2 && fMomentum.E() > 0.1) || fDebug > 10 )
-      printf("AliCaloTrackReader::FillInputEMCALAlgorithm() - Bad cluster E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
-             fMomentum.E(),fMomentum.Pt(),fMomentum.Phi()*TMath::RadToDeg(),fMomentum.Eta());
+    //if( (fDebug > 2 && fMomentum.E() > 0.1) || fDebug > 10 )
+    AliDebug(2,Form("Bad cluster E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f",
+                    fMomentum.E(),fMomentum.Pt(),fMomentum.Phi()*TMath::RadToDeg(),fMomentum.Eta()));
 
     return;
   }
@@ -1619,9 +1614,9 @@ void AliCaloTrackReader::FillInputEMCALAlgorithm(AliVCluster * clus, Int_t iclus
   {
     GetCaloUtils()->CorrectClusterEnergy(clus) ;
     
-    if( (fDebug > 5 && fMomentum.E() > 0.1) || fDebug > 10 )
-      printf("AliCaloTrackReader::FillInputEMCALAlgorithm() - Correct Non Lin: Old E %3.2f, New E %3.2f\n",
-             fMomentum.E(),clus->E());
+    //if( (fDebug > 5 && fMomentum.E() > 0.1) || fDebug > 10 )
+    AliDebug(5,Form("Correct Non Lin: Old E %3.2f, New E %3.2f",
+                    fMomentum.E(),clus->E()));
 
     // In case of MC analysis, to match resolution/calibration in real data
     // Not needed anymore, just leave for MC studies on systematics
@@ -1629,9 +1624,8 @@ void AliCaloTrackReader::FillInputEMCALAlgorithm(AliVCluster * clus, Int_t iclus
     {
       Float_t rdmEnergy = GetCaloUtils()->GetEMCALRecoUtils()->SmearClusterEnergy(clus);
       
-      if( (fDebug > 5 && fMomentum.E() > 0.1) || fDebug > 10 )
-        printf("AliCaloTrackReader::FillInputEMCALAlgorithm() - Smear energy: Old E %3.2f, New E %3.2f\n",
-               clus->E(),rdmEnergy);
+      //if( (fDebug > 5 && fMomentum.E() > 0.1) || fDebug > 10 )
+      AliDebug(5,Form("Smear energy: Old E %3.2f, New E %3.2f",clus->E(),rdmEnergy));
     
       clus->SetE(rdmEnergy);
     }
@@ -1658,9 +1652,9 @@ void AliCaloTrackReader::FillInputEMCALAlgorithm(AliVCluster * clus, Int_t iclus
   else
     fNNonPileUpClusters++;
   
-  if((fDebug > 2 && fMomentum.E() > 0.1) || fDebug > 10)
-    printf("AliCaloTrackReader::FillInputEMCALAlgorithm() - Selected clusters E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
-           fMomentum.E(),fMomentum.Pt(),fMomentum.Phi()*TMath::RadToDeg(),fMomentum.Eta());
+  //if((fDebug > 2 && fMomentum.E() > 0.1) || fDebug > 10)
+  AliDebug(2,Form("Selected clusters E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f",
+                  fMomentum.E(),fMomentum.Pt(),fMomentum.Phi()*TMath::RadToDeg(),fMomentum.Eta()));
   
   if (fMixedEvent)
     clus->SetID(iclus) ;
@@ -1690,7 +1684,7 @@ void AliCaloTrackReader::FillInputEMCAL()
 {
   //Return array with EMCAL clusters in aod format
   
-  if(fDebug > 2 ) printf("AliCaloTrackReader::FillInputEMCAL()\n");
+  AliDebug(1,"Begin");
   
   // First recalibrate cells, time or energy
   //  if(GetCaloUtils()->IsRecalibrationOn())
@@ -1741,7 +1735,7 @@ void AliCaloTrackReader::FillInputEMCAL()
     
     if(!clusterList)
     {
-      printf("AliCaloTrackReader::FillInputEMCAL() - Wrong name of list with clusters?  <%s>\n",fEMCALClustersListName.Data());
+      AliWarning(Form("Wrong name of list with clusters?  <%s>",fEMCALClustersListName.Data()));
       return;
     }
     
@@ -1751,7 +1745,7 @@ void AliCaloTrackReader::FillInputEMCAL()
       AliVCluster * clus = dynamic_cast<AliVCluster*> (clusterList->At(iclus));
       //printf("E %f\n",clus->E());
       if (clus) FillInputEMCALAlgorithm(clus, iclus);
-      else printf("AliCaloTrackReader::FillInputEMCAL() - Null cluster in list!\n");
+      else      AliWarning("Null cluster in list!");
     }// cluster loop
     
     // Recalculate the pile-up time, in case long time clusters removed during clusterization
@@ -1813,7 +1807,7 @@ void AliCaloTrackReader::FillInputEMCAL()
     
   }
   
-  if(fDebug > 1) printf("AliCaloTrackReader::FillInputEMCAL() - aod entries %d, n pile-up clusters %d, n non pile-up %d \n",  fEMCALClusters->GetEntriesFast(),fNPileUpClusters,fNNonPileUpClusters);
+  AliDebug(1,Form("AOD entries %d, n pile-up clusters %d, n non pile-up %d", fEMCALClusters->GetEntriesFast(),fNPileUpClusters,fNNonPileUpClusters));
   
 }
 
@@ -1822,7 +1816,7 @@ void AliCaloTrackReader::FillInputPHOS()
 {
   //Return array with PHOS clusters in aod format
   
-  if(fDebug > 2 ) printf("AliCaloTrackReader::FillInputPHOS()\n");
+  AliDebug(1,"Begin");
   
   //Loop to select clusters in fiducial cut and fill container with aodClusters
   Int_t nclusters = fInputEvent->GetNumberOfCaloClusters();
@@ -1860,10 +1854,9 @@ void AliCaloTrackReader::FillInputPHOS()
         
         if(fPHOSPtMin > fMomentum.E() || fPHOSPtMax < fMomentum.E())         continue;
         
-        if(fDebug > 2 && fMomentum.E() > 0.1)
-          printf("AliCaloTrackReader::FillInputPHOS() - Selected clusters E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
-                 fMomentum.E(),fMomentum.Pt(),fMomentum.Phi()*TMath::RadToDeg(),fMomentum.Eta());
-        
+        //if(fDebug > 2 && fMomentum.E() > 0.1)
+        AliDebug(2,Form("Selected clusters E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f",
+                        fMomentum.E(),fMomentum.Pt(),fMomentum.Phi()*TMath::RadToDeg(),fMomentum.Eta()));
         
         if (fMixedEvent)
         {
@@ -1876,7 +1869,7 @@ void AliCaloTrackReader::FillInputPHOS()
     }//cluster exists
   }//esd cluster loop
   
-  if(fDebug > 1) printf("AliCaloTrackReader::FillInputPHOS()  - aod entries %d\n",  fPHOSClusters->GetEntriesFast());
+ AliDebug(1,Form("AOD entries %d",fPHOSClusters->GetEntriesFast()));
   
 }
 
@@ -1919,13 +1912,12 @@ void AliCaloTrackReader::FillInputVZERO()
       fV0Mul[0] += (Int_t)v0->GetMultiplicityV0C(i);
       fV0Mul[1] += (Int_t)v0->GetMultiplicityV0A(i);
     }
-    if(fDebug > 0)
-      printf("AliCaloTrackReader::FillInputVZERO() - ADC (%d,%d), Multiplicity (%d,%d) \n",fV0ADC[0],fV0ADC[1],fV0Mul[0],fV0Mul[1]);
+    
+    AliDebug(1,Form("ADC (%d,%d), Multiplicity (%d,%d)",fV0ADC[0],fV0ADC[1],fV0Mul[0],fV0Mul[1]));
   }
   else
   {
-    if(fDebug > 0)
-      printf("AliCaloTrackReader::FillInputVZERO() - Cannot retrieve V0 ESD! Run w/ null V0 charges\n ");
+    AliDebug(1,"Cannot retrieve V0 ESD! Run w/ null V0 charges");
   }
 }
 
@@ -1937,14 +1929,14 @@ void AliCaloTrackReader::FillInputNonStandardJets()
   //
   // Adam T. Matyja
   
-  if(fDebug > 2 ) printf("AliCaloTrackReader::FillInputNonStandardJets()\n");
+  AliDebug(2,"Begin");
+  
   //
   //check if branch name is given
   if(!fInputNonStandardJetBranchName.Length())
   {
-    Printf("No non-standard jet branch name specified. Specify among existing ones.");
     fInputEvent->Print();
-    abort();
+    AliFatal("No non-standard jet branch name specified. Specify among existing ones.");
   }
   
   fNonStandardJets = dynamic_cast<TClonesArray*>(fInputEvent->FindListObject(fInputNonStandardJetBranchName.Data()));
@@ -1952,14 +1944,13 @@ void AliCaloTrackReader::FillInputNonStandardJets()
   if(!fNonStandardJets)
   {
     //check if jet branch exist; exit if not
-    Printf("%s:%d no reconstructed jet array with name %s in AOD", (char*)__FILE__,__LINE__,fInputNonStandardJetBranchName.Data());
     fInputEvent->Print();
-    abort();
+
+    AliFatal(Form("%s:%d no reconstructed jet array with name %s in AOD", (char*)__FILE__,__LINE__,fInputNonStandardJetBranchName.Data()));
   }
   else
   {
-    if(fDebug > 1)
-      printf("AliCaloTrackReader::FillInputNonStandardJets() - aod input jets %d\n", fNonStandardJets->GetEntriesFast() );
+    AliDebug(1,Form("AOD input jets %d", fNonStandardJets->GetEntriesFast()));
   }
   
 }
@@ -1972,14 +1963,14 @@ void AliCaloTrackReader::FillInputBackgroundJets()
   //
   // Adam T. Matyja
   
-  if(fDebug > 2 ) printf("AliCaloTrackReader::FillInputBackgroundJets()\n");
+  AliDebug(1,"Begin");
   //
   //check if branch name is given
   if(!fInputBackgroundJetBranchName.Length())
   {
-    Printf("No background jet branch name specified. Specify among existing ones.");
     fInputEvent->Print();
-    abort();
+    
+    AliFatal("No background jet branch name specified. Specify among existing ones.");
   }
   
   fBackgroundJets = (AliAODJetEventBackground*)(fInputEvent->FindListObject(fInputBackgroundJetBranchName.Data()));
@@ -1987,16 +1978,14 @@ void AliCaloTrackReader::FillInputBackgroundJets()
   if(!fBackgroundJets)
   {
     //check if jet branch exist; exit if not
-    Printf("%s:%d no reconstructed jet array with name %s in AOD", (char*)__FILE__,__LINE__,fInputBackgroundJetBranchName.Data());
     fInputEvent->Print();
-    abort();
+
+    AliFatal(Form("%s:%d no reconstructed jet array with name %s in AOD", (char*)__FILE__,__LINE__,fInputBackgroundJetBranchName.Data()));
   }
   else
   {
-    if(fDebug > 1){
-      printf("AliCaloTrackReader::FillInputBackgroundJets()\n");
-      fBackgroundJets->Print("");
-    }
+    AliDebug(1,"FillInputBackgroundJets");
+    fBackgroundJets->Print("");
   }
   
 }
@@ -2136,8 +2125,8 @@ TArrayI AliCaloTrackReader::GetTriggerPatches(Int_t tmin, Int_t tmax )
     } // trigger iterator
   } // go through triggers
   
-  if(patches.GetSize()<=0) printf("AliCaloTrackReader::GetTriggerPatches() - No patch found! for triggers: %s and selected <%s>\n",
-                                  GetFiredTriggerClasses().Data(),fFiredTriggerClassName.Data());
+  if(patches.GetSize()<=0) AliInfo(Form("No patch found! for triggers: %s and selected <%s>",
+                                        GetFiredTriggerClasses().Data(),fFiredTriggerClassName.Data()));
   //else                     printf(">>>>> N patches %d, test %d,first %d, last %d\n",patches.GetSize(), nPatch, patches.At(0), patches.At(patches.GetSize()-1));
                  
   return patches;
@@ -2525,9 +2514,8 @@ Bool_t  AliCaloTrackReader::RejectLEDEvents()
   
   if(ncellsSM3 >= ncellcut)
   {
-    if(fDebug > 0)
-      printf(" AliCaloTrackReader::FillInputEvent() - reject event with ncells in SM3 %d, cut %d, trig %s\n",
-             ncellsSM3,ncellcut,GetFiredTriggerClasses().Data());
+    AliDebug(1,Form("Reject event with ncells in SM3 %d, cut %d, trig %s",
+                    ncellsSM3,ncellcut,GetFiredTriggerClasses().Data()));
     return kTRUE;
   }
   
@@ -2610,9 +2598,8 @@ void AliCaloTrackReader::SetEventTriggerBit()
   fEventTrigEMCALL1Jet1   = kFALSE;
   fEventTrigEMCALL1Jet2   = kFALSE;
   
-  if(fDebug > 0)
-    printf("AliCaloTrackReader::SetEventTriggerBit() - Select trigger mask bit %d - Trigger Event %s - Select <%s>\n",
-           fEventTriggerMask,GetFiredTriggerClasses().Data(),fFiredTriggerClassName.Data());
+  AliDebug(1,Form("Select trigger mask bit %d - Trigger Event %s - Select <%s>",
+                  fEventTriggerMask,GetFiredTriggerClasses().Data(),fFiredTriggerClassName.Data()));
   
   if(fEventTriggerMask <=0 )// in case no mask set
   {
@@ -2729,11 +2716,10 @@ void AliCaloTrackReader::SetEventTriggerBit()
     }
        }
   
-  if(fDebug > 0 )
-    printf("AliCaloTrackReader::SetEventTriggerBit() - Event bits: \n \t MB   %d, Cen  %d, Sem  %d, L0   %d, L1G1 %d, L1G2 %d, L1J1 %d, L1J2 %d \n",
-           fEventTrigMinBias,      fEventTrigCentral,       fEventTrigSemiCentral,
-           fEventTrigEMCALL0 ,     fEventTrigEMCALL1Gamma1, fEventTrigEMCALL1Gamma2,
-           fEventTrigEMCALL1Jet1 , fEventTrigEMCALL1Jet2);
+  AliDebug(1,Form("Event bits: \n \t MB   %d, Cen  %d, Sem  %d, L0   %d, L1G1 %d, L1G2 %d, L1J1 %d, L1J2 %d",
+                  fEventTrigMinBias,      fEventTrigCentral,       fEventTrigSemiCentral,
+                  fEventTrigEMCALL0 ,     fEventTrigEMCALL1Gamma1, fEventTrigEMCALL1Gamma2,
+                  fEventTrigEMCALL1Jet1 , fEventTrigEMCALL1Jet2));
   
   if(fBitEGA == 0 && fBitEJE ==0)
   {
@@ -2766,8 +2752,8 @@ void AliCaloTrackReader::SetEventTriggerBit()
           fBitEGA = 6;
           fBitEJE = 8;
         }
-           }  else printf("AliCaloTrackReader()::SetEventTriggerBit() - Streamer info for trigger class not available, bit not changed\n");
-    } else printf("AliCaloTrackReader::SetEventTriggerBit() -  Streamer list not available!, bit not changed\n");
+           }  else AliInfo("AliCaloTrackReader()::SetEventTriggerBit() - Streamer info for trigger class not available, bit not changed");
+    } else AliInfo("AliCaloTrackReader::SetEventTriggerBit() -  Streamer list not available!, bit not changed");
     
   } // set once the EJE, EGA trigger bit
   
index ab95f2e..49b9898 100755 (executable)
@@ -43,6 +43,7 @@
 #include "AliOADBContainer.h"
 #include "AliAnalysisManager.h"
 #include "AliAODMCParticle.h"
+#include "AliLog.h"
 
 // --- Detector ---
 #include "AliEMCALGeometry.h"
@@ -123,7 +124,7 @@ void AliCalorimeterUtils::AccessOADB(AliVEvent* event)
   // EMCAL
   if(fOADBForEMCAL)
   {
-    printf("AliCalorimeterUtils::SetOADBParameters() - Get AODB parameters from EMCAL in %s for run %d, and <%s> \n",fOADBFilePathEMCAL.Data(),runnumber,pass.Data());
+    AliInfo(Form("Get AODB parameters from EMCAL in %s for run %d, and <%s>",fOADBFilePathEMCAL.Data(),runnumber,pass.Data()));
     
     Int_t nSM = fEMCALGeo->GetNumberOfSuperModules();
     
@@ -138,7 +139,7 @@ void AliCalorimeterUtils::AccessOADB(AliVEvent* event)
       if(arrayBC)
       {
         SwitchOnDistToBadChannelRecalculation();
-        printf("AliCalorimeterUtils::SetOADBParameters() - Remove EMCAL bad cells \n");
+        AliInfo("Remove EMCAL bad cells");
         
         for (Int_t i=0; i<nSM; ++i) 
         {
@@ -159,7 +160,7 @@ void AliCalorimeterUtils::AccessOADB(AliVEvent* event)
           SetEMCALChannelStatusMap(i,hbm);
           
         } // loop
-      } else printf("AliCalorimeterUtils::SetOADBParameters() - Do NOT remove EMCAL bad channels\n"); // run array
+      } else AliInfo("Do NOT remove EMCAL bad channels\n"); // run array
     }  // Remove bad
     
     // Energy Recalibration
@@ -181,8 +182,8 @@ void AliCalorimeterUtils::AccessOADB(AliVEvent* event)
           
           if(recalib)
           {
-            printf("AliCalorimeterUtils::SetOADBParameters() - Recalibrate EMCAL \n");
-            for (Int_t i=0; i<nSM; ++i) 
+            AliInfo("Recalibrate EMCAL");
+            for (Int_t i=0; i < nSM; ++i)
             {
               TH2F *h = GetEMCALChannelRecalibrationFactors(i);
               
@@ -201,9 +202,9 @@ void AliCalorimeterUtils::AccessOADB(AliVEvent* event)
               
               SetEMCALChannelRecalibrationFactors(i,h);
             } // SM loop
-          }else printf("AliCalorimeterUtils::SetOADBParameters() - Do NOT recalibrate EMCAL, no params object array\n"); // array ok
-        }else printf("AliCalorimeterUtils::SetOADBParameters() - Do NOT recalibrate EMCAL, no params for pass\n"); // array pass ok
-      }else printf("AliCalorimeterUtils::SetOADBParameters() - Do NOT recalibrate EMCAL, no params for run\n");  // run number array ok
+          } else AliInfo("Do NOT recalibrate EMCAL, no params object array"); // array ok
+        } else AliInfo("Do NOT recalibrate EMCAL, no params for pass"); // array pass ok
+      } else AliInfo("Do NOT recalibrate EMCAL, no params for run");  // run number array ok
       
       // once set, apply run dependent corrections if requested
       //fEMCALRecoUtils->SetRunDependentCorrections(runnumber);
@@ -245,7 +246,7 @@ void AliCalorimeterUtils::AccessOADB(AliVEvent* event)
       
       if(htd)
       {
-        printf("AliCalorimeterUtils::SetOADBParameters() - Recalibrate (Temperature) EMCAL \n");
+        AliInfo("Recalibrate (Temperature) EMCAL");
         
         for (Int_t ism=0; ism<nSM; ++ism) 
         {        
@@ -263,7 +264,7 @@ void AliCalorimeterUtils::AccessOADB(AliVEvent* event)
             } // columns
           } // rows 
         } // SM loop
-      }else printf("AliCalorimeterUtils::SetOADBParameters() - Do NOT recalibrate EMCAL with T variations, no params TH1 \n"); 
+      }else AliInfo("Do NOT recalibrate EMCAL with T variations, no params TH1");
     } // Run by Run T calibration    
     
     // Time Recalibration
@@ -283,7 +284,7 @@ void AliCalorimeterUtils::AccessOADB(AliVEvent* event)
         
         if(trecalpass)
         {
-          printf("AliCalorimeterUtils::SetOADBParameters() - Time Recalibrate EMCAL \n");
+          AliInfo("Time Recalibrate EMCAL");
           for (Int_t ibc = 0; ibc < 4; ++ibc) 
           {
             TH1F *h = GetEMCALChannelTimeRecalibrationFactors(ibc);
@@ -303,8 +304,8 @@ void AliCalorimeterUtils::AccessOADB(AliVEvent* event)
             
             SetEMCALChannelTimeRecalibrationFactors(ibc,h);
           } // bunch crossing loop
-        }else printf("AliCalorimeterUtils::SetOADBParameters() - Do NOT recalibrate time EMCAL, no params for pass\n"); // array pass ok
-      }else printf("AliCalorimeterUtils::SetOADBParameters() - Do NOT recalibrate time EMCAL, no params for run\n");  // run number array ok
+        } else AliInfo("Do NOT recalibrate time EMCAL, no params for pass"); // array pass ok
+      } else AliInfo("Do NOT recalibrate time EMCAL, no params for run");  // run number array ok
       
     } // Recalibration on    
     
@@ -313,7 +314,7 @@ void AliCalorimeterUtils::AccessOADB(AliVEvent* event)
   // PHOS
   if(fOADBForPHOS)
   {
-    printf("AliCalorimeterUtils::SetOADBParameters() - Get AODB parameters from PHOS in %s for run %d, and <%s> \n",fOADBFilePathPHOS.Data(),runnumber,pass.Data());
+    AliInfo(Form("Get AODB parameters from PHOS in %s for run %d, and <%s>",fOADBFilePathPHOS.Data(),runnumber,pass.Data()));
     
     // Bad map
     if(fRemoveBadChannels)
@@ -326,11 +327,11 @@ void AliCalorimeterUtils::AccessOADB(AliVEvent* event)
       TObjArray *maps = (TObjArray*)badmapContainer.GetObject(139000,"phosBadMap");
       if(!maps)
       {
-        printf("AliCalorimeterUtils::SetOADBParameters() - Can not read PHOS bad map for run %d.\n",runnumber) ;    
+        AliInfo(Form("Can not read PHOS bad map for run %d",runnumber)) ;
       }
       else
       {
-        printf("AliCalorimeterUtils::SetOADBParameters() - Setting PHOS bad map with name %s \n",maps->GetName()) ;
+        AliInfo(Form("Setting PHOS bad map with name %s",maps->GetName())) ;
         for(Int_t mod=1; mod<5;mod++)
         {
           TH2I *hbmPH = GetPHOSChannelStatusMap(mod);
@@ -369,7 +370,7 @@ void AliCalorimeterUtils::AccessGeometry(AliVEvent* inputEvent)
   {
     if(fLoadEMCALMatrices)
     {
-      printf("AliCalorimeterUtils::AccessGeometry() - Load user defined EMCAL geometry matrices\n");
+      AliInfo("Load user defined EMCAL geometry matrices");
       
       // OADB if available
       AliOADBContainer emcGeoMat("AliEMCALgeo");
@@ -380,9 +381,7 @@ void AliCalorimeterUtils::AccessGeometry(AliVEvent* inputEvent)
       {
         if (!fEMCALMatrix[mod]) // Get it from OADB
         {
-          if(fDebug > 1 ) 
-            printf("AliCalorimeterUtils::AccessGeometry() - EMCAL matrices SM %d, %p\n",
-                   mod,((TGeoHMatrix*) matEMCAL->At(mod)));
+          AliDebug(1,Form("EMCAL matrices SM %d, %p", mod,((TGeoHMatrix*) matEMCAL->At(mod))));
           //((TGeoHMatrix*) matEMCAL->At(mod))->Print();
           
           fEMCALMatrix[mod] = (TGeoHMatrix*) matEMCAL->At(mod) ;
@@ -402,8 +401,7 @@ void AliCalorimeterUtils::AccessGeometry(AliVEvent* inputEvent)
     }//Load matrices
     else if (!gGeoManager) 
     { 
-      if(fDebug > 1) 
-        printf(" AliCalorimeterUtils::AccessGeometry() - Load EMCAL misalignment matrices. \n");
+      AliDebug(1,"Load EMCAL misalignment matrices");
       if(!strcmp(inputEvent->GetName(),"AliESDEvent"))  
       {
         for(Int_t mod=0; mod < (fEMCALGeo->GetEMCGeometry())->GetNumberOfSuperModules(); mod++)
@@ -420,8 +418,7 @@ void AliCalorimeterUtils::AccessGeometry(AliVEvent* inputEvent)
       }//ESD as input
       else 
       {
-        if(fDebug > 1)
-          printf("AliCalorimeterUtils::SetGeometryTransformationMatrices() - Setting of EMCAL transformation matrixes for AODs not implemented yet. \n Import geometry.root file\n");
+        AliDebug(1,"Setting of EMCAL transformation matrixes for AODs not implemented yet. \n Import geometry.root file");
       }//AOD as input
     }//Get matrix from data
     else if(gGeoManager)
@@ -435,7 +432,7 @@ void AliCalorimeterUtils::AccessGeometry(AliVEvent* inputEvent)
   {
     if(fLoadPHOSMatrices)
     {
-      printf("AliCalorimeterUtils::SetGeometryTransformationMatrices() - Load user defined PHOS geometry matrices\n");
+      AliInfo("Load user defined PHOS geometry matrices");
       
       // OADB if available
       AliOADBContainer geomContainer("phosGeo");
@@ -446,9 +443,7 @@ void AliCalorimeterUtils::AccessGeometry(AliVEvent* inputEvent)
       {
         if (!fPHOSMatrix[mod]) // Get it from OADB
         {
-          if(fDebug > 1 ) 
-            printf("AliCalorimeterUtils::SetGeometryTransformationMatrices() - PHOS matrices module %d, %p\n",
-                   mod,((TGeoHMatrix*) matPHOS->At(mod)));
+          AliDebug(1,Form("PHOS matrices module %d, %p",mod,((TGeoHMatrix*) matPHOS->At(mod))));
           //((TGeoHMatrix*) matPHOS->At(mod))->Print();
           
           fPHOSMatrix[mod] = (TGeoHMatrix*) matPHOS->At(mod) ;
@@ -469,8 +464,7 @@ void AliCalorimeterUtils::AccessGeometry(AliVEvent* inputEvent)
     }//Load matrices
     else if (!gGeoManager) 
     { 
-      if(fDebug > 1) 
-        printf(" AliCalorimeterUtils::SetGeometryTransformationMatrices() - Load PHOS misalignment matrices. \n");
+      AliDebug(1,"Load PHOS misalignment matrices.");
                        if(!strcmp(inputEvent->GetName(),"AliESDEvent"))  
       {
                                for(Int_t mod = 0; mod < 5; mod++)
@@ -485,8 +479,7 @@ void AliCalorimeterUtils::AccessGeometry(AliVEvent* inputEvent)
                        }//ESD as input
                        else 
       {
-                               if(fDebug > 1) 
-                                       printf("AliCalorimeterUtils::SetGeometryTransformationMatrices() - Setting of EMCAL transformation matrixes for AODs not implemented yet. \n Import geometry.root file\n");
+                               AliDebug(1,"Setting of EMCAL transformation matrixes for AODs not implemented yet. \n Import geometry.root file");
       }//AOD as input
     }// get matrix from data
     else if(gGeoManager)
@@ -586,9 +579,8 @@ Bool_t AliCalorimeterUtils::CheckCellFiducialRegion(AliVCluster* cluster,
     }
   }
        
-       if(fDebug > 1)
-               printf("AliCalorimeterUtils::CheckCellFiducialRegion() - Cluster Max AbsId %d, Cell Energy %2.2f, Cluster Energy %2.2f\n", 
-           absIdMax, ampMax, cluster->E());
+       AliDebug(1,Form("Cluster Max AbsId %d, Cell Energy %2.2f, Cluster Energy %2.2f",
+                  absIdMax, ampMax, cluster->E()));
        
        if(absIdMax==-1) return kFALSE;
        
@@ -601,8 +593,9 @@ Bool_t AliCalorimeterUtils::CheckCellFiducialRegion(AliVCluster* cluster,
                Int_t iTower = -1, iIphi = -1, iIeta = -1, iphi = -1, ieta = -1, iSM = -1; 
                fEMCALGeo->GetCellIndex(absIdMax,iSM,iTower,iIphi,iIeta); 
                fEMCALGeo->GetCellPhiEtaIndexInSModule(iSM,iTower,iIphi, iIeta,iphi,ieta);
-               if(iSM < 0 || iphi < 0 || ieta < 0 ) {
-      Fatal("CheckCellFidutialRegion","Negative value for super module: %d, or cell ieta: %d, or cell iphi: %d, check EMCAL geometry name\n",iSM,ieta,iphi);
+               if(iSM < 0 || iphi < 0 || ieta < 0 )
+    {
+      AliFatal(Form("Negative value for super module: %d, or cell ieta: %d, or cell iphi: %d, check EMCAL geometry name",iSM,ieta,iphi));
     }
     
                //Check rows/phi
@@ -640,13 +633,9 @@ Bool_t AliCalorimeterUtils::CheckCellFiducialRegion(AliVCluster* cluster,
                        }
                }//eta 0 not checked
     
-               if(fDebug > 1)
-               {
-                       printf("AliCalorimeterUtils::CheckCellFiducialRegion() - EMCAL Cluster in %d cells fiducial volume: ieta %d, iphi %d, SM %d ?",
-             nborder, ieta, iphi, iSM);
-                       if (okcol && okrow ) printf(" YES \n");
-                       else  printf(" NO: column ok? %d, row ok? %d \n",okcol,okrow);
-               }
+    AliDebug(1,Form("EMCAL Cluster in %d cells fiducial volume: ieta %d, iphi %d, SM %d ? ok row %d, ok column %d",
+                    nborder, ieta, iphi, iSM,okrow,okcol));
+
        }//EMCAL
        else if ( cells->GetType() == AliVCaloCells::kPHOSCell )
   {
@@ -656,15 +645,12 @@ Bool_t AliCalorimeterUtils::CheckCellFiducialRegion(AliVCluster* cluster,
                irow = relId[2];
                icol = relId[3];
                //imod = relId[0]-1;
-               if(irow >= fNCellsFromPHOSBorder && irow < 64-fNCellsFromPHOSBorder) okrow =kTRUE; 
+       
+    if(irow >= fNCellsFromPHOSBorder && irow < 64-fNCellsFromPHOSBorder) okrow =kTRUE;
                if(icol >= fNCellsFromPHOSBorder && icol < 56-fNCellsFromPHOSBorder) okcol =kTRUE; 
-               if(fDebug > 1)
-               {
-                       printf("AliCalorimeterUtils::CheckCellFiducialRegion() - PHOS Cluster in %d cells fiducial volume: icol %d, irow %d, Module %d?",
-             fNCellsFromPHOSBorder, icol, irow, relId[0]-1);
-                       if (okcol && okrow ) printf(" YES \n");
-                       else  printf(" NO: column ok? %d, row ok? %d \n",okcol,okrow);
-               }
+
+    AliDebug(1,Form("PHOS Cluster in %d cells fiducial volume: ieta %d, iphi %d, SM %d ? ok row %d, ok column %d",
+                    fNCellsFromPHOSBorder, icol, irow, relId[0]-1,okrow,okcol));
        }//PHOS
        
        if (okcol && okrow) return kTRUE; 
@@ -794,7 +780,7 @@ AliVTrack * AliCalorimeterUtils::GetMatchedTrack(AliVCluster* cluster,
     
     if(trackIndex < 0 )
     { 
-      printf("AliCalorimeterUtils::GetMatchedTrack() - Wrong track index %d, from recalculation\n", trackIndex);
+      AliInfo(Form("Wrong track index %d, from recalculation", trackIndex));
     }
     else 
     {
@@ -813,7 +799,7 @@ AliVTrack * AliCalorimeterUtils::GetMatchedTrack(AliVCluster* cluster,
     
     if(iESDtrack < 0 )
     { 
-      printf("AliCalorimeterUtils::GetMatchedTrack() - Wrong track index %d\n", index);
+      AliWarning(Form("Wrong track index %d", index));
       return 0x0;
     }
     
@@ -835,7 +821,7 @@ Float_t AliCalorimeterUtils::GetMCECellClusFracCorrection(Float_t eCell, Float_t
   // Correction factor for cell energy in cluster to temptatively match Data and MC
   if( eCluster <= 0 || eCluster < eCell )
   {
-    printf("AliCalorimeterUtils::GetMCECellClusFracCorrection() - Bad values eCell=%f, eCluster %f\n",eCell,eCluster);
+    AliWarning(Form("Bad values eCell=%f, eCluster %f",eCell,eCluster));
     return 1;
   }
   
@@ -862,9 +848,8 @@ Int_t AliCalorimeterUtils::GetModuleNumber(AliAODPWG4Particle * particle, AliVEv
   {
                fEMCALGeo->GetAbsCellIdFromEtaPhi(particle->Eta(),particle->Phi(), absId);
     
-               if(fDebug > 2)
-                 printf("AliCalorimeterUtils::GetModuleNumber(PWG4AOD) - EMCAL: cluster eta %f, phi %f, absid %d, SuperModule %d\n",
-             particle->Eta(), particle->Phi()*TMath::RadToDeg(),absId, fEMCALGeo->GetSuperModuleNumber(absId));
+    AliDebug(2,Form("EMCAL: cluster eta %f, phi %f, absid %d, SuperModule %d",
+                    particle->Eta(), particle->Phi()*TMath::RadToDeg(),absId, fEMCALGeo->GetSuperModuleNumber(absId)));
     
                return fEMCALGeo->GetSuperModuleNumber(absId) ;
        }//EMCAL
@@ -884,7 +869,7 @@ Int_t AliCalorimeterUtils::GetModuleNumber(AliAODPWG4Particle * particle, AliVEv
       }
       else
       {
-        Fatal("GetModuleNumber(PWG4AOD)", "Stack not available, stop!");
+        AliFatal("Stack not available, stop!");
       }
       
       if(primary)
@@ -893,7 +878,7 @@ Int_t AliCalorimeterUtils::GetModuleNumber(AliAODPWG4Particle * particle, AliVEv
       }
       else
       {
-        Fatal("GetModuleNumber(PWG4AOD)", "Primary not available, stop!");
+        AliFatal("Primary not available, stop!");
       }
       return mod;
     }
@@ -918,7 +903,7 @@ Int_t AliCalorimeterUtils::GetModuleNumber(AliVCluster * cluster) const
   
   if(!cluster)
   {
-    if(fDebug > 1) printf("AliCalorimeterUtils::GetModuleNumber() - NUL Cluster, please check!!!");
+    AliDebug(1,"AliCalorimeterUtils::GetModuleNumber() - NUL Cluster, please check!!!");
     
     return -1;
   }
@@ -931,7 +916,7 @@ Int_t AliCalorimeterUtils::GetModuleNumber(AliVCluster * cluster) const
   
        if( cluster->IsEMCAL() )
   {
-               if(fDebug > 2) printf("AliCalorimeterUtils::GetModuleNumber() - EMCAL absid %d, SuperModule %d\n",absId, fEMCALGeo->GetSuperModuleNumber(absId));
+               AliDebug(2,Form("EMCAL absid %d, SuperModule %d",absId, fEMCALGeo->GetSuperModuleNumber(absId)));
     
                return fEMCALGeo->GetSuperModuleNumber(absId) ;
        }//EMCAL
@@ -940,7 +925,7 @@ Int_t AliCalorimeterUtils::GetModuleNumber(AliVCluster * cluster) const
                Int_t    relId[4];
     fPHOSGeo->AbsToRelNumbering(absId,relId);
     
-    if(fDebug > 2) printf("AliCalorimeterUtils::GetModuleNumber() - PHOS absid %d Module %d\n",absId, relId[0]-1);
+    AliDebug(2,Form("PHOS absid %d Module %d",absId, relId[0]-1));
     
     return relId[0]-1;
   }//PHOS
@@ -966,7 +951,7 @@ Int_t AliCalorimeterUtils::GetModuleNumberCellIndexes(Int_t absId, Int_t calo,
     
     if(imod < 0 || irow < 0 || icol < 0 )
     {
-      Fatal("GetModuleNumberCellIndexes()","Negative value for super module: %d, or cell icol: %d, or cell irow: %d, check EMCAL geometry name\n",imod,icol,irow);
+      AliFatal(Form("Negative value for super module: %d, or cell icol: %d, or cell irow: %d, check EMCAL geometry name",imod,icol,irow));
     }
     
     //RCU0
@@ -988,9 +973,9 @@ Int_t AliCalorimeterUtils::GetModuleNumberCellIndexes(Int_t absId, Int_t calo,
       iRCU = 0 ;
     }
     
-    if (iRCU<0)
+    if (iRCU < 0)
     {
-      Fatal("GetModuleNumberCellIndexes()","Wrong EMCAL RCU number = %d\n", iRCU);
+      AliFatal(Form("Wrong EMCAL RCU number = %d", iRCU));
     }
     
     return imod ;
@@ -1007,7 +992,7 @@ Int_t AliCalorimeterUtils::GetModuleNumberCellIndexes(Int_t absId, Int_t calo,
     //Int_t iBranch= (Int_t)(relid[3]-1)/28 ; //0 to 1
     if (iRCU >= 4)
     {
-      Fatal("GetModuleNumberCellIndexes()","Wrong PHOS RCU number = %d\n", iRCU);
+      AliFatal(Form("Wrong PHOS RCU number = %d", iRCU));
     }
     return imod;
   }//PHOS
@@ -1161,24 +1146,21 @@ Int_t AliCalorimeterUtils::GetNumberOfLocalMaxima(AliVCluster* cluster, AliVCalo
   
   if(iDigitN == 0)
   {
-    if(fDebug > 0) 
-      printf("AliCalorimeterUtils::GetNumberOfLocalMaxima() - No local maxima found, assign highest energy cell as maxima, id %d, en cell %2.2f, en cluster %2.2f\n",
-             idmax,emax,cluster->E());
+    AliDebug(1,Form("No local maxima found, assign highest energy cell as maxima, id %d, en cell %2.2f, en cluster %2.2f",
+                    idmax,emax,cluster->E()));
     iDigitN      = 1     ;
     maxEList[0]  = emax  ;
     absIdList[0] = idmax ; 
   }
   
-  if(fDebug > 0) 
-  {    
-    printf("AliCalorimeterUtils::GetNumberOfLocalMaxima() - In cluster E %2.2f (wth non lin. %2.2f), M02 %2.2f, M20 %2.2f, N maxima %d \n",
-           cluster->E(),eCluster, cluster->GetM02(),cluster->GetM20(), iDigitN);
   
-    if(fDebug > 1) for(Int_t imax = 0; imax < iDigitN; imax++) 
-    {
-      printf(" \t i %d, absId %d, Ecell %f\n",imax,absIdList[imax],maxEList[imax]);
-    }
-  }
+  AliDebug(1,Form("In cluster E %2.2f (wth non lin. %2.2f), M02 %2.2f, M20 %2.2f, N maxima %d",
+                  cluster->E(),eCluster, cluster->GetM02(),cluster->GetM20(), iDigitN));
+  
+//  if(fDebug > 1) for(Int_t imax = 0; imax < iDigitN; imax++)
+//  {
+//    printf(" \t i %d, absId %d, Ecell %f\n",imax,absIdList[imax],maxEList[imax]);
+//  }
   
   return iDigitN ;
   
@@ -1191,13 +1173,13 @@ TString AliCalorimeterUtils::GetPass()
     
   if (!AliAnalysisManager::GetAnalysisManager()->GetTree()) 
   {
-    AliError("AliCalorimeterUtils::GetPass() - Pointer to tree = 0, returning null\n");
+    AliError("AliCalorimeterUtils::GetPass() - Pointer to tree = 0, returning null");
     return TString("");
   }
   
   if (!AliAnalysisManager::GetAnalysisManager()->GetTree()->GetCurrentFile()) 
   {
-    AliError("AliCalorimeterUtils::GetPass() - Null pointer input file, returning null\n");
+    AliError("AliCalorimeterUtils::GetPass() - Null pointer input file, returning null");
     return TString("");
   }
   
@@ -1210,12 +1192,12 @@ TString AliCalorimeterUtils::GetPass()
   else if (pass.Contains("LHC11c") && pass.Contains("spc_calo") ) return TString("spc_calo");
   else if (pass.Contains("calo") || pass.Contains("high_lumi"))
   {
-    printf("AliCalorimeterUtils::GetPass() - Path contains <calo> or <high-lumi>, set as <pass1>\n");
+    AliInfo("Path contains <calo> or <high-lumi>, set as <pass1>");
     return TString("pass1");
   }
 
   // No condition fullfilled, give a default value
-  printf("AliCalorimeterUtils::GetPass() - Pass number string not found \n");
+  AliInfo("Pass number string not found");
   return TString("");            
   
 }
@@ -1269,7 +1251,8 @@ void AliCalorimeterUtils::InitParameters()
 void AliCalorimeterUtils::InitPHOSBadChannelStatusMap()
 {
   //Init PHOS bad channels map
-  if(fDebug > 0 )printf("AliCalorimeterUtils::InitPHOSBadChannelStatusMap()\n");
+  AliDebug(1,"Init bad channel map");
+  
   //In order to avoid rewriting the same histograms
   Bool_t oldStatus = TH1::AddDirectoryStatus();
   TH1::AddDirectory(kFALSE);
@@ -1288,8 +1271,9 @@ void AliCalorimeterUtils::InitPHOSBadChannelStatusMap()
 void AliCalorimeterUtils::InitPHOSRecalibrationFactors()
 {
        //Init EMCAL recalibration factors
-       if(fDebug > 0 )printf("AliCalorimeterUtils::InitPHOSRecalibrationFactors()\n");
-       //In order to avoid rewriting the same histograms
+  AliDebug(1,"Init recalibration map");
+
+  //In order to avoid rewriting the same histograms
        Bool_t oldStatus = TH1::AddDirectoryStatus();
        TH1::AddDirectory(kFALSE);
   
@@ -1314,47 +1298,44 @@ void AliCalorimeterUtils::InitPHOSRecalibrationFactors()
 //__________________________________________________________
 void AliCalorimeterUtils::InitEMCALGeometry(Int_t runnumber)
 {
-       //Initialize EMCAL geometry if it did not exist previously
+  //Initialize EMCAL geometry if it did not exist previously
+  
+  if (fEMCALGeo) return;
   
-       if (!fEMCALGeo)
+  AliDebug(1,Form(" for run=%d",runnumber));
+  
+  if(fEMCALGeoName=="")
   {
-    if(fEMCALGeoName=="")
-    {
-      if     (runnumber <  140000 && 
-              runnumber >= 100000)   fEMCALGeoName = "EMCAL_FIRSTYEARV1";
-      else if(runnumber >= 140000 &&
-              runnumber <  171000)   fEMCALGeoName = "EMCAL_COMPLETEV1";
-      else                           fEMCALGeoName = "EMCAL_COMPLETE12SMV1";  
-      printf("AliCalorimeterUtils::InitEMCALGeometry() - Set EMCAL geometry name to <%s> for run %d\n",fEMCALGeoName.Data(),runnumber);
-    }
-    
-               fEMCALGeo = AliEMCALGeometry::GetInstance(fEMCALGeoName);
-    
-    // Init geometry, I do not like much to do it like this ...
-    if(fImportGeometryFromFile && !gGeoManager)
+    if     (runnumber <  140000 &&
+            runnumber >= 100000)   fEMCALGeoName = "EMCAL_FIRSTYEARV1";
+    else if(runnumber >= 140000 &&
+            runnumber <  171000)   fEMCALGeoName = "EMCAL_COMPLETEV1";
+    else                           fEMCALGeoName = "EMCAL_COMPLETE12SMV1";
+    AliInfo(Form("Set EMCAL geometry name to <%s> for run %d",fEMCALGeoName.Data(),runnumber));
+  }
+  
+  fEMCALGeo = AliEMCALGeometry::GetInstance(fEMCALGeoName);
+  
+  // Init geometry, I do not like much to do it like this ...
+  if(fImportGeometryFromFile && !gGeoManager)
+  {
+    if(fImportGeometryFilePath=="") // If not specified, set location depending on run number
     {
-      if(fImportGeometryFilePath=="") // If not specified, set location depending on run number
-      {
-        // "$ALICE_ROOT/EVE/alice-data/default_geo.root"
-        if     (runnumber <  140000 &&
-                runnumber >= 100000) fImportGeometryFilePath = "$ALICE_ROOT/OADB/EMCAL/geometry_2010.root";
-        if     (runnumber >= 140000 &&
-                runnumber <  171000) fImportGeometryFilePath = "$ALICE_ROOT/OADB/EMCAL/geometry_2011.root";
-        else                         fImportGeometryFilePath = "$ALICE_ROOT/OADB/EMCAL/geometry_2012.root"; // 2012-2013
-
-      }
-      printf("AliCalorimeterUtils::InitEMCALGeometry() - Import %s\n",fImportGeometryFilePath.Data());
-      TGeoManager::Import(fImportGeometryFilePath) ; // default need file "geometry.root" in local dir!!!!
+      // "$ALICE_ROOT/EVE/alice-data/default_geo.root"
+      if     (runnumber <  140000 &&
+              runnumber >= 100000) fImportGeometryFilePath = "$ALICE_ROOT/OADB/EMCAL/geometry_2010.root";
+      if     (runnumber >= 140000 &&
+              runnumber <  171000) fImportGeometryFilePath = "$ALICE_ROOT/OADB/EMCAL/geometry_2011.root";
+      else                         fImportGeometryFilePath = "$ALICE_ROOT/OADB/EMCAL/geometry_2012.root"; // 2012-2013
+      
     }
-
     
-               if(fDebug > 0)
-    {
-                       printf("AliCalorimeterUtils::InitEMCALGeometry(run=%d)",runnumber);
-                       if (!gGeoManager) printf(" - Careful!, gGeoManager not loaded, load misalign matrices");
-                       printf("\n");
-               }
-       }
+    AliInfo(Form("Import %s",fImportGeometryFilePath.Data()));
+    TGeoManager::Import(fImportGeometryFilePath) ; // default need file "geometry.root" in local dir!!!!
+  }
+  else if (!gGeoManager) AliInfo("Careful!, gGeoManager not loaded, load misalign matrices");
+               
+  
 }
 
 //_________________________________________________________
@@ -1362,19 +1343,16 @@ void AliCalorimeterUtils::InitPHOSGeometry(Int_t runnumber)
 {
        //Initialize PHOS geometry if it did not exist previously
   
-       if (!fPHOSGeo)
-  {
-    if(fPHOSGeoName=="") fPHOSGeoName = "PHOSgeo";
-      
-               fPHOSGeo = new AliPHOSGeoUtils(fPHOSGeoName); 
-    
-               if(fDebug > 0)
-    {
-                       printf("AliCalorimeterUtils::InitPHOSGeometry(run=%d)",runnumber);
-                       if (!gGeoManager) printf(" - Careful!, gGeoManager not loaded, load misalign matrices");
-                       printf("\n");
-               }       
-       }       
+  if (fPHOSGeo) return;
+  
+  AliDebug(1,Form(" for run=%d",runnumber));
+  
+  if(fPHOSGeoName=="") fPHOSGeoName = "PHOSgeo";
+  
+  fPHOSGeo = new AliPHOSGeoUtils(fPHOSGeoName);
+  
+  //if (!gGeoManager) AliInfo("Careful!, gGeoManager not loaded, load misalign matrices");
+       
 }
 
 //_______________________________________________________________________________________________
@@ -1387,7 +1365,7 @@ Bool_t AliCalorimeterUtils::IsMCParticleInCalorimeterAcceptance(Int_t calo, TPar
   if( (!IsPHOSGeoMatrixSet () && calo == kPHOS ) ||
       (!IsEMCALGeoMatrixSet() && calo == kEMCAL)   )
   {
-    AliFatal(Form("Careful Geo Matrix for calo <%d> is not set, use AliFidutialCut instead \n",calo));
+    AliFatal(Form("Careful Geo Matrix for calo <%d> is not set, use AliFidutialCut instead",calo));
     return kFALSE ;
   }
 
@@ -1424,7 +1402,7 @@ Bool_t AliCalorimeterUtils::IsMCParticleInCalorimeterAcceptance(Int_t calo, AliA
   if( (!IsPHOSGeoMatrixSet () && calo == kPHOS ) ||
       (!IsEMCALGeoMatrixSet() && calo == kEMCAL)   )
   {
-    AliFatal(Form("Careful Geo Matrix for calo <%d> is not set, use AliFidutialCut instead \n",calo));
+    AliFatal(Form("Careful Geo Matrix for calo <%d> is not set, use AliFidutialCut instead",calo));
     return kFALSE ;
   }
 
@@ -1466,7 +1444,7 @@ Bool_t AliCalorimeterUtils::IsMCParticleInCalorimeterAcceptance(Int_t calo, Floa
   if( (!IsPHOSGeoMatrixSet () && calo == kPHOS ) ||
       (!IsEMCALGeoMatrixSet() && calo == kEMCAL)   )
   {
-    AliFatal(Form("Careful Geo Matrix for calo <%d> is not set, use AliFidutialCut instead \n",calo));
+    AliFatal(Form("Careful Geo Matrix for calo <%d> is not set, use AliFidutialCut instead",calo));
     return kFALSE ;
   }
 
@@ -1607,26 +1585,24 @@ Float_t AliCalorimeterUtils::RecalibrateClusterEnergy(AliVCluster * cluster,
       Float_t amp = cells->GetCellAmplitude(absId);
       RecalibrateCellAmplitude(amp,calo, absId);
       
-      if(fDebug>2)
-        printf("AliCalorimeterUtils::RecalibrateClusterEnergy() - recalibrate cell: calo <%d>, cell fraction %f, cell energy %f\n",
-               calo,frac,cells->GetCellAmplitude(absId));
+      AliDebug(2,Form("Recalibrate cell: calo <%d>, cell fraction %f, cell energy %f",
+                      calo,frac,cells->GetCellAmplitude(absId)));
       
       energy += amp*frac;
     }
     
-    if(fDebug>1)
-      printf("AliCalorimeterUtils::RecalibrateClusterEnergy() - Energy before %f, after %f\n",cluster->E(),energy);
+    AliDebug(1,Form("Energy before %f, after %f",cluster->E(),energy));
     
        }// cells available
   else
   {
-    Fatal("RecalibrateClusterEnergy()","Cells pointer does not exist!");
+    AliFatal("Cells pointer does not exist!");
   }
   
        return energy;
 }
 
-//__________________________________________________________________________
+//_______________________________________________________________________________________________________
 Float_t AliCalorimeterUtils::RecalibrateClusterEnergyWeightCell(AliVCluster * cluster,
                                                                 AliVCaloCells * cells, Float_t energyOrg)
 {
@@ -1649,8 +1625,8 @@ Float_t AliCalorimeterUtils::RecalibrateClusterEnergyWeightCell(AliVCluster * cl
     if(cluster->IsPHOS()) calo = kPHOS ;
     
     //Loop on the cells, get the cell amplitude and recalibration factor, multiply and and to the new energy
-    for(Int_t icell = 0; icell < ncells; icell++){
-      
+    for(Int_t icell = 0; icell < ncells; icell++)
+    {
       Int_t absId = index[icell];
       
       frac =  fraction[icell];
@@ -1661,20 +1637,18 @@ Float_t AliCalorimeterUtils::RecalibrateClusterEnergyWeightCell(AliVCluster * cl
       
       amp*=GetMCECellClusFracCorrection(amp,energyOrg);
       
-      if(fDebug>2)
-        printf("AliCalorimeterUtils::RecalibrateClusterEnergy() - recalibrate cell: calo <%d>, cell fraction %f, cell energy %f\n",
-               calo,frac,cells->GetCellAmplitude(absId));
+      AliDebug(2,Form("Recalibrate cell: calo <%d>, cell fraction %f, cell energy %f",
+                      calo,frac,cells->GetCellAmplitude(absId)));
       
       energy += amp*frac;
     }
     
-    if(fDebug>1)
-      printf("AliCalorimeterUtils::RecalibrateClusterEnergy() - Energy before %f, after %f\n",cluster->E(),energy);
+    AliDebug(1,Form("Energy before %f, after %f",cluster->E(),energy));
     
        }// cells available
   else
   {
-    Fatal("RecalibrateClusterEnergy()","Cells pointer does not exist!");
+    AliFatal("Cells pointer does not exist!");
   }
   
        return energy;
@@ -1746,7 +1720,7 @@ void AliCalorimeterUtils::SplitEnergy(Int_t absId1, Int_t absId2,
   if(cluster->IsPHOS())
   {
     calorimeter=kPHOS;
-    printf("AliCalorimeterUtils::SplitEnerg() Not supported for PHOS yet \n");
+    AliWarning("Not supported for PHOS yet");
     return;
   }
   
@@ -1865,9 +1839,9 @@ void AliCalorimeterUtils::SplitEnergy(Int_t absId1, Int_t absId2,
     
   }
   
-  if(GetDebug() > 1) printf("AliCalorimeterUtils::SplitEnergy() - n Local Max %d, Cluster energy  = %f, Ecell1 = %f, Ecell2 = %f, Enew1 = %f, Enew2 = %f, Remain %f, \n ncells %d, ncells1 %d, ncells2 %d, f1 %f, f2  %f, sum f12 = %f \n",
-                            nMax, eCluster,ecell1,ecell2,e1,e2,eCluster-e1-e2,ncells,ncells1,ncells2,shareFraction1,shareFraction2,shareFraction1+shareFraction2);
-  
+  AliDebug(1,Form("N Local Max %d, Cluster energy  = %f, Ecell1 = %f, Ecell2 = %f, Enew1 = %f, Enew2 = %f, Remain %f, \n ncells %d, ncells1 %d, ncells2 %d, f1 %f, f2  %f, sum f12 = %f",
+                  nMax, eCluster,ecell1,ecell2,e1,e2,eCluster-e1-e2,ncells,ncells1,ncells2,shareFraction1,shareFraction2,shareFraction1+shareFraction2));
+           
   cluster1->SetE(e1);
   cluster2->SetE(e2);  
   
index 6199173..f729379 100755 (executable)
@@ -42,6 +42,7 @@
 #include "AliCaloTrackReader.h"
 #include "AliMixedEvent.h"
 #include "AliCaloPID.h"
+#include "AliLog.h"
 
 ClassImp(AliIsolationCut)
   
@@ -56,7 +57,7 @@ fSumPtThresholdMax(10000.),
 fPtFraction(0.),
 fICMethod(0),
 fPartInCone(0),
-fDebug(-1),
+fDebug(0),
 fFracIsThresh(1),
 fMomentum(),
 fTrackVector()
@@ -168,8 +169,8 @@ Float_t AliIsolationCut::CalculateExcessAreaFraction(Float_t excess) const
   if(coneA > excessA) return coneA / (coneA-excessA);
   else
   {
-    printf("AliIsolationCut::CalculateExcessAreaFraction() - Please Check : Excess Track %2.3f, coneA %2.2f,  excessA %2.2f, angle %2.2f,factor %2.2f\n",
-           excess,coneA, excessA, angle*TMath::RadToDeg(), coneA / (coneA-excessA));
+    AliWarning(Form("Please Check : Excess Track %2.3f, coneA %2.2f,  excessA %2.2f, angle %2.2f,factor %2.2f",
+                    excess,coneA, excessA, angle*TMath::RadToDeg(), coneA / (coneA-excessA)));
     return  1;
   }
 }
@@ -245,8 +246,7 @@ Float_t AliIsolationCut::GetCellDensity(AliAODPWG4ParticleCorrelation * pCandida
         }
       }//end of cells loop
     }
-    
-    else if(fDebug>0) printf("cluster with bad (eta,phi) in EMCal for energy density calculation\n");
+    else AliWarning("Cluster with bad (eta,phi) in EMCal for energy density calculation");
     
     if (coneCells > 0.) 
     {
@@ -330,8 +330,7 @@ void AliIsolationCut::GetCoeffNormBadCell(AliAODPWG4ParticleCorrelation * pCandi
         }
       }//end of cells loop
     }
-    
-    else if(fDebug > 0) printf("cluster with bad (eta,phi) in EMCal for energy density coeff calculation\n");
+    else AliWarning("Cluster with bad (eta,phi) in EMCal for energy density coeff calculation");
     
     if (coneCells > 0.)
     {
@@ -437,16 +436,9 @@ void  AliIsolationCut::MakeIsolationCut(TObjArray * plCTS,
   nfrac     = 0 ;
   isolated  = kFALSE;
   
-  if(fDebug>0) 
-  {
-    printf("AliIsolationCut::MakeIsolationCut() - Cadidate pT %2.2f, eta %2.2f, phi %2.2f, cone %1.2f, thres %2.2f, Fill AOD? %d",
-           pCandidate->Pt(), pCandidate->Eta(), pCandidate->Phi()*TMath::RadToDeg(), fConeSize,fPtThreshold,bFillAOD);
-    if(plCTS) printf(", nTracks %d"  ,plCTS->GetEntriesFast());
-    if(plNe)  printf(", nClusters %d",plNe ->GetEntriesFast());
-    
-    printf("\n");
-  }
-  
+  AliDebug(1,Form("Candidate pT %2.2f, eta %2.2f, phi %2.2f, cone %1.2f, thres %2.2f, Fill AOD? %d",
+                  pCandidate->Pt(), pCandidate->Eta(), pCandidate->Phi()*TMath::RadToDeg(), fConeSize,fPtThreshold,bFillAOD));
+
   //Initialize the array with refrences
   TObjArray * refclusters  = 0x0;
   TObjArray * reftracks    = 0x0;
@@ -477,7 +469,7 @@ void  AliIsolationCut::MakeIsolationCut(TObjArray * plCTS,
         AliAODPWG4Particle * trackmix = dynamic_cast<AliAODPWG4Particle*>(plCTS->At(ipr)) ; 
         if(!trackmix)
         {
-          printf("AliIsolationCut::MakeIsolationCut() - Wrong track data type, continue\n");
+          AliWarning("Wrong track data type, continue");
           continue;
         }
         
@@ -526,12 +518,11 @@ void  AliIsolationCut::MakeIsolationCut(TObjArray * plCTS,
       // // Check if there is any particle inside cone with pt larger than  fPtThreshold
       // Check if the leading particule inside the cone has a ptLead larger than fPtThreshold
       
-      if( fDebug > 0 )
-        printf("\t track %d, pT %2.2f, eta %1.2f, phi %2.2f, R candidate %2.2f", ipr,pt,eta,phi,rad);
+      AliDebug(2,Form("\t Track %d, pT %2.2f, eta %1.2f, phi %2.2f, R candidate %2.2f", ipr,pt,eta,phi,rad));
       
       if(rad < fConeSize)
       {
-        if(fDebug > 0)  printf(" -  inside candidate cone");
+        AliDebug(2,"Inside candidate cone");
         
         if(bFillAOD)
         {
@@ -574,8 +565,6 @@ void  AliIsolationCut::MakeIsolationCut(TObjArray * plCTS,
 
       } // Inside cone
       
-      if( fDebug > 0 )  printf("\n");
-      
     }// charged particle loop
     
   }//Tracks
@@ -618,7 +607,7 @@ void  AliIsolationCut::MakeIsolationCut(TObjArray * plCTS,
         AliAODPWG4Particle * calomix = dynamic_cast<AliAODPWG4Particle*>(plNe->At(ipr)) ; 
         if(!calomix)
         {
-          printf("AliIsolationCut::MakeIsolationCut() - Wrong calo data type, continue\n");
+          AliWarning("Wrong calo data type, continue");
           continue;
         }
         
@@ -674,12 +663,11 @@ void  AliIsolationCut::MakeIsolationCut(TObjArray * plCTS,
       
       //Check if there is any particle inside cone with pt larger than  fPtThreshold
       
-      if(fDebug > 0 ) 
-        printf("\t cluster %d, pT %2.2f, eta %1.2f, phi %2.2f, R candidate %2.2f", ipr,pt,eta,phi,rad);
+      AliDebug(2,Form("\t Cluster %d, pT %2.2f, eta %1.2f, phi %2.2f, R candidate %2.2f", ipr,pt,eta,phi,rad));
       
       if(rad < fConeSize)
       {
-        if(fDebug > 0 )  printf(" - inside candidate cone");
+        AliDebug(2,"Inside candidate cone");
         
         if(bFillAOD) 
         {
@@ -722,8 +710,6 @@ void  AliIsolationCut::MakeIsolationCut(TObjArray * plCTS,
         
       }//in cone
       
-      if(fDebug > 0 )  printf("\n");
-    
     }// neutral particle loop
     
   }//neutrals
@@ -764,9 +750,8 @@ void  AliIsolationCut::MakeIsolationCut(TObjArray * plCTS,
   {
     if( n == 0 ) isolated = kTRUE ;
     
-    if(fDebug > 0 )
-      printf("pT Cand %2.2f, pT Lead %2.2f, %2.2f<pT Lead< %2.2f, isolated %d\n",
-             ptC,ptLead,fPtThreshold,fPtThresholdMax,isolated);
+    AliDebug(1,Form("pT Cand %2.2f, pT Lead %2.2f, %2.2f<pT Lead< %2.2f, isolated %d",
+                    ptC,ptLead,fPtThreshold,fPtThresholdMax,isolated));
   }
   else if( fICMethod == kSumPtIC )
   {
@@ -776,9 +761,8 @@ void  AliIsolationCut::MakeIsolationCut(TObjArray * plCTS,
     else
       isolated  =  kTRUE  ;
     
-    if(fDebug > 0 )
-      printf("pT Cand %2.2f, SumPt %2.2f, %2.2f<Sum pT< %2.2f, isolated %d\n",
-             ptC,ptLead,fSumPtThreshold,fSumPtThresholdMax,isolated);
+     AliDebug(1,Form("pT Cand %2.2f, SumPt %2.2f, %2.2f<Sum pT< %2.2f, isolated %d",
+                     ptC,ptLead,fSumPtThreshold,fSumPtThresholdMax,isolated));
   }
   else if( fICMethod == kPtFracIC )
   {
index 38c9bc4..c583857 100755 (executable)
@@ -86,6 +86,7 @@ class AliIsolationCut : public TObject {
   Float_t    GetPtFraction()          const { return fPtFraction     ; }
   Int_t      GetICMethod()            const { return fICMethod       ; }
   Int_t      GetParticleTypeInCone()  const { return fPartInCone     ; }
+  Int_t      GetDebug()               const { return fDebug          ; }
   Bool_t     GetFracIsThresh()        const { return fFracIsThresh   ; }
        
   void       SetConeSize(Float_t r)         { fConeSize       = r    ; }
index fc196b9..61b3702 100755 (executable)
@@ -37,6 +37,7 @@
 #include "AliStack.h"
 #include "AliGenPythiaEventHeader.h"
 #include "AliAODMCParticle.h"
+#include "AliLog.h"
 
 ClassImp(AliMCAnalysisUtils)
 
@@ -44,7 +45,7 @@ ClassImp(AliMCAnalysisUtils)
 AliMCAnalysisUtils::AliMCAnalysisUtils() : 
 TObject(), 
 fCurrentEvent(-1), 
-fDebug(-1), 
+fDebug(0),
 fJetsList(new TList), 
 fMCGenerator(kPythia),
 fMCGeneratorString("PYTHIA"),
@@ -158,8 +159,8 @@ Int_t AliMCAnalysisUtils::CheckCommonAncestor(Int_t index1, Int_t index2,
     }// Kine stack from ESDs
   }//First labels not the same
   
-  if((counter1==99 || counter2==99) && fDebug >=0)
-    printf("AliMCAnalysisUtils::CheckCommonAncestor() - Genealogy too large c1: %d, c2= %d\n", counter1, counter2);
+//  if((counter1==99 || counter2==99) && fDebug >=0)
+//    printf("AliMCAnalysisUtils::CheckCommonAncestor() - Genealogy too large c1: %d, c2= %d\n", counter1, counter2);
   //printf("CheckAncestor:\n");
   
   Int_t commonparents = 0;
@@ -225,8 +226,9 @@ Int_t AliMCAnalysisUtils::CheckOrigin(const Int_t * label, Int_t nlabels,
   
   Int_t tag = 0;
   
-  if(nlabels<=0) {
-    printf("AliMCAnalysisUtils::CheckOrigin(nlabel<=0) - No MC labels available, please check!!!\n");
+  if( nlabels <= 0 )
+  {
+    AliWarning("No MC labels available, please check!!!");
     return kMCBadLabel;
   }
   
@@ -252,8 +254,9 @@ Int_t AliMCAnalysisUtils::CheckOrigin(Int_t label, const AliCaloTrackReader* rea
   
   Int_t tag = 0;
   
-  if(label<0) {
-    printf("AliMCAnalysisUtils::CheckOrigin(label<0) - No MC labels available, please check!!!\n");
+  if( label < 0 )
+  {
+    AliWarning("No MC labels available, please check!!!");
     return kMCBadLabel;
   }
   
@@ -288,8 +291,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInStack(const Int_t *labels, Int_t nlabels,
   
   if(!stack)
   {
-    if (fDebug >=0) 
-      printf("AliMCAnalysisUtils::CheckOriginInStack() - Stack is not available, check analysis settings in configuration file, STOP!!\n");
+    AliDebug(1,"Stack is not available, check analysis settings in configuration file, STOP!!");
     return -1;
   }
   
@@ -306,7 +308,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInStack(const Int_t *labels, Int_t nlabels,
     Int_t mStatus  = mom->GetStatusCode() ;
     Int_t iParent  = mom->GetFirstMother() ;
     
-    if(fDebug > 0 && label < 8 && fMCGenerator != kBoxLike) printf("AliMCAnalysisUtils::CheckOriginInStack() - Mother is parton %d\n",iParent);
+    //if( label < 8 && fMCGenerator != kBoxLike ) AliDebug(1,Form("AliMCAnalysisUtils::CheckOriginInStack() - Mother is parton %d",iParent));
     
     //GrandParent of the entity
     TParticle * parent = NULL;
@@ -322,15 +324,12 @@ Int_t AliMCAnalysisUtils::CheckOriginInStack(const Int_t *labels, Int_t nlabels,
         pStatus = parent->GetStatusCode();  
       }
     }
-    else if(fDebug > 0 ) printf("AliMCAnalysisUtils::CheckOriginInStack() - Parent with label %d\n",iParent);
+    else AliDebug(1,Form("Parent with label %d",iParent));
+    
+    AliDebug(2,"Cluster most contributing mother and its parent:");
+    AliDebug(2,Form("\t Mother label %d, pdg %d, status %d",iMom, mPdg, mStatus));
+    AliDebug(2,Form("\t Parent label %d, pdg %d, status %d",iParent, pPdg, pStatus));
     
-    if(fDebug > 2 )
-    {
-      printf("AliMCAnalysisUtils::CheckOriginInStack() - Cluster most contributing mother and its parent: \n");
-      printf("\t Mother label %d, pdg %d, status %d\n",iMom, mPdg, mStatus);
-      printf("\t Parent label %d, pdg %d, status %d\n",iParent, pPdg, pStatus);
-    }
-         
     //Check if "mother" of entity is converted, if not, get the first non converted mother
     if((mPdg == 22 || mPdg == 11) && (pPdg == 22 || pPdg == 11) && mStatus == 0)
     {
@@ -347,7 +346,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInStack(const Int_t *labels, Int_t nlabels,
         mStatus  = mom->GetStatusCode() ;
         iParent  = mom->GetFirstMother() ;
        
-        if(fDebug > 0 && label < 8 ) printf("AliMCAnalysisUtils::CheckOriginInStack() - Mother is parton %d\n",iParent);
+        //if(label < 8 ) AliDebug(1,Form("AliMCAnalysisUtils::CheckOriginInStack() - Mother is parton %d\n",iParent));
         
         //GrandParent
         if(iParent >= 0)
@@ -366,12 +365,9 @@ Int_t AliMCAnalysisUtils::CheckOriginInStack(const Int_t *labels, Int_t nlabels,
         }
       }//while
       
-      if(fDebug > 2 )
-      {
-        printf("AliMCAnalysisUtils::CheckOriginInStack() - Converted photon/electron: \n");
-        printf("\t Mother label %d, pdg %d, status %d\n",iMom, mPdg, mStatus);
-        printf("\t Parent label %d, pdg %d, status %d\n",iParent, pPdg, pStatus);
-      }
+      AliDebug(2,"Converted photon/electron:");
+      AliDebug(2,Form("\t Mother label %d, pdg %d, status %d",iMom, mPdg, mStatus));
+      AliDebug(2,Form("\t Parent label %d, pdg %d, status %d",iParent, pPdg, pStatus));
       
     }//mother and parent are electron or photon and have status 0
     else if((mPdg == 22 || mPdg == 11) && mStatus == 0)
@@ -386,11 +382,9 @@ Int_t AliMCAnalysisUtils::CheckOriginInStack(const Int_t *labels, Int_t nlabels,
         mPdgSign = mom->GetPdgCode();
         mPdg     = TMath::Abs(mPdgSign);
         
-        if(fDebug > 2 )
-        {
-          printf("AliMCAnalysisUtils::CheckOriginInStack() - Converted hadron: \n");
-          printf("\t Mother label %d, pdg %d, status %d\n",iMom, mPdg, mStatus);
-        }
+        AliDebug(2,"Converted hadron:");
+        AliDebug(2,Form("\t Mother label %d, pdg %d, status %d",iMom, mPdg, mStatus));
+        
       }//hadron converted
       
       //Comment for the next lines, we do not check the parent of the hadron for the moment.
@@ -420,8 +414,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInStack(const Int_t *labels, Int_t nlabels,
       
       SetTagBit(tag,kMCPi0Decay);
       
-      if(fDebug > 2 )
-        printf("AliMCAnalysisUtils::CheckOriginInStack() - First mother is directly pi0, not decayed by generator \n");
+      AliDebug(2,"First mother is directly pi0, not decayed by generator");
       
       CheckOverlapped2GammaDecay(labels,nlabels, iMom, stack, tag); //set to kMCPi0 if 2 gammas in same cluster
       
@@ -430,8 +423,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInStack(const Int_t *labels, Int_t nlabels,
     {
       SetTagBit(tag,kMCEtaDecay);
       
-      if(fDebug > 2 )
-        printf("AliMCAnalysisUtils::CheckOriginInStack() - First mother is directly eta, not decayed by generator \n");
+      AliDebug(2,"First mother is directly eta, not decayed by generator");
       
       CheckOverlapped2GammaDecay(labels,nlabels, iMom, stack, tag); //set to kMCEta if 2 gammas in same cluster
     }
@@ -444,7 +436,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInStack(const Int_t *labels, Int_t nlabels,
       {
         SetTagBit(tag,kMCPi0Decay);
         
-        if(fDebug > 2 ) printf("AliMCAnalysisUtils::CheckOriginInStack() - PYTHIA pi0 decay photon,  parent pi0 with status 11 \n");
+        AliDebug(2,"PYTHIA pi0 decay photon,  parent pi0 with status 11");
         
         CheckOverlapped2GammaDecay(labels,nlabels, iParent, stack, tag); //set to kMCPi0 if 2 gammas in same cluster
         // In case it did not merge, check if the decay companion is lost
@@ -457,7 +449,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInStack(const Int_t *labels, Int_t nlabels,
       {
         SetTagBit(tag, kMCEtaDecay);
         
-        if(fDebug > 2 ) printf("AliMCAnalysisUtils::CheckOriginInStack() - PYTHIA eta decay photon,  parent pi0 with status 11 \n");
+        AliDebug(2,"PYTHIA eta decay photon,  parent pi0 with status 11");
         
         CheckOverlapped2GammaDecay(labels,nlabels, iParent, stack, tag);//set to kMCEta if 2 gammas in same cluster
         // In case it did not merge, check if the decay companion is lost
@@ -527,12 +519,13 @@ Int_t AliMCAnalysisUtils::CheckOriginInStack(const Int_t *labels, Int_t nlabels,
       
       SetTagBit(tag,kMCElectron);
       
-      if(fDebug > 0) printf("AliMCAnalysisUtils::CheckOriginInStack() - Checking ancestors of electrons\n");
+      AliDebug(1,"Checking ancestors of electrons");
      
       if      (pPdg == 111) { SetTagBit(tag,kMCPi0Decay); SetTagBit(tag, kMCDecayDalitz) ; } //Pi0 Dalitz decay
       else if (pPdg == 221) { SetTagBit(tag,kMCEtaDecay); SetTagBit(tag, kMCDecayDalitz) ; } //Eta Dalitz decay
       else if((499 < pPdg && pPdg < 600)||(4999 < pPdg && pPdg < 6000)) { SetTagBit(tag,kMCEFromB); } //b-->e decay
-      else if((399 < pPdg && pPdg < 500)||(3999 < pPdg && pPdg < 5000)) { //check charm decay
+      else if((399 < pPdg && pPdg < 500)||(3999 < pPdg && pPdg < 5000))
+      { //check charm decay
         if(parent)
         {
           Int_t iGrandma = parent->GetFirstMother();
@@ -541,7 +534,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInStack(const Int_t *labels, Int_t nlabels,
             TParticle* gma = (TParticle*)stack->Particle(iGrandma); //get mother of charm
             Int_t gPdg = TMath::Abs(gma->GetPdgCode());
             if((499 < gPdg && gPdg < 600)||(4999 < gPdg && gPdg < 6000)) SetTagBit(tag,kMCEFromCFromB); //b-->c-->e
-            else SetTagBit(tag,kMCEFromC); //c-->e 
+            else SetTagBit(tag,kMCEFromC); //c-->e
           }
           else SetTagBit(tag,kMCEFromC); //c-->e
         }//parent
@@ -552,27 +545,26 @@ Int_t AliMCAnalysisUtils::CheckOriginInStack(const Int_t *labels, Int_t nlabels,
         if(pPdg > 10000) SetTagBit(tag,kMCUnknown);
         
         else SetTagBit(tag,kMCOtherDecay);
-       
-        if(fDebug > 0 && parent) printf("AliMCAnalysisUtils::CheckOriginInStack() - Status %d Electron from other origin: %s (pPdg = %d) %s (mpdg = %d)\n",mStatus,parent->GetName(),pPdg,mom->GetName(),mPdg);
+        
+        //if(parent) AliDebug(1,Form("Status %d Electron from other origin: %s (pPdg = %d) %s (mpdg = %d)",mStatus,parent->GetName(),pPdg,mom->GetName(),mPdg));
       }
     }//electron check
     //Cluster was made by something else
     else
     {
-      if(fDebug > 0)
-        printf("AliMCAnalysisUtils::CheckOriginInStack() - \tSetting kMCUnknown for cluster from %s (pdg = %d, Parent pdg = %d)\n",
-               mom->GetName(),mPdg,pPdg);
+      AliDebug(2,Form("\t Setting kMCUnknown for cluster from %s (pdg = %d, Parent pdg = %d)",
+                      mom->GetName(),mPdg,pPdg));
       
       SetTagBit(tag,kMCUnknown);
     }
   }//Good label value
   else
   {// Bad label
-    if(label < 0 && (fDebug >= 0)) 
-      printf("AliMCAnalysisUtils::CheckOriginInStack() *** bad label or no stack ***:  label %d \n", label);
+    if(label < 0)
+      AliWarning(Form("*** bad label or no stack ***:  label %d ", label));
     
-    if(label >=  stack->GetNtrack() &&  (fDebug >= 0))
-      printf("AliMCAnalysisUtils::CheckOriginInStack() *** large label ***:  label %d, n tracks %d \n", label, stack->GetNtrack());
+    if(label >=  stack->GetNtrack())
+      AliWarning(Form("*** large label ***:  label %d, n tracks %d", label, stack->GetNtrack()));
     
     SetTagBit(tag,kMCUnknown);
   }//Bad label
@@ -591,8 +583,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInAOD(const Int_t *labels, Int_t nlabels,
   
   if(!mcparticles)
   {
-    if(fDebug >= 0)
-      printf("AliMCAnalysisUtils::CheckOriginInAOD() - AODMCParticles is not available, check analysis settings in configuration file!!\n");
+    AliDebug(1,"AODMCParticles is not available, check analysis settings in configuration file!!");
     return -1;
   }
        
@@ -609,7 +600,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInAOD(const Int_t *labels, Int_t nlabels,
     Int_t mPdg     = TMath::Abs(mPdgSign);
     Int_t iParent  = mom->GetMother() ;
     
-    if(fDebug > 0 && label < 8 && fMCGenerator != kBoxLike) printf("AliMCAnalysisUtils::CheckOriginInAOD() - Mother is parton %d\n",iParent);
+    //if(label < 8 && fMCGenerator != kBoxLike) AliDebug(1,Form("Mother is parton %d\n",iParent));
     
     //GrandParent
     AliAODMCParticle * parent = NULL ;
@@ -619,20 +610,12 @@ Int_t AliMCAnalysisUtils::CheckOriginInAOD(const Int_t *labels, Int_t nlabels,
       parent = (AliAODMCParticle *) mcparticles->At(iParent);
       pPdg = TMath::Abs(parent->GetPdgCode());
     }
-    else if(fDebug > 0 ) printf("AliMCAnalysisUtils::CheckOriginInAOD() - Parent with label %d\n",iParent);
+    else AliDebug(1,Form("Parent with label %d",iParent));
+    
+    AliDebug(2,"Cluster most contributing mother and its parent:");
+    AliDebug(2,Form("\t Mother label %d, pdg %d, Primary? %d, Physical Primary? %d",iMom, mPdg, mom->IsPrimary(), mom->IsPhysicalPrimary()));
+    AliDebug(2,Form("\t Parent label %d, pdg %d, Primary? %d, Physical Primary? %d",iParent, pPdg, parent->IsPrimary(), parent->IsPhysicalPrimary()));
     
-    if(fDebug > 2 )
-    {
-      printf("AliMCAnalysisUtils::CheckOriginInAOD() - Cluster most contributing mother and its parent: \n");
-      
-      printf("\t Mother label %d, pdg %d, Primary? %d, Physical Primary? %d\n",
-             iMom, mPdg, mom->IsPrimary(), mom->IsPhysicalPrimary());
-      
-      if(parent)
-        printf("\t Parent label %d, pdg %d, Primary? %d, Physical Primary? %d\n",
-               iParent, pPdg, parent->IsPrimary(), parent->IsPhysicalPrimary());
-    }
-         
     //Check if mother is converted, if not, get the first non converted mother
     if((mPdg == 22 || mPdg == 11) && (pPdg == 22 || pPdg == 11) && !mom->IsPrimary())
     {
@@ -647,7 +630,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInAOD(const Int_t *labels, Int_t nlabels,
         mPdgSign = mom->GetPdgCode();
         mPdg     = TMath::Abs(mPdgSign);
         iParent  = mom->GetMother() ;
-        if(fDebug > 0 && label < 8 ) printf("AliMCAnalysisUtils::CheckOriginInAOD() - Mother is parton %d\n",iParent);
+        //if(label < 8 ) AliDebug(1, Form("AliMCAnalysisUtils::CheckOriginInAOD() - Mother is parton %d\n",iParent));
         
         //GrandParent
         if(iParent >= 0 && parent)
@@ -660,17 +643,9 @@ Int_t AliMCAnalysisUtils::CheckOriginInAOD(const Int_t *labels, Int_t nlabels,
         
       }//while 
       
-      if(fDebug > 2 )
-      {
-        printf("AliMCAnalysisUtils::CheckOriginInAOD() - Converted photon/electron : \n");
-        
-        printf("\t Mother label %d, pdg %d, Primary? %d, Physical Primary? %d\n",
-               iMom, mPdg, mom->IsPrimary(), mom->IsPhysicalPrimary());
-       
-        if(parent)
-          printf("\t Parent label %d, pdg %d, Primary? %d, Physical Primary? %d\n",
-                 iParent, pPdg, parent->IsPrimary(), parent->IsPhysicalPrimary());
-      }
+      AliDebug(2,"AliMCAnalysisUtils::CheckOriginInAOD() - Converted photon/electron:");
+      AliDebug(2,Form("\t Mother label %d, pdg %d, Primary? %d, Physical Primary? %d",iMom, mPdg, mom->IsPrimary(), mom->IsPhysicalPrimary()));
+      AliDebug(2,Form("\t Parent label %d, pdg %d, Primary? %d, Physical Primary? %d",iParent, pPdg, parent->IsPrimary(), parent->IsPhysicalPrimary()));
       
     }//mother and parent are electron or photon and have status 0 and parent is photon or electron
     else if((mPdg == 22 || mPdg == 11) && !mom->IsPrimary())
@@ -685,11 +660,9 @@ Int_t AliMCAnalysisUtils::CheckOriginInAOD(const Int_t *labels, Int_t nlabels,
         mPdgSign = mom->GetPdgCode();
         mPdg     = TMath::Abs(mPdgSign);
         
-        if(fDebug > 2 )
-        {
-          printf("AliMCAnalysisUtils::CheckOriginInAOD() - Converted hadron : \n");
-          printf("\t Mother label %d, pdg %d, Primary? %d, Physical Primary? %d\n",iMom, mPdg, mom->IsPrimary(), mom->IsPhysicalPrimary());
-        }
+       AliDebug(2,"AliMCAnalysisUtils::CheckOriginInAOD() - Converted hadron:");
+       AliDebug(2,Form("\t Mother label %d, pdg %d, Primary? %d, Physical Primary? %d",iMom, mPdg, mom->IsPrimary(), mom->IsPhysicalPrimary()));
+      
       }//hadron converted
       
       //Comment for next lines, we do not check the parent of the hadron for the moment.
@@ -721,7 +694,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInAOD(const Int_t *labels, Int_t nlabels,
     {
       SetTagBit(tag,kMCPi0Decay);
       
-      if(fDebug > 2 ) printf("AliMCAnalysisUtils::CheckOriginInAOD() - First mother is directly pi0, not decayed by generator \n");
+      AliDebug(2,"First mother is directly pi0, not decayed by generator");
       
       CheckOverlapped2GammaDecay(labels,nlabels, iMom, mcparticles, tag); //set to kMCPi0 if 2 gammas in same cluster
     }
@@ -729,7 +702,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInAOD(const Int_t *labels, Int_t nlabels,
     {
       SetTagBit(tag,kMCEtaDecay);   
       
-      if(fDebug > 2 ) printf("AliMCAnalysisUtils::CheckOriginInAOD() - First mother is directly eta, not decayed by generator \n");
+      AliDebug(2,"First mother is directly eta, not decayed by generator");
       
       CheckOverlapped2GammaDecay(labels,nlabels, iMom, mcparticles, tag); //set to kMCEta if 2 gammas in same cluster
     }
@@ -742,7 +715,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInAOD(const Int_t *labels, Int_t nlabels,
       {
         SetTagBit(tag,kMCPi0Decay);
         
-        if(fDebug > 2 ) printf("AliMCAnalysisUtils::CheckOriginInAOD() - Generator pi0 decay photon \n");
+        AliDebug(2,"Generator pi0 decay photon");
         
         CheckOverlapped2GammaDecay(labels,nlabels, iParent, mcparticles, tag); //set to kMCPi0 if 2 gammas in same cluster
         // In case it did not merge, check if the decay companion is lost
@@ -757,7 +730,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInAOD(const Int_t *labels, Int_t nlabels,
       {
         SetTagBit(tag, kMCEtaDecay);
         
-        if(fDebug > 2 ) printf("AliMCAnalysisUtils::CheckOriginInAOD() - Generator eta decay photon \n");
+        AliDebug(2,"Generator eta decay photon");
         
         CheckOverlapped2GammaDecay(labels,nlabels, iParent, mcparticles, tag); //set to kMCEta if 2 gammas in same cluster
         // In case it did not merge, check if the decay companion is lost
@@ -799,7 +772,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInAOD(const Int_t *labels, Int_t nlabels,
       
       SetTagBit(tag,kMCElectron);
       
-      if(fDebug > 0) printf("AliMCAnalysisUtils::CheckOriginInAOD() - Checking ancestors of electrons");
+      AliDebug(1,"Checking ancestors of electrons");
       
       if      (pPdg == 111) { SetTagBit(tag,kMCPi0Decay); SetTagBit(tag,kMCDecayDalitz);} //Pi0 Dalitz decay
       else if (pPdg == 221) { SetTagBit(tag,kMCEtaDecay); SetTagBit(tag,kMCDecayDalitz);} //Eta Dalitz decay
@@ -823,8 +796,7 @@ Int_t AliMCAnalysisUtils::CheckOriginInAOD(const Int_t *labels, Int_t nlabels,
         TParticlePDG* foo = TDatabasePDG::Instance()->GetParticle(pPdg);
         TParticlePDG* foo1 = TDatabasePDG::Instance()->GetParticle(mPdg);
         
-        if(fDebug > 0) printf("AliMCAnalysisUtils::CheckOriginInAOD() - Electron from other origin: %s (pPdg = %d) %s (mPdg = %d)\n",
-                              foo->GetName(), pPdg,foo1->GetName(),mPdg);
+        AliDebug(1,Form("Electron from other origin: %s (pPdg = %d) %s (mPdg = %d)",foo->GetName(), pPdg,foo1->GetName(),mPdg));
         
         if(pPdg > 10000) SetTagBit(tag,kMCUnknown);
         else             SetTagBit(tag,kMCOtherDecay);
@@ -833,18 +805,17 @@ Int_t AliMCAnalysisUtils::CheckOriginInAOD(const Int_t *labels, Int_t nlabels,
     //cluster was made by something else
     else
     {
-      if(fDebug > 0) printf("AliMCAnalysisUtils::CheckOriginInAOD() - \tSetting kMCUnknown for cluster with pdg = %d, Parent pdg = %d\n",
-                            mPdg,pPdg);
+      AliDebug(1,Form("\t Setting kMCUnknown for cluster with pdg = %d, Parent pdg = %d",mPdg,pPdg));
       SetTagBit(tag,kMCUnknown);
     }
   }//Good label value
   else
   {//Bad label
-    if(label < 0 && (fDebug >= 0) ) 
-      printf("AliMCAnalysisUtils::CheckOriginInAOD() *** bad label or no mcparticles ***:  label %d \n", label);
+    if(label < 0)
+      AliWarning(Form("*** bad label or no mcparticles ***:  label %d", label));
     
-    if(label >=  mcparticles->GetEntriesFast() &&  (fDebug >= 0) )
-      printf("AliMCAnalysisUtils::CheckOriginInAOD() *** large label ***:  label %d, n tracks %d \n", label, mcparticles->GetEntriesFast());
+    if(label >=  mcparticles->GetEntriesFast())
+      AliWarning(Form("*** large label ***:  label %d, n tracks %d", label, mcparticles->GetEntriesFast()));
   
     SetTagBit(tag,kMCUnknown);
     
@@ -861,9 +832,9 @@ void AliMCAnalysisUtils::CheckOverlapped2GammaDecay(const Int_t *labels,    Int_
 {
   // Check if cluster is formed from the contribution of 2 decay photons from pi0 or eta. Input in stack.
   
-  if(labels[0] < 0 || labels[0] > stack->GetNtrack() || nlabels <= 1) {
-    if(fDebug > 2) printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(ESD) - Exit : label[0] %d, n primaries %d, nlabels %d \n",
-                          labels[0],stack->GetNtrack(), nlabels);
+  if(labels[0] < 0 || labels[0] > stack->GetNtrack() || nlabels <= 1)
+  {
+    AliDebug(2,Form("Exit : label[0] %d, n primaries %d, nlabels %d",labels[0],stack->GetNtrack(), nlabels));
     return;
   }
   
@@ -871,17 +842,16 @@ void AliMCAnalysisUtils::CheckOverlapped2GammaDecay(const Int_t *labels,    Int_
   Int_t mesonPdg    = meson->GetPdgCode();
   if(mesonPdg!=111 && mesonPdg!=221)
   {
-    printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(ESD) - Wrong pi0/eta PDG : %d \n",mesonPdg);
+    AliWarning(Form("Wrong pi0/eta PDG : %d",mesonPdg));
     return;
   }
   
-  if(fDebug > 2) printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(ESD) - %s, label %d\n",meson->GetName(), mesonIndex);
+  AliDebug(2,Form("%s, label %d",meson->GetName(), mesonIndex));
   
   //Check if meson decayed into 2 daughters or if both were kept.
   if(meson->GetNDaughters() != 2)
   {
-    if(fDebug > 2) 
-      printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(ESD) - Not overalapped. Number of daughters is %d, not 2 \n",meson->GetNDaughters());
+    AliDebug(2,Form("Not overalapped. Number of daughters is %d, not 2",meson->GetNDaughters()));
     return;
   }
   
@@ -894,25 +864,23 @@ void AliMCAnalysisUtils::CheckOverlapped2GammaDecay(const Int_t *labels,    Int_
   //Check if both daughters are photons
   if(photon0->GetPdgCode() != 22 || photon1->GetPdgCode()!=22)
   {
-    if(fDebug > 2) 
-      printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(ESD) - Not overalapped. PDG:  daughter 1 = %d, of daughter 2 = %d \n",photon0->GetPdgCode(),photon1->GetPdgCode());
+    AliDebug(2,Form("Not overalapped. PDG:  daughter 1 = %d, of daughter 2 = %d",photon0->GetPdgCode(),photon1->GetPdgCode()));
     return;
   }
   
-  if(fDebug > 2) 
-    printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(ESD) - Daughter labels : photon0 = %d, photon1 = %d \n",iPhoton0,iPhoton1);
+  AliDebug(2,Form("Daughter labels : photon0 = %d, photon1 = %d",iPhoton0,iPhoton1));
   
   //Check if both photons contribute to the cluster
   Bool_t okPhoton0 = kFALSE;
   Bool_t okPhoton1 = kFALSE;
   
-  if(fDebug > 3) printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(ESD) - Labels loop:\n");
+  AliDebug(3,"Labels loop:");
   
   Bool_t conversion = kFALSE;
   
   for(Int_t i = 0; i < nlabels; i++)
   {
-    if(fDebug > 3) printf("\t  at begin:label %d/%d: %d, ok? photon1 %d, photon2 %d\n", i+1, nlabels, labels[i], okPhoton0, okPhoton1);
+    AliDebug(3,Form("\t  at begin:label %d/%d: %d, ok? photon1 %d, photon2 %d", i+1, nlabels, labels[i], okPhoton0, okPhoton1));
     
     //If we already found both, break the loop
     if(okPhoton0 && okPhoton1) break;
@@ -933,18 +901,19 @@ void AliMCAnalysisUtils::CheckOverlapped2GammaDecay(const Int_t *labels,    Int_
     
     if(index >= stack->GetNtrack())
     {
-      printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(ESD) Particle index %d larger than size of list %d !!\n",
-             index,stack->GetNtrack());
+      AliWarning(Form("Particle index %d larger than size of list %d!!",index,stack->GetNtrack()));
       continue;
     }
     
     TParticle * daught = stack->Particle(index);
     Int_t tmpindex = daught->GetFirstMother();         
-    if(fDebug > 3) printf("\t Conversion? : mother %d\n",tmpindex);
+    
+    AliDebug(3,Form("\t Conversion? : mother %d",tmpindex));
+    
     while(tmpindex>=0)
     {
       //MC particle of interest is the mother
-      if(fDebug > 3) printf("\t \t parent index %d\n",tmpindex);
+      AliDebug(3,Form("\t \t parent index %d",tmpindex));
       daught   = stack->Particle(tmpindex);
       if      (iPhoton0 == tmpindex)
       {
@@ -963,16 +932,15 @@ void AliMCAnalysisUtils::CheckOverlapped2GammaDecay(const Int_t *labels,    Int_
       
     }//While to check if pi0/eta daughter was one of these contributors to the cluster
     
-    if(i == 0 && (!okPhoton0 && !okPhoton1) && fDebug>=0) 
-      printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(ESD) - Something happens, first label should be from a photon decay!\n");
+    //if(i == 0 && (!okPhoton0 && !okPhoton1))
+    //  AliDebug(1,Form("Something happens, first label should be from a photon decay!"));
     
   }//loop on list of labels
   
   //If both photons contribute tag as the corresponding meson.
   if(okPhoton0 && okPhoton1)
   {
-    if(fDebug > 2) 
-      printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(ESD) - %s OVERLAPPED DECAY \n", meson->GetName());
+    AliDebug(2,Form("%s OVERLAPPED DECAY", meson->GetName()));
     
     if(!CheckTagBit(tag,kMCConversion) && conversion) SetTagBit(tag,kMCConversion) ;
     
@@ -990,9 +958,7 @@ void AliMCAnalysisUtils::CheckOverlapped2GammaDecay(const Int_t *labels, Int_t n
   
   if(labels[0] < 0 || labels[0] > mcparticles->GetEntriesFast() || nlabels <= 1)
   {
-    if(fDebug > 2) 
-      printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(AOD) - Exit : label[0] %d, n primaries %d, nlabels %d \n",
-             labels[0],mcparticles->GetEntriesFast(), nlabels);
+    AliDebug(2,Form("Exit : label[0] %d, n primaries %d, nlabels %d",labels[0],mcparticles->GetEntriesFast(), nlabels));
     return;
   }
   
@@ -1000,18 +966,16 @@ void AliMCAnalysisUtils::CheckOverlapped2GammaDecay(const Int_t *labels, Int_t n
   Int_t mesonPdg = meson->GetPdgCode();
   if(mesonPdg != 111 && mesonPdg != 221)
   {
-    printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(AOD) - Wrong pi0/eta PDG : %d \n",mesonPdg);
+    AliWarning(Form("Wrong pi0/eta PDG : %d",mesonPdg));
     return;
   }
   
-  if(fDebug > 2) printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(AOD) - pdg %d, label %d, ndaughters %d\n", mesonPdg, mesonIndex, meson->GetNDaughters());
-  
+  AliDebug(2,Form("pdg %d, label %d, ndaughters %d", mesonPdg, mesonIndex, meson->GetNDaughters()));
   
   //Get the daughters
   if(meson->GetNDaughters() != 2)
   {
-    if(fDebug > 2) 
-      printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(ESD) - Not overalapped. Number of daughters is %d, not 2 \n",meson->GetNDaughters());
+    AliDebug(2,Form("Not overalapped. Number of daughters is %d, not 2",meson->GetNDaughters()));
     return;
   }
   
@@ -1028,26 +992,23 @@ void AliMCAnalysisUtils::CheckOverlapped2GammaDecay(const Int_t *labels, Int_t n
   //Check if both daughters are photons
   if(photon0->GetPdgCode() != 22 && photon1->GetPdgCode()!=22)
   {
-    printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(AOD) - Not overlapped. PDG:  daughter 1 = %d, of daughter 2 = %d \n",photon0->GetPdgCode(),photon1->GetPdgCode());
+    AliWarning(Form("Not overlapped. PDG:  daughter 1 = %d, of daughter 2 = %d",photon0->GetPdgCode(),photon1->GetPdgCode()));
     return;
   }
   
-  if(fDebug > 2) 
-    printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(AOD) - Daughter labels : photon0 = %d, photon1 = %d \n",iPhoton0,iPhoton1);
+  AliDebug(2,Form("Daughter labels : photon0 = %d, photon1 = %d",iPhoton0,iPhoton1));
   
   //Check if both photons contribute to the cluster
   Bool_t okPhoton0 = kFALSE;
   Bool_t okPhoton1 = kFALSE;
   
-  if(fDebug > 3) 
-    printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(AOD) - Labels loop:\n");
+  AliDebug(3,"Labels loop:");
   
   Bool_t conversion = kFALSE;
   
   for(Int_t i = 0; i < nlabels; i++)
   {
-    if(fDebug > 3)
-      printf("\t label %d/%d: %d, ok? %d, %d\n", i, nlabels, labels[i], okPhoton0, okPhoton1);
+    AliDebug(3, Form("\t label %d/%d: %d, ok? %d, %d", i, nlabels, labels[i], okPhoton0, okPhoton1));
     
     if(labels[i]<0) continue;
     
@@ -1070,20 +1031,18 @@ void AliMCAnalysisUtils::CheckOverlapped2GammaDecay(const Int_t *labels, Int_t n
     
     if(index >= mcparticles->GetEntriesFast())
     {
-      printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(AOD) Particle index %d larger than size of list %d !!\n",index,mcparticles->GetEntriesFast());
+      AliWarning(Form("Particle index %d larger than size of list %d!!",index,mcparticles->GetEntriesFast()));
       continue;
     }
     
     AliAODMCParticle * daught = (AliAODMCParticle*) mcparticles->At(index);
     Int_t tmpindex = daught->GetMother();
-    if(fDebug > 3) 
-      printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(AOD) - Conversion? : mother %d\n",tmpindex);
+    AliDebug(3,Form("Conversion? : mother %d",tmpindex));
     
     while(tmpindex>=0){
       
       //MC particle of interest is the mother
-      if(fDebug > 3) 
-        printf("\t parent index %d\n",tmpindex);
+      AliDebug(3,Form("\t parent index %d",tmpindex));
       daught   = (AliAODMCParticle*) mcparticles->At(tmpindex);
       //printf("tmpindex %d\n",tmpindex);
       if      (iPhoton0 == tmpindex)
@@ -1103,22 +1062,18 @@ void AliMCAnalysisUtils::CheckOverlapped2GammaDecay(const Int_t *labels, Int_t n
       
     }//While to check if pi0/eta daughter was one of these contributors to the cluster
     
-    if(i == 0 && (!okPhoton0 && !okPhoton1) && fDebug>=0 )
-      printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(AOD) - Something happens, first label should be from a photon decay!\n");
+    //if(i == 0 && (!okPhoton0 && !okPhoton1)) AliDebug(1,"Something happens, first label should be from a photon decay!");
     
   }//loop on list of labels
   
   //If both photons contribute tag as the corresponding meson.
   if(okPhoton0 && okPhoton1)
   {
-    if(fDebug > 2) 
-      printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(AOD) - %s OVERLAPPED DECAY \n",
-             (TDatabasePDG::Instance()->GetParticle(mesonPdg))->GetName());
+    AliDebug(2,Form("%s OVERLAPPED DECAY",(TDatabasePDG::Instance()->GetParticle(mesonPdg))->GetName()));
     
     if(!CheckTagBit(tag,kMCConversion) && conversion)
     {
-      if(fDebug > 2)
-        printf("AliMCAnalysisUtils::CheckOverlapped2GammaDecay(AOD) Second decay photon produced a conversion \n");
+      AliDebug(2,"Second decay photon produced a conversion");
       SetTagBit(tag,kMCConversion) ;
     }
     
@@ -1338,8 +1293,7 @@ void    AliMCAnalysisUtils::CheckLostDecayPair(const TObjArray   * arrayCluster,
 
   
   SetTagBit(tag, kMCDecayPairLost);
-
-
+  
 }
 
 //_____________________________________________________________________
@@ -1361,12 +1315,12 @@ TList * AliMCAnalysisUtils::GetJets(const AliCaloTrackReader * reader)
     if(stack->GetNtrack() < 8) return fJetsList;
     TParticle * parton1 =  stack->Particle(6);
     TParticle * parton2 =  stack->Particle(7);
-    if(fDebug > 2){
-      printf("AliMCAnalysisUtils::GetJets() - parton 6 : %s, pt %2.2f,E %2.2f, phi %2.2f, eta %2.2f \n",
-             parton1->GetName(),parton1->Pt(),parton1->Energy(),parton1->Phi()*TMath::RadToDeg(),parton1->Eta());
-      printf("AliMCAnalysisUtils::GetJets() - parton 7 : %s, pt %2.2f,E %2.2f, phi %2.2f, eta %2.2f \n",
-             parton2->GetName(),parton2->Pt(),parton2->Energy(),parton2->Phi()*TMath::RadToDeg(),parton2->Eta());
-               }
+    
+    AliDebug(2,Form("Parton 6 : %s, pt %2.2f,E %2.2f, phi %2.2f, eta %2.2f",
+                    parton1->GetName(),parton1->Pt(),parton1->Energy(),parton1->Phi()*TMath::RadToDeg(),parton1->Eta()));
+    AliDebug(2,Form("Parton 7 : %s, pt %2.2f,E %2.2f, phi %2.2f, eta %2.2f",
+                    parton2->GetName(),parton2->Pt(),parton2->Energy(),parton2->Phi()*TMath::RadToDeg(),parton2->Eta()));
+    
     //                 //Trace the jet from the mother parton
     //                 Float_t pt  = 0;
     //                 Float_t pt1 = 0;
@@ -1411,8 +1365,7 @@ TList * AliMCAnalysisUtils::GetJets(const AliCaloTrackReader * reader)
       TParticle * jet =  0x0;
       AliGenPythiaEventHeader* pygeh= (AliGenPythiaEventHeader*) geh;
       nTriggerJets =  pygeh->NTriggerJets();
-      if(fDebug > 1)
-        printf("AliMCAnalysisUtils::GetJets() - PythiaEventHeader: Njets: %d\n",nTriggerJets);
+      AliDebug(2,Form("PythiaEventHeader: Njets: %d",nTriggerJets));
       
       for(Int_t i = 0; i< nTriggerJets; i++)
       {
@@ -1424,9 +1377,8 @@ TList * AliMCAnalysisUtils::GetJets(const AliCaloTrackReader * reader)
         if(phidiff1 > phidiff2) jet->SetFirstMother(7);
         else  jet->SetFirstMother(6);
         //jet->Print();
-        if(fDebug > 1)
-          printf("AliMCAnalysisUtils::GetJets() - PYTHIA Jet %d: mother %d, pt %2.2f,E %2.2f, phi %2.2f, eta %2.2f \n",
-                 i, jet->GetFirstMother(),jet->Pt(),jet->Energy(),jet->Phi()*TMath::RadToDeg(),jet->Eta());
+        AliDebug(1,Form("PYTHIA Jet %d: mother %d, pt %2.2f,E %2.2f, phi %2.2f, eta %2.2f",
+                        i, jet->GetFirstMother(),jet->Pt(),jet->Energy(),jet->Phi()*TMath::RadToDeg(),jet->Eta()));
         fJetsList->Add(jet);                   
       }
     }//Pythia triggered jets
@@ -1452,9 +1404,8 @@ TList * AliMCAnalysisUtils::GetJets(const AliCaloTrackReader * reader)
         //tmp = stack->Particle(tmp->GetFirstDaughter());
         //tmp->Print();
         //jet1->Print();
-        if(fDebug > 1)                 
-          printf("AliMCAnalysisUtils::GetJets() - HERWIG Jet 1: mother %d, status %d, pt %2.2f,E %2.2f, phi %2.2f, eta %2.2f \n",
-                 jet1->GetFirstMother(),jet1->GetStatusCode(),jet1->Pt(),jet1->Energy(),jet1->Phi()*TMath::RadToDeg(),jet1->Eta());
+        AliDebug(1,Form("HERWIG Jet 1: mother %d, status %d, pt %2.2f,E %2.2f, phi %2.2f, eta %2.2f",
+                        jet1->GetFirstMother(),jet1->GetStatusCode(),jet1->Pt(),jet1->Energy(),jet1->Phi()*TMath::RadToDeg(),jet1->Eta()));
       }//not photon
       
       //Check parton 2
@@ -1474,9 +1425,8 @@ TList * AliMCAnalysisUtils::GetJets(const AliCaloTrackReader * reader)
         jet2->SetFirstMother(7);
         fJetsList->Add(jet2);
         //jet2->Print();
-        if(fDebug > 1)
-          printf("AliMCAnalysisUtils::GetJets() - HERWIG Jet 2: mother %d, status %d, pt %2.2f,E %2.2f, phi %2.2f, eta %2.2f \n",
-                 jet2->GetFirstMother(),jet2->GetStatusCode(),jet2->Pt(),jet2->Energy(),jet2->Phi()*TMath::RadToDeg(),jet2->Eta());
+        AliDebug(2,Form("HERWIG Jet 2: mother %d, status %d, pt %2.2f,E %2.2f, phi %2.2f, eta %2.2f",
+                        jet2->GetFirstMother(),jet2->GetStatusCode(),jet2->Pt(),jet2->Energy(),jet2->Phi()*TMath::RadToDeg(),jet2->Eta()));
         //Int_t first =  tmp->GetFirstDaughter();
         //Int_t last  =  tmp->GetLastDaughter();
         //printf("jet 2:  first daughter %d, last daughter %d, pdg %d\n",first, last, tmp->GetPdgCode());
@@ -1507,8 +1457,7 @@ TLorentzVector AliMCAnalysisUtils::GetDaughter(Int_t idaugh, Int_t label,
   {
     if(!reader->GetStack())
     {
-      if (fDebug >=0)
-        printf("AliMCAnalysisUtils::GetDaughter() - Stack is not available, check analysis settings in configuration file, STOP!!\n");
+      AliWarning("Stack is not available, check analysis settings in configuration file!!");
       
       ok=kFALSE;
       return fDaughMom;
@@ -1542,8 +1491,7 @@ TLorentzVector AliMCAnalysisUtils::GetDaughter(Int_t idaugh, Int_t label,
     TClonesArray* mcparticles = reader->GetAODMCParticles();
     if(!mcparticles)
     {
-      if(fDebug >= 0)
-        printf("AliMCAnalysisUtils::GetDaughter() - AODMCParticles is not available, check analysis settings in configuration file!!\n");
+      AliWarning("AODMCParticles is not available, check analysis settings in configuration file!!");
       
       ok=kFALSE;
       return fDaughMom;
@@ -1609,8 +1557,7 @@ TLorentzVector AliMCAnalysisUtils::GetMother(Int_t label, const AliCaloTrackRead
   {
     if(!reader->GetStack()) 
     {
-      if (fDebug >=0) 
-        printf("AliMCAnalysisUtils::GetMother() - Stack is not available, check analysis settings in configuration file, STOP!!\n");
+      AliWarning("Stack is not available, check analysis settings in configuration file, STOP!!");
      
       ok=kFALSE;
       return fMotherMom;
@@ -1634,8 +1581,7 @@ TLorentzVector AliMCAnalysisUtils::GetMother(Int_t label, const AliCaloTrackRead
     TClonesArray* mcparticles = reader->GetAODMCParticles();
     if(!mcparticles) 
     {
-      if(fDebug >= 0)
-        printf("AliMCAnalysisUtils::GetMother() - AODMCParticles is not available, check analysis settings in configuration file!!\n");
+      AliWarning("AODMCParticles is not available, check analysis settings in configuration file!!");
       
       ok=kFALSE;
       return fMotherMom;
@@ -1675,8 +1621,7 @@ TLorentzVector AliMCAnalysisUtils::GetMotherWithPDG(Int_t label, Int_t pdg,
   {
     if(!reader->GetStack())
     {
-      if (fDebug >=0) 
-        printf("AliMCAnalysisUtils::GetMotherWithPDG() - Stack is not available, check analysis settings in configuration file, STOP!!\n");
+     AliWarning("Stack is not available, check analysis settings in configuration file!!");
       
       ok = kFALSE;
       return fGMotherMom;
@@ -1704,7 +1649,7 @@ TLorentzVector AliMCAnalysisUtils::GetMotherWithPDG(Int_t label, Int_t pdg,
         
       }
       
-      if(grandmomPDG!=pdg) printf("AliMCAnalysisUtils::GetMotherWithPDG(ESD) - mother with PDG %d, not found! \n",pdg);
+      if(grandmomPDG!=pdg) AliInfo(Form("Mother with PDG %d, not found! \n",pdg));
     }
   }
   else if(reader->ReadAODMCParticles())
@@ -1712,8 +1657,7 @@ TLorentzVector AliMCAnalysisUtils::GetMotherWithPDG(Int_t label, Int_t pdg,
     TClonesArray* mcparticles = reader->GetAODMCParticles();
     if(!mcparticles) 
     {
-      if(fDebug >= 0)
-        printf("AliMCAnalysisUtils::GetMotherWithPDG() - AODMCParticles is not available, check analysis settings in configuration file!!\n");
+      AliWarning("AODMCParticles is not available, check analysis settings in configuration file!!");
       
       ok=kFALSE;
       return fGMotherMom;
@@ -1743,7 +1687,7 @@ TLorentzVector AliMCAnalysisUtils::GetMotherWithPDG(Int_t label, Int_t pdg,
         
       }
       
-      if(grandmomPDG!=pdg) printf("AliMCAnalysisUtils::GetMotherWithPDG(AOD) - mother with PDG %d, NOT found! \n",pdg);
+      if(grandmomPDG!=pdg) AliInfo(Form("Mother with PDG %d, NOT found!",pdg));
             
     }
   }
@@ -1766,8 +1710,7 @@ TLorentzVector AliMCAnalysisUtils::GetGrandMother(Int_t label, const AliCaloTrac
   {
     if(!reader->GetStack())
     {
-      if (fDebug >=0)
-        printf("AliMCAnalysisUtils::GetGrandMother() - Stack is not available, check analysis settings in configuration file, STOP!!\n");
+      AliWarning("Stack is not available, check analysis settings in configuration file, STOP!!");
       
       ok = kFALSE;
       return fGMotherMom;
@@ -1798,8 +1741,7 @@ TLorentzVector AliMCAnalysisUtils::GetGrandMother(Int_t label, const AliCaloTrac
     TClonesArray* mcparticles = reader->GetAODMCParticles();
     if(!mcparticles)
     {
-      if(fDebug >= 0)
-        printf("AliMCAnalysisUtils::GetGrandMother() - AODMCParticles is not available, check analysis settings in configuration file!!\n");
+      AliWarning("AODMCParticles is not available, check analysis settings in configuration file!!");
       
       ok=kFALSE;
       return fGMotherMom;
@@ -1842,8 +1784,7 @@ void AliMCAnalysisUtils::GetMCDecayAsymmetryAngleForPDG(Int_t label, Int_t pdg,
   {
     if(!reader->GetStack())
     {
-      if (fDebug >=0) 
-        printf("AliMCAnalysisUtils::GetMCDecayAsymmetryForPDG() - Stack is not available, check analysis settings in configuration file, STOP!!\n");
+      AliWarning("Stack is not available, check analysis settings in configuration file, STOP!!");
       
       ok = kFALSE;
     }
@@ -1880,7 +1821,7 @@ void AliMCAnalysisUtils::GetMCDecayAsymmetryAngleForPDG(Int_t label, Int_t pdg,
       else 
       {
         ok=kFALSE;
-        printf("AliMCAnalysisUtils::GetMCDecayAsymmetryForPDG(ESD) - mother with PDG %d, not found! \n",pdg);
+        AliInfo(Form("Mother with PDG %d, not found!",pdg));
       }
       
       } // good label
@@ -1890,8 +1831,7 @@ void AliMCAnalysisUtils::GetMCDecayAsymmetryAngleForPDG(Int_t label, Int_t pdg,
     TClonesArray* mcparticles = reader->GetAODMCParticles();
     if(!mcparticles) 
     {
-      if(fDebug >= 0)
-        printf("AliMCAnalysisUtils::GetMCDecayAsymmetryForPDG() - AODMCParticles is not available, check analysis settings in configuration file!!\n");
+      AliWarning("AODMCParticles is not available, check analysis settings in configuration file!!");
       
       ok=kFALSE;
       return;
@@ -1931,7 +1871,7 @@ void AliMCAnalysisUtils::GetMCDecayAsymmetryAngleForPDG(Int_t label, Int_t pdg,
       else 
       {
         ok=kFALSE;
-        printf("AliMCAnalysisUtils::GetMCDecayAsymmetryForPDG(AOD) - mother with PDG %d, not found! \n",pdg);
+        AliInfo(Form("Mother with PDG %d, not found! \n",pdg));
       }
       
     } // good label
@@ -1950,8 +1890,7 @@ Int_t AliMCAnalysisUtils::GetNDaughters(Int_t label, const AliCaloTrackReader* r
   {
     if(!reader->GetStack())
     {
-      if (fDebug >=0)
-        printf("AliMCAnalysisUtils::GetNDaughters() - Stack is not available, check analysis settings in configuration file, STOP!!\n");
+      AliWarning("Stack is not available, check analysis settings in configuration file, STOP!!");
       
       ok=kFALSE;
       return -1;
@@ -1973,8 +1912,7 @@ Int_t AliMCAnalysisUtils::GetNDaughters(Int_t label, const AliCaloTrackReader* r
     TClonesArray* mcparticles = reader->GetAODMCParticles();
     if(!mcparticles)
     {
-      if(fDebug >= 0)
-        printf("AliMCAnalysisUtils::GetNDaughters() - AODMCParticles is not available, check analysis settings in configuration file!!\n");
+      AliWarning("AODMCParticles is not available, check analysis settings in configuration file!!");
       
       ok=kFALSE;
       return -1;