]>
Commit | Line | Data |
---|---|---|
04366a57 | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-2003, 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 | **************************************************************************/ | |
d76cb3ad | 15 | |
16 | /* $Id$*/ | |
17 | ||
04366a57 | 18 | //////////////////////////////////////////////////////////////////////////// |
19 | // Implementation of the ITS clusterer V2 class // | |
20 | // // | |
21 | // Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch // | |
22 | // // | |
23 | /////////////////////////////////////////////////////////////////////////// | |
24 | ||
25 | ||
04366a57 | 26 | |
ae485645 | 27 | #include <TClonesArray.h> |
04366a57 | 28 | #include "AliITSClusterFinderV2SDD.h" |
00a7cc50 | 29 | #include "AliITSRecPoint.h" |
7d62fb64 | 30 | #include "AliITSDetTypeRec.h" |
04366a57 | 31 | #include "AliRawReader.h" |
0599a018 | 32 | #include "AliITSRawStreamSDD.h" |
2f8ed7ab | 33 | #include "AliITSRawStreamSDDCompressed.h" |
5bfe44ce | 34 | #include "AliITSCalibrationSDD.h" |
253e68a0 | 35 | #include "AliITSresponseSDD.h" |
5bfe44ce | 36 | #include "AliITSDetTypeRec.h" |
debf7b78 | 37 | #include "AliITSReconstructor.h" |
5bfe44ce | 38 | #include "AliITSsegmentationSDD.h" |
04366a57 | 39 | #include "AliITSdigitSDD.h" |
b4704be3 | 40 | #include "AliITSgeomTGeo.h" |
debf7b78 | 41 | |
04366a57 | 42 | ClassImp(AliITSClusterFinderV2SDD) |
43 | ||
4952f440 | 44 | AliITSClusterFinderV2SDD::AliITSClusterFinderV2SDD(AliITSDetTypeRec* dettyp):AliITSClusterFinderV2(dettyp) |
a66a0eb6 | 45 | { |
04366a57 | 46 | |
47 | //Default constructor | |
48 | ||
04366a57 | 49 | } |
50 | ||
51 | ||
52 | void AliITSClusterFinderV2SDD::FindRawClusters(Int_t mod){ | |
53 | ||
54 | //Find clusters V2 | |
55 | SetModule(mod); | |
56 | FindClustersSDD(fDigits); | |
57 | ||
58 | } | |
59 | ||
60 | void AliITSClusterFinderV2SDD::FindClustersSDD(TClonesArray *digits) { | |
61 | //------------------------------------------------------------ | |
62 | // Actual SDD cluster finder | |
63 | //------------------------------------------------------------ | |
4952f440 | 64 | Int_t nAnodes = GetSeg()->NpzHalf(); |
65 | Int_t nzBins = nAnodes+2; | |
66 | Int_t nTimeBins = GetSeg()->Npx(); | |
67 | Int_t nxBins = nTimeBins+2; | |
f6b6d58e | 68 | const Int_t kMaxBin=nzBins*nxBins; |
04366a57 | 69 | |
70 | AliBin *bins[2]; | |
4952f440 | 71 | bins[0]=new AliBin[kMaxBin]; |
72 | bins[1]=new AliBin[kMaxBin]; | |
404c1c29 | 73 | AliITSCalibrationSDD* cal = (AliITSCalibrationSDD*)GetResp(fModule); |
9c7f5947 | 74 | if(cal==0){ |
75 | AliError(Form("Calibration object not present for SDD module %d\n",fModule)); | |
76 | return; | |
77 | } | |
83ec5e27 | 78 | |
04366a57 | 79 | AliITSdigitSDD *d=0; |
80 | Int_t i, ndigits=digits->GetEntriesFast(); | |
81 | for (i=0; i<ndigits; i++) { | |
82 | d=(AliITSdigitSDD*)digits->UncheckedAt(i); | |
20f3f947 | 83 | Int_t ian=d->GetCoord1(); |
84 | Int_t itb=d->GetCoord2(); | |
85 | Int_t iSide=0; | |
86 | if (ian >= nAnodes) iSide=1; | |
87 | Float_t gain=cal->GetChannelGain(ian); | |
83ec5e27 | 88 | Float_t charge=d->GetSignal(); // returns expanded signal |
89 | // (10 bit, low threshold already added) | |
20f3f947 | 90 | Float_t baseline = cal->GetBaseline(ian); |
91 | if(charge>baseline) charge-=baseline; | |
92 | else charge=0; | |
5683bd96 | 93 | |
9c7f5947 | 94 | if(gain>0.){ // Bad channels have gain=0. |
95 | charge/=gain; | |
20f3f947 | 96 | if(charge<cal->GetThresholdAnode(ian)) continue; |
9c7f5947 | 97 | Int_t q=(Int_t)(charge+0.5); |
20f3f947 | 98 | Int_t y=itb+1; |
99 | Int_t z=ian+1; | |
9c7f5947 | 100 | if (z <= nAnodes){ |
101 | bins[0][y*nzBins+z].SetQ(q); | |
102 | bins[0][y*nzBins+z].SetMask(1); | |
103 | bins[0][y*nzBins+z].SetIndex(i); | |
104 | } else { | |
105 | z-=nAnodes; | |
106 | bins[1][y*nzBins+z].SetQ(q); | |
107 | bins[1][y*nzBins+z].SetMask(1); | |
108 | bins[1][y*nzBins+z].SetIndex(i); | |
109 | } | |
04366a57 | 110 | } |
111 | } | |
112 | ||
4952f440 | 113 | FindClustersSDD(bins, kMaxBin, nzBins, digits); |
04366a57 | 114 | |
115 | delete[] bins[0]; | |
116 | delete[] bins[1]; | |
117 | ||
118 | } | |
119 | ||
120 | void AliITSClusterFinderV2SDD:: | |
121 | FindClustersSDD(AliBin* bins[2], Int_t nMaxBin, Int_t nzBins, | |
122 | TClonesArray *digits, TClonesArray *clusters) { | |
123 | //------------------------------------------------------------ | |
124 | // Actual SDD cluster finder | |
125 | //------------------------------------------------------------ | |
b4704be3 | 126 | |
debf7b78 | 127 | static AliITSRecoParam *repa = NULL; |
128 | if(!repa){ | |
129 | repa = (AliITSRecoParam*) AliITSReconstructor::GetRecoParam(); | |
130 | if(!repa){ | |
ed446fa3 | 131 | repa = AliITSRecoParam::GetHighFluxParam(); |
debf7b78 | 132 | AliWarning("Using default AliITSRecoParam class"); |
133 | } | |
134 | } | |
b4704be3 | 135 | const TGeoHMatrix *mT2L=AliITSgeomTGeo::GetTracking2LocalMatrix(fModule); |
404c1c29 | 136 | AliITSCalibrationSDD* cal = (AliITSCalibrationSDD*)GetResp(fModule); |
9c7f5947 | 137 | if(cal==0){ |
138 | AliError(Form("Calibration object not present for SDD module %d\n",fModule)); | |
139 | return; | |
140 | } | |
04366a57 | 141 | Int_t ncl=0; |
142 | TClonesArray &cl=*clusters; | |
143 | for (Int_t s=0; s<2; s++) | |
144 | for (Int_t i=0; i<nMaxBin; i++) { | |
e3ee4984 | 145 | NoiseSuppress(i,s,nzBins,bins[s],cal); |
04366a57 | 146 | if (bins[s][i].IsUsed()) continue; |
147 | Int_t idx[32]; UInt_t msk[32]; Int_t npeaks=0; | |
148 | FindPeaks(i, nzBins, bins[s], idx, msk, npeaks); | |
149 | ||
150 | if (npeaks>30) continue; | |
151 | if (npeaks==0) continue; | |
152 | ||
153 | Int_t k,l; | |
154 | for (k=0; k<npeaks-1; k++){//mark adjacent peaks | |
155 | if (idx[k] < 0) continue; //this peak is already removed | |
156 | for (l=k+1; l<npeaks; l++) { | |
f6b6d58e | 157 | if (idx[l] < 0) continue; //this peak is already removed |
158 | Int_t ki=idx[k]/nzBins, kj=idx[k] - ki*nzBins; | |
159 | Int_t li=idx[l]/nzBins, lj=idx[l] - li*nzBins; | |
160 | Int_t di=TMath::Abs(ki - li); | |
161 | Int_t dj=TMath::Abs(kj - lj); | |
162 | if (di>1 || dj>1) continue; | |
163 | if (bins[s][idx[k]].GetQ() > bins[s][idx[l]].GetQ()) { | |
164 | msk[l]=msk[k]; | |
165 | idx[l]*=-1; | |
166 | } else { | |
167 | msk[k]=msk[l]; | |
168 | idx[k]*=-1; | |
169 | break; | |
170 | } | |
04366a57 | 171 | } |
172 | } | |
173 | ||
174 | for (k=0; k<npeaks; k++) { | |
debf7b78 | 175 | if(repa->GetUseUnfoldingInClusterFinderSDD()==kFALSE) msk[k]=msk[0]; |
04366a57 | 176 | MarkPeak(TMath::Abs(idx[k]), nzBins, bins[s], msk[k]); |
177 | } | |
178 | ||
179 | for (k=0; k<npeaks; k++) { | |
f6b6d58e | 180 | if (idx[k] < 0) continue; //removed peak |
181 | AliITSRecPoint c; | |
182 | MakeCluster(idx[k], nzBins, bins[s], msk[k], c); | |
183 | //mi change | |
184 | Int_t milab[10]; | |
185 | for (Int_t ilab=0;ilab<10;ilab++){ | |
186 | milab[ilab]=-2; | |
187 | } | |
188 | Int_t maxi=0,mini=0,maxj=0,minj=0; | |
189 | //AliBin *bmax=&bins[s][idx[k]]; | |
190 | //Float_t max = TMath::Max(TMath::Abs(bmax->GetQ())/5.,3.); | |
191 | ||
192 | for (Int_t di=-2; di<=2;di++){ | |
193 | for (Int_t dj=-3;dj<=3;dj++){ | |
194 | Int_t index = idx[k]+di+dj*nzBins; | |
195 | if (index<0) continue; | |
196 | if (index>=nMaxBin) continue; | |
197 | AliBin *b=&bins[s][index]; | |
198 | Int_t nAnode=index%nzBins-1; | |
199 | Int_t adcSignal=b->GetQ(); | |
200 | if(adcSignal>cal->GetThresholdAnode(nAnode)){ | |
201 | if (di>maxi) maxi=di; | |
202 | if (di<mini) mini=di; | |
203 | if (dj>maxj) maxj=dj; | |
204 | if (dj<minj) minj=dj; | |
205 | } | |
206 | // | |
207 | if(digits) { | |
208 | if (TMath::Abs(di)<2&&TMath::Abs(dj)<2){ | |
209 | AliITSdigitSDD* d=(AliITSdigitSDD*)digits->UncheckedAt(b->GetIndex()); | |
210 | for (Int_t itrack=0;itrack<10;itrack++){ | |
211 | Int_t track = (d->GetTracks())[itrack]; | |
212 | if (track>=0) { | |
213 | AddLabel(milab, track); | |
214 | } | |
215 | } | |
216 | } | |
217 | } | |
218 | } | |
219 | } | |
220 | ||
221 | ||
253e68a0 | 222 | AliITSresponseSDD* rsdd = fDetTypeRec->GetResponseSDD(); |
223 | Float_t y=c.GetY(),z=c.GetZ(), q=c.GetQ(); | |
f6b6d58e | 224 | y/=q; z/=q; |
225 | Float_t zAnode=z-0.5; // to have anode in range 0.-255. and centered on the mid of the pitch | |
226 | Float_t timebin=y-0.5; // to have time bin in range 0.-255. amd centered on the mid of the bin | |
227 | if(s==1) zAnode += GetSeg()->NpzHalf(); // right side has anodes from 256. to 511. | |
228 | Float_t zdet = GetSeg()->GetLocalZFromAnode(zAnode); | |
253e68a0 | 229 | Float_t driftTime = GetSeg()->GetDriftTimeFromTb(timebin) - rsdd->GetTimeOffset(); |
f6b6d58e | 230 | Float_t driftPathMicron = cal->GetDriftPath(driftTime,zAnode); |
231 | const Double_t kMicronTocm = 1.0e-4; | |
232 | Float_t xdet=(driftPathMicron-GetSeg()->Dx())*kMicronTocm; // xdet is negative | |
233 | if (s==0) xdet=-xdet; // left side has positive local x | |
c9a38d3d | 234 | |
235 | Float_t corrx=0, corrz=0; | |
236 | cal->GetCorrections(zdet,xdet,corrz,corrx,GetSeg()); | |
237 | zdet+=corrz; | |
238 | xdet+=corrx; | |
f6b6d58e | 239 | |
240 | Double_t loc[3]={xdet,0.,zdet},trk[3]={0.,0.,0.}; | |
241 | mT2L->MasterToLocal(loc,trk); | |
242 | y=trk[1]; | |
243 | z=trk[2]; | |
244 | ||
253e68a0 | 245 | q/=rsdd->GetADC2keV(); //to have MPV 1 MIP = 86.4 KeV |
f6b6d58e | 246 | Float_t hit[5] = {y, z, 0.0030*0.0030, 0.0020*0.0020, q}; |
247 | Int_t info[3] = {maxj-minj+1, maxi-mini+1, fNlayer[fModule]}; | |
248 | if (digits) { | |
249 | // AliBin *b=&bins[s][idx[k]]; | |
250 | // AliITSdigitSDD* d=(AliITSdigitSDD*)digits->UncheckedAt(b->GetIndex()); | |
251 | { | |
252 | //Int_t lab[3]; | |
253 | //lab[0]=(d->GetTracks())[0]; | |
254 | //lab[1]=(d->GetTracks())[1]; | |
255 | //lab[2]=(d->GetTracks())[2]; | |
256 | //CheckLabels(lab); | |
257 | CheckLabels2(milab); | |
258 | } | |
259 | } | |
260 | milab[3]=fNdet[fModule]; | |
261 | ||
262 | AliITSRecPoint cc(milab,hit,info); | |
263 | cc.SetType(npeaks); | |
264 | if(clusters) new (cl[ncl]) AliITSRecPoint(cc); | |
265 | else { | |
266 | fDetTypeRec->AddRecPoint(cc); | |
267 | } | |
268 | ncl++; | |
04366a57 | 269 | } |
270 | } | |
24c98369 | 271 | |
04366a57 | 272 | } |
f6b6d58e | 273 | //______________________________________________________________________ |
04366a57 | 274 | void AliITSClusterFinderV2SDD::RawdataToClusters(AliRawReader* rawReader,TClonesArray** clusters){ |
275 | //------------------------------------------------------------ | |
276 | // This function creates ITS clusters from raw data | |
277 | //------------------------------------------------------------ | |
278 | rawReader->Reset(); | |
2f8ed7ab | 279 | AliITSRawStream* inputSDD; |
280 | if(fDetTypeRec->IsHLTmodeC()==kTRUE){ | |
281 | inputSDD=new AliITSRawStreamSDDCompressed(rawReader); | |
282 | }else{ | |
283 | inputSDD=new AliITSRawStreamSDD(rawReader); | |
284 | } | |
285 | ||
979b5a5f | 286 | AliITSDDLModuleMapSDD *ddlmap=(AliITSDDLModuleMapSDD*)fDetTypeRec->GetDDLModuleMapSDD(); |
2f8ed7ab | 287 | inputSDD->SetDDLModuleMap(ddlmap); |
83ec5e27 | 288 | for(Int_t iddl=0; iddl<AliITSDDLModuleMapSDD::GetNDDLs(); iddl++){ |
289 | for(Int_t icar=0; icar<AliITSDDLModuleMapSDD::GetNModPerDDL();icar++){ | |
290 | Int_t iMod=ddlmap->GetModuleNumber(iddl,icar); | |
291 | if(iMod==-1) continue; | |
292 | AliITSCalibrationSDD* cal = (AliITSCalibrationSDD*)GetResp(iMod); | |
293 | if(cal==0){ | |
294 | AliError(Form("Calibration object not present for SDD module %d\n",iMod)); | |
295 | continue; | |
296 | } | |
253e68a0 | 297 | Bool_t isZeroSupp=cal->GetZeroSupp(); |
298 | if(isZeroSupp){ | |
2f8ed7ab | 299 | for(Int_t iSid=0; iSid<2; iSid++) inputSDD->SetZeroSuppLowThreshold(iMod-240,iSid,cal->GetZSLowThreshold(iSid)); |
83ec5e27 | 300 | }else{ |
2f8ed7ab | 301 | for(Int_t iSid=0; iSid<2; iSid++) inputSDD->SetZeroSuppLowThreshold(iMod-240,iSid,0); |
83ec5e27 | 302 | } |
303 | } | |
304 | } | |
2f8ed7ab | 305 | FindClustersSDD(inputSDD,clusters); |
306 | delete inputSDD; | |
04366a57 | 307 | } |
308 | ||
309 | void AliITSClusterFinderV2SDD::FindClustersSDD(AliITSRawStream* input, | |
310 | TClonesArray** clusters) | |
311 | { | |
312 | //------------------------------------------------------------ | |
313 | // Actual SDD cluster finder for raw data | |
314 | //------------------------------------------------------------ | |
315 | Int_t nClustersSDD = 0; | |
4952f440 | 316 | Int_t nAnodes = GetSeg()->NpzHalf(); |
317 | Int_t nzBins = nAnodes+2; | |
318 | Int_t nTimeBins = GetSeg()->Npx(); | |
319 | Int_t nxBins = nTimeBins+2; | |
f6b6d58e | 320 | const Int_t kMaxBin=nzBins*nxBins; |
ae485645 | 321 | AliBin *bins[2]; |
93a93288 | 322 | AliBin *ddlbins[kHybridsPerDDL]; // 12 modules (=24 hybrids) of 1 DDL read "in parallel" |
323 | for(Int_t iHyb=0;iHyb<kHybridsPerDDL;iHyb++) ddlbins[iHyb]=new AliBin[kMaxBin]; | |
04366a57 | 324 | // read raw data input stream |
14dceddf | 325 | while (input->Next()) { |
e4da2288 | 326 | Int_t iModule = input->GetModuleID(); |
93a93288 | 327 | if(iModule<0){ |
328 | AliWarning(Form("Invalid SDD module number %d\n", iModule)); | |
329 | continue; | |
330 | } | |
44c5b268 | 331 | |
6297cc01 | 332 | Int_t iCarlos =input->GetCarlosId(); |
333 | Int_t iSide = input->GetChannel(); | |
14dceddf | 334 | Int_t iHybrid=iCarlos*2+iSide; |
335 | if (input->IsCompletedModule()) { | |
04366a57 | 336 | // when all data from a module was read, search for clusters |
44c5b268 | 337 | if(iCarlos<0){ |
338 | AliWarning(Form("Invalid SDD carlos number %d on module %d\n", iCarlos,iModule)); | |
339 | continue; | |
340 | } | |
341 | clusters[iModule] = new TClonesArray("AliITSRecPoint"); | |
342 | fModule = iModule; | |
343 | bins[0]=ddlbins[iCarlos*2]; // first hybrid of the completed module | |
344 | bins[1]=ddlbins[iCarlos*2+1]; // second hybrid of the completed module | |
345 | FindClustersSDD(bins, kMaxBin, nzBins, NULL, clusters[iModule]); | |
346 | Int_t nClusters = clusters[iModule]->GetEntriesFast(); | |
347 | nClustersSDD += nClusters; | |
348 | for(Int_t iBin=0;iBin<kMaxBin; iBin++){ | |
349 | ddlbins[iCarlos*2][iBin].Reset(); | |
350 | ddlbins[iCarlos*2+1][iBin].Reset(); | |
351 | } | |
5dfa68c5 | 352 | }else{ |
04366a57 | 353 | // fill the current digit into the bins array |
44c5b268 | 354 | if(iHybrid<0 || iHybrid>=kHybridsPerDDL){ |
355 | AliWarning(Form("Invalid SDD hybrid number %d on module %d\n", iHybrid,iModule)); | |
356 | continue; | |
357 | } | |
e4da2288 | 358 | AliITSCalibrationSDD* cal = (AliITSCalibrationSDD*)GetResp(iModule); |
9c7f5947 | 359 | if(cal==0){ |
360 | AliError(Form("Calibration object not present for SDD module %d\n",iModule)); | |
361 | continue; | |
362 | } | |
5683bd96 | 363 | Float_t charge=input->GetSignal(); |
14dceddf | 364 | Int_t chan=input->GetCoord1()+nAnodes*iSide; |
365 | Float_t gain=cal->GetChannelGain(chan); | |
20f3f947 | 366 | Float_t baseline = cal->GetBaseline(chan); |
367 | if(charge>baseline) charge-=baseline; | |
368 | else charge=0; | |
9c7f5947 | 369 | if(gain>0.){ // Bad channels have gain=0 |
370 | charge/=gain; | |
371 | if(charge>=cal->GetThresholdAnode(chan)) { | |
372 | Int_t q=(Int_t)(charge+0.5); | |
373 | Int_t iz = input->GetCoord1(); | |
374 | Int_t itb = input->GetCoord2(); | |
375 | Int_t index = (itb+1) * nzBins + (iz+1); | |
376 | if(index<kMaxBin){ | |
377 | ddlbins[iHybrid][index].SetQ(q); | |
378 | ddlbins[iHybrid][index].SetMask(1); | |
379 | ddlbins[iHybrid][index].SetIndex(index); | |
380 | }else{ | |
381 | AliWarning(Form("Invalid SDD cell: Anode=%d TimeBin=%d",iz,itb)); | |
382 | } | |
93a93288 | 383 | } |
5dfa68c5 | 384 | } |
04366a57 | 385 | } |
386 | } | |
93a93288 | 387 | for(Int_t iHyb=0;iHyb<kHybridsPerDDL;iHyb++) delete [] ddlbins[iHyb]; |
04366a57 | 388 | Info("FindClustersSDD", "found clusters in ITS SDD: %d", nClustersSDD); |
04366a57 | 389 | } |
390 | ||
e3ee4984 | 391 | //______________________________________________________________________ |
392 | void AliITSClusterFinderV2SDD::NoiseSuppress(Int_t k, Int_t sid,Int_t nzBins, AliBin* bins, AliITSCalibrationSDD* cal) const { | |
393 | // applies zero suppression using the measured noise of each anode | |
394 | // threshold values from ALICE-INT-1999-28 V10 | |
83ec5e27 | 395 | Float_t xfactL=2.2; |
396 | Float_t xfactH=4.0; | |
e3ee4984 | 397 | // |
398 | ||
399 | Int_t iAn=(k%nzBins)-1; | |
400 | if(iAn<0 || iAn>255) return; | |
401 | if(sid==1) iAn+=256; | |
402 | Int_t nLow=0, nHigh=0; | |
403 | Float_t noise=cal->GetNoiseAfterElectronics(iAn); | |
404 | Float_t noisem1=noise; | |
405 | if(iAn>1) noisem1=cal->GetNoiseAfterElectronics(iAn-1); | |
406 | Float_t noisep1=noise; | |
407 | if(iAn<511) noisep1=cal->GetNoiseAfterElectronics(iAn+1); | |
83ec5e27 | 408 | Float_t tL=noise*xfactL; |
409 | Float_t tH=noise*xfactH; | |
410 | Float_t tLp1=noisep1*xfactL; | |
411 | Float_t tHp1=noisep1*xfactH; | |
412 | Float_t tLm1=noisem1*xfactL; | |
413 | Float_t tHm1=noisem1*xfactH; | |
e3ee4984 | 414 | Float_t cC=bins[k].GetQ(); |
415 | if(cC<=tL){ | |
416 | bins[k].SetQ(0); | |
417 | bins[k].SetMask(0xFFFFFFFE); | |
418 | return; | |
419 | } | |
420 | nLow++; // cC is greater than tL | |
421 | if(cC>tH) nHigh++; | |
422 | Int_t sS=bins[k-1].GetQ(); | |
423 | if(sS>tLm1) nLow++; | |
424 | if(sS>tHm1) nHigh++; | |
425 | Int_t nN=bins[k+1].GetQ(); | |
426 | if(nN>tLp1) nLow++; | |
427 | if(nN>tHp1) nHigh++; | |
428 | Int_t eE=bins[k-nzBins].GetQ(); | |
429 | if(eE>tL) nLow++; | |
430 | if(eE>tH) nHigh++; | |
431 | Int_t wW=bins[k+nzBins].GetQ(); | |
432 | if(wW>tL) nLow++; | |
433 | if(wW>tH) nHigh++; | |
434 | if(nLow<3 || nHigh<1){ | |
435 | bins[k].SetQ(0); | |
436 | bins[k].SetMask(0xFFFFFFFE); | |
437 | } | |
438 | } | |
439 | ||
04366a57 | 440 | |
d97eb3fe | 441 |