Updates for train running

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

#ifndef AliAnalysisTaskJetHadronCorrelation_cxx\r
#define AliAnalysisTaskJetHadronCorrelation_cxx\r
\r
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r
 * See cxx source for full Copyright notice                               */\r
\r
class AliJetHeader;\r
class AliESDEvent;\r
class AliAODEvent;\r
class AliAODJet;\r
class AliGenPythiaEventHeader;\r
class AliCFManager;\r
\r
class TList;\r
class TChain;\r
class TH2F;\r
class TH1F;\r
class TH3F;\r
class TProfile;\r
\r
\r
#include "AliAnalysisTaskSE.h"\r
#include  "THnSparse.h" // cannot forward declare ThnSparseF  \r
#include <TList.h>\r
#include <TH1F.h>\r
#include <TH2F.h>\r
\r
class AliAnalysisTaskJetHadronCorrelation : public AliAnalysisTaskSE \r
{\r
	public:\r
		AliAnalysisTaskJetHadronCorrelation();\r
		AliAnalysisTaskJetHadronCorrelation(const char *name);\r
		virtual ~AliAnalysisTaskJetHadronCorrelation() {;}\r
\r
		// Implementation of interface methods\r
		virtual void   UserCreateOutputObjects();\r
		virtual void   Init();\r
		virtual Bool_t Notify();\r
		virtual void   UserExec(Option_t *option);\r
		virtual void   Terminate(Option_t *);\r
		virtual void   SetDebug(Int_t debug = 0) {fDebug = debug;}\r
		virtual void   SetAlgorithm(TString jf="ANTIKT"){JFAlg=jf;}\r
		virtual void   SetRadius(Float_t radius=0.4){Radius=radius;}\r
		virtual void   SetFilterMask(UInt_t filter=256){Filtermask=filter;}\r
		virtual void   SetBackSubMode(Int_t backM=0){BackM=backM;}\r
		virtual void   SetTrackPtCut(Float_t tPtcut=0){TrackPtcut=tPtcut;}\r
		virtual void   SetSkipCone(Int_t skipCone=0){SkipCone=skipCone;}\r
		virtual void   SetMC(Bool_t ismc=true){IsMC=ismc;}\r
		virtual void   FinishTaskOutput();\r
\r
		enum {kNPTBINS=10};\r
\r
		// 0 all jets\r
		// 1 all jet in eta window\r
		// 2 all jets with partner\r
		// 3 all jets in eta window with partner\r
		// 4 all jets with partner in eta window\r
		enum {kStep0 = 0, kStep1, kStep2, kStep3, kStep4,kMaxStep};\r
\r
\r
\r
	private:\r
		AliAnalysisTaskJetHadronCorrelation(const AliAnalysisTaskJetHadronCorrelation &det); // not implemented\r
		AliAnalysisTaskJetHadronCorrelation& operator=(const AliAnalysisTaskJetHadronCorrelation &det); // not implemented\r
\r
		Bool_t   JetSelected(AliAODJet *jet);\r
		Double_t DeltaPhi(Double_t phi1,Double_t phi2);\r
\r
		Bool_t        fUseAODInput; // read jets from input AOD\r
		Bool_t        fFillAOD;     // option to fill AOD branch\r
		TString       fJetBranch;   // jet branch to read\r
		TString       fNonStdFile;\r
\r
		AliAODEvent   *fAODIn;       // AOD event\r
		AliAODEvent   *fAODOut;      // AOD event\r
		AliAODExtension *fAODExtension;\r
		TString         JFAlg;\r
		Float_t         Radius;\r
		UInt_t          Filtermask;\r
		Int_t           BackM;\r
		Float_t         TrackPtcut;\r
		Int_t           SkipCone;\r
		Bool_t          IsMC;\r
\r
\r
		Float_t       fxsec;\r
		Float_t       ftrial;\r
		Float_t       fJetRecEtaWindow;       // eta window for rec jets\r
		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\r
\r
		TList        *fHistList; // Output list\r
    Int_t        fIfiles;//!count no. of files\r
\r
\r
		TH1F         *fH1Events;\r
		TProfile     *fH1Xsec;\r
		TH1F         *fH1Trials;\r
\r
		TH1F         *fH1Track_pt          ;\r
		TH1F         *fH1Track_phi         ;\r
		TH1F         *fH1Track_eta         ;\r
		TH1F         *fH1Jet_pt            ;\r
		TH1F         *fH1Jet_phi           ;\r
		TH1F         *fH1Jet_eta           ;\r
		TH1F         *fH1leadJet_pt        ;\r
		TH1F         *fH1leadJet_pt_dijet  ;\r
		TH1F         *fH1subJet_pt_dijet   ;\r
		TH2F         *fH2JetsJet_dphi      ;\r
		TH2F         *fH2JetsJet_deta      ;\r
		TH2F         *fH2JetsJet_Aj        ;\r
		TH2F         *fH2JetsJet_pt        ;\r
		TH1F         *fH1JetMC_pt          ;\r
		TH1F         *fH1leadJetMC_pt      ;\r
		TH1F         *fH1leadJetMC_pt_dijet;\r
		TH1F         *fH1subJetMC_pt_dijet ;\r
		TH2F         *fH2JetsJetMC_dphi    ;\r
		TH2F         *fH2JetsJetMC_deta    ;\r
		TH2F         *fH2JetsJetMC_Aj      ;\r
		TH2F         *fH2JetsJetMC_pt      ;\r
\r
		TH2F         *fH2Mult_Mtrack      ;\r
		TH2F         *fH2Mult_Mlead       ;\r
		TH2F         *fH2Mult_Mjet        ;\r
		TH2F         *fH2Mult_Njet        ;\r
		TH2F         *fH2Mult_Aj          ;\r
		TH2F         *fH2Mlead_Aj         ;\r
		TH2F         *fH2Jet_pt_Mlead     ;\r
		TH2F         *fH2Jet_pt_Munder    ;\r
\r
		TH1F         *fH1ndiJ_ediv                     [5];\r
		TH1F         *fH1Aj                            [5];\r
		TH1F         *fH1Mlead                         [5];\r
\r
		TH1F         *fH1leadJetMC_dphiResolution      [5];\r
		TH1F         *fH1subJetMC_dphiResolution       [5];\r
		TH1F         *fH1leadJetMC_Efficiency          [5];\r
		TH1F         *fH1subJetMC_Efficiency           [5];\r
\r
		TH1F         *fH1JetHadron_dphi_ediv           [5][5];\r
		TH1F         *fH1JetHadron_dphi_tptweight_ediv [5][5];\r
		TH1F         *fH1JetHadron_dphi_tJptweight_ediv[5][5];\r
\r
		TH1F         *fH1ndiJ_2040Mlead                       [3];\r
		TH1F         *fH1ndiJ_2040Aj                          [3];\r
		TH1F         *fH1JetHadron_dphi_tptweight2040_Mleaddep[3][5];\r
		TH1F         *fH1JetHadron_dphi_tptweight2040_Ajdep   [3][5];\r
\r
		ClassDef(AliAnalysisTaskJetHadronCorrelation, 15); // Analysis task for JetHadronCorrelation\r
};\r
\r
#endif\r