]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG1/ITS/AliAnalysisTaskSDDRP.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Fix for savannah bug report 87728 (Laurent) + fix invalid read found with valgrind...
[u/mrichter/AliRoot.git] / PWG1 / ITS / AliAnalysisTaskSDDRP.cxx
CommitLineData
9bf660ff 1#include "AliAnalysisTask.h"
2#include "AliAnalysisManager.h"
3#include "AliAnalysisDataContainer.h"
4#include "AliITSRecPoint.h"
5#include "AliESDEvent.h"
6#include "AliTrackPointArray.h"
7#include "AliITSgeomTGeo.h"
8#include "AliESDfriend.h"
9#include "AliCDBManager.h"
10#include "AliCDBEntry.h"
11#include "AliITSCalibrationSDD.h"
12#include "AliITSresponseSDD.h"
13#include "AliGeomManager.h"
14#include <TSystem.h>
15#include <TTree.h>
16#include <TH1F.h>
c87da6ba 17#include <TH2F.h>
9bf660ff 18#include <TChain.h>
19#include <TGeoGlobalMagField.h>
20#include "AliESDInputHandlerRP.h"
c87da6ba 21/**************************************************************************
22 * Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
23 * *
24 * Author: The ALICE Off-line Project. *
25 * Contributors are mentioned in the code where appropriate. *
26 * *
27 * Permission to use, copy, modify and distribute this software and its *
28 * documentation strictly for non-commercial purposes is hereby granted *
29 * without fee, provided that the above copyright notice appears in all *
30 * copies and that both the copyright notice and this permission notice *
31 * appear in the supporting documentation. The authors make no claims *
32 * about the suitability of this software for any purpose. It is *
33 * provided "as is" without express or implied warranty. *
34 **************************************************************************/
35
36//*************************************************************************
37// Implementation of class AliAnalysiTaskSDDRP
38// AliAnalysisTaskSE to extract from ESD + ESDfreinds + ITS rec points
39// performance plots for SDD detector
40//
41// Author: F. Prino, prino@to.infn.it
42//*************************************************************************
43
44
9bf660ff 45#include "AliAnalysisTaskSDDRP.h"
46
47ClassImp(AliAnalysisTaskSDDRP)
48//______________________________________________________________________________
49AliAnalysisTaskSDDRP::AliAnalysisTaskSDDRP() : AliAnalysisTaskSE("SDD RecPoints"),
50 fOutput(0),
51 fHistNEvents(0),
307b389f 52 fHistCluInLay(0),
c87da6ba 53 fHistAllPMod(0),
54 fHistGoodPMod(0),
55 fHistBadRegMod(0),
56 fHistMissPMod(0),
57 fHistSkippedMod(0),
58 fHistOutAccMod(0),
59 fHistNoRefitMod(0),
307b389f 60 fHistAllPXloc(0),
61 fHistGoodPXloc(0),
62 fHistBadRegXloc(0),
63 fHistMissPXloc(0),
64 fHistAllPZloc(0),
65 fHistGoodPZloc(0),
66 fHistBadRegZloc(0),
67 fHistMissPZloc(0),
c87da6ba 68 fHistdEdxL3VsP(0),
69 fHistdEdxL4VsP(0),
70 fHistdEdxVsMod(0),
9bf660ff 71 fRecPMod(0),
72 fTrackPMod(0),
73 fGoodAnMod(0),
74 fRecPLadLay3(0),
75 fRecPLadLay4(0),
76 fTrackPLadLay3(0),
77 fTrackPLadLay4(0),
78 fGoodAnLadLay3(0),
79 fGoodAnLadLay4(0),
80 fDriftTimeRP(0),
c87da6ba 81 fDriftTimeTPAll(0),
82 fDriftTimeTPNoExtra(0),
83 fDriftTimeTPExtra(0),
307b389f 84 fCluSizAnVsTime(0),
85 fCluSizTbVsTime(0),
9bf660ff 86 fResp(0),
c87da6ba 87 fUseITSsaTracks(kFALSE),
9bf660ff 88 fMinITSpts(3),
c87da6ba 89 fMinTPCpts(70),
90 fMinPfordEdx(0.5),
9bf660ff 91 fOnlyCINT1BTrig(0),
70b1033a 92 fExcludeBadMod(kFALSE)
9bf660ff 93{
94 //
95 DefineOutput(1, TList::Class());
96}
97
98
99//___________________________________________________________________________
100AliAnalysisTaskSDDRP::~AliAnalysisTaskSDDRP(){
101 //
102 if (fOutput) {
103 delete fOutput;
104 fOutput = 0;
105 }
307b389f 106
9bf660ff 107}
108
109
110//___________________________________________________________________________
111
112void AliAnalysisTaskSDDRP::UserCreateOutputObjects() {
113
114 fOutput = new TList();
115 fOutput->SetOwner();
116 fOutput->SetName("OutputHistos");
117
118 fHistNEvents = new TH1F("hNEvents", "Number of processed events",3,-1.5,1.5);
119 fHistNEvents->Sumw2();
120 fHistNEvents->SetMinimum(0);
121 fOutput->Add(fHistNEvents);
122
307b389f 123 fHistCluInLay = new TH1F("hCluInLay","hCluInLay",7,-1.5,5.5);
124 fHistCluInLay->Sumw2();
125 fHistCluInLay->SetMinimum(0);
126 fOutput->Add(fHistCluInLay);
127
9bf660ff 128 // -- Module histos
129
c87da6ba 130 fHistAllPMod = new TH1F("hAllPmod","Crossing Tracks vs. Module",260,239.5,499.5);
131 fHistAllPMod->Sumw2();
132 fHistAllPMod->SetMinimum(0);
133 fOutput->Add(fHistAllPMod);
134
135 fHistGoodPMod = new TH1F("hGoodPmod","PointsAssocToTrack per Module",260,239.5,499.5);
136 fHistGoodPMod->Sumw2();
137 fHistGoodPMod->SetMinimum(0);
138 fOutput->Add(fHistGoodPMod);
139
140 fHistBadRegMod = new TH1F("hBadRegmod","Tracks in BadRegion per Module",260,239.5,499.5);
141 fHistBadRegMod->Sumw2();
142 fHistBadRegMod->SetMinimum(0);
143 fOutput->Add(fHistBadRegMod);
144
145 fHistMissPMod = new TH1F("hMissPmod","Missing Points per Module",260,239.5,499.5);
146 fHistMissPMod->Sumw2();
147 fHistMissPMod->SetMinimum(0);
148 fOutput->Add(fHistMissPMod);
149
150 fHistSkippedMod = new TH1F("hSkippedmod","Tracks in Skipped Module",260,239.5,499.5);
151 fHistSkippedMod->Sumw2();
152 fHistSkippedMod->SetMinimum(0);
153 fOutput->Add(fHistSkippedMod);
154
155 fHistOutAccMod = new TH1F("hOutAccmod","Tracks outside zAcc per Module",260,239.5,499.5);
156 fHistOutAccMod->Sumw2();
157 fHistOutAccMod->SetMinimum(0);
158 fOutput->Add(fHistOutAccMod);
159
160 fHistNoRefitMod = new TH1F("hNoRefitmod","Points rejected in refit per Module",260,239.5,499.5);
161 fHistNoRefitMod->Sumw2();
162 fHistNoRefitMod->SetMinimum(0);
163 fOutput->Add(fHistNoRefitMod);
164
c87da6ba 165
166
9bf660ff 167 fRecPMod = new TH1F("hRPMod","Rec Points per Module",260,239.5,499.5);
168 fRecPMod->Sumw2();
169 fRecPMod->SetMinimum(0);
170 fOutput->Add(fRecPMod);
171
172 fTrackPMod = new TH1F("hTPMod","Track Points per Module",260,239.5,499.5);
173 fTrackPMod->Sumw2();
174 fTrackPMod->SetMinimum(0);
175 fOutput->Add(fTrackPMod);
176
177 fGoodAnMod = new TH1F("hGAMod","Good Anodes per Module",260,239.5,499.5);
178 fOutput->Add(fGoodAnMod);
179
307b389f 180 // -- Local coordinates
181
182 fHistAllPXloc = new TH1F("hAllPxloc","Crossing Tracks vs. Xloc",75, -3.75, 3.75);
183 fHistAllPXloc->Sumw2();
184 fHistAllPXloc->SetMinimum(0);
185 fOutput->Add(fHistAllPXloc);
186
187 fHistGoodPXloc = new TH1F("hGoodPxloc","PointsAssocToTrack vs. Xloc",75, -3.75, 3.75);
188 fHistGoodPXloc->Sumw2();
189 fHistGoodPXloc->SetMinimum(0);
190 fOutput->Add(fHistGoodPXloc);
191
192 fHistBadRegXloc = new TH1F("hBadRegxloc","Tracks in BadRegion vs. Xloc",75, -3.75, 3.75);
193 fHistBadRegXloc->Sumw2();
194 fHistBadRegXloc->SetMinimum(0);
195 fOutput->Add(fHistBadRegXloc);
196
197 fHistMissPXloc = new TH1F("hMissPxloc","Missing Points vs. Xloc",75, -3.75, 3.75);
198 fHistMissPXloc->Sumw2();
199 fHistMissPXloc->SetMinimum(0);
200 fOutput->Add(fHistMissPXloc);
201
202 fHistAllPZloc = new TH1F("hAllPzloc","Crossing Tracks vs. Zloc",77, -3.85, 3.85);
203 fHistAllPZloc->Sumw2();
204 fHistAllPZloc->SetMinimum(0);
205 fOutput->Add(fHistAllPZloc);
206
207 fHistGoodPZloc = new TH1F("hGoodPzloc","PointsAssocToTrack vs. Zloc",77, -3.85, 3.85);
208 fHistGoodPZloc->Sumw2();
209 fHistGoodPZloc->SetMinimum(0);
210 fOutput->Add(fHistGoodPZloc);
211
212 fHistBadRegZloc = new TH1F("hBadRegzloc","Tracks in BadRegion vs. Zloc",77, -3.85, 3.85);
213 fHistBadRegZloc->Sumw2();
214 fHistBadRegZloc->SetMinimum(0);
215 fOutput->Add(fHistBadRegZloc);
216
217 fHistMissPZloc = new TH1F("hMissPzloc","Missing Points vs. Zloc",77, -3.85, 3.85);
218 fHistMissPZloc->Sumw2();
219 fHistMissPZloc->SetMinimum(0);
220 fOutput->Add(fHistMissPZloc);
221
9bf660ff 222 // -- Ladder histos
223
224 fRecPLadLay3 = new TH1F("hRPLad3","Rec Points per Ladder Layer 3",14,-0.5,13.5);
225 fRecPLadLay3->Sumw2();
226 fRecPLadLay3->SetMinimum(0);
227 fOutput->Add(fRecPLadLay3);
228
229 fRecPLadLay4 = new TH1F("hRPLad4","Rec Points per Ladder Layer 4",22,-0.5,21.5);
230 fRecPLadLay4->Sumw2();
231 fRecPLadLay4->SetMinimum(0);
232 fOutput->Add(fRecPLadLay4);
233
234 fTrackPLadLay3 = new TH1F("hTPLad3","Track Points per Ladder Layer 3",14,-0.5,13.5);
235 fTrackPLadLay3->Sumw2();
236 fTrackPLadLay3->SetMinimum(0);
237 fOutput->Add(fTrackPLadLay3);
238
239 fTrackPLadLay4 = new TH1F("hTPLad4","Track Points per Ladder Layer 4",22,-0.5,21.5);
240 fTrackPLadLay4->Sumw2();
241 fTrackPLadLay4->SetMinimum(0);
242 fOutput->Add(fTrackPLadLay4);
243
244 fGoodAnLadLay3 = new TH1F("hGALad3","Good Anodes per Ladder Layer 3",14,-0.5,13.5);
245 fOutput->Add(fGoodAnLadLay3);
246
247 fGoodAnLadLay4 = new TH1F("hGALad4","Good Anodes per Ladder Layer 4",22,-0.5,21.5);
248 fOutput->Add(fGoodAnLadLay4);
249
c87da6ba 250 fDriftTimeRP=new TH1F("hDrTimRP","Drift Time from Rec Points (ns)",640,0.,6400.);
9bf660ff 251 fDriftTimeRP->Sumw2();
252 fDriftTimeRP->SetMinimum(0.);
253 fOutput->Add(fDriftTimeRP);
254
c87da6ba 255 fDriftTimeTPAll=new TH1F("hDrTimTPAll","Drift Time from Track Points (ns)",640,0.,6400.);
256 fDriftTimeTPAll->Sumw2();
257 fDriftTimeTPAll->SetMinimum(0.);
258 fOutput->Add(fDriftTimeTPAll);
259
260 fDriftTimeTPNoExtra=new TH1F("hDrTimTPNoExtra","Drift Time from Track Points (ns)",640,0.,6400.);
261 fDriftTimeTPNoExtra->Sumw2();
262 fDriftTimeTPNoExtra->SetMinimum(0.);
263 fOutput->Add(fDriftTimeTPNoExtra);
264
265 fDriftTimeTPExtra=new TH1F("hDrTimTPExtra","Drift Time from Track Points (ns)",640,0.,6400.);
266 fDriftTimeTPExtra->Sumw2();
267 fDriftTimeTPExtra->SetMinimum(0.);
268 fOutput->Add(fDriftTimeTPExtra);
9bf660ff 269
307b389f 270 // dE/dx histos
271
272 fHistdEdxL3VsP=new TH2F("hdEdxL3VsP","dE/dx vs. p lay3",40,0.,2.,100,0.,500.);
273 fHistdEdxL3VsP->Sumw2();
274 fHistdEdxL3VsP->SetMinimum(0);
275 fOutput->Add(fHistdEdxL3VsP);
276
277 fHistdEdxL4VsP=new TH2F("hdEdxL4VsP","dE/dx vs. p lay4",40,0.,2.,100,0.,500);
278 fHistdEdxL4VsP->Sumw2();
279 fHistdEdxL4VsP->SetMinimum(0);
280 fOutput->Add(fHistdEdxL4VsP);
281
282 fHistdEdxVsMod=new TH2F("hdEdxVsMod","dE/dx vs. mod",260,239.5,499.5,100,0.,500.);
283 fHistdEdxVsMod->Sumw2();
284 fHistdEdxVsMod->SetMinimum(0);
285 fOutput->Add(fHistdEdxVsMod);
286
9bf660ff 287 for(Int_t it=0; it<8; it++){
288 fSignalTime[it]=new TH1F(Form("hSigTimeInt%d",it),Form("hSigTimeInt%d",it),100,0.,300.);
289 fSignalTime[it]->Sumw2();
290 fSignalTime[it]->SetMinimum(0);
291 fOutput->Add(fSignalTime[it]);
292 }
307b389f 293
294 // cluster size histos
295
296 fCluSizAnVsTime = new TH2F("hCluSizAn","hCluSizAn",40,0.,6400.,15,-0.5,14.5);
297 fCluSizAnVsTime->Sumw2();
298 fCluSizAnVsTime->SetMinimum(0);
299 fOutput->Add(fCluSizAnVsTime);
300
301 fCluSizTbVsTime = new TH2F("hCluSizTb","hCluSizTb",40,0.,6400.,15,-0.5,14.5);
302 fCluSizTbVsTime->Sumw2();
303 fCluSizTbVsTime->SetMinimum(0);
304 fOutput->Add(fCluSizTbVsTime);
305
306
307 // Read dead channels from OCDB
308
309 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
310 if (!mgr) AliFatal("No analysis manager");
311
312 AliCDBManager* man = AliCDBManager::Instance();
313 man->SetDefaultStorage("raw://");
314 Int_t nrun=mgr->GetRunFromPath();
315 man->SetRun(nrun);
316
317
318 AliCDBEntry* eR=(AliCDBEntry*)man->Get("ITS/Calib/RespSDD");
06482779 319 if (eR) {
320 eR->PrintId();
321 eR->PrintMetaData();
322 fResp=(AliITSresponseSDD*)eR->GetObject();
323 }
324
307b389f 325 AliCDBEntry* eC=(AliCDBEntry*)man->Get("ITS/Calib/CalibSDD");
b92e8782 326 Int_t countGood3[14];
327 Int_t countGood4[22];
328 Int_t countGoodMod[260];
06482779 329 if (eC) {
330 eC->PrintId();
331 eC->PrintMetaData();
332 TObjArray* calsdd=(TObjArray*)eC->GetObject();
06482779 333 for(Int_t ilad=0;ilad<14;ilad++) countGood3[ilad]=0;
334 for(Int_t ilad=0;ilad<22;ilad++) countGood4[ilad]=0;
335 for(Int_t imod=0;imod<260;imod++) countGoodMod[imod]=0;
336 for(Int_t imod=0;imod<260;imod++){
337 AliITSCalibrationSDD* cal=(AliITSCalibrationSDD*)calsdd->At(imod);
338 if(cal->IsBad()) continue;
339 Int_t modId=imod+AliITSgeomTGeo::GetModuleIndex(3,1,1);
340 Int_t lay,lad,det;
341 AliITSgeomTGeo::GetModuleId(modId,lay,lad,det);
342 if(fExcludeBadMod && !CheckModule(lay,lad,det)) continue;
343 for(Int_t ian=0; ian<512; ian++){
344 if(cal->IsBadChannel(ian)) continue;
345 countGoodMod[imod]++;
346 if(lay==3) countGood3[lad-1]++;
347 else if(lay==4) countGood4[lad-1]++;
348 }
349 }
307b389f 350 }
06482779 351
307b389f 352 for(Int_t imod=0;imod<260;imod++) fGoodAnMod->SetBinContent(imod+1,countGoodMod[imod]);
353 fGoodAnMod->SetMinimum(0);
354 for(Int_t ilad=0;ilad<14;ilad++) fGoodAnLadLay3->SetBinContent(ilad+1,countGood3[ilad]);
355 fGoodAnLadLay3->SetMinimum(0);
356 for(Int_t ilad=0;ilad<22;ilad++) fGoodAnLadLay4->SetBinContent(ilad+1,countGood4[ilad]);
357 fGoodAnLadLay4->SetMinimum(0);
358
359 PostData(1, fOutput);
360
9bf660ff 361}
362//______________________________________________________________________________
363void AliAnalysisTaskSDDRP::UserExec(Option_t *)
364{
365 //
307b389f 366 AliESDEvent *esd = (AliESDEvent*) (InputEvent());
9bf660ff 367
c87da6ba 368
307b389f 369 if(!esd) {
9bf660ff 370 printf("AliAnalysisTaskSDDRP::Exec(): bad ESD\n");
371 return;
372 }
9bf660ff 373
374
307b389f 375 if(!ESDfriend()) {
9bf660ff 376 printf("AliAnalysisTaskSDDRP::Exec(): bad ESDfriend\n");
377 return;
c87da6ba 378 }
379
380
9bf660ff 381 PostData(1, fOutput);
9bf660ff 382
383 fHistNEvents->Fill(0);
384 if(fOnlyCINT1BTrig){
307b389f 385 if(!esd->IsTriggerClassFired("CINT1B-ABCE-NOPF-ALL")) return;
9bf660ff 386 fHistNEvents->Fill(1);
387 }
388 const AliTrackPointArray *array = 0;
307b389f 389 Int_t ntracks = esd->GetNumberOfTracks();
9bf660ff 390 for (Int_t itrack=0; itrack < ntracks; itrack++) {
307b389f 391 AliESDtrack * track = esd->GetTrack(itrack);
9bf660ff 392 if (!track) continue;
c87da6ba 393
394 Bool_t accept=kTRUE;
395 if(fUseITSsaTracks){
396 if(track->GetNcls(1)>0) accept=kFALSE;
397 }else{
398 if(track->GetNcls(1)<fMinTPCpts) accept=kFALSE;
399 }
400 if(track->GetNcls(0) < fMinITSpts) accept=kFALSE;
401 Int_t trstatus=track->GetStatus();
402 if(!(trstatus&AliESDtrack::kITSrefit)) accept=kFALSE;
403 if(!accept) continue;
404
307b389f 405 fHistCluInLay->Fill(-1.); // bin -1 counts accepted tracks
406 UChar_t clumap=track->GetITSClusterMap();
407 for(Int_t iBit=0; iBit<6; iBit++){
408 if(clumap&(1<<iBit)) fHistCluInLay->Fill(iBit);
409 }
410
411
9bf660ff 412 Double_t dedx[4];
413 track->GetITSdEdxSamples(dedx);
c87da6ba 414 Float_t mom=track->P();
415 Int_t iMod,status;
416 Float_t xloc,zloc;
417 for(Int_t iLay=2; iLay<=3; iLay++){
418 Bool_t ok=track->GetITSModuleIndexInfo(iLay,iMod,status,xloc,zloc);
419 if(ok){
420 iMod+=240;
421 fHistAllPMod->Fill(iMod);
307b389f 422 fHistAllPXloc->Fill(xloc);
423 fHistAllPZloc->Fill(zloc);
c87da6ba 424 if(status==1){
425 fHistGoodPMod->Fill(iMod);
307b389f 426 fHistGoodPXloc->Fill(xloc);
427 fHistGoodPZloc->Fill(zloc);
c87da6ba 428 if(mom>fMinPfordEdx) fHistdEdxVsMod->Fill(iMod,dedx[iLay-2]);
429 if(iLay==2) fHistdEdxL3VsP->Fill(mom,dedx[0]);
430 else fHistdEdxL4VsP->Fill(mom,dedx[1]);
431 }
307b389f 432 else if(status==2){
433 fHistBadRegMod->Fill(iMod);
434 fHistBadRegXloc->Fill(xloc);
435 fHistBadRegZloc->Fill(zloc);
436 }
c87da6ba 437 else if(status==3) fHistSkippedMod->Fill(iMod);
438 else if(status==4) fHistOutAccMod->Fill(iMod);
307b389f 439 else if(status==5){
440 fHistMissPMod->Fill(iMod);
441 fHistMissPXloc->Fill(xloc);
442 fHistMissPZloc->Fill(zloc);
443 }
c87da6ba 444 else if(status==6) fHistNoRefitMod->Fill(iMod);
445 }
446 }
447
448
9bf660ff 449 array = track->GetTrackPointArray();
450 if(!array) continue;
451 for(Int_t ipt=0; ipt<array->GetNPoints(); ipt++) {
452 AliTrackPoint point;
453 Int_t modId;
454 array->GetPoint(point,ipt);
455 Int_t volId = point.GetVolumeID();
456 Int_t layerId = AliGeomManager::VolUIDToLayer(volId,modId);
457 if(layerId<3 || layerId>4) continue;
458 modId+=AliITSgeomTGeo::GetModuleIndex(layerId,1,1);
459 Int_t lay,lad,det;
460 AliITSgeomTGeo::GetModuleId(modId,lay,lad,det);
70b1033a 461 if(fExcludeBadMod && !CheckModule(lay,lad,det)) continue;
9bf660ff 462 fTrackPMod->Fill(modId);
c87da6ba 463 fDriftTimeTPAll->Fill(point.GetDriftTime());
464 if(point.IsExtra()) fDriftTimeTPExtra->Fill(point.GetDriftTime());
465 else fDriftTimeTPNoExtra->Fill(point.GetDriftTime());
9bf660ff 466 Float_t dtime=point.GetDriftTime()-fResp->GetTimeZero(modId);
307b389f 467 Int_t cluTyp=point.GetClusterType();
468 Int_t clSizAn=(cluTyp>>8)&0xFF;
469 Int_t clSizTb=cluTyp&0xFF;
470 fCluSizAnVsTime->Fill(dtime,clSizAn);
471 fCluSizTbVsTime->Fill(dtime,clSizTb);
9bf660ff 472 Int_t theBin=int(dtime/6500.*8.);
473 if(layerId==3){
474 fTrackPLadLay3->Fill(lad-1);
475 if(dedx[0]>0. && track->P()>fMinPfordEdx) fSignalTime[theBin]->Fill(dedx[0]);
476 }
477 if(layerId==4){
478 fTrackPLadLay4->Fill(lad-1);
c87da6ba 479 if(dedx[1]>0.&& track->P()>fMinPfordEdx) fSignalTime[theBin]->Fill(dedx[1]);
9bf660ff 480 }
481 }
482 }
483
484 AliESDInputHandlerRP *hand = dynamic_cast<AliESDInputHandlerRP*> (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
06482779 485 TTree* tR = 0;
486 if (hand) tR = hand->GetTreeR("ITS");
c87da6ba 487 if (tR){
488 TClonesArray *ITSrec= new TClonesArray("AliITSRecPoint");
489 TBranch *branch =tR->GetBranch("ITSRecPoints");
490 branch->SetAddress(&ITSrec);
491 for (Int_t modId=240; modId<500; modId++){
492 Int_t lay,lad,det;
493 AliITSgeomTGeo::GetModuleId(modId,lay,lad,det);
70b1033a 494 if(fExcludeBadMod && !CheckModule(lay,lad,det)) continue;
c87da6ba 495 branch->GetEvent(modId);
496 Int_t nrecp = ITSrec->GetEntries();
497 fRecPMod->Fill(modId,nrecp);
498 if(lay==3) fRecPLadLay3->Fill(lad-1,nrecp);
499 if(lay==4) fRecPLadLay4->Fill(lad-1,nrecp);
500 for(Int_t irec=0;irec<nrecp;irec++) {
501 AliITSRecPoint *recp = (AliITSRecPoint*)ITSrec->At(irec);
502 fDriftTimeRP->Fill(recp->GetDriftTime());
503 }
9bf660ff 504 }
c87da6ba 505 ITSrec->Delete();
506 delete ITSrec;
9bf660ff 507 }
9bf660ff 508 PostData(1,fOutput);
509
510}
511//______________________________________________________________________________
512Bool_t AliAnalysisTaskSDDRP::CheckModule(Int_t lay, Int_t lad, Int_t det) const{
513 //
514 if(lay==4){
515 if(lad==3 && det==5) return kFALSE; // 1500 V
516 if(lad==3 && det==6) return kFALSE; // 0 V
517 if(lad==3 && det==7) return kFALSE; // 1500 V
518 if(lad==4 && det==1) return kFALSE; // 0 V
519 if(lad==4 && det==2) return kFALSE; // 1500 V
520 if(lad==7 && det==5) return kFALSE; // 0 MV
521 if(lad==9 && det==3) return kFALSE; // 1500 V
522 if(lad==9 && det==4) return kFALSE; // 0 V
523 if(lad==9 && det==5) return kFALSE; // 1500 V
524 if(lad==11 && det==6) return kFALSE; // 1500 V
525 if(lad==11 && det==7) return kFALSE; // 0 V
526 if(lad==11 && det==8) return kFALSE; // 1500 V
527 if(lad==18 && det==5) return kFALSE; // 1500 V
528 if(lad==18 && det==6) return kFALSE; // 0 V
529 if(lad==18 && det==7) return kFALSE; // 1500 V
530 if(lad==22 && det==1) return kFALSE; // 0 V
531 if(lad==22 && det==2) return kFALSE; // 1500 V
532 }
533 if(lay==3){
534 if(lad==4 && det==4) return kFALSE; // 1500 V
535 if(lad==3) return kFALSE; // swapped in geometry
536 }
537 return kTRUE;
538}
539
540//______________________________________________________________________________
541void AliAnalysisTaskSDDRP::Terminate(Option_t */*option*/)
542{
543 // Terminate analysis
544 fOutput = dynamic_cast<TList*> (GetOutputData(1));
545 if (!fOutput) {
546 printf("ERROR: fOutput not available\n");
547 return;
548 }
549 fHistNEvents= dynamic_cast<TH1F*>(fOutput->FindObject("hNEvents"));
307b389f 550 fHistCluInLay= dynamic_cast<TH1F*>(fOutput->FindObject("hCluInLay"));
c87da6ba 551
552 fHistAllPMod= dynamic_cast<TH1F*>(fOutput->FindObject("hAllPMod"));
553 fHistGoodPMod= dynamic_cast<TH1F*>(fOutput->FindObject("hGoodPMod"));
554 fHistBadRegMod= dynamic_cast<TH1F*>(fOutput->FindObject("hBadRegMod"));
555 fHistMissPMod= dynamic_cast<TH1F*>(fOutput->FindObject("hMissPMod"));
556 fHistSkippedMod= dynamic_cast<TH1F*>(fOutput->FindObject("hSkippedMod"));
557 fHistOutAccMod= dynamic_cast<TH1F*>(fOutput->FindObject("hOutAccMod"));
558 fHistNoRefitMod= dynamic_cast<TH1F*>(fOutput->FindObject("hNoRefitMod"));
559
307b389f 560 fHistAllPXloc= dynamic_cast<TH1F*>(fOutput->FindObject("hAllPxloc"));
561 fHistGoodPXloc= dynamic_cast<TH1F*>(fOutput->FindObject("hGoodPxloc"));
562 fHistBadRegXloc= dynamic_cast<TH1F*>(fOutput->FindObject("hBadRegxloc"));
563 fHistMissPXloc= dynamic_cast<TH1F*>(fOutput->FindObject("hMissPxloc"));
564 fHistAllPZloc= dynamic_cast<TH1F*>(fOutput->FindObject("hAllPzloc"));
565 fHistGoodPZloc= dynamic_cast<TH1F*>(fOutput->FindObject("hGoodPzloc"));
566 fHistBadRegZloc= dynamic_cast<TH1F*>(fOutput->FindObject("fHistBadRegZloc"));
567 fHistMissPZloc= dynamic_cast<TH1F*>(fOutput->FindObject("hMissPzloc"));
568
c87da6ba 569 fHistdEdxL3VsP= dynamic_cast<TH2F*>(fOutput->FindObject("hdEdxL3VsP"));
570 fHistdEdxL4VsP= dynamic_cast<TH2F*>(fOutput->FindObject("hdEdxL4VsP"));
571 fHistdEdxVsMod= dynamic_cast<TH2F*>(fOutput->FindObject("hdEdxVsMod"));
572
9bf660ff 573 fRecPMod= dynamic_cast<TH1F*>(fOutput->FindObject("hRPMod"));
574 fTrackPMod= dynamic_cast<TH1F*>(fOutput->FindObject("hTPMod"));
575 fGoodAnMod= dynamic_cast<TH1F*>(fOutput->FindObject("hGAMod"));
576
577 fRecPLadLay3= dynamic_cast<TH1F*>(fOutput->FindObject("hRPLad3"));
578 fRecPLadLay4= dynamic_cast<TH1F*>(fOutput->FindObject("hRPLad4"));
579 fTrackPLadLay3= dynamic_cast<TH1F*>(fOutput->FindObject("hTPLad3"));
580 fTrackPLadLay4= dynamic_cast<TH1F*>(fOutput->FindObject("hTPLad4"));
581 fGoodAnLadLay3= dynamic_cast<TH1F*>(fOutput->FindObject("hGALad3"));
582 fGoodAnLadLay4= dynamic_cast<TH1F*>(fOutput->FindObject("hGALad4"));
583
584 fDriftTimeRP= dynamic_cast<TH1F*>(fOutput->FindObject("hDrTimRP"));
c87da6ba 585 fDriftTimeTPAll= dynamic_cast<TH1F*>(fOutput->FindObject("hDrTimTPAll"));
586 fDriftTimeTPNoExtra= dynamic_cast<TH1F*>(fOutput->FindObject("hDrTimTPNoExtra"));
587 fDriftTimeTPExtra= dynamic_cast<TH1F*>(fOutput->FindObject("hDrTimTPExtra"));
9bf660ff 588
589 for(Int_t it=0; it<8; it++){
590 fSignalTime[it]= dynamic_cast<TH1F*>(fOutput->FindObject(Form("hSigTimeInt%d",it)));
591 }
307b389f 592 fCluSizAnVsTime= dynamic_cast<TH2F*>(fOutput->FindObject("hCluSizAn"));
593 fCluSizTbVsTime= dynamic_cast<TH2F*>(fOutput->FindObject("hCluSizTb"));
9bf660ff 594
595 return;
596}
597
598
599
600
601