Add rho vs cent histogram
authorsaiola <salvatore.aiola@cern.ch>
Fri, 17 Oct 2014 21:33:27 +0000 (17:33 -0400)
committersaiola <salvatore.aiola@cern.ch>
Fri, 17 Oct 2014 21:33:27 +0000 (17:33 -0400)
PWGJE/EMCALJetTasks/AliAnalysisTaskDeltaPt.cxx
PWGJE/EMCALJetTasks/AliAnalysisTaskDeltaPt.h

index dbbc190..3a59483 100644 (file)
@@ -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
   // _________________________________
index 33502b1..5172d15 100644 (file)
@@ -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