Add rho vs cent histogram

diff --git a/PWGJE/EMCALJetTasks/AliAnalysisTaskDeltaPt.cxx b/PWGJE/EMCALJetTasks/AliAnalysisTaskDeltaPt.cxx
index dbbc190480a8bf0175e1453b757861152a4b87f2..3a5948324b10f0ccbc51ee96725e5108b8a72664 100644 (file)
--- a/PWGJE/EMCALJetTasks/AliAnalysisTaskDeltaPt.cxx
+++ b/PWGJE/EMCALJetTasks/AliAnalysisTaskDeltaPt.cxx
@@ -45,6 +45,7 @@ AliAnalysisTaskDeltaPt::AliAnalysisTaskDeltaPt() :
   fEmbCaloClustersCont(0),
   fRandTracksCont(0),
   fRandCaloClustersCont(0),
+  fHistRhovsCent(0),
   fHistRCPhiEta(0), 
   fHistRCPt(0),
   fHistRCPtExLJ(0),
@@ -114,6 +115,7 @@ AliAnalysisTaskDeltaPt::AliAnalysisTaskDeltaPt(const char *name) :
   fEmbCaloClustersCont(0),
   fRandTracksCont(0),
   fRandCaloClustersCont(0),
+  fHistRhovsCent(0),
   fHistRCPhiEta(0), 
   fHistRCPt(0),
   fHistRCPtExLJ(0),
@@ -219,6 +221,11 @@ void AliAnalysisTaskDeltaPt::UserCreateOutputObjects()
 
   AllocateHistogramArrays();
 
+  fHistRhovsCent = new TH2F("fHistRhovsCent", "fHistRhovsCent", 101, -1,  100, fNbins, 0, fMaxBinPt*2);
+  fHistRhovsCent->GetXaxis()->SetTitle("Centrality (%)");
+  fHistRhovsCent->GetYaxis()->SetTitle("#rho (GeV/c * rad^{-1})");
+  fOutput->Add(fHistRhovsCent);
+
   fJetsCont = GetJetContainer("Jets");
   fTracksCont = GetParticleContainer("Tracks");
   fCaloClustersCont = GetClusterContainer("CaloClusters");
@@ -432,6 +439,8 @@ Bool_t AliAnalysisTaskDeltaPt::FillHistograms()
 {
   // Fill histograms.
 
+  fHistRhovsCent->Fill(fCent, fRhoVal);
+
   // ************
   // Random cones
   // _________________________________

diff --git a/PWGJE/EMCALJetTasks/AliAnalysisTaskDeltaPt.h b/PWGJE/EMCALJetTasks/AliAnalysisTaskDeltaPt.h
index 33502b197307b3e908c662f8b3f7c7687cbbd1bb..5172d15365f0df3e6c8270e3d57cbbee6e936980 100644 (file)
--- a/PWGJE/EMCALJetTasks/AliAnalysisTaskDeltaPt.h
+++ b/PWGJE/EMCALJetTasks/AliAnalysisTaskDeltaPt.h
@@ -63,6 +63,9 @@ class AliAnalysisTaskDeltaPt : public AliAnalysisTaskEmcalJet {
   AliParticleContainer       *fRandTracksCont;             //!Randomized tracks
   AliClusterContainer        *fRandCaloClustersCont;       //!Randomized clusters
 
+  // General
+  TH2                        *fHistRhovsCent;              //!Rho vs. centrality
+
   // Random cones
   TH2                        *fHistRCPhiEta;               //!Phi-Eta distribution of random cones
   TH1                       **fHistRCPt;                   //!Random cone pt