removed duplications in AlidNdPtHelper and AliPWG0Helper and small changes
[u/mrichter/AliRoot.git] / PWG0 / dNdPt / AlidNdPtAnalysis.h
1 #ifndef ALIDEDPTANALYSIS_H
2 #define ALIDEDPTANALYSIS_H
3
4 //------------------------------------------------------------------------------
5 // AlidNdPtAnalysis class. 
6 // 
7 // a. functionality:
8 // - fills analysis control histograms
9 // - fills generic correction matrices 
10 // - generates correction matrices 
11 //
12 // b. data members:
13 // - generic correction matrices
14 // - control histograms
15 //
16 // Author: J.Otwinowski 04/11/2008 
17 //------------------------------------------------------------------------------
18
19 class TProfile;
20 class TFolder;
21 class TObjArray;
22 class TString;
23
24 class AliESDtrackCuts;
25 class AliVertexerTracks;
26 class AliESD;
27 class AliESDfriend;
28 class AliESDfriendTrack;
29
30 class AlidNdPtHelper;
31
32 #include "THnSparse.h"
33 #include "AlidNdPt.h"
34
35 class AlidNdPtAnalysis : public AlidNdPt {
36 public :
37   AlidNdPtAnalysis(); 
38   AlidNdPtAnalysis(Char_t* name, Char_t* title);
39   ~AlidNdPtAnalysis();
40
41   // Init data members
42   virtual void Init();
43
44   // Process events
45   virtual void Process(AliESDEvent *const esdEvent=0, AliMCEvent *const mcEvent=0);
46
47   // Merge output objects (needed by PROOF) 
48   virtual Long64_t Merge(TCollection* list);
49
50   // Analyse output histograms 
51   virtual void Analyse();
52
53   // Export objects to folder
54   virtual TFolder *ExportToFolder(TObjArray * array=0);
55
56   // Get analysis folder
57   TFolder* GetAnalysisFolder() {return fAnalysisFolder;}
58
59   // Fill control histograms
60   void SetHistogramsOn(Bool_t histOn=kTRUE) {fHistogramsOn = histOn;}
61   Bool_t IsHistogramsOn() {return fHistogramsOn;}
62
63   // Create folder for analysed histograms
64   TFolder *CreateFolder(TString folder = "folderdNdPtAnalysis",TString title = "Analysed dNdPt histograms");
65
66   // Fill histograms
67   void FillHistograms(AliESDtrack *const esdTrack, AliStack *const stack, AlidNdPtHelper::TrackObject trackObj);
68   void FillHistograms(AliStack *const stack, Int_t label, AlidNdPtHelper::TrackObject trackObj);
69   void FillHistograms(TObjArray *const allChargedTracks,Int_t *const labelsAll,Int_t multAll,Int_t *const labelsAcc,Int_t multAcc,Int_t *const labelsRec,Int_t multRec);
70
71   // Getters
72   THnSparseF *GetEventMultCorrelationMatrix() {return fEventMultCorrelationMatrix;}
73   THnSparseF *GetTrackPtCorrelationMatrix()   {return fTrackPtCorrelationMatrix;}
74
75   //
76   THnSparseF *GetGenEventMatrix() {return fGenEventMatrix;}
77   THnSparseF *GetGenSDEventMatrix() {return fGenSDEventMatrix;}
78   THnSparseF *GetGenDDEventMatrix() {return fGenDDEventMatrix;}
79   THnSparseF *GetGenNDEventMatrix() {return fGenNDEventMatrix;}
80   THnSparseF *GetGenNSDEventMatrix() {return fGenNSDEventMatrix;}
81
82   THnSparseF *GetTriggerEventMatrix() {return fTriggerEventMatrix;}
83   THnSparseF *GetTriggerSDEventMatrix() {return fTriggerSDEventMatrix;}
84   THnSparseF *GetTriggerDDEventMatrix() {return fTriggerDDEventMatrix;}
85   THnSparseF *GetTriggerNDEventMatrix() {return fTriggerNDEventMatrix;}
86   THnSparseF *GetTriggerNSDEventMatrix() {return fTriggerNSDEventMatrix;}
87
88   THnSparseF *GetRecEventMatrix() {return fRecEventMatrix;}
89   THnSparseF *GetRecSDEventMatrix() {return fRecSDEventMatrix;}
90   THnSparseF *GetRecDDEventMatrix() {return fRecDDEventMatrix;}
91   THnSparseF *GetRecNDEventMatrix() {return fRecNDEventMatrix;}
92   THnSparseF *GetRecNSDEventMatrix() {return fRecNSDEventMatrix;}
93
94   // 
95   THnSparseF *GetGenTrackEventMatrix() {return fGenTrackEventMatrix;}
96   THnSparseF *GetGenTrackSDEventMatrix() {return fGenTrackSDEventMatrix;}
97   THnSparseF *GetGenTrackDDEventMatrix() {return fGenTrackDDEventMatrix;}
98   THnSparseF *GetGenTrackNDEventMatrix() {return fGenTrackNDEventMatrix;}
99   THnSparseF *GetGenTrackNSDEventMatrix() {return fGenTrackNSDEventMatrix;}
100
101   THnSparseF *GetTriggerTrackEventMatrix() {return fTriggerTrackEventMatrix;}
102   THnSparseF *GetTriggerTrackSDEventMatrix() {return fTriggerTrackSDEventMatrix;}
103   THnSparseF *GetTriggerTrackDDEventMatrix() {return fTriggerTrackDDEventMatrix;}
104   THnSparseF *GetTriggerTrackNDEventMatrix() {return fTriggerTrackNDEventMatrix;}
105   THnSparseF *GetTriggerTrackNSDEventMatrix() {return fTriggerTrackNSDEventMatrix;}
106
107   THnSparseF *GetRecTrackEventMatrix() {return fRecTrackEventMatrix;}
108   THnSparseF *GetRecTrackSDEventMatrix() {return fRecTrackSDEventMatrix;}
109   THnSparseF *GetRecTrackDDEventMatrix() {return fRecTrackDDEventMatrix;}
110   THnSparseF *GetRecTrackNDEventMatrix() {return fRecTrackNDEventMatrix;}
111   THnSparseF *GetRecTrackNSDEventMatrix() {return fRecTrackNSDEventMatrix;}
112
113   //
114   THnSparseF *GetGenPrimTrackMatrix() {return fGenPrimTrackMatrix;}
115   THnSparseF *GetRecPrimTrackMatrix() {return fRecPrimTrackMatrix;}
116
117   THnSparseF *GetRecTrackMatrix() {return fRecTrackMatrix;}
118   THnSparseF *GetRecSecTrackMatrix() {return fRecSecTrackMatrix;}
119   THnSparseF *GetRecMultTrackMatrix() {return fRecMultTrackMatrix;}
120
121   //
122   // control histograms
123   //
124   THnSparseF *GetMCEventHist1() {return fMCEventHist1;}
125   THnSparseF *GetRecEventHist1() {return fRecEventHist1;}
126   THnSparseF *GetRecEventHist2() {return fRecEventHist2;}
127   THnSparseF *GetRecMCEventHist1() {return fRecMCEventHist1;}
128   THnSparseF *GetRecMCEventHist2() {return fRecMCEventHist2;}
129   THnSparseF *GetRecMCEventHist3() {return fRecMCEventHist3;}
130
131   THnSparseF *GetMCTrackHist1(Int_t i) {return fMCTrackHist1[i];}
132   THnSparseF *GetMCPrimTrackHist1(Int_t i) {return fMCPrimTrackHist1[i];}
133   THnSparseF *GetMCSecTrackHist1(Int_t i) {return fMCSecTrackHist1[i];}
134
135   THnSparseF *GetRecTrackHist1(Int_t i) {return fRecTrackHist1[i];}
136   THnSparseF *GetRecTrackMultHist1(Int_t i) {return fRecTrackMultHist1[i];}
137
138   THnSparseF *GetRecMCTrackHist1() {return fRecMCTrackHist1;}
139   THnSparseF *GetMCMultRecTrackHist1() {return fMCMultRecTrackHist1;}
140
141 private:
142
143   // analysis folder 
144   TFolder *fAnalysisFolder; // folder for analysed histograms
145
146   // switch on/off filling of control histograms
147   Bool_t fHistogramsOn; 
148
149   // 
150   // correlation matrices (histograms)
151   //
152
153   // event rec. track vs true track multiplicity correlation matrix 
154   THnSparseF *fEventMultCorrelationMatrix; //-> mult:mult_true_tracks
155
156   // rec. track pt vs true track pt correlation matrix for given eta
157   THnSparseF *fTrackPtCorrelationMatrix; //-> Pt:mcPt:mcEta
158
159   //
160   // event level correction 
161   //
162
163   // all genertated
164   THnSparseF *fGenEventMatrix; //-> mcZv:multMB (inelastic)
165   THnSparseF *fGenSDEventMatrix; //-> mcZv:multMB (single diffractive)
166   THnSparseF *fGenDDEventMatrix; //-> mcZv:multMB (single diffractive)
167   THnSparseF *fGenNDEventMatrix; //-> mcZv:multMB (non diffractive)
168   THnSparseF *fGenNSDEventMatrix; //-> mcZv:multMB (non single diffractive)
169
170   // trigger bias corrections (fTriggerEventMatrix / fGenEventMatrix)
171   THnSparseF *fTriggerEventMatrix; //-> mcZv:multMB
172   THnSparseF *fTriggerSDEventMatrix; //-> mcZv:multMB
173   THnSparseF *fTriggerDDEventMatrix; //-> mcZv:multMB
174   THnSparseF *fTriggerNDEventMatrix; //-> mcZv:multMB
175   THnSparseF *fTriggerNSDEventMatrix; //-> mcZv:multMB
176
177   // event vertex rec. eff correction (fRecEventMatrix / fTriggerEventMatrix)
178   THnSparseF *fRecEventMatrix; //-> mcZv:multMB 
179   THnSparseF *fRecSDEventMatrix; //-> mcZv:multMB
180   THnSparseF *fRecDDEventMatrix; //-> mcZv:multMB
181   THnSparseF *fRecNDEventMatrix; //-> mcZv:multMB
182   THnSparseF *fRecNSDEventMatrix; //-> mcZv:multMB
183
184   //
185   // track-event level correction 
186   //
187
188   THnSparseF *fGenTrackEventMatrix; //-> mcZv:mcPt:mcEta
189   THnSparseF *fGenTrackSDEventMatrix; //-> mcZv:mcPt:mcEta
190   THnSparseF *fGenTrackDDEventMatrix; //-> mcZv:mcPt:mcEta
191   THnSparseF *fGenTrackNDEventMatrix; //-> mcZv:mcPt:mcEta
192   THnSparseF *fGenTrackNSDEventMatrix; //-> mcZv:mcPt:mcEta
193
194   // trigger bias corrections (fTriggerTrackEventMatrix / fGenTrackEventMatrix)
195   THnSparseF *fTriggerTrackEventMatrix; //-> mcZv:mcPt:mcEta
196   THnSparseF *fTriggerTrackSDEventMatrix; //-> mcZv:mcPt:mcEta
197   THnSparseF *fTriggerTrackDDEventMatrix; //-> mcZv:mcPt:mcEta
198   THnSparseF *fTriggerTrackNDEventMatrix; //-> mcZv:mcPt:mcEta
199   THnSparseF *fTriggerTrackNSDEventMatrix; //-> mcZv:mcPt:mcEta
200
201   // event vertex rec. corrections (fRecTrackEventMatrix / fTriggerTrackEventMatrix)
202   THnSparseF *fRecTrackEventMatrix; //-> mcZv:Pt:mcEta
203   THnSparseF *fRecTrackSDEventMatrix; //-> mcZv:Pt:mcEta
204   THnSparseF *fRecTrackDDEventMatrix; //-> mcZv:Pt:mcEta
205   THnSparseF *fRecTrackNDEventMatrix; //-> mcZv:Pt:mcEta
206   THnSparseF *fRecTrackNSDEventMatrix; //-> mcZv:Pt:mcEta
207
208   //
209   // track level correction 
210   //
211
212   // track rec. efficiency correction (fRecPrimTrackMatrix / fGenPrimTrackMatrix)
213   THnSparseF *fGenPrimTrackMatrix; //-> mcZv:mcPt:mcEta
214   THnSparseF *fRecPrimTrackMatrix; //-> mcZv:mcPt:mcEta
215
216   // secondary track contamination correction (fRecSecTrackMatrix / fRecTrackMatrix)
217   THnSparseF *fRecTrackMatrix;    //-> mcZv:mcPt:mcEta
218   THnSparseF *fRecSecTrackMatrix; //-> mcZv:mcPt:mcEta
219
220   // multiple rec. track corrections (fRecMultTrackMatrix / fRecTrackMatrix)
221   THnSparseF *fRecMultTrackMatrix; //-> mcZv:Pt:mcEta
222
223   //
224   // ESD and MC control analysis histograms
225   //
226
227   // THnSparse event histograms
228   THnSparseF *fMCEventHist1;  //-> mcXv:mcYv:mcZv
229   THnSparseF *fRecEventHist1; //-> Xv:Yv:Zv
230   THnSparseF *fRecEventHist2; //-> Zv:multMB:mult
231   THnSparseF *fRecMCEventHist1; //-> Xv-mcXv:Yv-mcYv:Zv-mcZv
232   THnSparseF *fRecMCEventHist2; //-> Xv-mcXv:Zv-mcZv:mult
233   THnSparseF *fRecMCEventHist3; //-> mult:EventType (ND, DD, SD)
234
235   // THnSparse track histograms
236   // [0] - after charged track selection, [1] - after acceptance cuts, [2] - after esd track cuts
237
238   THnSparseF *fMCTrackHist1[AlidNdPtHelper::kCutSteps];     //-> mcPt:mcEta:mcPhi
239   THnSparseF *fMCPrimTrackHist1[AlidNdPtHelper::kCutSteps]; //-> mcPt:mcEta:pid:mech:mother
240   THnSparseF *fMCSecTrackHist1[AlidNdPtHelper::kCutSteps];  //-> mcPt:mcEta:pid:mech:mother
241
242   THnSparseF *fRecTrackHist1[AlidNdPtHelper::kCutSteps];     //-> Pt:Eta:Phi
243   THnSparseF *fRecTrackMultHist1[AlidNdPtHelper::kCutSteps]; //-> Pt:mult
244
245   THnSparseF *fRecMCTrackHist1; //-> mcPt:mcEta:(Pt-mcPt)/mcPt:(Eta-mcEta)
246
247   //multple reconstructed tracks
248   THnSparseF *fMCMultRecTrackHist1; //-> mcPt:mcEta:pid
249
250   AlidNdPtAnalysis(const AlidNdPtAnalysis&); // not implemented
251   AlidNdPtAnalysis& operator=(const AlidNdPtAnalysis&); // not implemented
252
253   ClassDef(AlidNdPtAnalysis,1);
254 };
255
256 #endif