updates concerning noise reduction
authormasera <massimo.masera@cern.ch>
Wed, 15 Oct 2014 09:56:19 +0000 (11:56 +0200)
committermasera <massimo.masera@cern.ch>
Wed, 15 Oct 2014 09:56:19 +0000 (11:56 +0200)
ITS/AliITSClusterFinderSDDfast.cxx

index 29e9714..67e822d 100644 (file)
@@ -338,6 +338,7 @@ void AliITSClusterFinderSDDfast::FindClustersSDD(std::vector<int>& bins0, std::v
       Float_t driftTime=driftTimeUncorr-rsdd->GetTimeZero(fModule);
 
       if(driftTime<fMaxDrTimeForTightCut && recp[i][10]<fCutOnPeakTight) continue;
+      if(recp[i][10]<fCutOnPeakLoose) continue;
 
       Float_t driftSpeed=cal->GetDriftSpeedAtAnode(zAnode) + rsdd->GetDeltaVDrift(fModule,zAnode>255);
       Float_t driftPathMicron=driftTime*driftSpeed;
@@ -521,15 +522,18 @@ void AliITSClusterFinderSDDfast::FindClustersSDD(AliITSRawStream* input){
           Int_t q=(Int_t)(charge+0.5);
           Int_t iz = input->GetCoord1();
           Int_t itb = input->GetCoord2();
-          Float_t noise=cal->GetNoiseAfterElectronics(iz)*2.2;  // applies zero suppression using the measured noise of each anode. Threshold values from ALICE-INT-1999-28 V10
+         Int_t ian=iz;
+         if(iSide==1)
+           ian+=256;
+          Float_t noise=cal->GetNoiseAfterElectronics(ian)*2.2;  // applies zero suppression using the measured noise of each anode. Threshold values from ALICE-INT-1999-28 V10
           Int_t index=(itb+1)*fNZbins+(iz+1);
           if((itb<fNTimeBins) && (iz<fNAnodes)){
             if(q<noise) continue;
-            fDDLBins[iHybrid].push_back(index);
-            fDDLBins[iHybrid].push_back(q);
-            fDDLBins[iHybrid].push_back(0);
-            fDDLBins[iHybrid].push_back(countRW);
-            mapsDDL[iHybrid][index/32]|=(1<<(index%32));
+           fDDLBins[iHybrid].push_back(index);
+           fDDLBins[iHybrid].push_back(q);
+           fDDLBins[iHybrid].push_back(0);
+           fDDLBins[iHybrid].push_back(countRW);
+           mapsDDL[iHybrid][index/32]|=(1<<(index%32));
           }
          else{
             AliWarning(Form("Invalid SDD cell: Anode=%d   TimeBin=%d",iz,itb));