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

// Jet v2 task using QA method, based on jet sample task (S.Aiola).
//
// Authors: Jason Mueller (CERN summer student 2014) & Alice Ohlson

#ifndef ALIANALYSISTASKEMCALJETVNQA_H
#define ALIANALYSISTASKEMCALJETVNQA_H

// $Id$

class TH1;
class TH2;
class TH3;
class AliJetContainer;
class AliParticleContainer;
class AliClusterContainer;

#include "AliAnalysisTaskEmcalJet.h"

class AliAnalysisTaskEmcalJetv2QA : public AliAnalysisTaskEmcalJet {
 public:

  AliAnalysisTaskEmcalJetv2QA();
  AliAnalysisTaskEmcalJetv2QA(const char *name);
  virtual ~AliAnalysisTaskEmcalJetv2QA();

  void                        UserCreateOutputObjects();
  void                        Terminate(Option_t *option);
  Double_t		      GetJetv2(){return fJetv2;};
  void			      SetJetv2(Double_t jetv2){fJetv2 = jetv2;};
  Bool_t		      GetDoPtWeight(){return doPtWeight;};
  void			      SetDoPtWeight(Bool_t ptWeight){doPtWeight = ptWeight;};

 protected:
  Int_t			      nCentBins;
  Int_t			      nJetPtBins;
  Double_t		      *centBins;
  Double_t		      *jetPtBins;
  void                        ExecOnce();
  Bool_t                      FillHistograms();
  Bool_t                      Run();
  Double_t		      fJetv2;
  Bool_t		      doPtWeight;

  // General histograms
  TH1F                       *fHistTracksPt;            //!Track pt spectrum
  TH1F                       *fHistClustersPt;          //!Cluster pt spectrum
  TH1F                       *fHistLeadingJetPt;        //!Leading jet pt spectrum
  TH1F                       *fHistLeadingJetPtCorr;    //!Leading jet pt spectrum, background subtracted
  TH2F                       *fHistJetsPhiEta;          //!Phi-Eta distribution of jets
  TH2F                       *fHistJetsPtArea;          //!Jet pt vs. area
  TH2F                       *fHistJetsPtLeadHad;       //!Jet pt vs. leading hadron
  TH2F                       *fHistJetsCorrPtArea;      //!Jet pt - bkg vs. area
  TH3F                       *fHistPtDEtaDPhiTrackClus; //!track pt, delta eta, delta phi to matched cluster
  TH3F                       *fHistPtDEtaDPhiClusTrack; //!cluster pt, delta eta, delta phi to matched track
  TH1F                       *fDPhiJet;                 //!
  TH1F                       *fDPhiJetPythia;           //!
  TH1F                       *fDPhiEP;                  //!
  TH2D 			     *hGx;		        //!
  TH2D                       *hGy2;                     //!
  TH2D                       *hGxGy2;                   //!
  TH2D                       *hGy4;                     //!
  TH2D                       *hGx2;                     //!
  TH2D                       *hGx2Gy2;                  //!
  TH2D                       *hGxGy4;                   //!
  TH2D                       *hGy6;                     //!
  TH2D                       *hGx2Gy4;                  //!
  TH2D                       *hGxGy6;                   //!
  TH2D                       *hGy8;                     //!
  TH2D                       *hGy;                      //!
  TH2D                       *hN;                       //!
  TH2D                       *htv2std;                  //!
  TH2D                       *htjv2std;                 //!
  TH2D                       *htj2v2std;                //!
  TH2D                       *hV0jv2std;                //!
  TH3D                       *htdPsi;                   //!
  TH3D                       *htjdPsi;                  //!
  TH3D                       *htj2dPsi;                 //!
  TH3D                       *hV0jdPsi;                 //!
  TH2D		 	     *hAx;		        //!
  TH2D			     *hAxDijet;		        //!
  TH2D			     *hQx;		        //!
  TH2D			     *hQy;		        //!
  TH1F			     *hEventData;	        //!
  TH2F			     *hNTracks;		        //!
  TH2D			     *hNTracksCent;		//!
  TH2D                       *hGxTracks;                //!
  TH2D                       *hGyTracks;                //!
  TH2D                       *hGy2Tracks;               //!
  TH2D                       *hGxGy2Tracks;             //!
  TH2D                       *hGy4Tracks;               //!
  TH2D                       *htEPRes;                  //!
  TH2D                       *htjEPRes;                 //!
  TH2D                       *htj2EPRes;                //!

  AliJetContainer            *fJetsCont;                   //!Jets
  AliParticleContainer       *fTracksCont;                 //!Tracks
  AliClusterContainer        *fCaloClustersCont;           //!Clusters


 private:
  AliAnalysisTaskEmcalJetv2QA(const AliAnalysisTaskEmcalJetv2QA&);            // not implemented
  AliAnalysisTaskEmcalJetv2QA &operator=(const AliAnalysisTaskEmcalJetv2QA&); // not implemented

  ClassDef(AliAnalysisTaskEmcalJetv2QA, 3) // jet v2QA analysis task
};
#endif
Commit	Line	Data
51a8c8fd	1	// Jet v2 task using QA method, based on jet sample task (S.Aiola).
	2	//
	3	// Authors: Jason Mueller (CERN summer student 2014) & Alice Ohlson
	4
	5	#ifndef ALIANALYSISTASKEMCALJETVNQA_H
	6	#define ALIANALYSISTASKEMCALJETVNQA_H
	7
	8	// $Id$
	9
	10	class TH1;
	11	class TH2;
	12	class TH3;
	13	class AliJetContainer;
	14	class AliParticleContainer;
	15	class AliClusterContainer;
	16
	17	#include "AliAnalysisTaskEmcalJet.h"
	18
	19	class AliAnalysisTaskEmcalJetv2QA : public AliAnalysisTaskEmcalJet {
	20	public:
	21
	22	AliAnalysisTaskEmcalJetv2QA();
	23	AliAnalysisTaskEmcalJetv2QA(const char *name);
	24	virtual ~AliAnalysisTaskEmcalJetv2QA();
	25
	26	void UserCreateOutputObjects();
	27	void Terminate(Option_t *option);
	28	Double_t GetJetv2(){return fJetv2;};
	29	void SetJetv2(Double_t jetv2){fJetv2 = jetv2;};
	30	Bool_t GetDoPtWeight(){return doPtWeight;};
	31	void SetDoPtWeight(Bool_t ptWeight){doPtWeight = ptWeight;};
	32
	33	protected:
	34	Int_t nCentBins;
	35	Int_t nJetPtBins;
	36	Double_t *centBins;
	37	Double_t *jetPtBins;
	38	void ExecOnce();
	39	Bool_t FillHistograms();
	40	Bool_t Run();
	41	Double_t fJetv2;
	42	Bool_t doPtWeight;
	43
	44	// General histograms
	45	TH1F *fHistTracksPt; //!Track pt spectrum
	46	TH1F *fHistClustersPt; //!Cluster pt spectrum
	47	TH1F *fHistLeadingJetPt; //!Leading jet pt spectrum
	48	TH1F *fHistLeadingJetPtCorr; //!Leading jet pt spectrum, background subtracted
	49	TH2F *fHistJetsPhiEta; //!Phi-Eta distribution of jets
	50	TH2F *fHistJetsPtArea; //!Jet pt vs. area
	51	TH2F *fHistJetsPtLeadHad; //!Jet pt vs. leading hadron
	52	TH2F *fHistJetsCorrPtArea; //!Jet pt - bkg vs. area
	53	TH3F *fHistPtDEtaDPhiTrackClus; //!track pt, delta eta, delta phi to matched cluster
	54	TH3F *fHistPtDEtaDPhiClusTrack; //!cluster pt, delta eta, delta phi to matched track
	55	TH1F *fDPhiJet; //!
	56	TH1F *fDPhiJetPythia; //!
	57	TH1F *fDPhiEP; //!
	58	TH2D *hGx; //!
	59	TH2D *hGy2; //!
	60	TH2D *hGxGy2; //!
	61	TH2D *hGy4; //!
	62	TH2D *hGx2; //!
	63	TH2D *hGx2Gy2; //!
	64	TH2D *hGxGy4; //!
65	TH2D *hGy6; //!
66	TH2D *hGx2Gy4; //!
67	TH2D *hGxGy6; //!
68	TH2D *hGy8; //!
69	TH2D *hGy; //!
70	TH2D *hN; //!
71	TH2D *htv2std; //!
72	TH2D *htjv2std; //!
73	TH2D *htj2v2std; //!
74	TH2D *hV0jv2std; //!
75	TH3D *htdPsi; //!
76	TH3D *htjdPsi; //!
77	TH3D *htj2dPsi; //!
78	TH3D *hV0jdPsi; //!
79	TH2D *hAx; //!
80	TH2D *hAxDijet; //!
81	TH2D *hQx; //!
82	TH2D *hQy; //!
83	TH1F *hEventData; //!
84	TH2F *hNTracks; //!
85	TH2D *hNTracksCent; //!
86	TH2D *hGxTracks; //!
87	TH2D *hGyTracks; //!
88	TH2D *hGy2Tracks; //!
89	TH2D *hGxGy2Tracks; //!
90	TH2D *hGy4Tracks; //!
91	TH2D *htEPRes; //!
92	TH2D *htjEPRes; //!
93	TH2D *htj2EPRes; //!
94
95	AliJetContainer *fJetsCont; //!Jets
96	AliParticleContainer *fTracksCont; //!Tracks
97	AliClusterContainer *fCaloClustersCont; //!Clusters
98
99
100
101	private:
102	AliAnalysisTaskEmcalJetv2QA(const AliAnalysisTaskEmcalJetv2QA&); // not implemented
103	AliAnalysisTaskEmcalJetv2QA &operator=(const AliAnalysisTaskEmcalJetv2QA&); // not implemented
104
105	ClassDef(AliAnalysisTaskEmcalJetv2QA, 3) // jet v2QA analysis task
106	};
107	#endif