From: gconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Fri, 25 May 2012 03:34:24 +0000 (+0000)
Subject: Add centrality control histogram
X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=e19ca838158ee8a42aae713cc3d9e578038006fa;p=u%2Fmrichter%2FAliRoot.git

Add centrality control histogram
---

diff --git a/PWG/CaloTrackCorrBase/AliAnalysisTaskCounter.cxx b/PWG/CaloTrackCorrBase/AliAnalysisTaskCounter.cxx
index 75564aab0ca..27f7aef8462 100644
--- a/PWG/CaloTrackCorrBase/AliAnalysisTaskCounter.cxx
+++ b/PWG/CaloTrackCorrBase/AliAnalysisTaskCounter.cxx
@@ -63,7 +63,8 @@ AliAnalysisTaskCounter::AliAnalysisTaskCounter(const char *name)
   fTriggerAnalysis (new AliTriggerAnalysis),
   fhNEvents(0),
   fhXVertex(0),    fhYVertex(0),    fhZVertex(0),
-  fhXGoodVertex(0),fhYGoodVertex(0),fhZGoodVertex(0)
+  fhXGoodVertex(0),fhYGoodVertex(0),fhZGoodVertex(0),
+  fhCentrality(0)
 {
   //ctor
   DefineOutput(1, TList::Class());
@@ -81,7 +82,8 @@ AliAnalysisTaskCounter::AliAnalysisTaskCounter()
     fTriggerAnalysis (new AliTriggerAnalysis),
     fhNEvents(0),
     fhXVertex(0),    fhYVertex(0),    fhZVertex(0),
-    fhXGoodVertex(0),fhYGoodVertex(0),fhZGoodVertex(0)
+    fhXGoodVertex(0),fhYGoodVertex(0),fhZGoodVertex(0),
+    fhCentrality(0)
 {
   // ctor
   DefineOutput(1, TList::Class());
@@ -133,6 +135,9 @@ void AliAnalysisTaskCounter::UserCreateOutputObjects()
   fhYGoodVertex->SetXTitle("v_{y} (cm)");
   fOutputContainer->Add(fhYGoodVertex);
   
+  fhCentrality   = new TH1F("hCentrality","Number of events in centrality bin, |vz|<10 cm, method <V0M> ",100,0.,100.) ;
+  fhCentrality->SetXTitle("Centrality bin");
+  fOutputContainer->Add(fhCentrality) ;  
   
   fhNEvents = new TH1I("hNEvents", "Number of analyzed events", 21, 0, 21) ;
   fhNEvents->SetXTitle("Selection");
@@ -337,6 +342,12 @@ void AliAnalysisTaskCounter::UserExec(Option_t *)
     if(bSelectVZ && bSelectTrack && bV0AND)    
                      fhNEvents->Fill(14.5); 
     if(!bPileup)     fhNEvents->Fill(15.5); 
+
+    if(TMath::Abs(v[2]) < 10 && InputEvent()->GetCentrality()) 
+    {
+      fhCentrality->Fill(InputEvent()->GetCentrality()->GetCentralityPercentile("V0M"));
+    }
+  
   }
 
   //printf("AliAnalysisTaskCounter::UserExec() : z vertex %d, good vertex %d, v0and %d, pile up %d, track mult %d\n ", bSelectVZ, bGoodV, bV0AND, bPileup, trackMult);
diff --git a/PWG/CaloTrackCorrBase/AliAnalysisTaskCounter.h b/PWG/CaloTrackCorrBase/AliAnalysisTaskCounter.h
index 5d735719c4b..127fa6c544f 100644
--- a/PWG/CaloTrackCorrBase/AliAnalysisTaskCounter.h
+++ b/PWG/CaloTrackCorrBase/AliAnalysisTaskCounter.h
@@ -57,11 +57,12 @@ class AliAnalysisTaskCounter : public AliAnalysisTaskSE {
   TH1F *  fhXGoodVertex;  //! X Vertex good distribution
   TH1F *  fhYGoodVertex;  //! Y Vertex good distribution
   TH1F *  fhZGoodVertex;  //! Z Vertex good distribution  
-
+  TH1F *  fhCentrality;   //! centrality
+  
   AliAnalysisTaskCounter(           const AliAnalysisTaskCounter&); // not implemented  
   AliAnalysisTaskCounter& operator=(const AliAnalysisTaskCounter&); // not implemented
   
-  ClassDef(AliAnalysisTaskCounter, 1);
+  ClassDef(AliAnalysisTaskCounter, 2);
 
 };