change threshold comparison from > to >= ..
authordsilverm <dsilverm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 21 May 2010 12:09:55 +0000 (12:09 +0000)
committerdsilverm <dsilverm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 21 May 2010 12:09:55 +0000 (12:09 +0000)
EMCAL/AliCaloRawAnalyzer.cxx
EMCAL/AliCaloRawAnalyzerFastFit.cxx
EMCAL/AliCaloRawAnalyzerLMS.cxx
EMCAL/AliCaloRawAnalyzerPeakFinder.cxx
EMCAL/AliEMCALRawUtils.cxx

index 6ce22c5..2c68461 100644 (file)
@@ -416,13 +416,13 @@ AliCaloRawAnalyzer::PreFitEvaluateSamples( const vector<AliCaloBunchInfo>  &bunc
   index = SelectBunch( bunchvector,  &maxampindex,  &maxamp ); // select the bunch with the highest amplitude unless any time constraints is set
 
   
-  if( index >= 0 && maxamp > fAmpCut) // something valid was found, and non-zero amplitude
+  if( index >= 0 && maxamp >= fAmpCut) // something valid was found, and non-zero amplitude
     {
       // use more convenient numbering and possibly subtract pedestal
       ped  = ReverseAndSubtractPed( &(bunchvector.at(index)),  altrocfg1, altrocfg2, fReversed  );
       maxf = TMath::MaxElement( bunchvector.at(index).GetLength(),  fReversed );
       
-      if ( maxf > fAmpCut ) // possibly significant signal
+      if ( maxf >= fAmpCut ) // possibly significant signal
        {
          // select array around max to possibly be used in fit
          maxrev = maxampindex - bunchvector.at(index).GetStartBin(); 
index 48e2b5e..4a7d1d8 100644 (file)
@@ -69,7 +69,7 @@ AliCaloRawAnalyzerFastFit::Evaluate( const vector<AliCaloBunchInfo> &bunchvector
       Float_t maxf = TMath::MaxElement( bunchvector.at(index).GetLength(),  fReversed );
       int maxrev =  maxampindex -  bunchvector.at(index).GetStartBin();
       short timebinOffset = maxampindex - (bunchvector.at(index).GetLength()-1);
-      if ( maxf > fAmpCut )
+      if ( maxf >= fAmpCut )
        {
          SelectSubarray( fReversed,  bunchvector.at(index).GetLength(), maxrev , &first, &last);
          int nsamples =  last - first + 1;
index ddeadb9..455ec98 100644 (file)
@@ -90,7 +90,7 @@ AliCaloRawAnalyzerLMS::Evaluate( const vector<AliCaloBunchInfo>  &bunchvector, c
       short maxrev = maxampindex  -  bunchvector.at(index).GetStartBin();
       // timebinOffset is timebin value at maximum (maxrev)
       short timebinOffset = maxampindex - (bunchvector.at(index).GetLength()-1);
-      if ( maxf > fAmpCut )
+      if ( maxf >= fAmpCut )
        {
          SelectSubarray( fReversed,  bunchvector.at(index).GetLength(),  maxrev, &first, &last);
          int nsamples =  last - first + 1;
index 44ab803..beb6b3a 100644 (file)
@@ -141,7 +141,7 @@ AliCaloRawAnalyzerPeakFinder::Evaluate( const vector<AliCaloBunchInfo> &bunchvec
       short maxrev = maxampindex  -  bunchvector.at(index).GetStartBin();
       short timebinOffset = maxampindex - (bunchvector.at( index ).GetLength()-1); 
             
-      if(  maxf < fAmpCut  ||  ( maxamp - ped) > 900  )         // (maxamp - ped) > 900 = Close to saturation (use low gain then)
+      if(  maxf <= fAmpCut  ||  ( maxamp - ped) > 900  )        // (maxamp - ped) > 900 = Close to saturation (use low gain then)
        {
          return  AliCaloFitResults( maxamp, ped, AliCaloFitResults::kCrude, maxf, timebinOffset);
        }
@@ -149,7 +149,7 @@ AliCaloRawAnalyzerPeakFinder::Evaluate( const vector<AliCaloBunchInfo> &bunchvec
       int first;
       int last;
       
-      if ( maxf > fAmpCut )
+      if ( maxf >= fAmpCut )
        {         
          SelectSubarray( fReversed,  bunchvector.at(index).GetLength(), maxrev, &first, &last);
          int nsamples =  last - first;
index 149e7c1..6d12115 100644 (file)
@@ -422,7 +422,7 @@ void AliEMCALRawUtils::Raw2Digits(AliRawReader* reader,TClonesArray *digitsArr,
        //
        int nsamples = fRawAnalyzer->PreFitEvaluateSamples( bunchlist, in.GetAltroCFG1(), in.GetAltroCFG2(), bunchIndex, ampEstimate, maxADC, timeEstimate, pedEstimate, first, last); 
        
-       if (ampEstimate > fNoiseThreshold) { // something worth looking at
+       if (ampEstimate >= fNoiseThreshold) { // something worth looking at
          
          time = timeEstimate; // maxrev in AliCaloRawAnalyzer speak; comes with an offset w.r.t. real timebin
          Int_t timebinOffset = bunchlist.at(bunchIndex).GetStartBin() - (bunchlist.at(bunchIndex).GetLength()-1); 
@@ -450,7 +450,7 @@ void AliEMCALRawUtils::Raw2Digits(AliRawReader* reader,TClonesArray *digitsArr,
        } 
       } // fitDone
     
-      if (amp > fNoiseThreshold  && amp<fgkRawSignalOverflow) { // something to be stored
+      if (amp >= fNoiseThreshold  && amp<fgkRawSignalOverflow) { // something to be stored
        if ( ! fitDone) { // smear ADC with +- 0.5 uniform (avoid discrete effects)
          amp += (0.5 - gRandom->Rndm()); // Rndm generates a number in ]0,1]
        }