avoid overflow signals/amplitudes
authordsilverm <dsilverm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 9 Mar 2010 17:47:22 +0000 (17:47 +0000)
committerdsilverm <dsilverm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 9 Mar 2010 17:47:22 +0000 (17:47 +0000)
EMCAL/AliEMCALRawUtils.cxx

index 1f69fef..537af71 100644 (file)
@@ -436,7 +436,7 @@ void AliEMCALRawUtils::Raw2Digits(AliRawReader* reader,TClonesArray *digitsArr,
        } // ampEstimate check
       } // method selection
     
-      if (amp > fNoiseThreshold) { // something to be stored
+      if (amp > fNoiseThreshold  && amp<fgkRawSignalOverflow) { // something to be stored
        Int_t id =  fGeom->GetAbsCellIdFromCellIndexes(in.GetModule(), in.GetRow(), in.GetColumn()) ;
        lowGain = in.IsLowGain();
 
@@ -680,6 +680,7 @@ void AliEMCALRawUtils::FitParabola(const TGraph *gSig, Float_t & amp) const
       Double_t a  = (sy-b*sx-c*sx2)/kN  ;
       Double_t xmax = -b/(2*c) ; 
       ymax = a + b*xmax + c*xmax*xmax ;//<========== This is the maximum amplitude
+      amp = ymax;
     }
   }