FastLinearCMNSubtractor.cc

Go to the documentation of this file.

#include "RecoLocalTracker/SiStripZeroSuppression/interface/FastLinearCMNSubtractor.h"

void FastLinearCMNSubtractor::subtract(uint32_t detId, uint16_t firstAPV, std::vector<int16_t>& digis){ subtract_(detId, firstAPV, digis);}
void FastLinearCMNSubtractor::subtract(uint32_t detId, uint16_t firstAPV, std::vector<float>& digis){ subtract_(detId, firstAPV, digis);}

template<typename T>
inline
void FastLinearCMNSubtractor::
subtract_(uint32_t detId, uint16_t firstAPV, std::vector<T>& digis){

  std::vector<T> tmp;  tmp.reserve(128);
  typename std::vector<T>::iterator
    strip( digis.begin() ),
    end(   digis.end()   ),
    endAPV, high, low;

  while( strip < end ) {
    endAPV = strip+128; tmp.clear();
    tmp.insert(tmp.end(),strip,endAPV);
    const float offset = median(tmp);

    low = strip;   high = strip+64;   tmp.clear();
    while( high < endAPV) tmp.push_back( *high++ - *low++ );
    const float slope = median(tmp)/64.;

    while (strip < endAPV) {
      *strip = static_cast<T>( *strip - (offset + slope*(65 - (endAPV-strip) ) ) );
      strip++;
    }

  }
}

// Details on http://abbaneo.web.cern.ch/abbaneo/cmode/cm.html