#include <IteratedMedianCMNSubtractor.h>

Inheritance diagram for IteratedMedianCMNSubtractor:

Public Member Functions
void	init (const edm::EventSetup &es)

void	subtract (const uint32_t &, const uint16_t &, std::vector< int16_t > &)

void	subtract (const uint32_t &, const uint16_t &, std::vector< float > &)

Public Member Functions inherited from SiStripCommonModeNoiseSubtractor
const std::vector< std::pair< short, float > > &	getAPVsCM ()

virtual	~SiStripCommonModeNoiseSubtractor ()

Private Member Functions
	IteratedMedianCMNSubtractor (double sigma, int iterations)

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

template<typename T >
void	subtract_ (const uint32_t &, const uint16_t &, std::vector< T > &)

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 26 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 )

virtual

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 90 of file IteratedMedianCMNSubtractor.cc.

References reco::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	(	const uint32_t &	detId,
		const uint16_t &	firstAPV,
		std::vector< int16_t > &	digis
	)

virtual

Implements SiStripCommonModeNoiseSubtractor.

Definition at line 23 of file IteratedMedianCMNSubtractor.cc.

References subtract_().

23 { subtract_(detId, firstAPV, digis);}

IteratedMedianCMNSubtractor::subtract_

void subtract_(const uint32_t &, const uint16_t &, std::vector< T > &)

Definition: IteratedMedianCMNSubtractor.cc:29

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

virtual

Implements SiStripCommonModeNoiseSubtractor.

Definition at line 24 of file IteratedMedianCMNSubtractor.cc.

References subtract_().

24 { subtract_(detId,firstAPV, digis);}

IteratedMedianCMNSubtractor::subtract_

void subtract_(const uint32_t &, const uint16_t &, std::vector< T > &)

Definition: IteratedMedianCMNSubtractor.cc:29

template<typename T >

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

inlineprivate

Definition at line 29 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, PFRecoTauDiscriminationByIsolation_cfi::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.size() != 0)
       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.size() == 0 ) 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++;
     }
 
   }
 }