1 #ifndef ALIANALYSISTASKJETCORRECTIONS_H
\r
2 #define ALIANALYSISTASKJETCORRECTIONS_H
\r
4 #include "AliAnalysisTaskSE.h"
\r
12 class AliGenPythiaEventHeader;
\r
26 class TLorentzVector;
\r
30 class AliAnalysisTaskJetCorrections : public AliAnalysisTaskSE
\r
33 AliAnalysisTaskJetCorrections();
\r
34 AliAnalysisTaskJetCorrections(const char * name);
\r
35 virtual ~AliAnalysisTaskJetCorrections() {;}
\r
37 //Implementation of interface methods
\r
38 virtual Bool_t Notify();
\r
39 virtual void UserCreateOutputObjects();
\r
40 virtual void Init();
\r
41 virtual void LocalInit() { Init(); };
\r
42 virtual void UserExec(Option_t * option);
\r
43 virtual void Terminate(Option_t * option);
\r
45 virtual void SetAODInput(Bool_t b){fUseAODInput = b;}
\r
47 virtual void SetBranchGen(char* c){fBranchGen = c;}
\r
48 virtual void SetBranchRec(char* c){fBranchRec = c;}
\r
50 virtual Double_t SetR(Double_t b){fR = b; return fR;}
\r
52 virtual void GetThrustAxis(TVector3 &n01, TVector3 * pTrack, Int_t &nTracks);
\r
54 AliAnalysisTaskJetCorrections(const AliAnalysisTaskJetCorrections&);
\r
55 AliAnalysisTaskJetCorrections& operator = (const AliAnalysisTaskJetCorrections&);
\r
57 enum {kMaxJets = 6};
\r
58 enum {kMaxEvents = 10};
\r
60 enum {kTracks = 1000};
\r
62 AliAODEvent *fAOD; // where we take the jets from can be input or output AOD
\r
64 TString fBranchRec; // AOD branch name for reconstructed
\r
65 TString fBranchGen; // AOD brnach for genereated
\r
67 Bool_t fUseAODInput;
\r
80 TH1F * fhEsumGenRest;
\r
81 TH1F * fhEsumRecRest;
\r
85 TH2F * fhE2E1vsEsumGen;
\r
86 TH2F * fhE2E1vsEsumRec;
\r
87 TH2F * fhE2E1vsE1Gen;
\r
88 TH2F * fhE2E1vsE1Rec;
\r
89 TH2F * fhE2E1vsdPhiGen;
\r
90 TH2F * fhE2E1vsdPhiRec;
\r
92 TH2F * fhTrackBalance2;
\r
93 TH2F * fhTrackBalance3;
\r
98 TProfile * fhECorrJet10[3];
\r
99 TProfile * fhECorrJet05[3];
\r
100 TProfile * fhECorrJet01[3];
\r
101 TProfile * fhECorrJet001[3];
\r
103 TProfile * fhdEvsErec10[3];
\r
104 TProfile * fhdEvsErec05[3];
\r
105 TProfile * fhdEvsErec01[3];
\r
106 TProfile * fhdEvsErec001[3];
\r
108 TH2F * fhdPhidEta10[3];
\r
109 TH2F * fhdPhidEta05[3];
\r
110 TH2F * fhdPhidEta01[3];
\r
111 TH2F * fhdPhidEta001[3];
\r
113 TH2F * fhdPhidEtaPt10[3];
\r
114 TH2F * fhdPhidEtaPt05[3];
\r
115 TH2F * fhdPhidEtaPt01[3];
\r
116 TH2F * fhdPhidEtaPt001[3];
\r
118 ClassDef(AliAnalysisTaskJetCorrections, 1)
\r