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

#ifndef AliAnalysisTaskEmcalJetSpectra_cxx
#define AliAnalysisTaskEmcalJetSpectra_cxx


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   SetTracksName(const char *n)             { fTracksName = n;     }
  virtual void   SetJetsName(const char *n)               { fJetsName =   n;     }
  virtual void   SetRhos1Name(const char *rn)             { fRhos1Name=   rn;    }
  virtual void   SetRhos2Name(const char *rn)             { fRhos2Name=   rn;    }
  virtual void   SetRhos3Name(const char *rn)             { fRhos3Name=   rn;    }
  virtual void   SetJetPhi(Double_t pmin, Double_t pmax)  { phimin =      pmin;  phimax = pmax; }
  virtual void   SetJetEta(Double_t emin, Double_t emax)  { etamin =      emin;  etamax = emax; }
  virtual void   SetAreaCut(Double_t a)                   { areacut =     a;     }
  
 protected:
  virtual Int_t GetCentBin(Double_t cent) const;
 
  
 private:
   AliESDEvent *fESD;    //! ESD object
  TList       *fOutputList; //! Output list
  TH1F        *fHistCentrality;
  TH1F        *fHistJetArea; 
  TH1F        *fHistJetMaxPt;
  TH1F        *fHistJetZ;
  TH1F        *fHistJetNEF;
  TH2F        *fHistJetPtvsCent;  
  TH2F        *fHistJetPtM3vsCent;  
  TH2F        *fHistLeadingJetPtvsCent;  
  TH2F        *fHistLeadingJetPtM3vsCent;  
  TH2F        *fHistJetAreavsCent;  
  TH2F        *fHistJetMaxPtvsCent;  
  TH2F        *fHistJetZvsCent;  
  TH2F        *fHistJetNEFvsCent;
  TH2F        *fHistNjetvsCent;
  
  TH2F        *fHistJetPtvsNtrack;
  TH2F        *fHistJetAreavsNtrack;
  TH2F        *fHistJetMaxPtvsNtrack;
  TH2F        *fHistJetZvsNtrack;
  TH2F        *fHistJetNEFvsNtrack;
  TH2F        *fHistNjetvsNtrack;

  TH2F       *fHistDeltaRho12vsCent;
  TH2F       *fHistDeltaRho13vsCent;
  TH2F       *fHistDeltaRho23vsCent;

  TH2F       *fHistDeltaJetPt12vsCent;
  TH2F       *fHistDeltaJetPt13vsCent;
  TH2F       *fHistDeltaJetPt23vsCent;
  
  TH2F       *fHistRho1vsCent;
  TH2F       *fHistRho2vsCent;
  TH2F       *fHistRho3vsCent;

  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 phimin;
  Double_t phimax;
  Double_t etamin;
  Double_t etamax;
  Double_t areacut;

  TH2F       *fHistNEFvsPt[6][4];
  TH2F       *fHistZvsPt[6][4];
  TH2F       *fHistZchvsPt[6][4];
  TH2F       *fHistZemvsPt[6][4];
  TH2F       *fHistAreavsPt[6][4];
  TH1F       *fHistJetPt[6][4];
  TH1F       *fHistJetPt5[6][4];
  TH1F       *fHistJetPt6[6][4];
  TH1F       *fHistJetPt7[6][4];
  TH1F       *fHistJetPt8[6][4];
  TH2F       *fHistNconsvsPt[6][4];
  TH1F       *fHistRawJetPt[6];
  TH2F       *fHistAreavsRawPt[6];


   AliAnalysisTaskEmcalJetSpectra(const AliAnalysisTaskEmcalJetSpectra&); // not implemented
  AliAnalysisTaskEmcalJetSpectra& operator=(const AliAnalysisTaskEmcalJetSpectra&); // not implemented
  
  ClassDef(AliAnalysisTaskEmcalJetSpectra, 1); // example of analysis
};

#endif
Commit	Line	Data
020052e4	1	#ifndef AliAnalysisTaskEmcalJetSpectra_cxx
	2	#define AliAnalysisTaskEmcalJetSpectra_cxx
	3
	4
	5	class TList;
	6	class TH1F;
	7	class TH2F;
	8	class AliESDEvent;
	9
	10	#include "AliAnalysisTaskSE.h"
	11
	12	class AliAnalysisTaskEmcalJetSpectra : public AliAnalysisTaskSE {
	13	public:
	14	AliAnalysisTaskEmcalJetSpectra();
	15	AliAnalysisTaskEmcalJetSpectra(const char *name);
	16	virtual ~AliAnalysisTaskEmcalJetSpectra() {}
	17
	18	virtual void UserCreateOutputObjects();
	19	virtual void UserExec(Option_t *option);
	20	virtual void Terminate(Option_t *);
6172360f	21	virtual void SetTracksName(const char *n) { fTracksName = n; }
	22	virtual void SetJetsName(const char *n) { fJetsName = n; }
	23	virtual void SetRhos1Name(const char *rn) { fRhos1Name= rn; }
	24	virtual void SetRhos2Name(const char *rn) { fRhos2Name= rn; }
	25	virtual void SetRhos3Name(const char *rn) { fRhos3Name= rn; }
	26	virtual void SetJetPhi(Double_t pmin, Double_t pmax) { phimin = pmin; phimax = pmax; }
	27	virtual void SetJetEta(Double_t emin, Double_t emax) { etamin = emin; etamax = emax; }
	28	virtual void SetAreaCut(Double_t a) { areacut = a; }
020052e4	29
	30	protected:
	31	virtual Int_t GetCentBin(Double_t cent) const;
	32
	33
	34	private:
74ec9ac9	35	AliESDEvent *fESD; //! ESD object
020052e4	36	TList *fOutputList; //! Output list
	37	TH1F *fHistCentrality;
	38	TH1F *fHistJetArea;
	39	TH1F *fHistJetMaxPt;
	40	TH1F *fHistJetZ;
	41	TH1F *fHistJetNEF;
	42	TH2F *fHistJetPtvsCent;
	43	TH2F *fHistJetPtM3vsCent;
	44	TH2F *fHistLeadingJetPtvsCent;
	45	TH2F *fHistLeadingJetPtM3vsCent;
	46	TH2F *fHistJetAreavsCent;
	47	TH2F *fHistJetMaxPtvsCent;
	48	TH2F *fHistJetZvsCent;
	49	TH2F *fHistJetNEFvsCent;
	50	TH2F *fHistNjetvsCent;
	51
	52	TH2F *fHistJetPtvsNtrack;
	53	TH2F *fHistJetAreavsNtrack;
	54	TH2F *fHistJetMaxPtvsNtrack;
	55	TH2F *fHistJetZvsNtrack;
	56	TH2F *fHistJetNEFvsNtrack;
	57	TH2F *fHistNjetvsNtrack;
	58
74ec9ac9	59	TH2F *fHistDeltaRho12vsCent;
	60	TH2F *fHistDeltaRho13vsCent;
	61	TH2F *fHistDeltaRho23vsCent;
6172360f	62
	63	TH2F *fHistDeltaJetPt12vsCent;
	64	TH2F *fHistDeltaJetPt13vsCent;
	65	TH2F *fHistDeltaJetPt23vsCent;
74ec9ac9	66
020052e4	67	TH2F *fHistRho1vsCent;
	68	TH2F *fHistRho2vsCent;
	69	TH2F *fHistRho3vsCent;
	70
74ec9ac9	71	TString fTracksName; // name of track collection
	72	TString fJetsName; // name of jet collection
	73	TString fClustersName; // name of clusters collection
	74	TString fRhos1Name; // name of Rho1 array output
	75	TString fRhos2Name; // name of Rho2 array output
	76	TString fRhos3Name; // name of Rho2 array output
020052e4	77
	78	Double_t phimin;
	79	Double_t phimax;
	80	Double_t etamin;
	81	Double_t etamax;
	82	Double_t areacut;
74ec9ac9	83
6172360f	84	TH2F *fHistNEFvsPt[6][4];
74ec9ac9	85	TH2F *fHistZvsPt[6][4];
	86	TH2F *fHistZchvsPt[6][4];
	87	TH2F *fHistZemvsPt[6][4];
	88	TH2F *fHistAreavsPt[6][4];
	89	TH1F *fHistJetPt[6][4];
6172360f	90	TH1F *fHistJetPt5[6][4];
	91	TH1F *fHistJetPt6[6][4];
	92	TH1F *fHistJetPt7[6][4];
	93	TH1F *fHistJetPt8[6][4];
74ec9ac9	94	TH2F *fHistNconsvsPt[6][4];
	95	TH1F *fHistRawJetPt[6];
	96	TH2F *fHistAreavsRawPt[6];
	97
020052e4	98
	99	AliAnalysisTaskEmcalJetSpectra(const AliAnalysisTaskEmcalJetSpectra&); // not implemented
	100	AliAnalysisTaskEmcalJetSpectra& operator=(const AliAnalysisTaskEmcalJetSpectra&); // not implemented
	101
	102	ClassDef(AliAnalysisTaskEmcalJetSpectra, 1); // example of analysis
	103	};
	104
	105	#endif