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

#ifndef ALIANALYSISTASKSAJF_H
#define ALIANALYSISTASKSAJF_H

class TH2;
class THnSparse;

#include <TH3F.h>

#include "AliAnalysisTaskEmcalJet.h"

class AliAnalysisTaskSAJF : public AliAnalysisTaskEmcalJet {
 public:

  AliAnalysisTaskSAJF();
  AliAnalysisTaskSAJF(const char *name);
  virtual ~AliAnalysisTaskSAJF() {;}

  void                        UserCreateOutputObjects();

  void                        SetHistoType(Int_t t) { fHistoType = t; }

 protected:
  void                        AllocateTHX();
  void                        AllocateTHnSparse();

  Bool_t                      FillHistograms();
  void                        FillJetHisto(Double_t cent, Double_t ep, Double_t eta, Double_t phi, Double_t pt, Double_t MCpt, Double_t corrpt, Double_t area, 
					   Double_t NEF, Double_t z, Int_t n, Double_t leadingpt);

  Int_t                       fHistoType;                      // histogram type (0=TH2, 1=THnSparse)

  TH2                       **fHistRejectionReason;         //!Rejection reason vs. jet pt
  TH2                       **fHistTracksJetPt;             //!Track pt vs. jet pt
  TH2                       **fHistClustersJetPt;           //!Cluster pt vs. jet pt
  TH2                       **fHistTracksPtDist;            //!Track pt vs. distance form jet axis
  TH2                       **fHistClustersPtDist;          //!Cluster pt vs. distance form jet axis

  // Inclusive jets histograms
  THnSparse                  *fHistJetObservables;          //!Jet observables

  // TH2/TH3 versions
  TH3                       **fHistJetPtEtaPhi;             //!Jet Pt vs. Eta vs. Phi
  TH2                       **fHistJetPtArea;               //!Jet Pt vs. Area
  TH2                       **fHistJetPtEP;                 //!Jet Pt vs. event plane
  TH2                       **fHistJetPtNEF;                //!Jet Pt vs. neutral energy fraction
  TH2                       **fHistJetPtZ;                  //!Jet Pt vs. z
  TH2                       **fHistJetPtLeadingPartPt;      //!Jet Pt vs. leading particle pt
  TH3                       **fHistJetCorrPtEtaPhi;         //!Jet corrPt vs. Eta vs. Phi
  TH2                       **fHistJetCorrPtArea;           //!Jet corrPt vs. Area
  TH2                       **fHistJetCorrPtEP;             //!Jet corrPt vs. event plane
  TH2                       **fHistJetCorrPtNEF;            //!Jet corrPt vs. neutral energy fraction
  TH2                       **fHistJetCorrPtZ;              //!Jet corrPt vs. z
  TH2                       **fHistJetCorrPtLeadingPartPt;  //!Jet corrPt vs. leading particle pt
  TH2                       **fHistJetPtCorrPt;             //!Jet Pt vs. corrPt
  TH2                       **fHistJetPtMCPt;               //!Jet Pt vs. MCPt
  TH2                       **fHistJetMCPtCorrPt;           //!Jet MCPt vs. corrPt

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

  ClassDef(AliAnalysisTaskSAJF, 17) // jet analysis task
};
#endif
Commit	Line	Data
00514d01	1	#ifndef ALIANALYSISTASKSAJF_H
00514d01	2	#define ALIANALYSISTASKSAJF_H
25283b37	3
43032ce2	4	class TH2;
	5	class THnSparse;
	6
85d48773	7	#include <TH3F.h>
85d48773	8
9239b066	9	#include "AliAnalysisTaskEmcalJet.h"
25283b37	10
9239b066	11	class AliAnalysisTaskSAJF : public AliAnalysisTaskEmcalJet {
25283b37	12	public:
91f4b7c5	13
00514d01	14	AliAnalysisTaskSAJF();
00514d01	15	AliAnalysisTaskSAJF(const char *name);
43032ce2	16	virtual ~AliAnalysisTaskSAJF() {;}
25283b37	17
25283b37	18	void UserCreateOutputObjects();
05077f28	19
85d48773	20	void SetHistoType(Int_t t) { fHistoType = t; }
85d48773	21
25283b37	22	protected:
85d48773	23	void AllocateTHX();
	24	void AllocateTHnSparse();
	25
	26	Bool_t FillHistograms();
3fe08cdb	27	void FillJetHisto(Double_t cent, Double_t ep, Double_t eta, Double_t phi, Double_t pt, Double_t MCpt, Double_t corrpt, Double_t area,
3fe08cdb	28	Double_t NEF, Double_t z, Int_t n, Double_t leadingpt);
a487deae	29
85d48773	30	Int_t fHistoType; // histogram type (0=TH2, 1=THnSparse)
85d48773	31
68e03fc3	32	TH2 **fHistRejectionReason; //!Rejection reason vs. jet pt
b15e47c9	33	TH2 **fHistTracksJetPt; //!Track pt vs. jet pt
	34	TH2 **fHistClustersJetPt; //!Cluster pt vs. jet pt
	35	TH2 **fHistTracksPtDist; //!Track pt vs. distance form jet axis
	36	TH2 **fHistClustersPtDist; //!Cluster pt vs. distance form jet axis
	37
a487deae	38	// Inclusive jets histograms
b15e47c9	39	THnSparse *fHistJetObservables; //!Jet observables
85d48773	40
85d48773	41	// TH2/TH3 versions
b15e47c9	42	TH3 **fHistJetPtEtaPhi; //!Jet Pt vs. Eta vs. Phi
	43	TH2 **fHistJetPtArea; //!Jet Pt vs. Area
	44	TH2 **fHistJetPtEP; //!Jet Pt vs. event plane
	45	TH2 **fHistJetPtNEF; //!Jet Pt vs. neutral energy fraction
	46	TH2 **fHistJetPtZ; //!Jet Pt vs. z
	47	TH2 **fHistJetPtLeadingPartPt; //!Jet Pt vs. leading particle pt
	48	TH3 **fHistJetCorrPtEtaPhi; //!Jet corrPt vs. Eta vs. Phi
	49	TH2 **fHistJetCorrPtArea; //!Jet corrPt vs. Area
	50	TH2 **fHistJetCorrPtEP; //!Jet corrPt vs. event plane
	51	TH2 **fHistJetCorrPtNEF; //!Jet corrPt vs. neutral energy fraction
	52	TH2 **fHistJetCorrPtZ; //!Jet corrPt vs. z
	53	TH2 **fHistJetCorrPtLeadingPartPt; //!Jet corrPt vs. leading particle pt
	54	TH2 **fHistJetPtCorrPt; //!Jet Pt vs. corrPt
	55	TH2 **fHistJetPtMCPt; //!Jet Pt vs. MCPt
	56	TH2 **fHistJetMCPtCorrPt; //!Jet MCPt vs. corrPt
25283b37	57
25283b37	58	private:
00514d01	59	AliAnalysisTaskSAJF(const AliAnalysisTaskSAJF&); // not implemented
00514d01	60	AliAnalysisTaskSAJF &operator=(const AliAnalysisTaskSAJF&); // not implemented
25283b37	61
7cd832c7	62	ClassDef(AliAnalysisTaskSAJF, 17) // jet analysis task
25283b37	63	};
25283b37	64	#endif