]>
Commit | Line | Data |
---|---|---|
c95f49f3 | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
16 | //============================================================================== | |
17 | // AliHMPIDAnalysysTask - Class representing a basic analysis tool of HMPID data | |
18 | // A set of histograms is created. | |
19 | //============================================================================== | |
20 | // | |
21 | // By means of AliHMPIDPerformanceTask.C macro it is possible to use this class | |
22 | // to perform the analysis on local data, on data on alien using local machine | |
23 | // and on CAF. | |
24 | ||
25 | #ifndef AliHMPIDPERFORMANCETASK_CXX | |
26 | #define AliHMPIDPERFORMANCETASK_CXX | |
27 | ||
28 | #include "TObject.h" | |
29 | #include "TH1.h" | |
30 | #include "TH2.h" | |
31 | #include "TFile.h" | |
32 | #include "TCanvas.h" | |
33 | #include "TTree.h" | |
34 | #include "TGraphErrors.h" | |
35 | #include "AliAnalysisManager.h" | |
36 | #include "AliESDInputHandler.h" | |
37 | #include "AliMCEventHandler.h" | |
38 | #include "AliMCEvent.h" | |
39 | #include "AliESDtrack.h" | |
40 | #include "AliESDEvent.h" | |
41 | #include "AliESDVertex.h" | |
42 | #include "AliESDfriend.h" | |
43 | #include "AliHMPIDCluster.h" | |
44 | #include "AliESDfriendTrack.h" | |
45 | #include "AliPID.h" | |
46 | #include "AliLog.h" | |
47 | #include "AliESDtrackCuts.h" | |
48 | #include "AliHMPIDPerformanceTask.h" | |
49 | #include "AliCentrality.h" | |
50 | ||
51 | ClassImp(AliHMPIDPerformanceTask) | |
52 | ||
53 | //__________________________________________________________________________ | |
54 | AliHMPIDPerformanceTask::AliHMPIDPerformanceTask() : | |
55 | fESD(0x0),fESDfriend(0x0),fEsdTrackCuts(0x0),fClu(0x0),fMC(0x0),fEsdVtx(0x0),fCentrality(0x0),fEsdTrack(0x0),fEsdFriendTrack(0x0), | |
56 | fCalibObject(0x0), fUseMC(kTRUE), fHmpHistList(0x0), fHmpNevents(0x0),fHmpNevPerTrigClass(0x0), | |
57 | fGlobalEventNumber(0), | |
58 | //tree stuff | |
59 | fvType(0), | |
60 | fvFiredTriggerClasses(0x0), fvRunNumber(0), fvBunchCrossNumber(0),fvOrbitNumber(0), fvPeriodNumber(0), fvMagField(0), | |
61 | fvVertexX(0),fvVertexY(0),fvVertexZ(0),fvVertexNContributors(0), fvPesd(0),fvPhmpMag(0),fvCentrality(0), | |
62 | fvAcceptedTracks(0),fvRefMultTpc(0), fvHmpChi2(0),fvHmpCluIndx(0),fvHmpMipX(0),fvHmpMipY(0),fvHmpMipQ(0),fvHmpMipNPhots(0),fvHmpSignal(0), | |
63 | fvHmpTrkX(0),fvHmpTrkY(0),fvHmpTrkTheta(0),fvHmpTrkPhi(0), fvHmpCluQ(0), fvHmpCluX(0),fvHmpCluY(0),fvHmpCluCh(0),fvHmpCluSize(0),fvHmpCluBox(0),fvHmpCluStatus(0), | |
64 | fvEsdTrackAccepted(0), fvKinkIndex(0),fvTofSignal(0), | |
65 | ||
66 | fTree(0x0) | |
67 | { | |
68 | //Default ctor | |
69 | // | |
70 | for (Int_t i=0; i<3; i++) fvPhmp[i]=0; | |
71 | for (Int_t i=0; i<5; i++) fvHmpPid[i] = 0; | |
72 | } | |
73 | ||
74 | //___________________________________________________________________________ | |
75 | AliHMPIDPerformanceTask::AliHMPIDPerformanceTask(const Char_t* name) : | |
76 | AliAnalysisTaskSE(name), | |
77 | fESD(0x0), fESDfriend(0x0), fEsdTrackCuts(0x0), fClu(0x0),fMC(0x0), fEsdVtx(0x0),fCentrality(0x0), | |
78 | fEsdTrack(0x0),fEsdFriendTrack(0x0), fCalibObject(0x0), fUseMC(kTRUE), | |
79 | fHmpHistList(0x0), fHmpNevents(0x0),fHmpNevPerTrigClass(0x0), | |
80 | fGlobalEventNumber(0), fvType(0), fvFiredTriggerClasses(0x0), fvRunNumber(0), fvBunchCrossNumber(0),fvOrbitNumber(0), | |
81 | fvPeriodNumber(0),fvMagField(0), fvVertexX(0),fvVertexY(0),fvVertexZ(0),fvVertexNContributors(0), | |
82 | fvPesd(0),fvPhmpMag(0),fvCentrality(0), fvAcceptedTracks(0),fvRefMultTpc(0), | |
83 | fvHmpChi2(0),fvHmpCluIndx(0),fvHmpMipX(0),fvHmpMipY(0),fvHmpMipQ(0),fvHmpMipNPhots(0),fvHmpSignal(0), | |
84 | fvHmpTrkX(0),fvHmpTrkY(0),fvHmpTrkTheta(0),fvHmpTrkPhi(0), fvHmpCluQ(0), fvHmpCluX(0),fvHmpCluY(0),fvHmpCluCh(0),fvHmpCluSize(0),fvHmpCluBox(0),fvHmpCluStatus(0), | |
85 | fvEsdTrackAccepted(0), fvKinkIndex(0),fvTofSignal(0), | |
86 | ||
87 | fTree(0x0) | |
88 | { | |
89 | // | |
90 | // Constructor. Initialization of Inputs and Outputs | |
91 | // | |
92 | for (Int_t i=0; i<3; i++) fvPhmp[i]=0; | |
93 | for (Int_t i=0; i<5; i++) fvHmpPid[i] = 0; | |
94 | DefineOutput(1,TList::Class()); | |
95 | DefineOutput(2,TTree::Class()); | |
96 | } | |
97 | ||
98 | //___________________________________________________________________________ | |
99 | AliHMPIDPerformanceTask& AliHMPIDPerformanceTask::operator=(const AliHMPIDPerformanceTask& c) | |
100 | { | |
101 | // | |
102 | // Assignment operator | |
103 | // | |
104 | if (this!=&c) { | |
105 | AliAnalysisTaskSE::operator=(c); | |
106 | fESD = c.fESD; | |
107 | fESDfriend = c.fESDfriend; | |
108 | fEsdTrackCuts = c.fEsdTrackCuts; | |
109 | fClu = c.fClu; | |
110 | fMC = c.fMC; | |
111 | fEsdVtx = c.fEsdVtx; | |
112 | fCentrality = c.fCentrality; | |
113 | fEsdTrack = c.fEsdTrack; | |
114 | fEsdFriendTrack =c.fEsdFriendTrack; | |
115 | fCalibObject = c.fCalibObject; | |
116 | ||
117 | fUseMC = c.fUseMC; | |
118 | fHmpHistList = c.fHmpHistList; | |
119 | fHmpNevents = c.fHmpNevents; | |
120 | fHmpNevPerTrigClass = c.fHmpNevPerTrigClass; | |
121 | fGlobalEventNumber = c.fGlobalEventNumber; | |
122 | ||
123 | fvFiredTriggerClasses =c.fvFiredTriggerClasses; | |
124 | fvRunNumber = c.fvRunNumber; | |
125 | fvBunchCrossNumber = c.fvBunchCrossNumber; | |
126 | fvOrbitNumber = c.fvOrbitNumber; | |
127 | fvPeriodNumber = c.fvPeriodNumber; | |
128 | fvMagField =c.fvMagField; | |
129 | fvVertexX = c.fvVertexX; | |
130 | fvVertexY = c.fvVertexY; | |
131 | fvVertexZ = c.fvVertexZ; | |
132 | fvVertexNContributors = c.fvVertexNContributors; | |
133 | fvPesd = c.fvPesd; | |
134 | fvPhmpMag = c.fvPhmpMag; | |
135 | fvCentrality = c.fvCentrality; | |
136 | fvAcceptedTracks = c.fvAcceptedTracks; | |
137 | fvRefMultTpc = c.fvRefMultTpc; | |
138 | ||
139 | fvHmpChi2 = c.fvHmpChi2; | |
140 | fvHmpCluIndx = c.fvHmpCluIndx; | |
141 | fvHmpMipX = c.fvHmpMipX; | |
142 | fvHmpMipY = c.fvHmpMipY; | |
143 | fvHmpMipQ = c.fvHmpMipQ; | |
144 | fvHmpMipNPhots = c.fvHmpMipNPhots; | |
145 | fvHmpSignal = c.fvHmpSignal; | |
146 | fvHmpTrkX = c.fvHmpTrkX; | |
147 | fvHmpTrkY = c.fvHmpTrkY; | |
148 | fvHmpTrkTheta = c.fvHmpTrkTheta; | |
149 | fvHmpTrkPhi = c.fvHmpTrkPhi; | |
150 | ||
151 | fvHmpCluQ = c.fvHmpCluQ; | |
152 | fvHmpCluX = c.fvHmpCluX; | |
153 | fvHmpCluY = c.fvHmpCluY; | |
154 | fvHmpCluCh = c.fvHmpCluCh; | |
155 | fvHmpCluSize = c.fvHmpCluSize; | |
156 | fvHmpCluBox = c.fvHmpCluBox; | |
157 | fvHmpCluStatus = c.fvHmpCluStatus; | |
158 | fvEsdTrackAccepted = c.fvEsdTrackAccepted; | |
159 | ||
160 | fvKinkIndex = c.fvKinkIndex; | |
161 | fvTofSignal = c.fvTofSignal; | |
162 | ||
163 | fTree = c.fTree; | |
164 | ||
165 | for (Int_t i=0; i<3; i++) fvPhmp[i] = c.fvPhmp[i]; | |
166 | for (Int_t i=0; i<5; i++) fvHmpPid[i] = c.fvHmpPid[i]; | |
167 | } | |
168 | return *this; | |
169 | } | |
170 | ||
171 | //___________________________________________________________________________ | |
172 | AliHMPIDPerformanceTask::AliHMPIDPerformanceTask(const AliHMPIDPerformanceTask& c) : | |
173 | AliAnalysisTaskSE(c), | |
174 | fESD(c.fESD),fESDfriend(c.fESDfriend),fEsdTrackCuts(c.fEsdTrackCuts),fClu(c.fClu), | |
175 | fMC(c.fMC), fEsdVtx(c.fEsdVtx), fCentrality(c.fCentrality),fEsdTrack(c.fEsdTrack), | |
176 | fEsdFriendTrack(c.fEsdFriendTrack), | |
177 | fCalibObject(c.fCalibObject), | |
178 | ||
179 | fUseMC(c.fUseMC), | |
180 | fHmpHistList(c.fHmpHistList), | |
181 | fHmpNevents(c.fHmpNevents), | |
182 | fHmpNevPerTrigClass(c.fHmpNevPerTrigClass), | |
183 | fGlobalEventNumber(c.fGlobalEventNumber), | |
184 | fvFiredTriggerClasses(c.fvFiredTriggerClasses), | |
185 | fvRunNumber(c.fvRunNumber), | |
186 | fvBunchCrossNumber(c.fvBunchCrossNumber), | |
187 | fvOrbitNumber(c.fvOrbitNumber), | |
188 | fvPeriodNumber(c.fvPeriodNumber), | |
189 | fvMagField(c.fvMagField), | |
190 | fvVertexX(c.fvVertexX), | |
191 | fvVertexY(c.fvVertexY), | |
192 | fvVertexZ(c.fvVertexZ), | |
193 | fvVertexNContributors(c.fvVertexNContributors), | |
194 | fvPesd(c.fvPesd), | |
195 | fvPhmpMag(c.fvPhmpMag), | |
196 | fvCentrality(c.fvCentrality), | |
197 | fvAcceptedTracks(c.fvAcceptedTracks), | |
198 | fvRefMultTpc(c.fvRefMultTpc), | |
199 | fvHmpChi2 ( c.fvHmpChi2), | |
200 | fvHmpCluIndx ( c.fvHmpCluIndx), | |
201 | fvHmpMipX ( c.fvHmpMipX), | |
202 | fvHmpMipY ( c.fvHmpMipY), | |
203 | fvHmpMipQ ( c.fvHmpMipQ), | |
204 | fvHmpMipNPhots ( c.fvHmpMipNPhots), | |
205 | fvHmpSignal ( c.fvHmpSignal), | |
206 | fvHmpTrkX ( c.fvHmpTrkX), | |
207 | fvHmpTrkY ( c.fvHmpTrkY), | |
208 | fvHmpTrkTheta ( c.fvHmpTrkTheta), | |
209 | fvHmpTrkPhi ( c.fvHmpTrkPhi), | |
210 | ||
211 | fvHmpCluQ ( c.fvHmpCluQ), | |
212 | fvHmpCluX ( c.fvHmpCluX), | |
213 | fvHmpCluY ( c.fvHmpCluY), | |
214 | fvHmpCluCh ( c.fvHmpCluCh), | |
215 | fvHmpCluSize ( c.fvHmpCluSize), | |
216 | fvHmpCluBox ( c.fvHmpCluBox), | |
217 | fvHmpCluStatus ( c.fvHmpCluStatus), | |
218 | ||
219 | fvEsdTrackAccepted(c.fvEsdTrackAccepted), | |
220 | fvKinkIndex(c.fvKinkIndex), | |
221 | fvTofSignal(c.fvTofSignal), | |
222 | ||
223 | fTree(c.fTree) | |
224 | { | |
225 | // | |
226 | // Copy Constructor | |
227 | // | |
228 | for (Int_t i=0; i<3; i++) fvPhmp[i] = c.fvPhmp[i]; | |
229 | for (Int_t i=0; i<5; i++) fvHmpPid[i] = c.fvHmpPid[i]; | |
230 | } | |
231 | ||
232 | //___________________________________________________________________________ | |
233 | AliHMPIDPerformanceTask::~AliHMPIDPerformanceTask() { | |
234 | // | |
235 | //destructor | |
236 | // | |
237 | Info("~AliHMPIDPerformanceTask","Calling Destructor"); | |
238 | if (fHmpHistList && !AliAnalysisManager::GetAnalysisManager()->IsProofMode()) delete fHmpHistList; | |
239 | } | |
240 | ||
241 | //___________________________________________________________________________ | |
242 | void AliHMPIDPerformanceTask::ConnectInputData(Option_t *) | |
243 | { | |
244 | // Connect ESD here | |
245 | ||
246 | AliESDInputHandler *esdH = dynamic_cast<AliESDInputHandler*> (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()); | |
247 | if (!esdH) { | |
248 | AliDebug(2,Form("ERROR: Could not get ESDInputHandler")); | |
249 | } else | |
250 | fESD = esdH->GetEvent(); | |
251 | ||
252 | if (fUseMC){ | |
253 | // Connect MC | |
254 | AliMCEventHandler *mcH = dynamic_cast<AliMCEventHandler*> (AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()); | |
255 | if (!mcH) { | |
256 | AliDebug(2,Form("ERROR: Could not get MCEventHandler")); | |
257 | fUseMC = kFALSE; | |
258 | } else | |
259 | fMC = mcH->MCEvent(); | |
260 | if (!fMC) AliDebug(2,Form("ERROR: Could not get MCEvent")); | |
261 | } | |
262 | ||
263 | //Init ESD track cuts | |
264 | fEsdTrackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kTRUE); | |
265 | ||
266 | } | |
267 | //___________________________________________________________________________ | |
268 | void AliHMPIDPerformanceTask::UserExec(Option_t *) | |
269 | { | |
270 | ||
271 | if (!fESD) { | |
272 | Printf("ERROR: fESD not available"); | |
273 | return; | |
274 | } | |
275 | ||
276 | fHmpNevents->Fill(0); | |
277 | ||
278 | //fEsdVtx = fESD->GetPrimaryVertex(); | |
279 | ||
280 | /* check which trigger is fired */ fHmpNevPerTrigClass->Fill(0); | |
281 | if(fESD->IsTriggerClassFired("kMB")) fHmpNevPerTrigClass->Fill(1); | |
282 | if(fESD->IsTriggerClassFired("kINT7")) fHmpNevPerTrigClass->Fill(2); | |
283 | if(fESD->IsTriggerClassFired("kCINT5")) fHmpNevPerTrigClass->Fill(3); | |
284 | if(fESD->IsTriggerClassFired("kFastOnly")) fHmpNevPerTrigClass->Fill(4); | |
285 | if(fESD->IsTriggerClassFired("kHighMult")) fHmpNevPerTrigClass->Fill(5); | |
286 | if(fESD->IsTriggerClassFired("kAnyINT")) fHmpNevPerTrigClass->Fill(6); | |
287 | ||
288 | fvFiredTriggerClasses = Form("%s",(fESD->GetFiredTriggerClasses()).Data()); | |
289 | ||
290 | //___ check vertex | |
291 | if(!fESD->GetPrimaryVertex() || !fESD->GetPrimaryVertex() || fESD->GetPrimaryVertex()->GetNContributors()<1) { // SPD vertex | |
292 | if(!fESD->GetPrimaryVertexSPD()) return; | |
293 | if(!fESD->GetPrimaryVertexSPD()->GetStatus()) return; | |
294 | if(fESD->GetPrimaryVertexSPD()->GetNContributors()<1) return; // no good vertex, skip event | |
295 | } | |
296 | ||
297 | fHmpNevents->Fill(1); | |
298 | if(fESD->GetPrimaryVertex()) { | |
299 | fvVertexX = fESD->GetPrimaryVertex()->GetX(); | |
300 | fvVertexY = fESD->GetPrimaryVertex()->GetY(); | |
301 | fvVertexZ = fESD->GetPrimaryVertex()->GetZ(); | |
302 | fvVertexNContributors = fESD->GetPrimaryVertex()->GetNContributors(); | |
303 | fvMagField = fESD->GetMagneticField(); | |
304 | } | |
305 | ||
306 | fCentrality = fESD->GetCentrality(); | |
307 | if(fCentrality) | |
308 | { | |
309 | fvCentrality = fCentrality->GetCentralityPercentile("V0M"); | |
310 | } | |
311 | ||
312 | // Exec function | |
313 | // Loop over tracks and call Process function | |
314 | ||
315 | ||
316 | if(fEsdTrackCuts) | |
317 | { | |
318 | fvAcceptedTracks = fEsdTrackCuts->CountAcceptedTracks(fESD); | |
319 | fvRefMultTpc = fEsdTrackCuts->GetReferenceMultiplicity(fESD, kTRUE); | |
320 | } | |
321 | ||
322 | ||
323 | ||
324 | fESDfriend=(AliESDfriend*)(((AliESDEvent*)fESD)->FindListObject("AliESDfriend")); | |
325 | ||
326 | // if (!fESDfriend) { Printf("ERROR: fESDfriend not available"); } | |
327 | ||
328 | Double_t ktol = 0.001; | |
329 | ||
330 | fvRunNumber = fESD->GetRunNumber(); | |
331 | fvBunchCrossNumber = fESD->GetBunchCrossNumber(); | |
332 | fvOrbitNumber = fESD->GetOrbitNumber(); | |
333 | fvPeriodNumber = fESD->GetPeriodNumber(); | |
334 | ||
335 | fGlobalEventNumber = (ULong64_t)fESD->GetBunchCrossNumber() + (ULong64_t)fESD->GetOrbitNumber()*3564+ (ULong64_t)fESD->GetPeriodNumber()*16777215*3564; | |
336 | ||
337 | fEsdTrackCuts -> CountAcceptedTracks(fESD); | |
338 | ||
339 | ||
340 | for (Int_t iTrack=0; iTrack<fESD->GetNumberOfTracks(); iTrack++) { // tracks loop | |
341 | ||
342 | fEsdTrack=fESD->GetTrack(iTrack); | |
343 | if(!fEsdTrack) continue; | |
344 | ||
345 | if(fEsdTrackCuts && fEsdTrackCuts->AcceptTrack(fEsdTrack) == kTRUE ) fvEsdTrackAccepted = 1; | |
346 | else fvEsdTrackAccepted = 0; | |
347 | ||
348 | ||
349 | if(Equal(fEsdTrack->GetHMPIDsignal(),-20.,ktol)) continue; | |
350 | if(fEsdTrack->GetHMPIDcluIdx() < 0) continue; | |
351 | ||
352 | fvPesd = fEsdTrack->P(); | |
353 | ||
354 | //fEsdTrack->GetOuterHmpPxPyPz(fvPhmp); | |
355 | fvPhmpMag = TMath::Sqrt(fvPhmp[0]*fvPhmp[0] + fvPhmp[1]*fvPhmp[1] + fvPhmp[2]*fvPhmp[2] ); | |
356 | ||
357 | fvTofSignal = fEsdTrack->GetTOFsignal(); | |
358 | ||
359 | fEsdFriendTrack=fESDfriend->GetTrack(iTrack); | |
360 | ||
361 | //___ fill HMP stuff | |
362 | fvHmpChi2 = fEsdTrack -> GetHMPIDchi2(); | |
363 | fvHmpCluIndx = fEsdTrack -> GetHMPIDcluIdx(); | |
364 | fEsdTrack->GetHMPIDmip(fvHmpMipX,fvHmpMipY,fvHmpMipQ,fvHmpMipNPhots); | |
365 | ||
366 | //oid GetHMPIDpid(Double_t* p) const | |
367 | fvHmpSignal= fEsdTrack->GetHMPIDsignal(); | |
368 | fEsdTrack->GetHMPIDtrk(fvHmpTrkX,fvHmpTrkY,fvHmpTrkTheta, fvHmpTrkPhi); | |
369 | fEsdTrack->GetHMPIDpid(fvHmpPid); | |
370 | ||
371 | fvKinkIndex = fEsdTrack->GetKinkIndex(0); | |
372 | ||
373 | fvType = 0; fTree->Fill(); | |
374 | if(fEsdFriendTrack) { | |
375 | for(Int_t j=0;( fCalibObject=fEsdFriendTrack->GetCalibObject(j));++j){ // calib object loop | |
376 | if((fClu=dynamic_cast<AliHMPIDCluster*>(fCalibObject))){ // if HMPID calib object | |
377 | fvHmpCluQ = fClu->Q(); | |
378 | fvHmpCluX = fClu->X(); | |
379 | fvHmpCluY = fClu->Y(); | |
380 | fvHmpCluCh = fClu->Ch(); | |
381 | fvHmpCluSize = fClu->Size(); | |
382 | fvHmpCluBox = fClu->Box(); | |
383 | fvHmpCluStatus = fClu->Status(); | |
384 | ||
385 | fvType = 1; fTree->Fill(); | |
386 | ||
387 | } // if HMPID calib object | |
388 | } // calib object loop | |
389 | }//there is fEsdFriendTrack | |
390 | ||
391 | ||
392 | ||
393 | } // tracks loop | |
394 | ||
395 | /* PostData(0) is taken care of by AliAnalysisTaskSE */ | |
396 | PostData(1,fHmpHistList); | |
397 | PostData(2,fTree); | |
398 | } | |
399 | //___________________________________________________________________________ | |
400 | void AliHMPIDPerformanceTask::Terminate(Option_t*) | |
401 | { | |
402 | // The Terminate() function is the last function to be called during | |
403 | // a query. It always runs on the client, it can be used to present | |
404 | // the results graphically or save the results to file. | |
405 | ||
406 | Info("Terminate"," "); | |
407 | ||
408 | if (!fUseMC) return; | |
409 | ||
410 | fHmpHistList = dynamic_cast<TList*> (GetOutputData(1)); | |
411 | ||
412 | if (!fHmpHistList) { | |
413 | AliError("Histogram List is not available"); | |
414 | return; | |
415 | } | |
416 | ||
417 | AliAnalysisTaskSE::Terminate(); | |
418 | ||
419 | } | |
420 | ||
421 | //___________________________________________________________________________ | |
422 | void AliHMPIDPerformanceTask::UserCreateOutputObjects() { | |
423 | // | |
424 | //HERE ONE CAN CREATE OUTPUT OBJECTS | |
425 | //TO BE SET BEFORE THE EXECUTION OF THE TASK | |
426 | // | |
427 | ||
428 | //slot #1 | |
429 | // OpenFile(1); | |
430 | fHmpHistList = new TList(); | |
431 | fHmpHistList->SetOwner(); | |
432 | ||
433 | fHmpNevents = new TH1F("fHmpNevents","Number of events",2,0,2); | |
434 | fHmpHistList->Add(fHmpNevents); | |
435 | ||
436 | fHmpNevPerTrigClass = new TH1I("fHmpNevPerTrigClass","Number of event per trigger class;Trigger class;Events",10,0,10); | |
437 | fHmpNevPerTrigClass->GetXaxis()->SetBinLabel(1,"All coll."); | |
438 | fHmpNevPerTrigClass->GetXaxis()->SetBinLabel(2,"MB"); | |
439 | fHmpNevPerTrigClass->GetXaxis()->SetBinLabel(3,"INT7"); | |
440 | fHmpNevPerTrigClass->GetXaxis()->SetBinLabel(4,"CINT5"); | |
441 | fHmpNevPerTrigClass->GetXaxis()->SetBinLabel(5,"FastOnly"); | |
442 | fHmpNevPerTrigClass->GetXaxis()->SetBinLabel(6,"HighMult"); | |
443 | fHmpNevPerTrigClass->GetXaxis()->SetBinLabel(7,"AnyINT"); | |
444 | fHmpHistList->Add(fHmpNevPerTrigClass); | |
445 | ||
446 | ||
447 | // OpenFile(2); | |
448 | fTree = new TTree("TreePhot","Tree with photon clusters information"); | |
449 | fTree->Branch("fvFiredTriggerClasses",fvFiredTriggerClasses,"fvFiredTriggerClasses/C"); | |
450 | fTree->Branch("fvRunNumber",&fvRunNumber,"fvRunNumber/I"); | |
451 | fTree->Branch("fvBunchCrossNumber",&fvBunchCrossNumber,"fvBunchCrossNumber/I"); | |
452 | fTree->Branch("fvOrbitNumber",&fvOrbitNumber,"fvOrbitNumber/I"); | |
453 | fTree->Branch("fvPeriodNumber",&fvPeriodNumber,"fvPeriodNumber/I"); | |
454 | fTree->Branch("fvMagField",&fvMagField,"fvMagField/D"); | |
455 | fTree->Branch("fvVertexX",&fvVertexX,"fvVertexX/D"); | |
456 | fTree->Branch("fvVertexY",&fvVertexY,"fvVertexY/D"); | |
457 | fTree->Branch("fvVertexZ",&fvVertexZ,"fvVertexZ/D"); | |
458 | fTree->Branch("fvVertexNContributors",&fvVertexNContributors,"fvVertexNContributors/I"); | |
459 | fTree->Branch("fvCentrality",&fvCentrality,"fvCentrality/D"); | |
460 | fTree->Branch("fvAcceptedTracks",&fvAcceptedTracks,"fvAcceptedTracks/I"); | |
461 | fTree->Branch("fvRefMultTpc",&fvRefMultTpc,"fvRefMultTpc/I"); | |
462 | fTree->Branch("fvHmpChi2",&fvHmpChi2,"fvHmpChi2/D"); | |
463 | fTree->Branch("fvHmpCluIndx",&fvHmpCluIndx,"fvHmpCluIndx/I"); | |
464 | fTree->Branch("fvHmpMipX",&fvHmpMipX,"fvHmpMipX/F"); | |
465 | fTree->Branch("fvHmpMipY",&fvHmpMipY,"fvHmpMipY/F"); | |
466 | fTree->Branch("fvHmpMipQ",&fvHmpMipQ,"fvHmpMipQ/I"); | |
467 | fTree->Branch("fvHmpMipNPhots",&fvHmpMipNPhots,"fvHmpMipNPhots/I"); | |
468 | fTree->Branch("fvHmpSignal",&fvHmpSignal,"fvHmpSignal/D"); | |
469 | fTree->Branch("fvHmpTrkX",&fvHmpTrkX,"fvHmpTrkX/F"); | |
470 | fTree->Branch("fvHmpTrkY",&fvHmpTrkY,"fvHmpTrkY/F"); | |
471 | fTree->Branch("fvHmpTrkTheta",&fvHmpTrkTheta,"fvHmpTrkTheta/F"); | |
472 | fTree->Branch("fvHmpTrkPhi",&fvHmpTrkPhi,"fvHmpTrkPhi/F"); | |
473 | fTree->Branch("fvHmpPid",fvHmpPid,"fvHmpPid/D"); | |
474 | fTree->Branch("fvType",&fvType,"fvType/I"); | |
475 | fTree->Branch("fvHmpCluQ",&fvHmpCluQ,"fvHmpCluQ/D"); | |
476 | fTree->Branch("fvHmpCluX",&fvHmpCluX,"fvHmpCluX/D"); | |
477 | fTree->Branch("fvHmpCluY",&fvHmpCluY,"fvHmpCluY/D"); | |
478 | fTree->Branch("fvHmpCluCh",&fvHmpCluCh,"fvHmpCluCh/I"); | |
479 | fTree->Branch("fvHmpCluSize",&fvHmpCluSize,"fvHmpCluSize/I"); | |
480 | fTree->Branch("fvHmpCluBox",&fvHmpCluBox,"fvHmpCluBox/I"); | |
481 | fTree->Branch("fvHmpCluStatus",&fvHmpCluStatus,"fvHmpCluStatus/I"); | |
482 | fTree->Branch("fvEsdTrackAccepted",&fvEsdTrackAccepted,"fvEsdTrackAccepted/I"); | |
483 | fTree->Branch("fvKinkIndex",&fvKinkIndex,"fvKinkIndex/I"); | |
484 | fTree->Branch("fvTofSignal",&fvTofSignal,"fvTofSignal/D"); | |
485 | ||
486 | PostData(1,fHmpHistList); | |
487 | PostData(2,fTree); | |
488 | } | |
489 | //____________________________________________________________________________________________________________________________________ | |
490 | Bool_t AliHMPIDPerformanceTask::Equal(Double_t x, Double_t y, Double_t tolerance) | |
491 | { | |
492 | return abs(x - y) <= tolerance ; | |
493 | } | |
494 | //____________________________________________________________________________________________________________________________________ | |
495 | ||
496 | #endif |