]>
Commit | Line | Data |
---|---|---|
0f67a57c | 1 | /* $Id$ */ |
2 | ||
3 | #ifndef AlidNdEtaCorrectionTask_H | |
4 | #define AlidNdEtaCorrectionTask_H | |
5 | ||
6 | #include "AliAnalysisTask.h" | |
7 | #include <TString.h> | |
770a1f1d | 8 | #include "AliPWG0Helper.h" |
0f67a57c | 9 | |
10 | class AliESDtrackCuts; | |
11 | class dNdEtaAnalysis; | |
12 | class AlidNdEtaCorrection; | |
1c15d51a | 13 | class TH1; |
0f67a57c | 14 | class TH1F; |
15 | class AliESDEvent; | |
16 | class TParticlePDG; | |
54b096ef | 17 | class TH2F; |
18 | class TProfile; | |
0f67a57c | 19 | |
20 | class AlidNdEtaCorrectionTask : public AliAnalysisTask { | |
21 | public: | |
1c15d51a | 22 | AlidNdEtaCorrectionTask(); |
23 | AlidNdEtaCorrectionTask(const char* opt); | |
0f67a57c | 24 | virtual ~AlidNdEtaCorrectionTask(); |
25 | ||
26 | virtual void ConnectInputData(Option_t *); | |
27 | virtual void CreateOutputObjects(); | |
28 | virtual void Exec(Option_t*); | |
29 | virtual void Terminate(Option_t *); | |
30 | ||
31 | void SetTrackCuts(AliESDtrackCuts* cuts) { fEsdTrackCuts = cuts; } | |
770a1f1d | 32 | void SetAnalysisMode(AliPWG0Helper::AnalysisMode mode) { fAnalysisMode = mode; } |
3d7758c1 | 33 | void SetOnlyPrimaries(Bool_t flag = kTRUE) { fOnlyPrimaries = flag; } |
0fc41645 | 34 | void SetTrigger(AliPWG0Helper::Trigger trigger) { fTrigger = trigger; } |
69b09e3b | 35 | void SetFillPhi(Bool_t flag = kTRUE) { fFillPhi = flag; } |
36 | void SetDeltaPhiCut(Float_t cut) { fDeltaPhiCut = cut; } | |
0fc41645 | 37 | |
38 | void SetOption(const char* opt) { fOption = opt; } | |
0f67a57c | 39 | |
40 | protected: | |
41 | Bool_t SignOK(TParticlePDG* particle); | |
42 | ||
43 | AliESDEvent *fESD; //! ESD object | |
44 | TList* fOutput; //! list send on output slot 0 | |
45 | ||
46 | TString fOption; // option string | |
770a1f1d | 47 | AliPWG0Helper::AnalysisMode fAnalysisMode; // detector that is used for analysis |
0fc41645 | 48 | AliPWG0Helper::Trigger fTrigger; // trigger used in the analysis |
69b09e3b | 49 | Bool_t fFillPhi; // if true phi is filled as 3rd coordinate in all maps |
50 | Float_t fDeltaPhiCut; // cut in delta phi (only SPD) | |
51 | ||
0f67a57c | 52 | Int_t fSignMode; // if 0 process all particles, if +-1 process only particles with that sign |
3d7758c1 | 53 | Bool_t fOnlyPrimaries; // only process primaries (syst. studies) |
69b09e3b | 54 | Int_t fStatError; // statistical error evaluation: if set to 1 we only count unique primaries (binomial errors are valid), for 2 all the rest |
0f67a57c | 55 | |
56 | AliESDtrackCuts* fEsdTrackCuts; // Object containing the parameters of the esd track cuts | |
57 | ||
58 | AlidNdEtaCorrection* fdNdEtaCorrection; //! contains the intermediate histograms (on each slave) | |
59 | ||
60 | dNdEtaAnalysis* fdNdEtaAnalysisMC; //! analysis from MC (only triggered, vertex events) | |
61 | dNdEtaAnalysis* fdNdEtaAnalysisESD; //! analysis from ESD (not yet corrected!) | |
62 | ||
63 | // control histograms | |
64 | TH1F* fPIDParticles; //! pid of primary particles | |
65 | TH1F* fPIDTracks; //! pid of reconstructed tracks | |
1c15d51a | 66 | |
54b096ef | 67 | TH2F* fVertexCorrelation; //! ESD z-vtx vs MC z-vtx |
51f6de65 | 68 | TH2F* fVertexCorrelationShift; //! (MC z-vtx - ESD z-vtx) vs MC z-vtx |
54b096ef | 69 | TProfile* fVertexProfile; //! Profile of MC z-vtx - ESD z-vtx vs. MC z-vtx |
1c15d51a | 70 | TH1F* fVertexShift; //! (MC z-vtx - ESD z-vtx) in +- 10 cm |
54b096ef | 71 | TH1F* fVertexShiftNorm; //! (MC z-vtx - ESD z-vtx) / (sigma_ESD-z-vtx) histogrammed |
0f67a57c | 72 | |
3d7758c1 | 73 | TH2F* fEtaCorrelation; //! ESD eta vs MC eta |
51f6de65 | 74 | TH2F* fEtaCorrelationShift; //! (MC eta - ESD eta) vs MC eta |
3d7758c1 | 75 | TProfile* fEtaProfile; //! Profile of MC eta - ESD eta vs. MC eta |
1c15d51a | 76 | TH1F* fEtaResolution; //! MC eta - ESD eta in |eta| < 1 |
69b09e3b | 77 | TH2F* fDeltaPhiCorrelation; //! delta phi ESD vs. MC |
1c15d51a | 78 | |
69b09e3b | 79 | TH2F* fpTResolution; //! (MC pT - ESD pT) / MC pT vs. MC pT in |eta| < 0.9 |
1c15d51a | 80 | |
81 | AliESDtrackCuts* fEsdTrackCutsPrim; //! control histograms for primaries | |
82 | AliESDtrackCuts* fEsdTrackCutsSec; //! control histograms for secondaries | |
3d7758c1 | 83 | |
96fcc8a7 | 84 | // histograms for systematic studies (must be enabled with option) |
85 | ||
1c15d51a | 86 | TH1* fTemp1; //! temp histogram for quick study of variables |
87 | TH1* fTemp2; //! temp histogram for quick study of variables | |
96fcc8a7 | 88 | |
4ef701a0 | 89 | TH1F* fMultAll; //! primary particles in |eta| < 1 and pT > 0.2 in all events |
90 | TH1F* fMultTr; //! primary particles in |eta| < 1 and pT > 0.2 in triggered events | |
91 | TH1F* fMultVtx; //! primary particles in |eta| < 1 and pT > 0.2 in triggered events with vertex | |
92 | ||
69b09e3b | 93 | TH2* fDeltaPhi[8]; //! delta phi of primaries, secondaries, other (= unclear cases) |
4ef701a0 | 94 | |
69b09e3b | 95 | TH2F* fEventStats; //! some stats on number of events, see CreateOutputObjects for a detailed definition |
50ec344d | 96 | |
69b09e3b | 97 | AlidNdEtaCorrection* fdNdEtaCorrectionSpecial[4]; //! correction maps used for systematic studies, may contain: |
98 | // for specific process type (ND, SD, DD), enable with option: process-types | |
99 | // for particle species (pi, K, p, rest), enable with: particle-species | |
96fcc8a7 | 100 | |
0f67a57c | 101 | private: |
102 | AlidNdEtaCorrectionTask(const AlidNdEtaCorrectionTask&); | |
103 | AlidNdEtaCorrectionTask& operator=(const AlidNdEtaCorrectionTask&); | |
104 | ||
105 | ClassDef(AlidNdEtaCorrectionTask, 1); | |
106 | }; | |
107 | ||
108 | #endif |