]>
Commit | Line | Data |
---|---|---|
2a436fa1 | 1 | #ifndef ALIANALYSISTASKJETRESPONSE_H |
2 | #define ALIANALYSISTASKJETRESPONSE_H | |
3 | ||
7a88ef30 | 4 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * |
5 | * See cxx source for full Copyright notice */ | |
6 | ||
7 | // | |
8 | // task compares jets in two branches, | |
9 | // written for analysis of jet embedding in HI events | |
10 | // | |
11 | // newer class: AliAnalysisTaskJetResponseV2 | |
12 | // | |
13 | ||
2a436fa1 | 14 | class TH1F; |
15 | class TH2F; | |
46465e39 | 16 | class TH3F; |
2a436fa1 | 17 | class AliESDEvent; |
18 | class AliAODEvent; | |
19 | ||
20 | #include "AliAnalysisTaskSE.h" | |
21 | #include "AliVEvent.h" | |
22 | ||
23 | class AliAnalysisTaskJetResponse : public AliAnalysisTaskSE { | |
24 | public: | |
25 | AliAnalysisTaskJetResponse(); | |
26 | AliAnalysisTaskJetResponse(const char *name); | |
27 | virtual ~AliAnalysisTaskJetResponse(); | |
28 | ||
29 | virtual void LocalInit() {Init();} | |
30 | virtual void Init(); | |
31 | virtual void UserCreateOutputObjects(); | |
32 | virtual void UserExec(Option_t *option); | |
33 | virtual void Terminate(const Option_t*); | |
34 | ||
35 | virtual AliVEvent::EOfflineTriggerTypes GetOfflineTrgMask() const { return fOfflineTrgMask; } | |
36 | virtual void GetBranchNames(TString &branch1, TString &branch2) const { branch1 = fJetBranchName[0]; branch2 = fJetBranchName[1]; } | |
46465e39 | 37 | virtual Int_t GetMinContribVtx() const { return fMinContribVtx; }; |
2a436fa1 | 38 | virtual Float_t GetVtxZMin() const { return fVtxZMin; } |
39 | virtual Float_t GetVtxZMax() const { return fVtxZMax; } | |
40 | virtual Int_t GetEvtClassMin() const { return fEvtClassMin; } | |
41 | virtual Int_t GetEvtClassMax() const { return fEvtClassMax; } | |
42 | virtual Float_t GetCentMin() const { return fCentMin; } | |
43 | virtual Float_t GetCentMax() const { return fCentMax; } | |
3c6a60f7 | 44 | virtual Int_t GetNInputTracksMin() const { return fNInputTracksMin; } |
45 | virtual Int_t GetNInputTracksMax() const { return fNInputTracksMax; } | |
46 | virtual Float_t GetReactionPlaneBin() const { return fReactionPlaneBin; } | |
2a436fa1 | 47 | virtual Float_t GetJetEtaMin() const { return fJetEtaMin; } |
48 | virtual Float_t GetJetEtaMax() const { return fJetEtaMax; } | |
46465e39 | 49 | virtual Float_t GetJetPtMin() const { return fJetPtMin; } |
50 | virtual Float_t GetJetPtFractionMin() const { return fJetPtFractionMin; } | |
51 | virtual Int_t GetNMatchJets() const { return fNMatchJets; } | |
3c6a60f7 | 52 | virtual Int_t GetEventClassMode() const { return fEvtClassMode; } |
2a436fa1 | 53 | |
54 | virtual void SetBranchNames(const TString &branch1, const TString &branch2); | |
55 | virtual void SetOfflineTrgMask(AliVEvent::EOfflineTriggerTypes mask) { fOfflineTrgMask = mask; } | |
46465e39 | 56 | virtual void SetMinContribVtx(Int_t n) { fMinContribVtx = n; } |
2a436fa1 | 57 | virtual void SetVtxZMin(Float_t z) { fVtxZMin = z; } |
58 | virtual void SetVtxZMax(Float_t z) { fVtxZMax = z; } | |
59 | virtual void SetEvtClassMin(Int_t evtClass) { fEvtClassMin = evtClass; } | |
60 | virtual void SetEvtClassMax(Int_t evtClass) { fEvtClassMax = evtClass; } | |
61 | virtual void SetCentMin(Float_t cent) { fCentMin = cent; } | |
62 | virtual void SetCentMax(Float_t cent) { fCentMax = cent; } | |
3c6a60f7 | 63 | virtual void SetNInputTracksMin(Int_t nTr) { fNInputTracksMin = nTr; } |
64 | virtual void SetNInputTracksMax(Int_t nTr) { fNInputTracksMax = nTr; } | |
65 | virtual void SetReactionPlaneBin(Int_t rpBin) { fReactionPlaneBin = rpBin; } | |
2a436fa1 | 66 | virtual void SetJetEtaMin(Float_t eta) { fJetEtaMin = eta; } |
67 | virtual void SetJetEtaMax(Float_t eta) { fJetEtaMax = eta; } | |
46465e39 | 68 | virtual void SetJetPtMin(Float_t pt) { fJetPtMin = pt; } |
69 | virtual void SetJetPtFractionMin(Float_t pt) { fJetPtFractionMin = pt; } | |
70 | virtual void SetNMatchJets(Int_t n) { fNMatchJets = n; } | |
3c6a60f7 | 71 | virtual void SetEventClassMode(Int_t mode) { fEvtClassMode = mode; } |
2a436fa1 | 72 | |
73 | private: | |
74 | // ESD/AOD events | |
75 | AliESDEvent *fESD; //! ESD object | |
76 | AliAODEvent *fAOD; //! AOD event | |
77 | ||
78 | // jets to compare | |
79 | TString fJetBranchName[2]; // name of jet branches to compare | |
80 | TList *fListJets[2]; //! jet lists | |
81 | ||
82 | // event selection | |
83 | AliVEvent::EOfflineTriggerTypes fOfflineTrgMask; // mask of offline triggers to accept | |
84 | Int_t fMinContribVtx; // minimum number of track contributors for primary vertex | |
85 | Float_t fVtxZMin; // lower bound on vertex z | |
86 | Float_t fVtxZMax; // upper bound on vertex z | |
87 | Int_t fEvtClassMin; // lower bound on event class | |
88 | Int_t fEvtClassMax; // upper bound on event class | |
89 | Float_t fCentMin; // lower bound on centrality | |
90 | Float_t fCentMax; // upper bound on centrality | |
3c6a60f7 | 91 | Int_t fNInputTracksMin; // lower bound of nb. of input tracks |
92 | Int_t fNInputTracksMax; // upper bound of nb. of input tracks | |
93 | Float_t fReactionPlaneBin; // reaction plane bin | |
2a436fa1 | 94 | Float_t fJetEtaMin; // lower bound on eta for found jets |
95 | Float_t fJetEtaMax; // upper bound on eta for found jets | |
46465e39 | 96 | Float_t fJetPtMin; // minimum jet pT |
97 | Float_t fJetPtFractionMin; // minimum fraction for positiv match of jets | |
98 | Int_t fNMatchJets; // maximal nb. of jets taken for matching | |
3c6a60f7 | 99 | |
100 | Int_t fEvtClassMode; // event class mode; 0: centrality (default), 1: multiplicity | |
2a436fa1 | 101 | |
102 | // output objects | |
103 | const Int_t fkNbranches; //! number of branches to be read | |
104 | const Int_t fkEvtClasses; //! number of event classes | |
105 | TList *fOutputList; //! output data container | |
46465e39 | 106 | TH1I *fHistEvtSelection; //! event selection statistic |
107 | TH1I *fHistEvtClass; //! event classes (from helper task) | |
108 | TH1F *fHistCentrality; //! centrality of the event | |
3c6a60f7 | 109 | TH1F *fHistNInputTracks; //! nb. of input tracks in the event |
110 | TH2F *fHistCentVsTracks; //! centrality vs. nb. of input tracks of the event | |
111 | TH1F *fHistReactionPlane; //! reaction plane of the event | |
112 | TH1F *fHistReactionPlaneWrtJets; //! reaction plane of the event wrt the jet | |
46465e39 | 113 | TH1F **fHistPtJet; //! pt distribution of jets |
114 | TH2F **fHistEtaPhiJet; //! eta-phi distribution of jets (before acceptance cuts) | |
115 | TH2F **fHistEtaPhiJetCut; //! eta-phi distribution of jets in eta acceptace per event class | |
116 | TH2F **fHistDeltaEtaDeltaPhiJet; //! delta eta vs. delta phi of matched jets (before acceptance cuts) | |
117 | TH2F **fHistDeltaEtaDeltaPhiJetCut; //! delta eta vs. delta phi of matched jets | |
3c6a60f7 | 118 | //TH2F **fHistDeltaEtaDeltaPhiJetNOMatching; //! delta eta vs. delta phi of jets which do not match |
46465e39 | 119 | TH2F **fHistDeltaEtaEtaJet; //! delta eta vs. eta of matched jets per event class |
120 | TH2F **fHistDeltaPtEtaJet; //! delta eta vs. eta of matched jets per event class | |
121 | TH2F **fHistPtFraction; //! fraction from embedded jet in reconstructed jet per event class | |
46465e39 | 122 | TH2F **fHistPtResponse; //! jet pt response per event class |
a11972e9 | 123 | TH2F **fHistPtSmearing; //! emb-jet pt vs (emb+UE - emb) pt |
46465e39 | 124 | TH2F **fHistDeltaR; //! shift dR of jets vs (emb+UE - emb) pt |
3c6a60f7 | 125 | TH3F **fHistArea; //! area of jets vs (emb+UE - emb) pt |
126 | //TH2F **fHistJetPtArea; //! area of jets vs (emb+UE - emb) pt | |
46465e39 | 127 | TH2F **fHistDeltaArea; //! delta area of jets vs (emb+UE - emb) pt |
3c6a60f7 | 128 | TH2F **fHistJetPtDeltaArea; //! delta area of jets vs emb pt |
2a436fa1 | 129 | |
130 | AliAnalysisTaskJetResponse(const AliAnalysisTaskJetResponse&); // not implemented | |
131 | AliAnalysisTaskJetResponse& operator=(const AliAnalysisTaskJetResponse&); // not implemented | |
132 | ||
3c6a60f7 | 133 | ClassDef(AliAnalysisTaskJetResponse, 3); |
2a436fa1 | 134 | }; |
135 | ||
136 | #endif |