From 68cd26247ea27f25f1053a434e5c6a4008d58565 Mon Sep 17 00:00:00 2001
From: gconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Tue, 19 Jun 2012 11:59:49 +0000
Subject: [PATCH] Add histograms to check the event centrality and event plane
 of the leading particle

---
 .../AliAnaParticleHadronCorrelation.cxx       | 39 ++++++++++++++++---
 .../AliAnaParticleHadronCorrelation.h         |  6 ++-
 2 files changed, 38 insertions(+), 7 deletions(-)

diff --git a/PWGGA/CaloTrackCorrelations/AliAnaParticleHadronCorrelation.cxx b/PWGGA/CaloTrackCorrelations/AliAnaParticleHadronCorrelation.cxx
index 085f2262c10..41e862201a6 100755
--- a/PWGGA/CaloTrackCorrelations/AliAnaParticleHadronCorrelation.cxx
+++ b/PWGGA/CaloTrackCorrelations/AliAnaParticleHadronCorrelation.cxx
@@ -73,7 +73,9 @@ ClassImp(AliAnaParticleHadronCorrelation)
     fListMixEvents(),               fUseMixStoredInReader(0),
     //Histograms
     fhPtLeading(0),                 fhPhiLeading(0),       
-    fhEtaLeading(0),                fhDeltaPhiDeltaEtaCharged(0),
+    fhEtaLeading(0),                
+    fhPtLeadingCentrality(0),       fhPtLeadingEventPlane(0), 
+    fhLeadingEventPlaneCentrality(0),fhDeltaPhiDeltaEtaCharged(0),
     fhPhiCharged(0),                fhEtaCharged(0), 
     fhDeltaPhiCharged(0),           fhDeltaEtaCharged(0), 
     fhDeltaPhiChargedPt(0),         fhDeltaPhiUeChargedPt(0), 
@@ -127,7 +129,7 @@ ClassImp(AliAnaParticleHadronCorrelation)
     fhNtracksMB(0),               
     fhMixDeltaPhiCharged(0),        fhMixDeltaPhiDeltaEtaCharged(0),
     fhMixDeltaPhiChargedAssocPtBin(),fhMixDeltaPhiDeltaEtaChargedAssocPtBin(0),
-    fhEventBin(0),                   fhEventMixBin(0)
+    fhEventBin(0),                  fhEventMixBin(0)
 { 
   //Default Ctor
     
@@ -696,6 +698,21 @@ TList *  AliAnaParticleHadronCorrelation::GetCreateOutputObjects()
   outputContainer->Add(fhPhiLeading);
   outputContainer->Add(fhEtaLeading);
   
+  fhPtLeadingCentrality   = new TH2F("hPtLeadingCentrality","Leading particle p_{T} vs centrality",nptbins,ptmin,ptmax,100,0.,100) ;
+  fhPtLeadingCentrality->SetXTitle("p_{T}^{trig} (GeV/c)");
+  fhPtLeadingCentrality->SetYTitle("Centrality (%)");
+  outputContainer->Add(fhPtLeadingCentrality) ;  
+  
+  fhPtLeadingEventPlane  = new TH2F("hPtLeadingEventPlane","Leading particle p_{T} vs event plane angle",nptbins,ptmin,ptmax, 100,0.,TMath::Pi()) ;
+  fhPtLeadingEventPlane->SetXTitle("p_{T}^{trig} (GeV/c)");
+  fhPtLeadingEventPlane->SetXTitle("EP angle (rad)");
+  outputContainer->Add(fhPtLeadingEventPlane) ;
+  
+  fhLeadingEventPlaneCentrality  = new TH2F("hLeadingEventPlane","Leading particle centrality vs event plane angle",100,0.,100,100,0.,TMath::Pi()) ;
+  fhLeadingEventPlaneCentrality->SetXTitle("Centrality (%)");
+  fhLeadingEventPlaneCentrality->SetYTitle("EP angle (rad)");
+  outputContainer->Add(fhLeadingEventPlaneCentrality) ;
+  
   //Correlation with charged hadrons
   if(GetReader()->IsCTSSwitchedOn()) 
   {
@@ -1781,13 +1798,23 @@ void  AliAnaParticleHadronCorrelation::MakeAnalysisFillHistograms()
     // no problem was found, like not absolute leading, or bad vertex in mixing.
     if(okcharged && okneutral)
     {
-      fhPtLeading->Fill(particle->Pt());
+      Float_t pt = particle->Pt();
+      fhPtLeading->Fill(pt);
+      
       Float_t phi = particle->Phi();
       if(phi<0)phi+=TMath::TwoPi();
-      fhPhiLeading->Fill(particle->Pt(), phi);
-      fhEtaLeading->Fill(particle->Pt(), particle->Eta());
+      fhPhiLeading->Fill(pt, phi);
+      
+      fhEtaLeading->Fill(pt, particle->Eta());
       //printf("AliAnaParticleHadronCorrelation::MakeAnalysisFillHistograms() - Leading particle : pt %f, eta %f, phi %f\n",particle->Pt(),particle->Eta(),phi);
-
+      
+      Float_t cen = GetEventCentrality();
+      Float_t ep  = GetEventPlaneAngle();
+      
+      fhPtLeadingCentrality        ->Fill(pt,cen);
+      fhPtLeadingEventPlane        ->Fill(pt,ep);
+      fhLeadingEventPlaneCentrality->Fill(cen,ep);
+      
     }//ok charged && neutral
   }//Aod branch loop
   
diff --git a/PWGGA/CaloTrackCorrelations/AliAnaParticleHadronCorrelation.h b/PWGGA/CaloTrackCorrelations/AliAnaParticleHadronCorrelation.h
index 2572b2a6aeb..89210d1f484 100755
--- a/PWGGA/CaloTrackCorrelations/AliAnaParticleHadronCorrelation.h
+++ b/PWGGA/CaloTrackCorrelations/AliAnaParticleHadronCorrelation.h
@@ -200,6 +200,10 @@ class AliAnaParticleHadronCorrelation : public AliAnaCaloTrackCorrBaseClass {
   TH2F *       fhPhiLeading;                   //! phi distribution vs pT of leading particles
   TH2F *       fhEtaLeading;                   //! eta distribution vs pT of leading particles
   
+  TH2F *       fhPtLeadingCentrality;          //! pT distribution of leading particles vs centrality
+  TH2F *       fhPtLeadingEventPlane;          //! pT distribution of leading particles vs centrality
+  TH2F *       fhLeadingEventPlaneCentrality;  //! event plane vs centrality for leading particles
+  
   //trigger-charged histograms
   TH2F *       fhDeltaPhiDeltaEtaCharged ;     //! differences of eta and phi between trigger and charged hadrons
   TH2F *       fhPhiCharged  ;                 //! Phi distribution of charged particles
@@ -328,7 +332,7 @@ class AliAnaParticleHadronCorrelation : public AliAnaCaloTrackCorrBaseClass {
   AliAnaParticleHadronCorrelation(              const AliAnaParticleHadronCorrelation & ph) ; // cpy ctor
   AliAnaParticleHadronCorrelation & operator = (const AliAnaParticleHadronCorrelation & ph) ; // cpy assignment
 	
-  ClassDef(AliAnaParticleHadronCorrelation,15)
+  ClassDef(AliAnaParticleHadronCorrelation,16)
 } ;
  
 
-- 
2.43.0