[u/mrichter/AliRoot.git] / PWGJE / EMCALJetTasks / UserTasks / AliAnalysisTaskRhoMass.h

#ifndef ALIANALYSISTASKRHOMASS_H
#define ALIANALYSISTASKRHOMASS_H

// $Id$

#include "AliAnalysisTaskRhoMassBase.h"

class AliAnalysisTaskRhoMass : public AliAnalysisTaskRhoMassBase {

 public:
  AliAnalysisTaskRhoMass();
  AliAnalysisTaskRhoMass(const char *name, Bool_t histo=kFALSE);
  virtual ~AliAnalysisTaskRhoMass() {}

  enum JetRhoMassType {
    kMd     = 0,            //rho_m from arXiv:1211.2811
    kMdP    = 1,            //rho_m using P instead of pT
    kMd4    = 2             //rho_m using addition of 4-vectors
  };

  void             UserCreateOutputObjects();

  void             SetExcludeLeadJets(UInt_t n)     { fNExclLeadJets  = n   ; }
  void             SetRhoMassType(JetRhoMassType t) { fJetRhoMassType = t   ; }

 protected:
  Bool_t           Run();

  Double_t         GetSumMConstituents(AliEmcalJet *jet);
  Double_t         GetSumPtConstituents(AliEmcalJet *jet);
  Double_t         GetMd(AliEmcalJet *jet);

  UInt_t           fNExclLeadJets;                 // number of leading jets to be excluded from the median calculation
  JetRhoMassType   fJetRhoMassType;                // method for rho_m calculation

  TH2F            *fHistMdAreavsCent;              //! Md/Area vs cent for all kt clusters

  AliAnalysisTaskRhoMass(const AliAnalysisTaskRhoMass&);             // not implemented
  AliAnalysisTaskRhoMass& operator=(const AliAnalysisTaskRhoMass&);  // not implemented
  
  ClassDef(AliAnalysisTaskRhoMass, 1); // Rho_m task
};
#endif
Commit	Line	Data
4d3b366f	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
	26	protected:
	27	Bool_t Run();
	28
	29	Double_t GetSumMConstituents(AliEmcalJet *jet);
	30	Double_t GetSumPtConstituents(AliEmcalJet *jet);
	31	Double_t GetMd(AliEmcalJet *jet);
	32
	33	UInt_t fNExclLeadJets; // number of leading jets to be excluded from the median calculation
	34	JetRhoMassType fJetRhoMassType; // method for rho_m calculation
	35
	36	TH2F *fHistMdAreavsCent; //! Md/Area vs cent for all kt clusters
	37
	38	AliAnalysisTaskRhoMass(const AliAnalysisTaskRhoMass&); // not implemented
	39	AliAnalysisTaskRhoMass& operator=(const AliAnalysisTaskRhoMass&); // not implemented
	40
	41	ClassDef(AliAnalysisTaskRhoMass, 1); // Rho_m task
	42	};
	43	#endif