PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskRhoMass.h

   1 #ifndef ALIANALYSISTASKRHOMASS_H
   2 #define ALIANALYSISTASKRHOMASS_H
   3
   4 // $Id$
   5
   6 #include "AliAnalysisTaskRhoMassBase.h"
   7
   8 class AliAnalysisTaskRhoMass : public AliAnalysisTaskRhoMassBase {
   9
  10  public:
  11   AliAnalysisTaskRhoMass();
  12   AliAnalysisTaskRhoMass(const char *name, Bool_t histo=kFALSE);
  13   virtual ~AliAnalysisTaskRhoMass() {}
  14
  15   enum JetRhoMassType {
  16     kMd     = 0,            //rho_m from arXiv:1211.2811
  17     kMdP    = 1,            //rho_m using P instead of pT
  18     kMd4    = 2             //rho_m using addition of 4-vectors
  19   };
  20
  21   void             UserCreateOutputObjects();
  22
  23   void             SetExcludeLeadJets(UInt_t n)     { fNExclLeadJets  = n   ; }
  24   void             SetRhoMassType(JetRhoMassType t) { fJetRhoMassType = t   ; }
  25   void             SetPionMassForClusters(Bool_t b) { fPionMassClusters = b ; }
  26
  27  protected:
  28   Bool_t           Run();
  29
  30   Double_t         GetSumMConstituents(AliEmcalJet *jet);
  31   Double_t         GetSumPtConstituents(AliEmcalJet *jet);
  32   Double_t         GetMd(AliEmcalJet *jet);
  33
  34   UInt_t           fNExclLeadJets;                 // number of leading jets to be excluded from the median calculation
  35   JetRhoMassType   fJetRhoMassType;                // method for rho_m calculation
  36   Bool_t           fPionMassClusters;              // assume pion mass for clusters
  37
  38   TH2F            *fHistMdAreavsCent;              //! Md/Area vs cent for all kt clusters
  39
  40   AliAnalysisTaskRhoMass(const AliAnalysisTaskRhoMass&);             // not implemented
  41   AliAnalysisTaskRhoMass& operator=(const AliAnalysisTaskRhoMass&);  // not implemented
  42
  43   ClassDef(AliAnalysisTaskRhoMass, 2); // Rho_m task
  44 };
  45 #endif