From c3a893125772effda93a60720ecef4c0cecc95b7 Mon Sep 17 00:00:00 2001 From: mverweij Date: Mon, 10 Nov 2014 10:42:08 +0100 Subject: [PATCH] add correation with V0mult to Rho task --- .../EMCALJetTasks/AliAnalysisTaskRhoBase.cxx | 36 +++++++++++-------- PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.h | 7 ++-- 2 files changed, 26 insertions(+), 17 deletions(-) diff --git a/PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.cxx b/PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.cxx index 891e57378de..e18eac75b88 100644 --- a/PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.cxx +++ b/PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.cxx @@ -8,6 +8,7 @@ #include #include #include +#include #include #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 4b32efea224..1daf5a290a6 100644 --- 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 -- 2.43.0