- */
-
- Float_t y=c.GetY(),z=c.GetZ(), q=c.GetQ();
- y/=q; z/=q;
- //
- //Float_t s2 = c.GetSigmaY2()/c.GetQ() - y*y;
- // c.SetSigmaY2(s2);
- //s2 = c.GetSigmaZ2()/c.GetQ() - z*z;
- //c.SetSigmaZ2(s2);
- //
-
- Float_t yyyy = y;
- //y=(y-0.5)*fYpitchSDD;
- //y-=fHwSDD;
- //y-=fYoffSDD; //delay ?
- //if (s) y=-y;
-
- z=(z-0.5)*fZpitchSDD;
- z-=fHlSDD;
- Float_t zdet = z;
- // y=-(-y+fYshift[fModule]);
- // z= -z+fZshift[fModule];
- // c.SetY(y);
- // c.SetZ(z);
- Float_t xdet = cal->GetDriftPath((yyyy-0.5)*25,0);
- xdet=xdet/10000.-fHwSDD-fYoffSDD;
- if (s) xdet=-xdet;
-
- CorrectPosition(zdet,xdet);
-
- y=-(-xdet+fYshift[fModule]);
- z= -zdet+fZshift[fModule];
-
- q/=12.7; //this WAS consistent with SSD. To be reassessed
- // 23-MAR-2007
- Float_t hit[5] = {y, z, 0.0030*0.0030, 0.0020*0.0020, q};
- Int_t info[3] = {maxj-minj+1, maxi-mini+1, fNlayer[fModule]};
-
- //if (c.GetQ() < 20.) continue; //noise cluster
-
- if (digits) {
- // AliBin *b=&bins[s][idx[k]];
- // AliITSdigitSDD* d=(AliITSdigitSDD*)digits->UncheckedAt(b->GetIndex());
- {
- //Int_t lab[3];
- //lab[0]=(d->GetTracks())[0];
- //lab[1]=(d->GetTracks())[1];
- //lab[2]=(d->GetTracks())[2];
- //CheckLabels(lab);
- CheckLabels2(milab);
- }
- }
- milab[3]=fNdet[fModule];
-
- AliITSRecPoint cc(milab,hit,info);
- cc.SetType(npeaks);
-
- if(clusters) new (cl[ncl]) AliITSRecPoint(cc);
- else {
- fDetTypeRec->AddRecPoint(cc);
- }
- ncl++;
+ }
+ else { // raw data
+ if (fRawID2ClusID) milab[0] = fNClusters+1; // RS: store clID+1 as a reference to the cluster
+ }
+
+
+ Int_t clSizAnode=fZmax-fZmin+1;
+ Int_t clSizTb=fXmax-fXmin+1;
+ if(repa->GetUseSDDClusterSizeSelection()){
+ if(clSizTb==1) continue; // cut common mode noise spikes
+ if(clSizAnode>5) continue; // cut common mode noise spikes
+ if(clSizTb>10) continue; // cut clusters on noisy anodes
+ if(cal-> IsAMAt20MHz() && clSizTb>8) continue; // cut clusters on noisy anodes
+ }
+
+ AliITSresponseSDD* rsdd = fDetTypeRec->GetResponseSDD();
+ Float_t y=c.GetY(),z=c.GetZ(), q=c.GetQ();
+ y/=q; z/=q;
+ Float_t zAnode=z-0.5; // to have anode in range 0.-255. and centered on the mid of the pitch
+ Float_t timebin=y-0.5; // to have time bin in range 0.-255. amd centered on the mid of the bin
+ if(s==1) zAnode += GetSeg()->NpzHalf(); // right side has anodes from 256. to 511.
+ Float_t zdet = GetSeg()->GetLocalZFromAnode(zAnode);
+ Float_t driftTimeUncorr = GetSeg()->GetDriftTimeFromTb(timebin)+jitter*rsdd->GetCarlosRXClockPeriod();
+ Float_t driftTime=driftTimeUncorr-rsdd->GetTimeZero(fModule);
+ if(driftTime<fMaxDrTimeForTightCut && maxADC<fCutOnPeakTight) continue;
+
+ Float_t driftSpeed = cal->GetDriftSpeedAtAnode(zAnode) + rsdd->GetDeltaVDrift(fModule,zAnode>255);
+ Float_t driftPathMicron = driftTime*driftSpeed;
+ const Double_t kMicronTocm = 1.0e-4;
+ Float_t xdet=(driftPathMicron-GetSeg()->Dx())*kMicronTocm; // xdet is negative
+ if (s==0) xdet=-xdet; // left side has positive local x
+
+ if(repa->GetUseSDDCorrectionMaps()){
+ Float_t corrx=0, corrz=0;
+ cal->GetCorrections(zdet,xdet,corrz,corrx,GetSeg());
+ zdet+=corrz;
+ xdet+=corrx;
+ }
+
+ Double_t loc[3]={xdet,0.,zdet},trk[3]={0.,0.,0.};
+ mT2L->MasterToLocal(loc,trk);
+ y=trk[1];
+ z=trk[2];
+
+ q+=(driftTime*rsdd->GetADCvsDriftTime(fModule)); // correction for zero supp.
+ q/=rsdd->GetADCtokeV(fModule);
+ if(cal-> IsAMAt20MHz()) q*=2.; // account for 1/2 sampling freq.
+ if(q<repa->GetMinClusterChargeSDD()) continue; // remove noise clusters
+
+ Float_t hit[6] = {y, z, 0.0030*0.0030, 0.0020*0.0020, q, 0.};
+ Int_t info[3] = {clSizTb, clSizAnode, fNlayer[fModule]};
+ if (digits) CheckLabels2(milab);
+ milab[3]=fNdet[fModule];
+ AliITSRecPoint cc(milab,hit,info);
+ cc.SetType(nClust*100+npeaks);
+ cc.SetDriftTime(driftTimeUncorr);
+ cc.SetDriftSide(s);
+ cc.SetChargeRatio(maxADC);
+ if(clusters) new (cl[ncl]) AliITSRecPoint(cc);
+ else {
+ fDetTypeRec->AddRecPoint(cc);
+ }
+ fNClusters++; // RS
+ ncl++;
+ }