Changes to reduce noise in SDD reconstruction of cosmic ray data.
authormasera <masera@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 12 Feb 2009 14:40:46 +0000 (14:40 +0000)
committermasera <masera@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 12 Feb 2009 14:40:46 +0000 (14:40 +0000)
This noise reduction must be activated in the recoparam object.

Settings for best cosmic data reconstruction are (in rec.C):
itsRecoParam->SetUseSDDClusterSizeSelection();
itsRecoParam->SetMinClusterChargeSDD(40.);

F. Prino

ITS/AliITSClusterFinderV2SDD.cxx
ITS/AliITSRecoParam.cxx
ITS/AliITSRecoParam.h

index 84fcb96..9e6fbf9 100644 (file)
@@ -187,8 +187,8 @@ FindClustersSDD(AliBin* bins[2], Int_t nMaxBin, Int_t nzBins,
        }
        Int_t maxi=0,mini=0,maxj=0,minj=0;
 
-       for (Int_t di=-2; di<=2;di++){
-         for (Int_t dj=-3;dj<=3;dj++){
+       for (Int_t di=-5; di<=5;di++){
+         for (Int_t dj=-10;dj<=10;dj++){
            Int_t index = idx[k]+di+dj*nzBins;
            if (index<0) continue;
            if (index>=nMaxBin) continue;
@@ -215,6 +215,13 @@ FindClustersSDD(AliBin* bins[2], Int_t nMaxBin, Int_t nzBins,
            }
          }
        }
+       Int_t clSizAnode=maxi-mini+1;
+       Int_t clSizTb=maxj-minj+1;
+       if(repa->GetUseSDDClusterSizeSelection()){
+         if(clSizTb==1) continue; // cut common mode noise spikes
+         if(clSizAnode>3)  continue; // cut common mode noise spikes
+         if(clSizTb>10)  continue; // cut clusters on noisy anodes
+       }
 
        AliITSresponseSDD* rsdd = fDetTypeRec->GetResponseSDD();
        Float_t y=c.GetY(),z=c.GetZ(), q=c.GetQ();
@@ -242,8 +249,10 @@ FindClustersSDD(AliBin* bins[2], Int_t nMaxBin, Int_t nzBins,
 
        q/=rsdd->GetADC2keV();  //to have MPV 1 MIP = 86.4 KeV
        if(cal-> IsAMAt20MHz()) q*=2.; // account for 1/2 sampling freq.
+       if(q<repa->GetMinClusterChargeSDD()) continue; // remove noise clusters
+
        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]};
+       Int_t  info[3] = {clSizTb, clSizAnode, fNlayer[fModule]};
        if (digits) {     
          //       AliBin *b=&bins[s][idx[k]];
          //       AliITSdigitSDD* d=(AliITSdigitSDD*)digits->UncheckedAt(b->GetIndex());
@@ -257,7 +266,6 @@ FindClustersSDD(AliBin* bins[2], Int_t nMaxBin, Int_t nzBins,
          }
        }
        milab[3]=fNdet[fModule];
-
        AliITSRecPoint cc(milab,hit,info);
        cc.SetType(npeaks);
        cc.SetDriftTime(driftTimeUncorr);
index 176bbde..451607b 100644 (file)
@@ -139,6 +139,8 @@ fFindV0s(kTRUE),
 fUseUnfoldingInClusterFinderSPD(kFALSE),
 fUseUnfoldingInClusterFinderSDD(kTRUE),
 fUseUnfoldingInClusterFinderSSD(kTRUE),
+fUseSDDClusterSizeSelection(kFALSE),
+fMinClusterChargeSDD(0.),
 fUseChargeMatchingInClusterFinderSSD(kTRUE),
 fTrackleterOnlyOneTrackletPerC2(kTRUE),
 fTrackleterPhiWindow(0.08),
index 0c3c562..9e330eb 100644 (file)
@@ -269,6 +269,11 @@ class AliITSRecoParam : public AliDetectorRecoParam
   void   SetUseUnfoldingInClusterFinderSSD(Bool_t use=kTRUE) { fUseUnfoldingInClusterFinderSSD=use; return; }
   Bool_t GetUseUnfoldingInClusterFinderSSD() const { return fUseUnfoldingInClusterFinderSSD; }
 
+  void   SetUseSDDClusterSizeSelection(Bool_t use=kTRUE) {fUseSDDClusterSizeSelection=use;}
+  Bool_t GetUseSDDClusterSizeSelection() const {return fUseSDDClusterSizeSelection;}
+  void   SetMinClusterChargeSDD(Float_t qcut=0.){fMinClusterChargeSDD=qcut;}
+  Float_t GetMinClusterChargeSDD() const {return fMinClusterChargeSDD;}
+
   void   SetUseChargeMatchingInClusterFinderSSD(Bool_t use=kTRUE) { fUseChargeMatchingInClusterFinderSSD=use; return; }
   Bool_t GetUseChargeMatchingInClusterFinderSSD() const { return fUseChargeMatchingInClusterFinderSSD; }
 
@@ -452,6 +457,9 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Bool_t fUseUnfoldingInClusterFinderSDD; // SDD
   Bool_t fUseUnfoldingInClusterFinderSSD; // SSD
 
+  Bool_t  fUseSDDClusterSizeSelection; // cut on SDD cluster size
+  Float_t fMinClusterChargeSDD; // cut on SDD cluster charge
+
   Bool_t fUseChargeMatchingInClusterFinderSSD; // SSD
 
   // SPD Tracklets (D. Elia)
@@ -463,7 +471,8 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Float_t fTrackleterZetaOverlapCut;               // Fiducial window in eta for overlap cut
   Bool_t fUseCosmicRunShiftsSSD; // SSD time shifts for cosmic run 2007/2008 (use for data taken up to 18 sept 2008)
 
-  ClassDef(AliITSRecoParam,14) // ITS reco parameters
+
+  ClassDef(AliITSRecoParam,15) // ITS reco parameters
 };
 
 #endif