added emc area cut

diff --git a/PWGGA/EMCALJetTasks/AliAnalysisTaskRho.cxx b/PWGGA/EMCALJetTasks/AliAnalysisTaskRho.cxx
index 74fd12446cf9ea4cc65624b4a3e93a922a75aa4f..54895db9f9d5c217004657039c93c32ecb650ebd 100644 (file)
--- 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()<fAreaEmcCut)
+      continue;
     if ((jet->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 f000cc53d2ea791ff7b52e66fc0356fd5d61e05f..81182b725ecf432c254064c450986c0855408f14 100644 (file)
--- 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