TT6CMNSubtractor Class Reference

#include <TT6CMNSubtractor.h>

Inheritance diagram for TT6CMNSubtractor:

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

Private Member Functions
template<typename T >
void	subtract_ (uint32_t detId, uint16_t firstAPV, std::vector< T > &digis)
	TT6CMNSubtractor (double in)

Private Attributes
double	cut_to_avoid_signal_
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 9 of file TT6CMNSubtractor.h.

Constructor & Destructor Documentation

TT6CMNSubtractor()

TT6CMNSubtractor::TT6CMNSubtractor ( double  in )

inlineprivate

Definition at line 21 of file TT6CMNSubtractor.h.

: cut_to_avoid_signal_(in), noise_cache_id(0), quality_cache_id(0){};

Member Function Documentation

init()

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

overridevirtual

Reimplemented from SiStripCommonModeNoiseSubtractor.

Definition at line 9 of file TT6CMNSubtractor.cc.

                                                   {
  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;
  }
}

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

subtract() [1/2]

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

overridevirtual

Implements SiStripCommonModeNoiseSubtractor.

Definition at line 26 of file TT6CMNSubtractor.cc.

                                                                                          {
  subtract_(detId, firstAPV, digis);
}

References subtract_().

subtract() [2/2]

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

overridevirtual

Implements SiStripCommonModeNoiseSubtractor.

Definition at line 23 of file TT6CMNSubtractor.cc.

                                                                                            {
  subtract_(detId, firstAPV, digis);
}

References subtract_().

subtract_()

template<typename T >

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

inlineprivate

Definition at line 31 of file TT6CMNSubtractor.cc.

                                                                                              {
  short FixedBias = 128;
  SiStripNoises::Range detNoiseRange = noiseHandle->getRange(detId);
  SiStripQuality::Range detQualityRange = qualityHandle->getRange(detId);
 
  typename std::vector<T>::iterator fs, ls;
  double sumVal = 0.0;
  double sumWt = 0.0;
 
  for (uint16_t istrip = 0; istrip < digis.size(); ++istrip) {
    if (!qualityHandle->IsStripBad(detQualityRange, istrip + firstAPV * 128)) {
      float stripNoise = noiseHandle->getNoiseFast(istrip + firstAPV * 128, detNoiseRange);
 
      if (std::abs(digis[istrip] - FixedBias) < cut_to_avoid_signal_ * stripNoise) {
        double nWeight = 1 / (stripNoise * stripNoise);
        sumVal += (digis[istrip] - FixedBias) * nWeight;
        sumWt += nWeight;
      }
    }
 
    if (istrip % 128 == 127) {
      double CM = (sumWt) ? sumVal / sumWt : 0.0;
 
      fs = digis.begin() + istrip - 127;
      ls = digis.begin() + istrip + 1;
 
      while (fs < ls) {
        *fs = static_cast<T>(*fs - FixedBias - CM);
        fs++;
      }
 
      sumVal = 0.0;
      sumWt = 0.0;
    }
  }
}

References funct::abs(), cut_to_avoid_signal_, SiStripNoises::getNoiseFast(), SiStripNoises::getRange(), SiStripBadStrip::getRange(), SiStripQuality::IsStripBad(), eostools::ls(), noiseHandle, and qualityHandle.

Referenced by subtract().

Friends And Related Function Documentation

SiStripRawProcessingFactory

friend class SiStripRawProcessingFactory

friend

Definition at line 10 of file TT6CMNSubtractor.h.

Member Data Documentation

cut_to_avoid_signal_

double TT6CMNSubtractor::cut_to_avoid_signal_

private

Definition at line 21 of file TT6CMNSubtractor.h.

Referenced by subtract_().

noise_cache_id

uint32_t TT6CMNSubtractor::noise_cache_id

private

Definition at line 25 of file TT6CMNSubtractor.h.

Referenced by init().

noiseHandle

edm::ESHandle<SiStripNoises> TT6CMNSubtractor::noiseHandle

private

Definition at line 23 of file TT6CMNSubtractor.h.

Referenced by init(), and subtract_().

quality_cache_id

uint32_t TT6CMNSubtractor::quality_cache_id

private

Definition at line 25 of file TT6CMNSubtractor.h.

Referenced by init().

qualityHandle

edm::ESHandle<SiStripQuality> TT6CMNSubtractor::qualityHandle

private

Definition at line 24 of file TT6CMNSubtractor.h.

Referenced by init(), and subtract_().