]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ITS/AliITSChannelStatus.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Macro to display the output of SDD task in QA-train
[u/mrichter/AliRoot.git] / ITS / AliITSChannelStatus.cxx
CommitLineData
4a60ccc4 1/**************************************************************************
2 * Copyright(c) 2007-2009, 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
6a0fdcbd 16/* $Id$ */
4a60ccc4 17
18//////////////////////////////////////////////////////////////////////////
19// //
20// Implementation of the class for bad channel treatment in the tracker //
21// Stores 1 status bit for each SPD pixel and SDD anode: //
22// 0 = bad channel //
23// 1 = good channel //
23197852 24// Dead and noisy channels are read from AliITSCalibration objects //
4a60ccc4 25// Origin: F.Prino, Torino, prino@to.infn.it //
26// //
27//////////////////////////////////////////////////////////////////////////
28
29#include "AliITSChannelStatus.h"
30#include "AliITSCalibrationSPD.h"
31#include "AliITSCalibrationSDD.h"
6a0fdcbd 32#include "AliITSCalibrationSSD.h"
b0702f40 33#include "AliITSsegmentationSPD.h"
34#include "AliITSsegmentationSDD.h"
6a0fdcbd 35#include "AliITSsegmentationSSD.h"
4a60ccc4 36#include "AliCDBEntry.h"
b0702f40 37#include "TMath.h"
4a60ccc4 38#include "AliLog.h"
39
40ClassImp(AliITSChannelStatus)
41
42
43//______________________________________________________________________
44AliITSChannelStatus::AliITSChannelStatus():
45TObject(),
46fSPDChannelStatus(0),
6a0fdcbd 47fSDDChannelStatus(0),
48fSSDChannelStatus(0)
4a60ccc4 49{
50 // default constructor
51 UInt_t nSPDchan=kSPDModules*kSPDNpxPerModule*kSPDNpzPerModule;
52 fSPDChannelStatus=new TBits(nSPDchan);
53 UInt_t nSDDchan=kSDDModules*kSDDAnodesPerModule;
54 fSDDChannelStatus=new TBits(nSDDchan);
6a0fdcbd 55 UInt_t nSSDchan=kSSDModules*kSSDStripsPerModule;
56 fSSDChannelStatus=new TBits(nSSDchan);
4a60ccc4 57 InitDefaults();
58}
59//______________________________________________________________________
60AliITSChannelStatus::AliITSChannelStatus(AliCDBManager *cdb):
61TObject(),
62fSPDChannelStatus(0),
6a0fdcbd 63fSDDChannelStatus(0),
64fSSDChannelStatus(0)
4a60ccc4 65{
66 AliCDBEntry* spdEntryD = cdb->Get("ITS/Calib/SPDDead");
67 if (!spdEntryD) AliFatal("Cannot get CDB entry for SPDDead");
68 TObjArray* deadArrSPD = (TObjArray*)spdEntryD->GetObject();
69 if (!deadArrSPD) AliFatal("No object found in SPDDead file");
70
71 AliCDBEntry* spdEntryN = cdb->Get("ITS/Calib/SPDNoisy");
72 if (!spdEntryN) AliFatal("Cannot get CDB entry for SPDNoisy");
73 TObjArray* noisArrSPD = (TObjArray*)spdEntryN->GetObject();
74 if (!noisArrSPD) AliFatal("No object found in SPDNoisy file");
75
76 AliCDBEntry* sddEntry = cdb->Get("ITS/Calib/CalibSDD");
77 if (!sddEntry) AliFatal("Cannot get CDB entry for CalibSDD");
78 TObjArray* calArrSDD = (TObjArray*)sddEntry->GetObject();
79 if (!calArrSDD) AliFatal("No object found in CalibSDD file");
80
6a0fdcbd 81 AliCDBEntry* ssdEntry = cdb->Get("ITS/Calib/CalibSSD");
82 if (!ssdEntry) AliFatal("Cannot get CDB entry for CalibSSD");
83 TObjArray* calArrSSD = (TObjArray*)ssdEntry->GetObject();
84 if (!calArrSSD) AliFatal("No object found in CalibSSD file");
85
4a60ccc4 86 UInt_t nSPDchan=kSPDModules*kSPDNpxPerModule*kSPDNpzPerModule;
87 fSPDChannelStatus=new TBits(nSPDchan);
88 UInt_t nSDDchan=kSDDModules*kSDDAnodesPerModule;
89 fSDDChannelStatus=new TBits(nSDDchan);
6a0fdcbd 90 UInt_t nSSDchan=kSSDModules*kSSDStripsPerModule;
91 fSSDChannelStatus=new TBits(nSSDchan);
92 InitFromOCDB(deadArrSPD,noisArrSPD,calArrSDD,calArrSSD);
4a60ccc4 93}
94//______________________________________________________________________
04c81913 95AliITSChannelStatus::AliITSChannelStatus(const AliITSDetTypeRec *dtrec):
23197852 96TObject(),
97fSPDChannelStatus(0),
6a0fdcbd 98fSDDChannelStatus(0),
99fSSDChannelStatus(0)
23197852 100{
101 UInt_t nSPDchan=kSPDModules*kSPDNpxPerModule*kSPDNpzPerModule;
102 fSPDChannelStatus=new TBits(nSPDchan);
103
104 UInt_t nSDDchan=kSDDModules*kSDDAnodesPerModule;
105 fSDDChannelStatus=new TBits(nSDDchan);
6a0fdcbd 106
107 UInt_t nSSDchan=kSSDModules*kSSDStripsPerModule;
108 fSSDChannelStatus=new TBits(nSSDchan);
23197852 109
110 // SPD modules
111 for(Int_t imod=0; imod<kSPDModules; imod++){
112 for(Int_t ix=0; ix<kSPDNpxPerModule; ix++){
113 for(Int_t iz=0; iz<kSPDNpzPerModule; iz++){
114 Int_t index=imod*kSPDNpxPerModule*kSPDNpzPerModule+ix*kSPDNpzPerModule+iz;
115 fSPDChannelStatus->SetBitNumber(index,kTRUE);
116 }
117 }
118 Int_t ix,iz;
119
120 // Mask SPD dead pixels
121 AliITSCalibrationSPD* deadspd=(AliITSCalibrationSPD*)dtrec->GetSPDDeadModel(imod);
122 for(Int_t ipix=0; ipix<deadspd->GetNrBad();ipix++){
123 deadspd->GetBadPixel(ipix,ix,iz);
124 Int_t index=imod*kSPDNpxPerModule*kSPDNpzPerModule+ix*kSPDNpzPerModule+iz;
125 fSPDChannelStatus->SetBitNumber(index,kFALSE);
126 }
127 // Mask SPD noisy pixels
128 AliITSCalibrationSPD* noisspd=(AliITSCalibrationSPD*)dtrec->GetCalibrationModel(imod);
129 for(Int_t ipix=0; ipix<noisspd->GetNrBad();ipix++){
130 noisspd->GetBadPixel(ipix,ix,iz);
131 Int_t index=imod*kSPDNpxPerModule*kSPDNpzPerModule+ix*kSPDNpzPerModule+iz;
132 fSPDChannelStatus->SetBitNumber(index,kFALSE);
133 }
134 }
135
136 // SDD modules
137 for(Int_t imod=0; imod<kSDDModules; imod++){
138 AliITSCalibrationSDD* calsdd=(AliITSCalibrationSDD*)dtrec->GetCalibrationModel(imod+kSPDModules);
139 for(Int_t ian=0; ian<kSDDAnodesPerModule; ian++){
140 Bool_t cstatus=kTRUE;
141 if(calsdd->IsBadChannel(ian)) cstatus=kFALSE;
142 Int_t index=imod*kSDDAnodesPerModule+ian;
143 fSDDChannelStatus->SetBitNumber(index,cstatus);
144 }
145 }
6a0fdcbd 146
147 // SSD modules
148 for (Int_t imod = 0; imod < kSSDModules; imod++) {
149 AliITSCalibrationSSD* calssd=(AliITSCalibrationSSD*)dtrec->GetCalibrationModel(kSSDFirstModule+imod);
150 for(Int_t ip=0; ip<kSSDStripsPerModule; ip++) {
151 Int_t index=imod*kSSDStripsPerModule+ip;
152 Bool_t cstatus = kTRUE;
153 if (ip < 768 && calssd->IsPChannelBad(ip))
154 cstatus = kFALSE;
155 if (ip >= 768 && calssd->IsNChannelBad(ip-768))
156 cstatus = kFALSE;
157
158 fSSDChannelStatus->SetBitNumber(index,cstatus);
159 }
160 }
23197852 161}
162//______________________________________________________________________
4a60ccc4 163void AliITSChannelStatus::InitDefaults(){
164 // fill bitmaps setting all channels as good
165 for(Int_t imod=0; imod<kSPDModules; imod++){
166 for(Int_t ix=0; ix<kSPDNpxPerModule; ix++){
167 for(Int_t iz=0; iz<kSPDNpzPerModule; iz++){
168 Int_t index=imod*kSPDNpxPerModule*kSPDNpzPerModule+ix*kSPDNpzPerModule+iz;
169 fSPDChannelStatus->SetBitNumber(index,kTRUE);
170 }
171 }
172 }
173 for(Int_t imod=0; imod<kSDDModules; imod++){
174 for(Int_t ian=0; ian<kSDDAnodesPerModule; ian++){
175 Int_t index=imod*kSDDAnodesPerModule+ian;
176 fSDDChannelStatus->SetBitNumber(index,kTRUE);
177 }
178 }
6a0fdcbd 179 for(Int_t imod=0; imod<kSSDModules; imod++){
180 for(Int_t is=0; is<kSSDStripsPerModule; is++){
181 Int_t index=imod*kSSDStripsPerModule+is;
182 fSSDChannelStatus->SetBitNumber(index,kTRUE);
183 }
184 }
4a60ccc4 185}
186//______________________________________________________________________
6a0fdcbd 187void AliITSChannelStatus::InitFromOCDB(TObjArray* deadArrSPD, TObjArray* noisArrSPD, TObjArray* calArrSDD, TObjArray *calArrSSD){
4a60ccc4 188// fills bitmaps from arrays of AliITSCalibrationSXD objects
189
190 // SPD modules
191 for(Int_t imod=0; imod<kSPDModules; imod++){
192 for(Int_t ix=0; ix<kSPDNpxPerModule; ix++){
193 for(Int_t iz=0; iz<kSPDNpzPerModule; iz++){
194 Int_t index=imod*kSPDNpxPerModule*kSPDNpzPerModule+ix*kSPDNpzPerModule+iz;
195 fSPDChannelStatus->SetBitNumber(index,kTRUE);
196 }
197 }
198 Int_t ix,iz;
199
200 // Mask SPD dead pixels
201 AliITSCalibrationSPD* deadspd=(AliITSCalibrationSPD*)deadArrSPD->At(imod);
202 for(Int_t ipix=0; ipix<deadspd->GetNrBad();ipix++){
203 deadspd->GetBadPixel(ipix,ix,iz);
204 Int_t index=imod*kSPDNpxPerModule*kSPDNpzPerModule+ix*kSPDNpzPerModule+iz;
205 fSPDChannelStatus->SetBitNumber(index,kFALSE);
206 }
207
208 // Mask SPD noisy pixels
209 AliITSCalibrationSPD* noisspd=(AliITSCalibrationSPD*)noisArrSPD->At(imod);
210 for(Int_t ipix=0; ipix<noisspd->GetNrBad();ipix++){
211 noisspd->GetBadPixel(ipix,ix,iz);
212 Int_t index=imod*kSPDNpxPerModule*kSPDNpzPerModule+ix*kSPDNpzPerModule+iz;
213 fSPDChannelStatus->SetBitNumber(index,kFALSE);
214 }
215 }
216
217 // SDD modules
218 for(Int_t imod=0; imod<kSDDModules; imod++){
219 AliITSCalibrationSDD* calsdd=(AliITSCalibrationSDD*)calArrSDD->At(imod);
220 for(Int_t ian=0; ian<kSDDAnodesPerModule; ian++){
221 Bool_t cstatus=kTRUE;
222 if(calsdd->IsBadChannel(ian)) cstatus=kFALSE;
223 Int_t index=imod*kSDDAnodesPerModule+ian;
224 fSDDChannelStatus->SetBitNumber(index,cstatus);
225 }
226 }
6a0fdcbd 227
228 // SSD modules
229 for (Int_t imod = 0; imod < kSSDModules; imod++) {
230 AliITSCalibrationSSD* calssd=(AliITSCalibrationSSD*)calArrSSD->At(imod);
231 for(Int_t ip=0; ip<kSSDStripsPerModule; ip++) {
232 Int_t index=imod*kSSDStripsPerModule+ip;
233 Bool_t cstatus = kTRUE;
234 if (ip < 768 && calssd->IsPChannelBad(ip))
235 cstatus = kFALSE;
236 if (ip >= 768 && calssd->IsNChannelBad(ip-768))
237 cstatus = kFALSE;
238
239 fSSDChannelStatus->SetBitNumber(index,cstatus);
240 }
241 }
4a60ccc4 242}
243//______________________________________________________________________
244AliITSChannelStatus::AliITSChannelStatus(const AliITSChannelStatus& cstatus):
245TObject(),
246fSPDChannelStatus(cstatus.fSPDChannelStatus),
6a0fdcbd 247fSDDChannelStatus(cstatus.fSDDChannelStatus),
248fSSDChannelStatus(cstatus.fSSDChannelStatus)
4a60ccc4 249{
250 // copy constructor
251}
252//______________________________________________________________________
253AliITSChannelStatus& AliITSChannelStatus::operator=(const AliITSChannelStatus& cstatus)
254{
255 // assignment operator
256 this->~AliITSChannelStatus();
257 new(this) AliITSChannelStatus(cstatus);
258 return *this;
259}
260
261//______________________________________________________________________
262AliITSChannelStatus::~AliITSChannelStatus(){
263 // destructor
264 if(fSPDChannelStatus) delete fSPDChannelStatus;
265 if(fSDDChannelStatus) delete fSDDChannelStatus;
6a0fdcbd 266 if(fSSDChannelStatus) delete fSSDChannelStatus;
4a60ccc4 267}
268
269//______________________________________________________________________
270Bool_t AliITSChannelStatus::CheckBounds(Int_t imod, Int_t iz, Int_t ix) const {
271 // check for out of bounds
6a0fdcbd 272 if(imod<0 || imod>=kSPDModules+kSDDModules+kSSDModules){
4a60ccc4 273 AliError(Form("Module number out of range 0-%d",kSPDModules+kSDDModules));
274 return kFALSE;
275 }
276 if(imod<kSPDModules){
277 if(ix<0 || ix>=kSPDNpxPerModule || iz<0 || iz>=kSPDNpzPerModule){
278 AliError("SPD: Pixel number out of range");
279 return kFALSE;
280 }
6a0fdcbd 281 }else if (imod<kSSDFirstModule){
4a60ccc4 282 if(iz<0 || iz>=kSDDAnodesPerModule){
283 AliError("SDD: anode number out of range");
284 return kFALSE;
285 }
286 }
6a0fdcbd 287 else {
288 if(iz<0 || iz>=kSSDStripsPerModule){
289 AliError("SSD: strip number out of range");
290 return kFALSE;
291 }
292 }
4a60ccc4 293 return kTRUE;
294}
295//______________________________________________________________________
296Bool_t AliITSChannelStatus::GetChannelStatus(Int_t imod, Int_t iz, Int_t ix) const {
297 // return status of inquired channel
298 if(CheckBounds(imod,iz,ix)==kFALSE) return kFALSE;
299 if(imod<kSPDModules){
300 Int_t index=imod*kSPDNpxPerModule*kSPDNpzPerModule+ix*kSPDNpzPerModule+iz;
301 return fSPDChannelStatus->TestBitNumber(index);
6a0fdcbd 302 } else if (imod < kSSDFirstModule) {
4a60ccc4 303 imod-=kSPDModules;
304 Int_t index=imod*kSDDAnodesPerModule+iz;
305 return fSDDChannelStatus->TestBitNumber(index);
306 }
6a0fdcbd 307 else { // SSD: iz is strip number 0-767 P-side, 768 - 1535 N-side
308 imod-=kSSDFirstModule;
309 Int_t index=imod*kSSDStripsPerModule+iz;
310 return fSSDChannelStatus->TestBitNumber(index);
311 }
4a60ccc4 312}
313//______________________________________________________________________
314void AliITSChannelStatus::SetChannelStatus(Bool_t cstatus, Int_t imod, Int_t iz, Int_t ix){
315 // set status for given channel
316 if(CheckBounds(imod,iz,ix)==kFALSE) return;
317 if(imod<kSPDModules){
318 Int_t index=imod*kSPDNpxPerModule*kSPDNpzPerModule+ix*kSPDNpzPerModule+iz;
319 fSPDChannelStatus->SetBitNumber(index,cstatus);
6a0fdcbd 320 }else if (imod < kSSDFirstModule) {
4a60ccc4 321 imod-=kSPDModules;
322 Int_t index=imod*kSDDAnodesPerModule+iz;
323 fSDDChannelStatus->SetBitNumber(index,cstatus);
324 }
6a0fdcbd 325 else { // SSD: iz is strip number 0-767 P-side, 768 - 1535 N-side
326 imod-=kSSDFirstModule;
327 Int_t index=imod*kSSDStripsPerModule+iz;
328 return fSSDChannelStatus->SetBitNumber(index,cstatus);
329 }
4a60ccc4 330}
b0702f40 331//______________________________________________________________________
332Bool_t AliITSChannelStatus::GetSDDLimits(Float_t zlocmin, Float_t zlocmax, Float_t xlocmin, Float_t xlocmax, Int_t& izmin, Int_t& izmax, Int_t& izmin2, Int_t& izmax2) const {
333 // Returns min. and max. anode numbers from local coordindate
334 AliITSsegmentationSDD *seg=new AliITSsegmentationSDD();
335 Float_t dummySpeed=6.5; // to avoid warnings in SDD segmentation
336 Float_t tolerance=0.9999;
337 seg->SetDriftSpeed(dummySpeed);
338 Float_t zHalfSize=0.5*seg->Dz()/10000.;
339 zHalfSize*=tolerance;
340 if(zlocmin<-zHalfSize) zlocmin=-zHalfSize;
341 if(zlocmax>zHalfSize) zlocmax=zHalfSize;
342 if(zlocmax<-zHalfSize || zlocmin>zHalfSize){
343 AliWarning("Search region completely outside module");
344 return kFALSE;
345 }
346 Float_t xHalfSize=seg->Dx()/10000.;
347 xHalfSize*=tolerance;
348 if(xlocmin<-xHalfSize) xlocmin=-xHalfSize;
349 if(xlocmax>xHalfSize) xlocmax=xHalfSize;
350 if(xlocmax<-xHalfSize || xlocmin>xHalfSize){
351 AliWarning("Search region completely outside module");
352 return kFALSE;
353 }
354
355 Int_t iSid1=seg->GetSideFromLocalX(xlocmin);
356 Int_t iSid2=seg->GetSideFromLocalX(xlocmax);
357 Int_t iz1,iz2,ixdummy;
358 seg->LocalToDet(xlocmin,zlocmin,ixdummy,iz1);
359 seg->LocalToDet(xlocmin,zlocmax,ixdummy,iz2);
360 izmin=TMath::Min(iz1,iz2);
361 izmax=TMath::Max(iz1,iz2);
362 if(iSid1==iSid2){
363 izmax2=izmin2=-1;
364 }else{
365 seg->LocalToDet(xlocmax,zlocmin,ixdummy,iz1);
366 seg->LocalToDet(xlocmax,zlocmax,ixdummy,iz2);
367 izmin2=TMath::Min(iz1,iz2);
368 izmax2=TMath::Max(iz1,iz2);
369 }
370 delete seg;
371 return kTRUE;
372}
373//______________________________________________________________________
374Bool_t AliITSChannelStatus::GetSPDLimits(Float_t zlocmin, Float_t zlocmax, Float_t xlocmin, Float_t xlocmax, Int_t& izmin, Int_t& izmax, Int_t& ixmin, Int_t& ixmax) const {
375 // Returns min. and max. pixel numbers from local coordindate
376 Float_t tolerance=0.9999;
377 AliITSsegmentationSPD *seg=new AliITSsegmentationSPD();
378 Float_t zHalfSize=0.5*seg->Dz()/10000.;
379 zHalfSize*=tolerance;
380 if(zlocmin<-zHalfSize) zlocmin=-zHalfSize;
381 if(zlocmax>zHalfSize) zlocmax=zHalfSize;
382 if(zlocmax<-zHalfSize || zlocmin>zHalfSize){
383 AliWarning("Search region completely outside module");
384 return kFALSE;
385 }
386 Float_t xHalfSize=0.5*seg->Dx()/10000.;
387 xHalfSize*=tolerance;
388 if(xlocmin<-xHalfSize) xlocmin=-xHalfSize;
389 if(xlocmax>xHalfSize) xlocmax=xHalfSize;
390 if(xlocmax<-xHalfSize || xlocmin>xHalfSize){
391 AliWarning("Search region completely outside module");
392 return kFALSE;
393 }
394
395 Int_t iz1,ix1,iz2,ix2;
396 seg->LocalToDet(xlocmin,zlocmin,ix1,iz1);
397 seg->LocalToDet(xlocmax,zlocmax,ix2,iz2);
398 izmin=TMath::Min(iz1,iz2);
399 izmax=TMath::Max(iz1,iz2);
400 ixmin=TMath::Min(ix1,ix2);
401 ixmax=TMath::Max(ix1,ix2);
402 delete seg;
403 return kTRUE;
404}
405//______________________________________________________________________
6a0fdcbd 406Bool_t AliITSChannelStatus::GetSSDLimits(Int_t layer, Float_t zlocmin, Float_t zlocmax, Float_t xlocmin, Float_t xlocmax, Int_t& iPmin, Int_t& iPmax, Int_t& iNmin, Int_t& iNmax) const {
407 // Returns min, max strip for SSD, given a search window
408 static AliITSsegmentationSSD seg;
409
410 AliDebug(2,Form("xmin %f zmin %f xmax %f zmax %f\n",xlocmin,zlocmin,xlocmax,zlocmax));
411
412 Int_t p,n;
413 seg.SetLayer(layer);
414 seg.GetPadIxz(xlocmin,zlocmin,p,n);
415 iPmin = iPmax = p;
416 iNmin = iNmax = n;
417 AliDebug(5,Form("lay %d xmin, zmin p %d n %d\n",layer,p,n));
418 seg.GetPadIxz(xlocmax,zlocmin,p,n);
419 iPmin = TMath::Min(iPmin,p);
7ef72c1e 420 iPmax = TMath::Max(iPmax,p);
6a0fdcbd 421 iNmin = TMath::Min(iNmin,n);
7ef72c1e 422 iNmax = TMath::Max(iNmax,n);
6a0fdcbd 423 AliDebug(5,Form("lay %d xmax, zmin p %d n %d\n",layer,p,n));
424 seg.GetPadIxz(xlocmax,zlocmax,p,n);
425 iPmin = TMath::Min(iPmin,p);
7ef72c1e 426 iPmax = TMath::Max(iPmax,p);
6a0fdcbd 427 iNmin = TMath::Min(iNmin,n);
7ef72c1e 428 iNmax = TMath::Max(iNmax,n);
6a0fdcbd 429 AliDebug(5,Form("lay %d xmax, zmax p %d n %d\n",layer,p,n));
430 seg.GetPadIxz(xlocmin,zlocmax,p,n);
431 iPmin = TMath::Min(iPmin,p);
7ef72c1e 432 iPmax = TMath::Max(iPmax,p);
6a0fdcbd 433 iNmin = TMath::Min(iNmin,n);
7ef72c1e 434 iNmax = TMath::Max(iNmax,n);
6a0fdcbd 435 AliDebug(5,Form("lay %d xmin, zmax p %d n %d\n",layer,p,n));
436
437 if (iPmin < 0)
438 iPmin = 0;
439 if (iNmin < 0)
440 iNmin = 0;
441 if (iPmax >= kSSDStripsPerSide)
442 iPmax = kSSDStripsPerSide-1;
443 if (iNmax >= kSSDStripsPerSide)
444 iNmax = kSSDStripsPerSide-1;
445 AliDebug(2,Form("lay %d p %d %d n %d %d\n",layer,iPmin,iPmax,iNmin,iNmax));
446 return kTRUE;
447}
448//______________________________________________________________________
b0702f40 449Bool_t AliITSChannelStatus::AnyBadInRoad(Int_t imod, Float_t zlocmin, Float_t zlocmax, Float_t xlocmin, Float_t xlocmax) const{
450 // Checks if there is at least one bad channel in the search road
6a0fdcbd 451 // For SSD: return kTRUE if there is at least one bad strip on both sides
452 // if only bad strips on one side 1D clusters can be used
453 AliDebug(2,Form("checking for module %d",imod));
b0702f40 454 if(imod<kSPDModules){
455 Int_t izmin,izmax,ixmin,ixmax;
456 Bool_t retcode=GetSPDLimits(zlocmin,zlocmax,xlocmin,xlocmax,izmin,izmax,ixmin,ixmax);
457 if(!retcode) return kFALSE;
458 for(Int_t iz=izmin; iz<=izmax;iz++){
459 for(Int_t ix=ixmin; ix<=ixmax;ix++){
460 if(GetChannelStatus(imod,iz,ix)==kFALSE) return kTRUE;
461 }
462 }
6a0fdcbd 463 }else if (imod < kSSDFirstModule) {
b0702f40 464 Int_t izmin,izmax,izmin2,izmax2;
465 Bool_t retcode=GetSDDLimits(zlocmin,zlocmax,xlocmin,xlocmax,izmin,izmax,izmin2,izmax2);
466 if(!retcode) return kFALSE;
467 for(Int_t iz=izmin; iz<=izmax;iz++){
6a0fdcbd 468 if(GetChannelStatus(imod,iz,0)==kFALSE) { return kTRUE;}
b0702f40 469 }
470 if(izmin2!=-1 && izmax2!=-1){
471 for(Int_t iz=izmin2; iz<=izmax2;iz++){
6a0fdcbd 472 if(GetChannelStatus(imod,iz,0)==kFALSE) { return kTRUE;}
b0702f40 473 }
474 }
475 }
6a0fdcbd 476 else {
477 Int_t layer = 5;
478 if (imod > kSSDMaxModLay5)
479 layer = 6;
480 Int_t iPmin,iPmax,iNmin,iNmax;
481 Bool_t retcode=GetSSDLimits(layer,zlocmin,zlocmax,xlocmin,xlocmax,iPmin,iPmax,iNmin,iNmax);
482 if(!retcode) return kFALSE;
483 Int_t nPbad = 0;
484 for(Int_t iP=iPmin; iP<=iPmax;iP++){
485 if(GetChannelStatus(imod,iP,0)==kFALSE) nPbad++;
486 }
487 if (nPbad == 0)
488 return kFALSE;
489 for(Int_t iN=iNmin; iN<=iNmax;iN++){
490 if(GetChannelStatus(imod,iN+768,0)==kFALSE) { return kTRUE; }
491 }
492 }
b0702f40 493 return kFALSE;
494}
495//______________________________________________________________________
496Float_t AliITSChannelStatus::FractionOfBadInRoad(Int_t imod, Float_t zlocmin, Float_t zlocmax, Float_t xlocmin, Float_t xlocmax) const{
497 // Calculate the fraction of bad channels in the road
6a0fdcbd 498 // Note: SSD returns fraction of dead strips on 'best' side. This
499 // is not proportional to dead surface area.
b0702f40 500 Float_t totChan=0.;
501 Float_t badChan=0.;
502 if(imod<kSPDModules){
503 Int_t izmin,izmax,ixmin,ixmax;
504 Bool_t retcode=GetSPDLimits(zlocmin,zlocmax,xlocmin,xlocmax,izmin,izmax,ixmin,ixmax);
505 if(!retcode) return 0.;
506 for(Int_t iz=izmin; iz<=izmax;iz++){
507 for(Int_t ix=ixmin; ix<=ixmax;ix++){
508 totChan+=1;
509 if(GetChannelStatus(imod,iz,ix)==kFALSE) badChan+=1.;
510 }
511 }
6a0fdcbd 512 }else if (imod < kSSDFirstModule) {
b0702f40 513 Int_t izmin,izmax,izmin2,izmax2;
514 Bool_t retcode=GetSDDLimits(zlocmin,zlocmax,xlocmin,xlocmax,izmin,izmax,izmin2,izmax2);
515 if(!retcode) return 0.;
516 for(Int_t iz=izmin; iz<=izmax;iz++){
517 totChan+=1;
518 if(GetChannelStatus(imod,iz,0)==kFALSE) badChan+=1.;
519 }
520 if(izmin2!=-1 && izmax2!=-1){
521 for(Int_t iz=izmin2; iz<=izmax2;iz++){
522 totChan+=1;
523 if(GetChannelStatus(imod,iz,0)==kFALSE) badChan+=1.;
524 }
525 }
526 }
6a0fdcbd 527 else {
528 Int_t layer = 5;
529 if (imod > kSSDMaxModLay5)
530 layer = 6;
531 Int_t iPmin,iPmax,iNmin,iNmax;
532 Bool_t retcode=GetSSDLimits(layer,zlocmin,zlocmax,xlocmin,xlocmax,iPmin,iPmax,iNmin,iNmax);
533 if(!retcode) return kFALSE;
534 Int_t nPbad = 0;
535 for(Int_t iP=iPmin; iP<=iPmax;iP++){
536 if(GetChannelStatus(imod,iP,0)==kFALSE) nPbad++;
537 }
538 Float_t fracP = (Float_t) nPbad / (iPmax-iPmin+1);
539 if (nPbad == 0)
540 return 0;
541 Int_t nNbad = 0;
542 for(Int_t iN=iNmin; iN<=iNmax;iN++){
543 if(GetChannelStatus(imod,iN+768,0)==kFALSE) nNbad++;
544 }
545 Float_t fracN = (Float_t) nNbad / (iPmax-iPmin+1);
546 return TMath::Min(fracP,fracN);
547 }
b0702f40 548 if(totChan==0.) return 0.;
549 else return badChan/totChan;
550}
6a0fdcbd 551