[u/mrichter/AliRoot.git] / PWGGA / EMCALJetTasks / AliAnalysisTaskEmcalJetSpectra.h

#ifndef AliAnalysisTaskEmcalJetSpectra_h
#define AliAnalysisTaskEmcalJetSpectra_h

// $Id$

class TList;
class TH1F;
class TH2F;
class AliESDEvent;

#include "AliAnalysisTaskSE.h"

class AliAnalysisTaskEmcalJetSpectra : public AliAnalysisTaskSE {
 public:
  AliAnalysisTaskEmcalJetSpectra();
  AliAnalysisTaskEmcalJetSpectra(const char *name);
  virtual ~AliAnalysisTaskEmcalJetSpectra() {}
  
  virtual void           UserCreateOutputObjects();
  virtual void           UserExec(Option_t *option);
  virtual void           Terminate(Option_t *);

  virtual void           SetAreaCut(Double_t a)                   { fAreacut    = a; }
  virtual void           SetJetEta(Double_t emin, Double_t emax)  { fEtamin = emin; fEtamax = emax; }
  virtual void           SetJetPhi(Double_t pmin, Double_t pmax)  { fPhimin = pmin; fPhimax = pmax; }
  virtual void           SetJetsName(const char *n)               { fJetsName   = n; }
  virtual void           SetRhos1Name(const char *n)              { fRhos1Name  = n; }
  virtual void           SetRhos2Name(const char *n)              { fRhos2Name  = n; }
  virtual void           SetRhos3Name(const char *n)              { fRhos3Name  = n; }
  virtual void           SetTracksName(const char *n)             { fTracksName = n; }
  
 protected:
  virtual Int_t          GetCentBin(Double_t cent) const;
   
 private:
  TString                fTracksName;              // name of track collection
  TString                fJetsName;                // name of jet collection
  TString                fClustersName;            // name of clusters collection
  TString                fRhos1Name;               // name of Rho1 array output
  TString                fRhos2Name;               // name of Rho2 array output
  TString                fRhos3Name;               // name of Rho2 array output
  Double_t               fPhimin;                  // phi min
  Double_t               fPhimax;                  // phi max
  Double_t               fEtamin;                  // eta min
  Double_t               fEtamax;                  // eta max
  Double_t               fAreacut;                 // area cut

  AliESDEvent           *fESD;                     //!esd event
  TList                 *fOutputList;              //!output list
  TH1F                  *fHistCentrality;          //!centrality
  TH2F                  *fHistDeltaRho12vsCent;    //!delta rho1 and rho2 vs centrality
  TH2F                  *fHistDeltaRho13vsCent;    //!delta rho1 and rho3 vs centrality
  TH2F                  *fHistDeltaRho23vsCent;    //!delta rho2 and rho3 vs centrality
  TH2F                  *fHistDeltaJetPt12vsCent;  //!delta jet pt1 and pt2 vs centrality 
  TH2F                  *fHistDeltaJetPt13vsCent;  //!delta jet pt1 and pt3 vs centrality 
  TH2F                  *fHistDeltaJetPt23vsCent;  //!delta jet pt2 and pt3 vs centrality 
  TH2F                  *fHistRho1vsCent;          //!rho1 vs centrality
  TH2F                  *fHistRho2vsCent;          //!rho2 vs centrality
  TH2F                  *fHistRho3vsCent;          //!rho3 vs centrality
  TH2F                  *fHistNEFvsPt[6][4];       //!neutral energy fraction vs pt
  TH2F                  *fHistZvsPt[6][4];         //!z all vs pt
  TH2F                  *fHistZchvsPt[6][4];       //!z charged vs pt
  TH2F                  *fHistZemvsPt[6][4];       //!z neutral vs pt
  TH1F                  *fHistJetPt[6][4];         //!jet pt
  TH1F                  *fHistJetPt5[6][4];        //!jet pt>5
  TH1F                  *fHistJetPt6[6][4];        //!jet pt>6
  TH1F                  *fHistJetPt7[6][4];        //!jet pt>7
  TH1F                  *fHistJetPt8[6][4];        //!jet pt>8
  TH2F                  *fHistNconsvsPt[6][4];     //!constituents vs pt
  TH1F                  *fHistRawJetPt[6];         //!raw jet pt
  TH2F                  *fHistAreavsRawPt[6];      //!area vs raw pt

  AliAnalysisTaskEmcalJetSpectra(const AliAnalysisTaskEmcalJetSpectra&); // not implemented
  AliAnalysisTaskEmcalJetSpectra& operator=(const AliAnalysisTaskEmcalJetSpectra&); // not implemented
  
  ClassDef(AliAnalysisTaskEmcalJetSpectra, 2); // Emcal jet spectra task
};
#endif
Commit	Line	Data
9993af2d	1	#ifndef AliAnalysisTaskEmcalJetSpectra_h
9993af2d	2	#define AliAnalysisTaskEmcalJetSpectra_h
020052e4	3
9993af2d	4	// $Id$
020052e4	5
	6	class TList;
	7	class TH1F;
	8	class TH2F;
	9	class AliESDEvent;
	10
	11	#include "AliAnalysisTaskSE.h"
	12
	13	class AliAnalysisTaskEmcalJetSpectra : public AliAnalysisTaskSE {
	14	public:
	15	AliAnalysisTaskEmcalJetSpectra();
	16	AliAnalysisTaskEmcalJetSpectra(const char *name);
	17	virtual ~AliAnalysisTaskEmcalJetSpectra() {}
	18
9993af2d	19	virtual void UserCreateOutputObjects();
	20	virtual void UserExec(Option_t *option);
	21	virtual void Terminate(Option_t *);
	22
	23	virtual void SetAreaCut(Double_t a) { fAreacut = a; }
	24	virtual void SetJetEta(Double_t emin, Double_t emax) { fEtamin = emin; fEtamax = emax; }
	25	virtual void SetJetPhi(Double_t pmin, Double_t pmax) { fPhimin = pmin; fPhimax = pmax; }
	26	virtual void SetJetsName(const char *n) { fJetsName = n; }
	27	virtual void SetRhos1Name(const char *n) { fRhos1Name = n; }
	28	virtual void SetRhos2Name(const char *n) { fRhos2Name = n; }
	29	virtual void SetRhos3Name(const char *n) { fRhos3Name = n; }
	30	virtual void SetTracksName(const char *n) { fTracksName = n; }
020052e4	31
020052e4	32	protected:
9993af2d	33	virtual Int_t GetCentBin(Double_t cent) const;
9993af2d	34
020052e4	35	private:
9993af2d	36	TString fTracksName; // name of track collection
	37	TString fJetsName; // name of jet collection
	38	TString fClustersName; // name of clusters collection
	39	TString fRhos1Name; // name of Rho1 array output
	40	TString fRhos2Name; // name of Rho2 array output
	41	TString fRhos3Name; // name of Rho2 array output
	42	Double_t fPhimin; // phi min
	43	Double_t fPhimax; // phi max
	44	Double_t fEtamin; // eta min
	45	Double_t fEtamax; // eta max
	46	Double_t fAreacut; // area cut
74ec9ac9	47
9993af2d	48	AliESDEvent *fESD; //!esd event
	49	TList *fOutputList; //!output list
	50	TH1F *fHistCentrality; //!centrality
	51	TH2F *fHistDeltaRho12vsCent; //!delta rho1 and rho2 vs centrality
	52	TH2F *fHistDeltaRho13vsCent; //!delta rho1 and rho3 vs centrality
	53	TH2F *fHistDeltaRho23vsCent; //!delta rho2 and rho3 vs centrality
	54	TH2F *fHistDeltaJetPt12vsCent; //!delta jet pt1 and pt2 vs centrality
	55	TH2F *fHistDeltaJetPt13vsCent; //!delta jet pt1 and pt3 vs centrality
	56	TH2F *fHistDeltaJetPt23vsCent; //!delta jet pt2 and pt3 vs centrality
	57	TH2F *fHistRho1vsCent; //!rho1 vs centrality
	58	TH2F *fHistRho2vsCent; //!rho2 vs centrality
	59	TH2F *fHistRho3vsCent; //!rho3 vs centrality
	60	TH2F *fHistNEFvsPt[6][4]; //!neutral energy fraction vs pt
	61	TH2F *fHistZvsPt[6][4]; //!z all vs pt
	62	TH2F *fHistZchvsPt[6][4]; //!z charged vs pt
	63	TH2F *fHistZemvsPt[6][4]; //!z neutral vs pt
	64	TH1F *fHistJetPt[6][4]; //!jet pt
	65	TH1F *fHistJetPt5[6][4]; //!jet pt>5
	66	TH1F *fHistJetPt6[6][4]; //!jet pt>6
	67	TH1F *fHistJetPt7[6][4]; //!jet pt>7
	68	TH1F *fHistJetPt8[6][4]; //!jet pt>8
	69	TH2F *fHistNconsvsPt[6][4]; //!constituents vs pt
	70	TH1F *fHistRawJetPt[6]; //!raw jet pt
	71	TH2F *fHistAreavsRawPt[6]; //!area vs raw pt
74ec9ac9	72
9993af2d	73	AliAnalysisTaskEmcalJetSpectra(const AliAnalysisTaskEmcalJetSpectra&); // not implemented
020052e4	74	AliAnalysisTaskEmcalJetSpectra& operator=(const AliAnalysisTaskEmcalJetSpectra&); // not implemented
020052e4	75
9993af2d	76	ClassDef(AliAnalysisTaskEmcalJetSpectra, 2); // Emcal jet spectra task
020052e4	77	};
020052e4	78	#endif