]>
Commit | Line | Data |
---|---|---|
95392a57 | 1 | #ifndef ALIANALYSISTASKJETCORRECTIONS_H\r |
2 | #define ALIANALYSISTASKJETCORRECTIONS_H\r | |
3 | \r | |
4 | #include "AliAnalysisTaskSE.h"\r | |
47dcf196 | 5 | \r |
6 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r | |
7 | * See cxx source for full Copyright notice */\r | |
393afcc4 | 8 | \r |
47dcf196 | 9 | // ***\r |
10 | // corrections to jet energy by sona.pochybova@cern.ch \r | |
11 | // ***\r | |
95392a57 | 12 | \r |
13 | class AliJetFinder;\r | |
14 | class AliESDEvent;\r | |
15 | class AliAODEvent;\r | |
16 | class AliAODJet;\r | |
17 | class AliGenPythiaEventHeader;\r | |
18 | class AliAODPid;\r | |
19 | \r | |
20 | class TList;\r | |
21 | class TArrayD;\r | |
22 | class TChain;\r | |
23 | class TH1;\r | |
24 | class TH2;\r | |
25 | class TH1F;\r | |
26 | class TH2F;\r | |
27 | class TH2I;\r | |
28 | class TH3D;\r | |
29 | class TTree;\r | |
30 | class TProfile;\r | |
31 | class TLorentzVector;\r | |
32 | class TVector3;\r | |
33 | class TParticle;\r | |
34 | \r | |
35 | class AliAnalysisTaskJetCorrections : public AliAnalysisTaskSE\r | |
36 | {\r | |
37 | public:\r | |
38 | AliAnalysisTaskJetCorrections();\r | |
39 | AliAnalysisTaskJetCorrections(const char * name);\r | |
40 | virtual ~AliAnalysisTaskJetCorrections() {;}\r | |
41 | \r | |
42 | //Implementation of interface methods\r | |
43 | virtual Bool_t Notify(); \r | |
44 | virtual void UserCreateOutputObjects();\r | |
45 | virtual void Init();\r | |
46 | virtual void LocalInit() { Init(); };\r | |
47 | virtual void UserExec(Option_t * option);\r | |
48 | virtual void Terminate(Option_t * option);\r | |
49 | \r | |
50 | virtual void SetAODInput(Bool_t b){fUseAODInput = b;}\r | |
51 | \r | |
393afcc4 | 52 | virtual void SetBranchGen(const char* c){fBranchGen = c;}\r |
53 | virtual void SetBranchRec(const char* c){fBranchRec = c;}\r | |
95392a57 | 54 | \r |
55 | virtual Double_t SetR(Double_t b){fR = b; return fR;} \r | |
56 | \r | |
95392a57 | 57 | private:\r |
58 | AliAnalysisTaskJetCorrections(const AliAnalysisTaskJetCorrections&);\r | |
59 | AliAnalysisTaskJetCorrections& operator = (const AliAnalysisTaskJetCorrections&);\r | |
60 | \r | |
61 | enum {kMaxJets = 6};\r | |
62 | enum {kMaxEvents = 10};\r | |
63 | enum {kJets = 3};\r | |
64 | enum {kTracks = 1000};\r | |
65 | \r | |
66 | AliAODEvent *fAOD; // where we take the jets from can be input or output AOD\r | |
67 | \r | |
68 | TString fBranchRec; // AOD branch name for reconstructed\r | |
69 | TString fBranchGen; // AOD brnach for genereated\r | |
70 | \r | |
393afcc4 | 71 | Bool_t fUseAODInput; // use aod input not output\r |
72 | Double_t fR; // radius\r | |
73 | TList * fList; // output list\r | |
74 | \r | |
6f3f79de | 75 | Int_t fGlobVar; // switch enabling checking out just one event\r |
76 | Double_t fXsection; // cross-section \r | |
77 | \r | |
78 | \r | |
79 | TH1F * fhEGen; // generated energy histogram\r | |
80 | TH1F * fhERec; // reconstructed energy histogram\r | |
81 | TH1F * fhEGenRest; // generated energy in the rest frame of three-jet events\r | |
82 | TH1F * fhERecRest; // generated energy in the rest frame of three-jet events\r | |
83 | TH1F * fhEsumGenRest; //generated summed energy in the rest frame of three-jet events\r | |
84 | TH1F * fhEsumRecRest; // reconstructed summed energy in the rest frame of three-jet events\r | |
85 | \r | |
86 | TH2F * fhE2vsE1Gen; // leading vs next-to leading jet energy, generated\r | |
87 | TH2F * fhE2vsE1Rec; // leading vs next-to leading jet energy, reconstruted\r | |
88 | TH2F * fhE2E1vsEsumGen; // E2,E1 diference vs. summed energy, generated\r | |
89 | TH2F * fhE2E1vsEsumRec; // E2, E1 difference vs summed jet energy, reconstructed\r | |
90 | TH2F * fhE2E1vsE1Gen; // E2, E1 difference vs E1, generated\r | |
91 | TH2F * fhE2E1vsE1Rec; // E2, E1 diference vs E1, reeconstructed\r | |
92 | TH2F * fhE2E1vsdPhiGen; // E2, E1 difference vs dPhi, generated\r | |
93 | TH2F * fhE2E1vsdPhiRec; // E2, E1 difference vs dPhi, reconstrted\r | |
94 | \r | |
95 | TH2F * fhTrackBalance2; // track balance in 2-jet events\r | |
96 | TH2F * fhTrackBalance3; // track balance in 3-jet events\r | |
97 | \r | |
98 | TH2F * fhEt1Et22; // Et1,2 in back-to-back cones in 2-jet events\r | |
99 | TH2F * fhEt1Et23; // Et1,2 in back-to-back cones in 3-jet events\r | |
100 | \r | |
101 | TProfile * fhECorrJet10[3]; // corr. factor for jets matched within dR = 1.\r | |
102 | TProfile * fhECorrJet05[3]; // corr. factor for jets matched within dR = .5\r | |
103 | TProfile * fhECorrJet01[3]; // corr. factor for jets matched within dR = .1\r | |
104 | TProfile * fhECorrJet001[3]; // corr. factor for jets matched within dR = .01\r | |
95392a57 | 105 | \r |
6f3f79de | 106 | TProfile * fhdEvsErec10[3]; // energy difference vs rec. energy, dR = 1.\r |
107 | TProfile * fhdEvsErec05[3]; // energy difference vs rec. energy, dR = .5\r | |
108 | TProfile * fhdEvsErec01[3]; // energy difference vs rec. energy, dR = .1\r | |
109 | TProfile * fhdEvsErec001[3]; // energy difference vs rec. energy, dR = .01\r | |
110 | \r | |
111 | TH2F * fhdPhidEta10[3]; // dPhi vs dEta of particles, dR = 1.\r | |
112 | TH2F * fhdPhidEta05[3]; // dPhi vs dEta of particles, dR = .5\r | |
113 | TH2F * fhdPhidEta01[3]; // dPhi vs dEta of particles, dR = .1\r | |
114 | TH2F * fhdPhidEta001[3]; // dPhi vs dEta of particles, dR = .01\r | |
115 | \r | |
116 | TH2F * fhdPhidEtaPt10[3]; // Pt weighted dPhi vs dEta of particles, dR = 1.\r | |
117 | TH2F * fhdPhidEtaPt05[3]; // Pt weighted dPhi vs dEta of particles, dR = .5\r | |
118 | TH2F * fhdPhidEtaPt01[3]; // Pt weighted dPhi vs dEta of particles, dR = .1\r | |
119 | TH2F * fhdPhidEtaPt001[3]; // Pt weighted dPhi vs dEta of particles, dR = .01\r | |
95392a57 | 120 | ClassDef(AliAnalysisTaskJetCorrections, 1)\r |
121 | };\r | |
122 | \r | |
123 | #endif\r |