IteratedMedianCMNSubtractor Class Reference

#include <IteratedMedianCMNSubtractor.h>

Inheritance diagram for IteratedMedianCMNSubtractor:

Public Member Functions
void	init (const edm::EventSetup &es) override
void	subtract (uint32_t detId, uint16_t firstAPV, std::vector< int16_t > &digis) override
void	subtract (uint32_t detId, uint16_t firstAPV, std::vector< float > &digis) override
Public Member Functions inherited from SiStripCommonModeNoiseSubtractor
const std::vector< std::pair< short, float > > &	getAPVsCM () const
virtual	~SiStripCommonModeNoiseSubtractor ()

Private Member Functions
	IteratedMedianCMNSubtractor (double sigma, int iterations)
float	pairMedian (std::vector< std::pair< float, float > > &sample)
template<typename T >
void	subtract_ (uint32_t detId, uint16_t firstAPV, std::vector< T > &digis)

Private Attributes
double	cut_to_avoid_signal_
int	iterations_
uint32_t	noise_cache_id
edm::ESHandle< SiStripNoises >	noiseHandle
uint32_t	quality_cache_id
edm::ESHandle< SiStripQuality >	qualityHandle

Friends
class	SiStripRawProcessingFactory

Additional Inherited Members
Protected Member Functions inherited from SiStripCommonModeNoiseSubtractor
template<typename T >
float	median (std::vector< T > &)
	SiStripCommonModeNoiseSubtractor ()
Protected Attributes inherited from SiStripCommonModeNoiseSubtractor
std::vector< std::pair< short, float > >	_vmedians

Detailed Description

Definition at line 11 of file IteratedMedianCMNSubtractor.h.

Constructor & Destructor Documentation

IteratedMedianCMNSubtractor::IteratedMedianCMNSubtractor ( double  sigma, int  iterations  )

inlineprivate

Definition at line 24 of file IteratedMedianCMNSubtractor.h.

      : cut_to_avoid_signal_(sigma), iterations_(iterations), noise_cache_id(0), quality_cache_id(0){};

Member Function Documentation

void IteratedMedianCMNSubtractor::init ( const edm::EventSetup &  es )

overridevirtual

Reimplemented from SiStripCommonModeNoiseSubtractor.

Definition at line 9 of file IteratedMedianCMNSubtractor.cc.

References edm::EventSetup::get(), noise_cache_id, noiseHandle, quality_cache_id, and qualityHandle.

                                                              {
  uint32_t n_cache_id = es.get<SiStripNoisesRcd>().cacheIdentifier();
  uint32_t q_cache_id = es.get<SiStripQualityRcd>().cacheIdentifier();

  if (n_cache_id != noise_cache_id) {
    es.get<SiStripNoisesRcd>().get(noiseHandle);
    noise_cache_id = n_cache_id;
  }
  if (q_cache_id != quality_cache_id) {
    es.get<SiStripQualityRcd>().get(qualityHandle);
    quality_cache_id = q_cache_id;
  }
}

float IteratedMedianCMNSubtractor::pairMedian ( std::vector< std::pair< float, float > > &  sample )

inlineprivate

Definition at line 85 of file IteratedMedianCMNSubtractor.cc.

References mathSSE::return(), and simplePhotonAnalyzer_cfi::sample.

Referenced by subtract_().

                                                                                              {
  std::vector<std::pair<float, float> >::iterator mid = sample.begin() + sample.size() / 2;
  std::nth_element(sample.begin(), mid, sample.end());
  if (sample.size() & 1)  //odd size
    return (*mid).first;
  return ((*std::max_element(sample.begin(), mid)).first + (*mid).first) / 2.;
}

void IteratedMedianCMNSubtractor::subtract ( uint32_t  detId, uint16_t  firstAPV, std::vector< int16_t > &  digis  )

overridevirtual

Implements SiStripCommonModeNoiseSubtractor.

Definition at line 23 of file IteratedMedianCMNSubtractor.cc.

References subtract_().

                                                                                                       {
  subtract_(detId, firstAPV, digis);
}

void IteratedMedianCMNSubtractor::subtract ( uint32_t  detId, uint16_t  firstAPV, std::vector< float > &  digis  )

overridevirtual

Implements SiStripCommonModeNoiseSubtractor.

Definition at line 26 of file IteratedMedianCMNSubtractor.cc.

References subtract_().

                                                                                                     {
  subtract_(detId, firstAPV, digis);
}

template<typename T >

void IteratedMedianCMNSubtractor::subtract_ ( uint32_t  detId, uint16_t  firstAPV, std::vector< T > &  digis  )

inlineprivate

Definition at line 31 of file IteratedMedianCMNSubtractor.cc.

References SiStripCommonModeNoiseSubtractor::_vmedians, sistrip::APV, cut_to_avoid_signal_, SiStripNoises::getNoiseFast(), SiStripNoises::getRange(), SiStripBadStrip::getRange(), cuy::ii, SiStripQuality::IsStripBad(), iterations_, eostools::ls(), noiseHandle, hltrates_dqm_sourceclient-live_cfg::offset, pairMedian(), and qualityHandle.

Referenced by subtract().

                                                                                                         {
  SiStripNoises::Range detNoiseRange = noiseHandle->getRange(detId);
  SiStripQuality::Range detQualityRange = qualityHandle->getRange(detId);

  typename std::vector<T>::iterator fs, ls;
  float offset = 0;
  std::vector<std::pair<float, float> > subset;
  subset.reserve(128);

  _vmedians.clear();

  uint16_t APV = firstAPV;
  for (; APV < digis.size() / 128 + firstAPV; ++APV) {
    subset.clear();
    // fill subset vector with all good strips and their noises
    for (uint16_t istrip = APV * 128; istrip < (APV + 1) * 128; ++istrip) {
      if (!qualityHandle->IsStripBad(detQualityRange, istrip)) {
        std::pair<float, float> pin((float)digis[istrip - firstAPV * 128],
                                    (float)noiseHandle->getNoiseFast(istrip, detNoiseRange));
        subset.push_back(pin);
      }
    }

    // caluate offset for all good strips (first iteration)
    if (!subset.empty())
      offset = pairMedian(subset);

    // for second, third... iterations, remove strips over threshold
    // and recalculate offset on remaining strips
    for (int ii = 0; ii < iterations_ - 1; ++ii) {
      std::vector<std::pair<float, float> >::iterator si = subset.begin();
      while (si != subset.end()) {
        if (si->first - offset > cut_to_avoid_signal_ * si->second)
          si = subset.erase(si);
        else
          ++si;
      }
      if (subset.empty())
        break;
      offset = pairMedian(subset);
    }

    _vmedians.push_back(std::pair<short, float>(APV, offset));

    // remove offset
    fs = digis.begin() + (APV - firstAPV) * 128;
    ls = digis.begin() + (APV - firstAPV + 1) * 128;
    while (fs < ls) {
      *fs = static_cast<T>(*fs - offset);
      fs++;
    }
  }
}

Friends And Related Function Documentation

friend class SiStripRawProcessingFactory

friend

Definition at line 12 of file IteratedMedianCMNSubtractor.h.

Member Data Documentation

double IteratedMedianCMNSubtractor::cut_to_avoid_signal_

private

Definition at line 25 of file IteratedMedianCMNSubtractor.h.

Referenced by subtract_().

int IteratedMedianCMNSubtractor::iterations_

private

Definition at line 27 of file IteratedMedianCMNSubtractor.h.

Referenced by subtract_().

uint32_t IteratedMedianCMNSubtractor::noise_cache_id

private

Definition at line 30 of file IteratedMedianCMNSubtractor.h.

Referenced by init().

edm::ESHandle<SiStripNoises> IteratedMedianCMNSubtractor::noiseHandle

private

Definition at line 28 of file IteratedMedianCMNSubtractor.h.

Referenced by init(), and subtract_().

uint32_t IteratedMedianCMNSubtractor::quality_cache_id

private

Definition at line 30 of file IteratedMedianCMNSubtractor.h.

Referenced by init().

edm::ESHandle<SiStripQuality> IteratedMedianCMNSubtractor::qualityHandle

private

Definition at line 29 of file IteratedMedianCMNSubtractor.h.

Referenced by init(), and subtract_().