add correation with V0mult to Rho task

diff --git a/PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.cxx b/PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.cxx
index 891e57378def6fd9e3d7128e5d7a13802c7aeefb..e18eac75b8896dd8e91854d1e7848f8c7c606855 100644 (file)
--- a/PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.cxx
+++ b/PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.cxx
@@ -8,6 +8,7 @@
 #include <TF1.h>
 #include <TH1F.h>
 #include <TH2F.h>
+#include <TH3F.h>
 #include <TClonesArray.h>
 
 #include "AliLog.h"
@@ -15,6 +16,7 @@
 #include "AliEmcalJet.h"
 #include "AliParticleContainer.h"
 #include "AliClusterContainer.h"
+#include "AliVVZERO.h"
 
 #include "AliAnalysisTaskRhoBase.h"
 
@@ -46,8 +48,8 @@ AliAnalysisTaskRhoBase::AliAnalysisTaskRhoBase() :
   fHistRhoScaledvsCent(0),
   fHistDeltaRhovsCent(0),
   fHistDeltaRhoScalevsCent(0),
-  fHistRhovsNtrack(0),
-  fHistRhoScaledvsNtrack(0),
+  fHistRhovsNtrackvsV0Mult(0),
+  fHistRhoScaledvsNtrackvsV0Mult(0),
   fHistDeltaRhovsNtrack(0),
   fHistDeltaRhoScalevsNtrack(0),
   fHistRhovsNcluster(0),
@@ -90,8 +92,8 @@ AliAnalysisTaskRhoBase::AliAnalysisTaskRhoBase(const char *name, Bool_t histo) :
   fHistRhoScaledvsCent(0),
   fHistDeltaRhovsCent(0),
   fHistDeltaRhoScalevsCent(0),
-  fHistRhovsNtrack(0),
-  fHistRhoScaledvsNtrack(0),
+  fHistRhovsNtrackvsV0Mult(0),
+  fHistRhoScaledvsNtrackvsV0Mult(0),
   fHistDeltaRhovsNtrack(0),
   fHistDeltaRhoScalevsNtrack(0),
   fHistRhovsNcluster(0),
@@ -125,10 +127,11 @@ void AliAnalysisTaskRhoBase::UserCreateOutputObjects()
   fOutput->Add(fHistRhovsCent);
 
   if (fParticleCollArray.GetEntriesFast()>0) {
-    fHistRhovsNtrack = new TH2F("fHistRhovsNtrack", "fHistRhovsNtrack", 150, 0, 6000, fNbins, fMinBinPt, fMaxBinPt*2);
-    fHistRhovsNtrack->GetXaxis()->SetTitle("No. of tracks");
-    fHistRhovsNtrack->GetYaxis()->SetTitle("#rho (GeV/c * rad^{-1})");
-    fOutput->Add(fHistRhovsNtrack);
+    fHistRhovsNtrackvsV0Mult = new TH3F("fHistRhovsNtrackvsV0Mult", "fHistRhovsNtrackvsV0Mult", 150, 0, 6000, fNbins, fMinBinPt, fMaxBinPt*2,100,0.,20000.);
+    fHistRhovsNtrackvsV0Mult->GetXaxis()->SetTitle("No. of tracks");
+    fHistRhovsNtrackvsV0Mult->GetYaxis()->SetTitle("#rho (GeV/c * rad^{-1})");
+    fHistRhovsNtrackvsV0Mult->GetZaxis()->SetTitle("V0 mult");
+    fOutput->Add(fHistRhovsNtrackvsV0Mult);
   }
 
   if (fClusterCollArray.GetEntriesFast()>0) {
@@ -223,10 +226,11 @@ void AliAnalysisTaskRhoBase::UserCreateOutputObjects()
     fOutput->Add(fHistRhoScaledvsCent);
 
     if (fParticleCollArray.GetEntriesFast()>0) {
-      fHistRhoScaledvsNtrack = new TH2F("fHistRhoScaledvsNtrack", "fHistRhoScaledvsNtrack", 150, 0, 6000, fNbins, fMinBinPt, fMaxBinPt*2);
-      fHistRhoScaledvsNtrack->GetXaxis()->SetTitle("No. of tracks");
-      fHistRhoScaledvsNtrack->GetYaxis()->SetTitle("#rho_{scaled} (GeV/c * rad^{-1})");
-      fOutput->Add(fHistRhoScaledvsNtrack);
+      fHistRhoScaledvsNtrackvsV0Mult = new TH3F("fHistRhoScaledvsNtrackvsV0Mult", "fHistRhoScaledvsNtrackvsV0Mult", 150, 0, 6000, fNbins, fMinBinPt, fMaxBinPt*2,100,0.,20000.);
+      fHistRhoScaledvsNtrackvsV0Mult->GetXaxis()->SetTitle("No. of tracks");
+      fHistRhoScaledvsNtrackvsV0Mult->GetYaxis()->SetTitle("#rho (GeV/c * rad^{-1})");
+      fHistRhoScaledvsNtrackvsV0Mult->GetZaxis()->SetTitle("V0 mult");
+      fOutput->Add(fHistRhoScaledvsNtrackvsV0Mult);
     }
 
     if (fClusterCollArray.GetEntriesFast()>0) {
@@ -276,6 +280,10 @@ Bool_t AliAnalysisTaskRhoBase::FillHistograms()
   Int_t Ntracks   = 0;
   Int_t Nclusters = 0;
 
+  AliVVZERO* vV0 = InputEvent()->GetVZEROData();
+  Float_t multV0A = vV0->GetMTotV0A();
+  Float_t multV0C = vV0->GetMTotV0C();
+
   if (GetParticleContainer(0))
     Ntracks = GetParticleContainer(0)->GetNAcceptedParticles();
   if (GetClusterContainer(0))
@@ -340,7 +348,7 @@ Bool_t AliAnalysisTaskRhoBase::FillHistograms()
   fHistRhovsCent->Fill(fCent, fOutRho->GetVal());
 
   if (fTracks)
-    fHistRhovsNtrack->Fill(Ntracks, fOutRho->GetVal());
+    fHistRhovsNtrackvsV0Mult->Fill(Ntracks, fOutRho->GetVal(),multV0A+multV0C);
   if (fCaloClusters)
     fHistRhovsNcluster->Fill(Nclusters, fOutRho->GetVal());
   if (fCompareRho) {
@@ -352,7 +360,7 @@ Bool_t AliAnalysisTaskRhoBase::FillHistograms()
   if (fOutRhoScaled) {
     fHistRhoScaledvsCent->Fill(fCent, fOutRhoScaled->GetVal());
     if (fTracks)
-      fHistRhoScaledvsNtrack->Fill(Ntracks, fOutRhoScaled->GetVal());
+      fHistRhoScaledvsNtrackvsV0Mult->Fill(Ntracks, fOutRhoScaled->GetVal(),multV0A+multV0C);
     if (fCaloClusters)
       fHistRhoScaledvsNcluster->Fill(Nclusters,  fOutRhoScaled->GetVal());
     if (fCompareRhoScaled) {

diff --git a/PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.h b/PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.h
index 4b32efea224d618cb3eedb5034d6afded685d1dd..1daf5a290a6507374311cd736fa43fb3d48aa19f 100644 (file)
--- a/PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.h
+++ b/PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.h
@@ -7,6 +7,7 @@ class TString;
 class TF1;
 class TH1F;
 class TH2F;
+class TH3F;
 class AliRhoParameter;
 
 #include "AliAnalysisTaskEmcalJet.h"
@@ -68,8 +69,8 @@ class AliAnalysisTaskRhoBase : public AliAnalysisTaskEmcalJet {
   TH2F                  *fHistDeltaRhovsCent;            //!delta rho vs. centrality
   TH2F                  *fHistDeltaRhoScalevsCent;       //!delta rhoscaled vs. centrality
 
-  TH2F                  *fHistRhovsNtrack;               //!rho vs. no. of tracks
-  TH2F                  *fHistRhoScaledvsNtrack;         //!rhoscaled vs. no. of tracks
+  TH3F                  *fHistRhovsNtrackvsV0Mult;       //!rho vs. no. of tracks vs V0mult
+  TH3F                  *fHistRhoScaledvsNtrackvsV0Mult; //!rhoscaled vs. no. of tracks vs V0mult
   TH2F                  *fHistDeltaRhovsNtrack;          //!delta rho vs. no. of tracks
   TH2F                  *fHistDeltaRhoScalevsNtrack;     //!delta rho scaled vs. no. of tracks
  
@@ -79,6 +80,6 @@ class AliAnalysisTaskRhoBase : public AliAnalysisTaskEmcalJet {
   AliAnalysisTaskRhoBase(const AliAnalysisTaskRhoBase&);             // not implemented
   AliAnalysisTaskRhoBase& operator=(const AliAnalysisTaskRhoBase&);  // not implemented
   
-  ClassDef(AliAnalysisTaskRhoBase, 10); // Rho base task
+  ClassDef(AliAnalysisTaskRhoBase, 11); // Rho base task
 };
 #endif