samll fixes (Marta)
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 7 Jun 2011 09:09:47 +0000 (09:09 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 7 Jun 2011 09:09:47 +0000 (09:09 +0000)
PWG4/JetTasks/AliPWG4HighPtTrackQA.cxx
PWG4/JetTasks/AliPWG4HighPtTrackQA.h

index ba26738..da114d0 100644 (file)
@@ -697,6 +697,10 @@ Bool_t AliPWG4HighPtTrackQA::SelectEvent() {
 
 //________________________________________________________________________
 Int_t AliPWG4HighPtTrackQA::CalculateCentrality(AliVEvent *ev){
+  //
+  // Get centrality from ESD or AOD
+  //
+
   if(fDataType==kESD)
     return CalculateCentrality(dynamic_cast<AliESDEvent*>(ev));
   else if(fDataType==kAOD)
@@ -707,35 +711,48 @@ Int_t AliPWG4HighPtTrackQA::CalculateCentrality(AliVEvent *ev){
 
 //________________________________________________________________________
 Int_t AliPWG4HighPtTrackQA::CalculateCentrality(AliESDEvent *esd){
+  //
+  // Get centrality from ESD
+  //
 
-
-  Float_t cent = 999;
+  Float_t cent = -1;
 
   if(esd){
     if(esd->GetCentrality()){
       cent = esd->GetCentrality()->GetCentralityPercentile("V0M");
-      if(fDebug>3) cout << "centrality: " << cent << endl;
+      if(fDebug>3) printf("centrality: %f\n",cent);
     }
   }
 
-  if(cent>80)return 4;
-  if(cent>50)return 3;
-  if(cent>30)return 2;
-  if(cent>10)return 1;
-  return 0;
+  return GetCentralityClass(cent);
 
 }
 
 //________________________________________________________________________
 Int_t AliPWG4HighPtTrackQA::CalculateCentrality(AliAODEvent *aod){
+  //
+  // Get centrality from AOD
+  //
 
-  if(!aod)return 4;
+  if(!aod) return 5;
   Float_t cent = aod->GetHeader()->GetCentrality();
-  cout << "centrality: " << cent << endl;
-  if(cent>80)return 4;
-  if(cent>50)return 3;
-  if(cent>30)return 2;
-  if(cent>10)return 1;
+  if(fDebug>3) printf("centrality: %f\n",cent);
+
+  return GetCentralityClass(cent);
+
+}
+
+//________________________________________________________________________
+Int_t AliPWG4HighPtTrackQA::GetCentralityClass(Float_t cent) {
+  //
+  // Get centrality class
+  //
+
+  if(cent<0)  return 5; // OB - cent sometimes negative
+  if(cent>80) return 4;
+  if(cent>50) return 3;
+  if(cent>30) return 2;
+  if(cent>10) return 1;
   return 0;
 
 }
index b490d61..4dced3c 100644 (file)
@@ -62,6 +62,7 @@ class AliPWG4HighPtTrackQA: public AliAnalysisTaskSE {
   Int_t CalculateCentrality(AliVEvent *ev);
   Int_t CalculateCentrality(AliESDEvent *esd);
   Int_t CalculateCentrality(AliAODEvent *aod);
+  Int_t GetCentralityClass(Float_t cent=-1.);
   void DoAnalysisESD();
   void DoAnalysisAOD();
   void FillHistograms();
@@ -107,7 +108,7 @@ class AliPWG4HighPtTrackQA: public AliAnalysisTaskSE {
   Float_t fPtMax;                 // Maximum pT for histograms
   Float_t fPtBinEdges[3][2];      // 3 regions total with different binning for pT axis of histos
 
-  Bool_t   fIsPbPb;               //  kTRUE if PbPb
+  Bool_t   fIsPbPb;               // kTRUE if PbPb
   Int_t fCentClass;               // Select only events from predefined centrality class
 
   /*