From 6a12db523773edaedcf4238f10ac50f1dfcc91b3 Mon Sep 17 00:00:00 2001 From: loizides Date: Thu, 21 Jun 2012 04:33:04 +0000 Subject: [PATCH] added emc area cut --- PWGGA/EMCALJetTasks/AliAnalysisTaskRho.cxx | 13 ++++++---- PWGGA/EMCALJetTasks/AliAnalysisTaskRho.h | 28 ++++++++++++---------- 2 files changed, 23 insertions(+), 18 deletions(-) diff --git a/PWGGA/EMCALJetTasks/AliAnalysisTaskRho.cxx b/PWGGA/EMCALJetTasks/AliAnalysisTaskRho.cxx index 74fd12446cf..54895db9f9d 100644 --- a/PWGGA/EMCALJetTasks/AliAnalysisTaskRho.cxx +++ b/PWGGA/EMCALJetTasks/AliAnalysisTaskRho.cxx @@ -28,14 +28,15 @@ ClassImp(AliAnalysisTaskRho) //________________________________________________________________________ AliAnalysisTaskRho::AliAnalysisTaskRho() : AliAnalysisTaskRhoBase(), - fTracksName("tracks"), - fJetsName("KtJets"), + fTracksName(), + fJetsName(), fRhoScaledName(""), fPhiMin(0), fPhiMax(0), fEtaMin(0), fEtaMax(0), fAreaCut(0), + fAreaEmcCut(0), fNExclLeadJets(0), fScaleFunction(0), fCreateHisto(kFALSE), @@ -59,7 +60,7 @@ AliAnalysisTaskRho::AliAnalysisTaskRho() : fHistNjetvsNtrack(0), fRhoScaled(0) { - // Constructor + // Constructor. } //________________________________________________________________________ @@ -73,6 +74,7 @@ AliAnalysisTaskRho::AliAnalysisTaskRho(const char *name, Bool_t histo) : fEtaMin(-0.5), fEtaMax(+0.5), fAreaCut(0.01), + fAreaEmcCut(0), fNExclLeadJets(1), fScaleFunction(0), fCreateHisto(histo), @@ -96,7 +98,7 @@ AliAnalysisTaskRho::AliAnalysisTaskRho(const char *name, Bool_t histo) : fHistNjetvsNtrack(0), fRhoScaled(0) { - // Constructor + // Constructor. if (fCreateHisto) DefineOutput(1, TList::Class()); @@ -224,6 +226,8 @@ void AliAnalysisTaskRho::UserExec(Option_t *) // applying some other cuts if (jet->Area() < fAreaCut) continue; + if (jet->AreaEmc()Phi() < fPhiMin) || (jet->Phi() > fPhiMax)) continue; if ((jet->Eta() < fEtaMin) || (jet->Eta() > fEtaMax)) @@ -249,7 +253,6 @@ void AliAnalysisTaskRho::UserExec(Option_t *) //find median value Double_t rho0 = TMath::Median(NjetAcc, rhovec); fRho->SetVal(rho0); - Double_t rhoScaled = rho0; if (fRhoScaled) { rhoScaled *= GetScaleFactor(fCent); diff --git a/PWGGA/EMCALJetTasks/AliAnalysisTaskRho.h b/PWGGA/EMCALJetTasks/AliAnalysisTaskRho.h index f000cc53d2e..81182b725ec 100644 --- a/PWGGA/EMCALJetTasks/AliAnalysisTaskRho.h +++ b/PWGGA/EMCALJetTasks/AliAnalysisTaskRho.h @@ -3,12 +3,12 @@ // $Id$ -class TList; +class TClonesArray; +class TF1; class TH1F; class TH2F; -class TClonesArray; +class TList; class TString; -class TF1; class AliRhoParameter; #include "AliAnalysisTaskRhoBase.h" @@ -23,17 +23,18 @@ class AliAnalysisTaskRho : public AliAnalysisTaskRhoBase { void UserCreateOutputObjects(); void UserExec(Option_t*); - const char *GetRhoScaled() const { return fRhoScaledName; } - void SetAreaCut(Double_t a = 0.0) { fAreaCut = a ; } - void SetExcludeLeadJets(UInt_t n) { fNExclLeadJets = n ; } + const char *GetRhoScaled() const { return fRhoScaledName ; } + void SetAreaCut(Double_t a = 0.05) { fAreaCut = a ; } + void SetAreaEmcCut(Double_t a = 0.99) { fAreaEmcCut = a ; } + void SetExcludeLeadJets(UInt_t n) { fNExclLeadJets = n ; } void SetJetEta(Double_t emin, Double_t emax) { fEtaMin = emin ; fEtaMax = emax ; } void SetJetPhi(Double_t pmin, Double_t pmax) { fPhiMin = pmin ; fPhiMax = pmax ; } - void SetJetsName(const char *n) { fJetsName = n ; } - void SetScaleFunction(TF1* sf) { fScaleFunction = sf ; } - void SetTracksName(const char *n) { fTracksName = n ; } + void SetJetsName(const char *n) { fJetsName = n ; } void SetRhoName(const char *name) { fRhoName = name ; - fRhoScaledName = name; - fRhoScaledName += "_Scaled"; } + fRhoScaledName = name ; + fRhoScaledName += "_Scaled"; } + void SetScaleFunction(TF1* sf) { fScaleFunction = sf ; } + void SetTracksName(const char *n) { fTracksName = n ; } protected: virtual void ExecOnce(); @@ -46,7 +47,8 @@ class AliAnalysisTaskRho : public AliAnalysisTaskRhoBase { Double_t fPhiMax; // maximum phi Double_t fEtaMin; // minimum eta Double_t fEtaMax; // maximum eta - Double_t fAreaCut; // cut on jet area + Double_t fAreaCut; // minimum cut on jet area + Double_t fAreaEmcCut; // minimum cut on jet emcal area UInt_t fNExclLeadJets; // number of leading jets to be excluded from the median calculation TF1 *fScaleFunction; // pre-computed scale factor as a function of centrality Bool_t fCreateHisto; // whether or not create histograms @@ -73,6 +75,6 @@ class AliAnalysisTaskRho : public AliAnalysisTaskRhoBase { AliAnalysisTaskRho(const AliAnalysisTaskRho&); // not implemented AliAnalysisTaskRho& operator=(const AliAnalysisTaskRho&); // not implemented - ClassDef(AliAnalysisTaskRho, 6); // Rho task + ClassDef(AliAnalysisTaskRho, 7); // Rho task }; #endif -- 2.43.0