[u/mrichter/AliRoot.git] / PWGJE / UserTasks / AliAnalysisTaskCheckSingleTrackJetRejection.h

#ifndef AliAnalysisTaskCheckSingleTrackJetRejection_cxx
#define AliAnalysisTaskCheckSingleTrackJetRejection_cxx

/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

class AliJetHeader;
class AliESDEvent;
class AliAODEvent;
class AliAODJet;
class AliGenPythiaEventHeader;
class AliCFManager;

class TList;
class TChain;
class TH2F;
class TH1F;
class TH3F;
class TProfile;


#include "AliAnalysisTaskSE.h"
#include  "THnSparse.h" // cannot forward declare ThnSparseF  
#include <TList.h>
#include <TH1F.h>
#include <TH2F.h>

class AliAnalysisTaskCheckSingleTrackJetRejection : public AliAnalysisTaskSE 
{
	public:
		AliAnalysisTaskCheckSingleTrackJetRejection();
		AliAnalysisTaskCheckSingleTrackJetRejection(const char *name);
		virtual ~AliAnalysisTaskCheckSingleTrackJetRejection() {;}

		// Implementation of interface methods
		virtual void    UserCreateOutputObjects();
		virtual void    Init();
		virtual Bool_t  Notify();
		virtual void    UserExec(Option_t *option);
		virtual void    Terminate(Option_t *);
		virtual void    SetDebug(Int_t debug = 0) {fDebug = debug;}
		virtual void    SetAlgorithm(TString jf="ANTIKT"){JFAlg=jf;}
		virtual void    SetRadius(Float_t radius=0.4){Radius=radius;}
		virtual void    SetFilterMask(UInt_t filter=256){Filtermask=filter;}
		virtual void    SetBackSubMode(Int_t backM=0){BackM=backM;}
		virtual void    SetTrackPtCut(Float_t tPtcut=0){TrackPtcut=tPtcut;}
		virtual void    SetSkipCone(Int_t skipCone=0){SkipCone=skipCone;}
		virtual void    SetMC(Bool_t ismc=kFALSE){IsMC=ismc;}
    virtual void    FinishTaskOutput();
    virtual Bool_t  PythiaInfoFromFile(const char* currFile,Float_t &fXsec,Float_t &fTrials);
    //virtual Float_t GetTotalEvents(const char* currFile);
    virtual Double_t DeltaPhi(Double_t phi1,Double_t phi2);
    virtual Double_t DeltaR(Double_t phi1,Double_t phi2,Double_t eta1,Double_t eta2);


		enum {kNPTBINS=10};

		// 0 all jets
		// 1 all jet in eta window
		// 2 all jets with partner
		// 3 all jets in eta window with partner
		// 4 all jets with partner in eta window
		enum {kStep0 = 0, kStep1, kStep2, kStep3, kStep4,kMaxStep};


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

    Bool_t  JetSelected(AliAODJet *jet);

		Bool_t        fUseAODInput; // read jets from input AOD
		Bool_t        fFillAOD;     // option to fill AOD branch
		TString       fNonStdFile;
		TString       fJetBranch;   // jet branch to read

		AliAODEvent     *fAODIn;      // AOD event
		AliAODEvent     *fAODOut;      // AOD event
		AliAODExtension *fAODExtension;
		TString         JFAlg;
		Float_t         Radius;
		UInt_t          Filtermask;
		Int_t           BackM;
		Float_t         TrackPtcut;
		Int_t           SkipCone;
		Bool_t          IsMC;

		TList        *fHistList; // Output list
		Float_t       fxsec;
		Float_t       ftrial;
		Float_t       fJetRecEtaWindow;       // eta window for rec jets
		Float_t       fMinJetPt;              // limits the jet p_T in addition to what already is done in the jet finder, this is important for jet matching for JF with lo threshold

		TProfile     *fH1Xsec;
		TH1F         *fH1Trials;
		TH1F         *fH1Events; 	

		TH1F         *fH1jetMCAKT04_pt         [6];
		TH2F         *fH2jetMCAKT04_Jetpt_maxpt;
		TH1F         *fH1jetAKT04_pt           [6];     
		TH2F         *fH2jetAKT04_Jetpt_maxpt  ;

		TH2F         *fH2jetMCAKT04_Eratio       [6];
		TH1F         *fH1jetMCAKT04_ptmatch        [6];


		ClassDef(AliAnalysisTaskCheckSingleTrackJetRejection, 13); // Analysis task for standard dijet analysis
};

#endif
Commit	Line	Data
7fac8669	1	#ifndef AliAnalysisTaskCheckSingleTrackJetRejection_cxx
	2	#define AliAnalysisTaskCheckSingleTrackJetRejection_cxx
	3
	4	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	5	* See cxx source for full Copyright notice */
	6
	7	class AliJetHeader;
	8	class AliESDEvent;
	9	class AliAODEvent;
	10	class AliAODJet;
	11	class AliGenPythiaEventHeader;
	12	class AliCFManager;
	13
	14	class TList;
	15	class TChain;
	16	class TH2F;
	17	class TH1F;
	18	class TH3F;
	19	class TProfile;
	20
	21
	22	#include "AliAnalysisTaskSE.h"
	23	#include "THnSparse.h" // cannot forward declare ThnSparseF
	24	#include <TList.h>
	25	#include <TH1F.h>
	26	#include <TH2F.h>
	27
	28	class AliAnalysisTaskCheckSingleTrackJetRejection : public AliAnalysisTaskSE
	29	{
	30	public:
	31	AliAnalysisTaskCheckSingleTrackJetRejection();
	32	AliAnalysisTaskCheckSingleTrackJetRejection(const char *name);
	33	virtual ~AliAnalysisTaskCheckSingleTrackJetRejection() {;}
	34
	35	// Implementation of interface methods
	36	virtual void UserCreateOutputObjects();
	37	virtual void Init();
	38	virtual Bool_t Notify();
	39	virtual void UserExec(Option_t *option);
	40	virtual void Terminate(Option_t *);
	41	virtual void SetDebug(Int_t debug = 0) {fDebug = debug;}
	42	virtual void SetAlgorithm(TString jf="ANTIKT"){JFAlg=jf;}
	43	virtual void SetRadius(Float_t radius=0.4){Radius=radius;}
	44	virtual void SetFilterMask(UInt_t filter=256){Filtermask=filter;}
	45	virtual void SetBackSubMode(Int_t backM=0){BackM=backM;}
	46	virtual void SetTrackPtCut(Float_t tPtcut=0){TrackPtcut=tPtcut;}
	47	virtual void SetSkipCone(Int_t skipCone=0){SkipCone=skipCone;}
	48	virtual void SetMC(Bool_t ismc=kFALSE){IsMC=ismc;}
	49	virtual void FinishTaskOutput();
	50	virtual Bool_t PythiaInfoFromFile(const char* currFile,Float_t &fXsec,Float_t &fTrials);
	51	//virtual Float_t GetTotalEvents(const char* currFile);
	52	virtual Double_t DeltaPhi(Double_t phi1,Double_t phi2);
	53	virtual Double_t DeltaR(Double_t phi1,Double_t phi2,Double_t eta1,Double_t eta2);
	54
	55
	56	enum {kNPTBINS=10};
	57
	58	// 0 all jets
	59	// 1 all jet in eta window
	60	// 2 all jets with partner
	61	// 3 all jets in eta window with partner
	62	// 4 all jets with partner in eta window
	63	enum {kStep0 = 0, kStep1, kStep2, kStep3, kStep4,kMaxStep};
	64
65
66
67	private:
68	AliAnalysisTaskCheckSingleTrackJetRejection(const AliAnalysisTaskCheckSingleTrackJetRejection &det); // not implemented
69	AliAnalysisTaskCheckSingleTrackJetRejection& operator=(const AliAnalysisTaskCheckSingleTrackJetRejection &det); // not implemented
70
71	Bool_t JetSelected(AliAODJet *jet);
72
73	Bool_t fUseAODInput; // read jets from input AOD
74	Bool_t fFillAOD; // option to fill AOD branch
75	TString fNonStdFile;
76	TString fJetBranch; // jet branch to read
77
78	AliAODEvent *fAODIn; // AOD event
79	AliAODEvent *fAODOut; // AOD event
80	AliAODExtension *fAODExtension;
81	TString JFAlg;
82	Float_t Radius;
83	UInt_t Filtermask;
84	Int_t BackM;
85	Float_t TrackPtcut;
86	Int_t SkipCone;
87	Bool_t IsMC;
88
89	TList *fHistList; // Output list
90	Float_t fxsec;
91	Float_t ftrial;
92	Float_t fJetRecEtaWindow; // eta window for rec jets
93	Float_t fMinJetPt; // limits the jet p_T in addition to what already is done in the jet finder, this is important for jet matching for JF with lo threshold
94
95	TProfile *fH1Xsec;
96	TH1F *fH1Trials;
97	TH1F *fH1Events;
98
99	TH1F *fH1jetMCAKT04_pt [6];
100	TH2F *fH2jetMCAKT04_Jetpt_maxpt;
101	TH1F *fH1jetAKT04_pt [6];
102	TH2F *fH2jetAKT04_Jetpt_maxpt ;
103
104	TH2F *fH2jetMCAKT04_Eratio [6];
105	TH1F *fH1jetMCAKT04_ptmatch [6];
106
107
108
109	ClassDef(AliAnalysisTaskCheckSingleTrackJetRejection, 13); // Analysis task for standard dijet analysis
110	};
111
112	#endif