---

SiStripMedianCommonModeNoiseSubtraction.cc

Go to the documentation of this file.

#include "CommonTools/SiStripZeroSuppression/interface/SiStripMedianCommonModeNoiseSubtraction.h"
#include <sstream>

void SiStripMedianCommonModeNoiseSubtraction::subtract(const uint32_t& detId,std::vector<int16_t>& digis){
  
#ifdef DEBUG_SiStripZeroSuppression_
  LogDebug("SiStripZeroSuppression") << "[SiStripMedianCommonModeNoiseSubtraction::subtract] digis.size()= " << digis.size();
#endif
  std::vector<int16_t> APVdigis;
  std::vector<int16_t>::iterator fs;
  std::vector<int16_t>::iterator ls;
  float CM;
  APVdigis.reserve(128);
  int nAPV = digis.size()/128;
  std::vector<int16_t>::iterator mid = APVdigis.begin() + 64;
  
#ifdef DEBUG_SiStripZeroSuppression_
  if (edm::isDebugEnabled())
    LogDebug("SiStripZeroSuppression") << "[SiStripMedianCommonModeNoiseSubtraction::subtract] DetId " << detId << " number of apvs: nAPV= " << nAPV;
#endif
  for (int iAPV=0; iAPV<nAPV; iAPV++){
    APVdigis.clear(); //added verify
    fs = digis.begin()+iAPV*128;
    ls = digis.begin()+(iAPV+1)*128;
    APVdigis.insert(APVdigis.end(), fs, ls );

    // std::sort(APVdigis.begin(),APVdigis.end());
    // CM = (APVdigis[63]+APVdigis[64])/2.;
    std::nth_element(APVdigis.begin(), mid, APVdigis.end());
    std::vector<int16_t>::iterator i63 = std::max_element(APVdigis.begin(), mid);
    CM = (*mid + *i63)/2.;

#ifdef DEBUG_SiStripZeroSuppression_
    std::stringstream ss;
    if (edm::isDebugEnabled())
      ss << "[SiStripMedianCommonModeNoiseSubtraction::subtract] DetId " << detId << " iApv= " <<iAPV << " CM= " << CM << std::endl;
#endif    

    while (fs < ls) {
      *fs = (int16_t) (*fs-CM);
#ifdef DEBUG_SiStripZeroSuppression_
      if (edm::isDebugEnabled())
        ss << "[SiStripMedianCommonModeNoiseSubtraction::subtract] DetId " << detId << " strip " << fs-digis.begin() << " adc CM subtr " << *fs << std::endl;
#endif
      fs++;
    }
#ifdef DEBUG_SiStripZeroSuppression_
    LogDebug("SiStripZeroSuppression") << ss.str().c_str();
#endif
  }
}

---