[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);

  void                        SetCentBins(Int_t n, Double_t* bins);
  void                        SetJetPtBins(Int_t n, Double_t* bins);

  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                       nCentBins1;
  Double_t*		      centBins; //[nCentBins1]
  Int_t			      nJetPtBins;
  Int_t			      nJetPtBins1;
  Double_t*		      jetPtBins; //[nJetPtBins1]
  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);
de2a3fac	28
	29	void SetCentBins(Int_t n, Double_t* bins);
	30	void SetJetPtBins(Int_t n, Double_t* bins);
	31
51a8c8fd	32	Double_t GetJetv2(){return fJetv2;};
	33	void SetJetv2(Double_t jetv2){fJetv2 = jetv2;};
	34	Bool_t GetDoPtWeight(){return doPtWeight;};
	35	void SetDoPtWeight(Bool_t ptWeight){doPtWeight = ptWeight;};
	36
	37	protected:
	38	Int_t nCentBins;
de2a3fac	39	Int_t nCentBins1;
de2a3fac	40	Double_t* centBins; //[nCentBins1]
51a8c8fd	41	Int_t nJetPtBins;
de2a3fac	42	Int_t nJetPtBins1;
de2a3fac	43	Double_t* jetPtBins; //[nJetPtBins1]
51a8c8fd	44	void ExecOnce();
	45	Bool_t FillHistograms();
	46	Bool_t Run();
	47	Double_t fJetv2;
	48	Bool_t doPtWeight;
	49
	50	// General histograms
	51	TH1F *fHistTracksPt; //!Track pt spectrum
	52	TH1F *fHistClustersPt; //!Cluster pt spectrum
	53	TH1F *fHistLeadingJetPt; //!Leading jet pt spectrum
	54	TH1F *fHistLeadingJetPtCorr; //!Leading jet pt spectrum, background subtracted
	55	TH2F *fHistJetsPhiEta; //!Phi-Eta distribution of jets
	56	TH2F *fHistJetsPtArea; //!Jet pt vs. area
	57	TH2F *fHistJetsPtLeadHad; //!Jet pt vs. leading hadron
	58	TH2F *fHistJetsCorrPtArea; //!Jet pt - bkg vs. area
	59	TH3F *fHistPtDEtaDPhiTrackClus; //!track pt, delta eta, delta phi to matched cluster
	60	TH3F *fHistPtDEtaDPhiClusTrack; //!cluster pt, delta eta, delta phi to matched track
	61	TH1F *fDPhiJet; //!
	62	TH1F *fDPhiJetPythia; //!
	63	TH1F *fDPhiEP; //!
	64	TH2D *hGx; //!
	65	TH2D *hGy2; //!
	66	TH2D *hGxGy2; //!
	67	TH2D *hGy4; //!
	68	TH2D *hGx2; //!
	69	TH2D *hGx2Gy2; //!
	70	TH2D *hGxGy4; //!
	71	TH2D *hGy6; //!
	72	TH2D *hGx2Gy4; //!
	73	TH2D *hGxGy6; //!
	74	TH2D *hGy8; //!
	75	TH2D *hGy; //!
	76	TH2D *hN; //!
	77	TH2D *htv2std; //!
	78	TH2D *htjv2std; //!
	79	TH2D *htj2v2std; //!
	80	TH2D *hV0jv2std; //!
	81	TH3D *htdPsi; //!
	82	TH3D *htjdPsi; //!
	83	TH3D *htj2dPsi; //!
	84	TH3D *hV0jdPsi; //!
	85	TH2D *hAx; //!
	86	TH2D *hAxDijet; //!
	87	TH2D *hQx; //!
	88	TH2D *hQy; //!
	89	TH1F *hEventData; //!
	90	TH2F *hNTracks; //!
	91	TH2D *hNTracksCent; //!
	92	TH2D *hGxTracks; //!
	93	TH2D *hGyTracks; //!
	94	TH2D *hGy2Tracks; //!
	95	TH2D *hGxGy2Tracks; //!
	96	TH2D *hGy4Tracks; //!
	97	TH2D *htEPRes; //!
	98	TH2D *htjEPRes; //!
	99	TH2D *htj2EPRes; //!
	100
	101	AliJetContainer *fJetsCont; //!Jets
	102	AliParticleContainer *fTracksCont; //!Tracks
	103	AliClusterContainer *fCaloClustersCont; //!Clusters
	104
	105
	106
	107	private:
108	AliAnalysisTaskEmcalJetv2QA(const AliAnalysisTaskEmcalJetv2QA&); // not implemented
109	AliAnalysisTaskEmcalJetv2QA &operator=(const AliAnalysisTaskEmcalJetv2QA&); // not implemented
110
111	ClassDef(AliAnalysisTaskEmcalJetv2QA, 3) // jet v2QA analysis task
112	};
113	#endif