#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, edm::ConsumesCollector iC)

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_

const SiStripNoises *	noiseHandle

edm::ESGetToken< SiStripNoises, SiStripNoisesRcd >	noiseToken_

edm::ESWatcher< SiStripNoisesRcd >	noiseWatcher_

const SiStripQuality *	qualityHandle

edm::ESGetToken< SiStripQuality, SiStripQualityRcd >	qualityToken_

edm::ESWatcher< SiStripQualityRcd >	qualityWatcher_

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 14 of file IteratedMedianCMNSubtractor.h.

Constructor & Destructor Documentation

◆ IteratedMedianCMNSubtractor()

IteratedMedianCMNSubtractor::IteratedMedianCMNSubtractor	(	double	sigma,
		int	iterations,
		edm::ConsumesCollector	iC
	)

inlineprivate

Definition at line 27 of file IteratedMedianCMNSubtractor.h.

       : cut_to_avoid_signal_(sigma),
         iterations_(iterations),
         noiseToken_(iC.esConsumes<SiStripNoises, SiStripNoisesRcd>()),
         qualityToken_(iC.esConsumes<SiStripQuality, SiStripQualityRcd>()) {}

Member Function Documentation

◆ init()

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

overridevirtual

Reimplemented from SiStripCommonModeNoiseSubtractor.

Definition at line 7 of file IteratedMedianCMNSubtractor.cc.

References edm::ESWatcher< T >::check(), edm::EventSetup::getData(), noiseHandle, noiseToken_, noiseWatcher_, qualityHandle, qualityToken_, and qualityWatcher_.

                                                               {
   if (noiseWatcher_.check(es)) {
     noiseHandle = &es.getData(noiseToken_);
   }
   if (qualityWatcher_.check(es)) {
     qualityHandle = &es.getData(qualityToken_);
   }
 }

◆ pairMedian()

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

inlineprivate

Definition at line 78 of file IteratedMedianCMNSubtractor.cc.

References mathSSE::return(), and ecalGpuTask_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.;
 }

◆ subtract() [1/2]

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

overridevirtual

Implements SiStripCommonModeNoiseSubtractor.

Definition at line 16 of file IteratedMedianCMNSubtractor.cc.

References hcalRecHitTable_cff::detId, and subtract_().

                                                                                                        {
   subtract_(detId, firstAPV, digis);
 }

◆ subtract() [2/2]

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

overridevirtual

Implements SiStripCommonModeNoiseSubtractor.

Definition at line 19 of file IteratedMedianCMNSubtractor.cc.

References hcalRecHitTable_cff::detId, and subtract_().

                                                                                                      {
   subtract_(detId, firstAPV, digis);
 }

◆ subtract_()

template<typename T >

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

inlineprivate

Definition at line 24 of file IteratedMedianCMNSubtractor.cc.

References SiStripCommonModeNoiseSubtractor::_vmedians, sistrip::APV, cut_to_avoid_signal_, hcalRecHitTable_cff::detId, make_classfiles::fs, SiStripNoises::getNoiseFast(), SiStripNoises::getRange(), SiStripBadStrip::getRange(), cuy::ii, SiStripQuality::IsStripBad(), iterations_, relval_nano::ls, noiseHandle, HLT_IsoTrack_cff::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++;
     }
   }
 }