3 #ifndef AlidNdEtaCorrectionTask_H
4 #define AlidNdEtaCorrectionTask_H
6 #include "AliAnalysisTask.h"
8 #include "AliPWG0Helper.h"
10 class AliESDtrackCuts;
12 class AlidNdEtaCorrection;
20 class AlidNdEtaCorrectionTask : public AliAnalysisTask {
22 AlidNdEtaCorrectionTask();
23 AlidNdEtaCorrectionTask(const char* opt);
24 virtual ~AlidNdEtaCorrectionTask();
26 virtual void ConnectInputData(Option_t *);
27 virtual void CreateOutputObjects();
28 virtual void Exec(Option_t*);
29 virtual void Terminate(Option_t *);
31 void SetTrackCuts(AliESDtrackCuts* cuts) { fEsdTrackCuts = cuts; }
32 void SetAnalysisMode(AliPWG0Helper::AnalysisMode mode) { fAnalysisMode = mode; }
33 void SetOnlyPrimaries(Bool_t flag = kTRUE) { fOnlyPrimaries = flag; }
34 void SetTrigger(AliPWG0Helper::Trigger trigger) { fTrigger = trigger; }
35 void SetFillPhi(Bool_t flag = kTRUE) { fFillPhi = flag; }
36 void SetDeltaPhiCut(Float_t cut) { fDeltaPhiCut = cut; }
38 void SetOption(const char* opt) { fOption = opt; }
41 Bool_t SignOK(TParticlePDG* particle);
43 AliESDEvent *fESD; //! ESD object
44 TList* fOutput; //! list send on output slot 0
46 TString fOption; // option string
47 AliPWG0Helper::AnalysisMode fAnalysisMode; // detector that is used for analysis
48 AliPWG0Helper::Trigger fTrigger; // trigger used in the analysis
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)
52 Int_t fSignMode; // if 0 process all particles, if +-1 process only particles with that sign
53 Bool_t fOnlyPrimaries; // only process primaries (syst. studies)
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
56 AliESDtrackCuts* fEsdTrackCuts; // Object containing the parameters of the esd track cuts
58 AlidNdEtaCorrection* fdNdEtaCorrection; //! contains the intermediate histograms (on each slave)
60 dNdEtaAnalysis* fdNdEtaAnalysisMC; //! analysis from MC (only triggered, vertex events)
61 dNdEtaAnalysis* fdNdEtaAnalysisESD; //! analysis from ESD (not yet corrected!)
64 TH1F* fPIDParticles; //! pid of primary particles
65 TH1F* fPIDTracks; //! pid of reconstructed tracks
67 TH2F* fVertexCorrelation; //! ESD z-vtx vs MC z-vtx
68 TH2F* fVertexCorrelationShift; //! (MC z-vtx - ESD z-vtx) vs MC z-vtx
69 TProfile* fVertexProfile; //! Profile of MC z-vtx - ESD z-vtx vs. MC z-vtx
70 TH1F* fVertexShift; //! (MC z-vtx - ESD z-vtx) in +- 10 cm
71 TH1F* fVertexShiftNorm; //! (MC z-vtx - ESD z-vtx) / (sigma_ESD-z-vtx) histogrammed
73 TH2F* fEtaCorrelation; //! ESD eta vs MC eta
74 TH2F* fEtaCorrelationShift; //! (MC eta - ESD eta) vs MC eta
75 TProfile* fEtaProfile; //! Profile of MC eta - ESD eta vs. MC eta
76 TH1F* fEtaResolution; //! MC eta - ESD eta in |eta| < 1
77 TH2F* fDeltaPhiCorrelation; //! delta phi ESD vs. MC
79 TH2F* fpTResolution; //! (MC pT - ESD pT) / MC pT vs. MC pT in |eta| < 0.9
81 AliESDtrackCuts* fEsdTrackCutsPrim; //! control histograms for primaries
82 AliESDtrackCuts* fEsdTrackCutsSec; //! control histograms for secondaries
84 // histograms for systematic studies (must be enabled with option)
86 TH1* fTemp1; //! temp histogram for quick study of variables
87 TH1* fTemp2; //! temp histogram for quick study of variables
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
93 TH2* fDeltaPhi[8]; //! delta phi of primaries, secondaries, other (= unclear cases)
95 TH2F* fEventStats; //! some stats on number of events, see CreateOutputObjects for a detailed definition
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
102 AlidNdEtaCorrectionTask(const AlidNdEtaCorrectionTask&);
103 AlidNdEtaCorrectionTask& operator=(const AlidNdEtaCorrectionTask&);
105 ClassDef(AlidNdEtaCorrectionTask, 1);