]>
Commit | Line | Data |
---|---|---|
2a436fa1 | 1 | #ifndef ALIANALYSISTASKJETRESPONSE_H |
2 | #define ALIANALYSISTASKJETRESPONSE_H | |
3 | ||
4 | class TH1F; | |
5 | class TH2F; | |
46465e39 | 6 | class TH3F; |
2a436fa1 | 7 | class AliESDEvent; |
8 | class AliAODEvent; | |
9 | ||
10 | #include "AliAnalysisTaskSE.h" | |
11 | #include "AliVEvent.h" | |
12 | ||
13 | class AliAnalysisTaskJetResponse : public AliAnalysisTaskSE { | |
14 | public: | |
15 | AliAnalysisTaskJetResponse(); | |
16 | AliAnalysisTaskJetResponse(const char *name); | |
17 | virtual ~AliAnalysisTaskJetResponse(); | |
18 | ||
19 | virtual void LocalInit() {Init();} | |
20 | virtual void Init(); | |
21 | virtual void UserCreateOutputObjects(); | |
22 | virtual void UserExec(Option_t *option); | |
23 | virtual void Terminate(const Option_t*); | |
24 | ||
25 | virtual AliVEvent::EOfflineTriggerTypes GetOfflineTrgMask() const { return fOfflineTrgMask; } | |
26 | virtual void GetBranchNames(TString &branch1, TString &branch2) const { branch1 = fJetBranchName[0]; branch2 = fJetBranchName[1]; } | |
46465e39 | 27 | virtual Int_t GetMinContribVtx() const { return fMinContribVtx; }; |
2a436fa1 | 28 | virtual Float_t GetVtxZMin() const { return fVtxZMin; } |
29 | virtual Float_t GetVtxZMax() const { return fVtxZMax; } | |
30 | virtual Int_t GetEvtClassMin() const { return fEvtClassMin; } | |
31 | virtual Int_t GetEvtClassMax() const { return fEvtClassMax; } | |
32 | virtual Float_t GetCentMin() const { return fCentMin; } | |
33 | virtual Float_t GetCentMax() const { return fCentMax; } | |
34 | virtual Float_t GetJetEtaMin() const { return fJetEtaMin; } | |
35 | virtual Float_t GetJetEtaMax() const { return fJetEtaMax; } | |
46465e39 | 36 | virtual Float_t GetJetPtMin() const { return fJetPtMin; } |
37 | virtual Float_t GetJetPtFractionMin() const { return fJetPtFractionMin; } | |
38 | virtual Int_t GetNMatchJets() const { return fNMatchJets; } | |
39 | //virtual Float_t GetJetDeltaEta() const { return fJetDeltaEta; } | |
40 | //virtual Float_t GetJetDeltaPhi() const { return fJetDeltaPhi; } | |
2a436fa1 | 41 | |
42 | virtual void SetBranchNames(const TString &branch1, const TString &branch2); | |
43 | virtual void SetOfflineTrgMask(AliVEvent::EOfflineTriggerTypes mask) { fOfflineTrgMask = mask; } | |
46465e39 | 44 | virtual void SetMinContribVtx(Int_t n) { fMinContribVtx = n; } |
2a436fa1 | 45 | virtual void SetVtxZMin(Float_t z) { fVtxZMin = z; } |
46 | virtual void SetVtxZMax(Float_t z) { fVtxZMax = z; } | |
47 | virtual void SetEvtClassMin(Int_t evtClass) { fEvtClassMin = evtClass; } | |
48 | virtual void SetEvtClassMax(Int_t evtClass) { fEvtClassMax = evtClass; } | |
49 | virtual void SetCentMin(Float_t cent) { fCentMin = cent; } | |
50 | virtual void SetCentMax(Float_t cent) { fCentMax = cent; } | |
51 | virtual void SetJetEtaMin(Float_t eta) { fJetEtaMin = eta; } | |
52 | virtual void SetJetEtaMax(Float_t eta) { fJetEtaMax = eta; } | |
46465e39 | 53 | virtual void SetJetPtMin(Float_t pt) { fJetPtMin = pt; } |
54 | virtual void SetJetPtFractionMin(Float_t pt) { fJetPtFractionMin = pt; } | |
55 | virtual void SetNMatchJets(Int_t n) { fNMatchJets = n; } | |
56 | //virtual void SetJetDeltaEta(Float_t eta) { fJetDeltaEta = eta; } | |
57 | //virtual void SetJetDeltaPhi(Float_t phi) { fJetDeltaPhi = phi; } | |
2a436fa1 | 58 | |
59 | private: | |
60 | // ESD/AOD events | |
61 | AliESDEvent *fESD; //! ESD object | |
62 | AliAODEvent *fAOD; //! AOD event | |
63 | ||
64 | // jets to compare | |
65 | TString fJetBranchName[2]; // name of jet branches to compare | |
66 | TList *fListJets[2]; //! jet lists | |
67 | ||
68 | // event selection | |
69 | AliVEvent::EOfflineTriggerTypes fOfflineTrgMask; // mask of offline triggers to accept | |
70 | Int_t fMinContribVtx; // minimum number of track contributors for primary vertex | |
71 | Float_t fVtxZMin; // lower bound on vertex z | |
72 | Float_t fVtxZMax; // upper bound on vertex z | |
73 | Int_t fEvtClassMin; // lower bound on event class | |
74 | Int_t fEvtClassMax; // upper bound on event class | |
75 | Float_t fCentMin; // lower bound on centrality | |
76 | Float_t fCentMax; // upper bound on centrality | |
77 | Float_t fJetEtaMin; // lower bound on eta for found jets | |
78 | Float_t fJetEtaMax; // upper bound on eta for found jets | |
46465e39 | 79 | Float_t fJetPtMin; // minimum jet pT |
80 | Float_t fJetPtFractionMin; // minimum fraction for positiv match of jets | |
81 | Int_t fNMatchJets; // maximal nb. of jets taken for matching | |
82 | //Float_t fJetDeltaEta; // max difference in eta to match leading jets | |
83 | //Float_t fJetDeltaPhi; // max difference in phi to match leading jets | |
2a436fa1 | 84 | |
85 | // output objects | |
86 | const Int_t fkNbranches; //! number of branches to be read | |
87 | const Int_t fkEvtClasses; //! number of event classes | |
88 | TList *fOutputList; //! output data container | |
46465e39 | 89 | TH1I *fHistEvtSelection; //! event selection statistic |
90 | TH1I *fHistEvtClass; //! event classes (from helper task) | |
91 | TH1F *fHistCentrality; //! centrality of the event | |
92 | TH1F **fHistPtJet; //! pt distribution of jets | |
93 | TH2F **fHistEtaPhiJet; //! eta-phi distribution of jets (before acceptance cuts) | |
94 | TH2F **fHistEtaPhiJetCut; //! eta-phi distribution of jets in eta acceptace per event class | |
95 | TH2F **fHistDeltaEtaDeltaPhiJet; //! delta eta vs. delta phi of matched jets (before acceptance cuts) | |
96 | TH2F **fHistDeltaEtaDeltaPhiJetCut; //! delta eta vs. delta phi of matched jets | |
97 | TH2F **fHistDeltaEtaDeltaPhiJetNOMatching; //! delta eta vs. delta phi of jets which do not match | |
98 | TH2F **fHistDeltaEtaEtaJet; //! delta eta vs. eta of matched jets per event class | |
99 | TH2F **fHistDeltaPtEtaJet; //! delta eta vs. eta of matched jets per event class | |
100 | TH2F **fHistPtFraction; //! fraction from embedded jet in reconstructed jet per event class | |
2a436fa1 | 101 | TH2F *fHistPtPtExtra; //! jet pt response |
46465e39 | 102 | TH2F **fHistPtResponse; //! jet pt response per event class |
a11972e9 | 103 | TH2F **fHistPtSmearing; //! emb-jet pt vs (emb+UE - emb) pt |
46465e39 | 104 | TH2F **fHistDeltaR; //! shift dR of jets vs (emb+UE - emb) pt |
a11972e9 | 105 | TH2F **fHistArea; //! area of jets vs (emb+UE - emb) pt |
46465e39 | 106 | TH2F **fHistDeltaArea; //! delta area of jets vs (emb+UE - emb) pt |
2a436fa1 | 107 | |
108 | AliAnalysisTaskJetResponse(const AliAnalysisTaskJetResponse&); // not implemented | |
109 | AliAnalysisTaskJetResponse& operator=(const AliAnalysisTaskJetResponse&); // not implemented | |
110 | ||
46465e39 | 111 | ClassDef(AliAnalysisTaskJetResponse, 2); |
2a436fa1 | 112 | }; |
113 | ||
114 | #endif |