SiStripFedZeroSuppression Class Reference

#include <SiStripFedZeroSuppression.h>

Public Member Functions
void	init (const edm::EventSetup &es)
	SiStripFedZeroSuppression (uint16_t fedalgo, edm::ConsumesCollector *iC=nullptr, bool trunc=true, bool trunc10bits=false)
void	suppress (const std::vector< SiStripDigi > &in, std::vector< SiStripDigi > &selectedSignal, uint32_t detId, const SiStripNoises &, const SiStripThreshold &)
void	suppress (const std::vector< SiStripDigi > &in, std::vector< SiStripDigi > &selectedSignal, uint32_t detId)
void	suppress (const edm::DetSet< SiStripRawDigi > &in, edm::DetSet< SiStripDigi > &out)
void	suppress (const std::vector< int16_t > &in, uint16_t firstAPV, edm::DetSet< SiStripDigi > &out)
uint16_t	truncate (int16_t adc) const
	~SiStripFedZeroSuppression ()

Private Member Functions
void	fillThresholds_ (const uint32_t detID, size_t size)
bool	isAValidDigi ()

Private Attributes
int16_t	adc
int16_t	adcMaxNeigh
int16_t	adcNext
int16_t	adcNext2
int16_t	adcPrev
int16_t	adcPrev2
bool	doTruncate
bool	doTruncate10bits
std::vector< int16_t >	highThr_
std::vector< float >	highThrSN_
std::vector< int16_t >	lowThr_
std::vector< float >	lowThrSN_
const SiStripNoises *	noise_
std::vector< float >	noises_
edm::ESGetToken< SiStripNoises, SiStripNoisesRcd >	noiseToken_
edm::ESWatcher< SiStripNoisesRcd >	noiseWatcher_
uint16_t	theFEDalgorithm
int16_t	theFEDhighThresh
int16_t	theFEDlowThresh
int16_t	theNeighFEDhighThresh
int16_t	theNeighFEDlowThresh
int16_t	theNext2FEDlowThresh
int16_t	theNextFEDhighThresh
int16_t	theNextFEDlowThresh
int16_t	thePrev2FEDlowThresh
int16_t	thePrevFEDhighThresh
int16_t	thePrevFEDlowThresh
const SiStripThreshold *	threshold_
edm::ESGetToken< SiStripThreshold, SiStripThresholdRcd >	thresholdToken_
edm::ESWatcher< SiStripThresholdRcd >	thresholdWatcher_

Friends
class	SiStripRawProcessingFactory

Detailed Description

Definition at line 20 of file SiStripFedZeroSuppression.h.

Constructor & Destructor Documentation

SiStripFedZeroSuppression()

SiStripFedZeroSuppression::SiStripFedZeroSuppression ( uint16_t  fedalgo, edm::ConsumesCollector *  iC = nullptr, bool  trunc = true, bool  trunc10bits = false  )

inline

Definition at line 24 of file SiStripFedZeroSuppression.h.

References noiseToken_.

      : noiseToken_{iC ? decltype(noiseToken_){iC->esConsumes<SiStripNoises, SiStripNoisesRcd>()}
                       : decltype(noiseToken_){}},
        thresholdToken_{iC ? decltype(thresholdToken_){iC->esConsumes<SiStripThreshold, SiStripThresholdRcd>()}
                           : decltype(thresholdToken_){}},
        theFEDalgorithm(fedalgo),
        doTruncate(trunc),
        doTruncate10bits(trunc10bits) {}

~SiStripFedZeroSuppression()

SiStripFedZeroSuppression::~SiStripFedZeroSuppression ( )

inline

Definition at line 35 of file SiStripFedZeroSuppression.h.

{};

Member Function Documentation

fillThresholds_()

void SiStripFedZeroSuppression::fillThresholds_ ( const uint32_t  detID, size_t  size  )

private

Definition at line 227 of file SiStripFedZeroSuppression.cc.

References MillePedeFileConverter_cfg::e, hgchebackDigitizer_cfi::noise, sistrip::extrainfo::noise_, findQualityFiles::size, and digitizers_cfi::strip.

                                                                                 {
  SiStripNoises::Range detNoiseRange = noise_->getRange(detID);
  SiStripThreshold::Range detThRange = threshold_->getRange(detID);

  if (highThr_.size() != size) {
    highThr_.resize(size);
    lowThr_.resize(size);
    noises_.resize(size);
    highThrSN_.resize(size);
    lowThrSN_.resize(size);
  }

  noise_->allNoises(noises_, detNoiseRange);
  threshold_->allThresholds(lowThrSN_, highThrSN_, detThRange);  // thresholds as S/N
  for (size_t strip = 0; strip < size; ++strip) {
    float noise = noises_[strip];
    //  uncomment line below to check bluk noise decoding
    //assert( noise == noiseHandle->getNoiseFast(strip,detNoiseRange) );
    highThr_[strip] = static_cast<int16_t>(highThrSN_[strip] * noise + 0.5 + 1e-6);
    lowThr_[strip] = static_cast<int16_t>(lowThrSN_[strip] * noise + 0.5 + 1e-6);
    // Note: it's a bit wierd, but there are some cases for which 'highThrSN_[strip]*noise' is an exact integer
    //   but due to roundoffs it gets rounded to the integer below if.
    //   Apparently the optimized code inlines differently and this changes the roundoff.
    //   The +1e-6 fixes the problem.   [GPetruc]
  }
}

init()

void SiStripFedZeroSuppression::init

(

const edm::EventSetup &

es

)

Definition at line 12 of file SiStripFedZeroSuppression.cc.

References edm::EventSetup::getData(), and sistrip::extrainfo::noise_.

                                                            {
  if (noiseWatcher_.check(es)) {
    noise_ = &es.getData(noiseToken_);
  }
  if (thresholdWatcher_.check(es)) {
    threshold_ = &es.getData(thresholdToken_);
  }
}

isAValidDigi()

bool SiStripFedZeroSuppression::isAValidDigi ( )

private

Definition at line 363 of file SiStripFedZeroSuppression.cc.

References accept(), gpuClustering::adc, edm::isDebugEnabled(), and LogTrace.

                                             {
#ifdef DEBUG_SiStripZeroSuppression_

  if (edm::isDebugEnabled()) {
    LogTrace("SiStripZeroSuppression") << "[SiStripFedZeroSuppression::suppress] "
                                       << "\n\t adc " << adc << "\n\t adcPrev " << adcPrev << "\n\t adcNext " << adcNext
                                       << "\n\t adcMaxNeigh " << adcMaxNeigh << "\n\t adcPrev2 " << adcPrev2
                                       << "\n\t adcNext2 " << adcNext2 << std::endl;

    LogTrace("SiStripZeroSuppression") << "[SiStripFedZeroSuppression::suppress] "
                                       << "\n\t theFEDlowThresh " << theFEDlowThresh << "\n\t theFEDhighThresh "
                                       << theFEDhighThresh << "\n\t thePrevFEDlowThresh " << thePrevFEDlowThresh
                                       << "\n\t thePrevFEDhighThresh " << thePrevFEDhighThresh
                                       << "\n\t theNextFEDlowThresh " << theNextFEDlowThresh
                                       << "\n\t theNextFEDhighThresh " << theNextFEDhighThresh
                                       << "\n\t theNeighFEDlowThresh " << theNeighFEDlowThresh
                                       << "\n\t theNeighFEDhighThresh " << theNeighFEDhighThresh
                                       << "\n\t thePrev2FEDlowThresh " << thePrev2FEDlowThresh
                                       << "\n\t theNext2FEDlowThresh " << theNext2FEDlowThresh << std::endl;
  }
#endif
  // Decide if this strip should be accepted.
  bool accept = false;
  switch (theFEDalgorithm) {
    case 1:
      accept = (adc >= theFEDlowThresh);
      break;
    case 2:
      accept = (adc >= theFEDhighThresh || (adc >= theFEDlowThresh && adcMaxNeigh >= theNeighFEDlowThresh));
      break;
    case 3:
      accept = (adc >= theFEDhighThresh || (adc >= theFEDlowThresh && adcMaxNeigh >= theNeighFEDhighThresh));
      break;
    case 4:
      accept =
          ((adc >= theFEDhighThresh)     //Test for adc>highThresh (same as algorithm 2)
           || ((adc >= theFEDlowThresh)  //Test for adc>lowThresh, with neighbour adc>lowThresh (same as algorithm 2)
               && (adcMaxNeigh >= theNeighFEDlowThresh)) ||
           ((adc < theFEDlowThresh)                 //Test for adc<lowThresh
            && (((adcPrev >= thePrevFEDhighThresh)  //with both neighbours>highThresh
                 && (adcNext >= theNextFEDhighThresh)) ||
                ((adcPrev >= thePrevFEDhighThresh)    //OR with previous neighbour>highThresh and
                 && (adcNext >= theNextFEDlowThresh)  //both the next neighbours>lowThresh
                 && (adcNext2 >= theNext2FEDlowThresh)) ||
                ((adcNext >= theNextFEDhighThresh)    //OR with next neighbour>highThresh and
                 && (adcPrev >= thePrevFEDlowThresh)  //both the previous neighbours>lowThresh
                 && (adcPrev2 >= thePrev2FEDlowThresh)) ||
                ((adcNext >= theNextFEDlowThresh)       //OR with both next neighbours>lowThresh and
                 && (adcNext2 >= theNext2FEDlowThresh)  //both the previous neighbours>lowThresh
                 && (adcPrev >= thePrevFEDlowThresh) && (adcPrev2 >= thePrev2FEDlowThresh)))));
      break;
    case 5:
      accept = adc > 0;
      break;
  }
  return accept;
}

suppress() [1/4]

void SiStripFedZeroSuppression::suppress	(	const std::vector< SiStripDigi > &	in,
		std::vector< SiStripDigi > &	selectedSignal,
		uint32_t	detId,
		const SiStripNoises &	noise,
		const SiStripThreshold &	threshold
	)

Definition at line 27 of file SiStripFedZeroSuppression.cc.

References gpuClustering::adc, SiStripThreshold::Data::getHth(), SiStripThreshold::Data::getLth(), mps_fire::i, recoMuon::in, hgchebackDigitizer_cfi::noise, digitizers_cfi::strip, remoteMonitoring_LASER_era2018_cfg::threshold, and particleFlowZeroSuppressionECAL_cff::thresholds.

Referenced by DigiSimLinkAlgorithm::run().

                                                                            {
  int i;
  int inSize = in.size();
  SiStripNoises::Range detNoiseRange = noise.getRange(detID);
  SiStripThreshold::Range detThRange = threshold.getRange(detID);

  // reserving more than needed, but quicker than one at a time
  selectedSignal.clear();
  selectedSignal.reserve(inSize);
  for (i = 0; i < inSize; i++) {
    //Find adc values for neighbouring strips
    const uint32_t strip = (uint32_t)in[i].strip();

    adc = in[i].adc();

    SiStripThreshold::Data thresholds = threshold.getData(strip, detThRange);
    theFEDlowThresh = static_cast<int16_t>(thresholds.getLth() * noise.getNoiseFast(strip, detNoiseRange) + 0.5);
    theFEDhighThresh = static_cast<int16_t>(thresholds.getHth() * noise.getNoiseFast(strip, detNoiseRange) + 0.5);

    adcPrev = -9999;
    adcNext = -9999;
    adcPrev2 = -9999;
    adcNext2 = -9999;

    /*
      Since we are not running on 
      Raw data we need to initialize
      all the FED threshold
    */

    theNextFEDlowThresh = theFEDlowThresh;
    theNext2FEDlowThresh = theFEDlowThresh;
    thePrevFEDlowThresh = theFEDlowThresh;
    thePrev2FEDlowThresh = theFEDlowThresh;
    theNeighFEDlowThresh = theFEDlowThresh;

    theNextFEDhighThresh = theFEDhighThresh;
    thePrevFEDhighThresh = theFEDhighThresh;
    theNeighFEDhighThresh = theFEDhighThresh;

    if (((strip) % 128) == 127) {
      adcNext = 0;
      theNextFEDlowThresh = 9999;
      theNextFEDhighThresh = 9999;
    } else if (i + 1 < inSize && in[i + 1].strip() == strip + 1) {
      adcNext = in[i + 1].adc();
      SiStripThreshold::Data thresholds_1 = threshold.getData(strip + 1, detThRange);
      theNextFEDlowThresh =
          static_cast<int16_t>(thresholds_1.getLth() * noise.getNoiseFast(strip + 1, detNoiseRange) + 0.5);
      theNextFEDhighThresh =
          static_cast<int16_t>(thresholds_1.getHth() * noise.getNoiseFast(strip + 1, detNoiseRange) + 0.5);
      if (((strip) % 128) == 126) {
        adcNext2 = 0;
        theNext2FEDlowThresh = 9999;
      } else if (i + 2 < inSize && in[i + 2].strip() == strip + 2) {
        adcNext2 = in[i + 2].adc();
        theNext2FEDlowThresh = static_cast<int16_t>(
            threshold.getData(strip + 2, detThRange).getLth() * noise.getNoiseFast(strip + 2, detNoiseRange) + 0.5);
      }
    }

    if (((strip) % 128) == 0) {
      adcPrev = 0;
      thePrevFEDlowThresh = 9999;
      thePrevFEDhighThresh = 9999;
    } else if (i - 1 >= 0 && in[i - 1].strip() == strip - 1) {
      adcPrev = in[i - 1].adc();
      SiStripThreshold::Data thresholds_1 = threshold.getData(strip - 1, detThRange);
      thePrevFEDlowThresh =
          static_cast<int16_t>(thresholds_1.getLth() * noise.getNoiseFast(strip - 1, detNoiseRange) + 0.5);
      thePrevFEDhighThresh =
          static_cast<int16_t>(thresholds_1.getHth() * noise.getNoiseFast(strip - 1, detNoiseRange) + 0.5);
      if (((strip) % 128) == 1) {
        adcPrev2 = 0;
        thePrev2FEDlowThresh = 9999;
      } else if (i - 2 >= 0 && in[i - 2].strip() == strip - 2) {
        adcPrev2 = in[i - 2].adc();
        thePrev2FEDlowThresh = static_cast<int16_t>(
            threshold.getData(strip - 2, detThRange).getLth() * noise.getNoiseFast(strip - 2, detNoiseRange) + 0.5);
      }
    }

    if (adcNext <= adcPrev) {
      adcMaxNeigh = adcPrev;
      theNeighFEDlowThresh = thePrevFEDlowThresh;
      theNeighFEDhighThresh = thePrevFEDhighThresh;
    } else {
      adcMaxNeigh = adcNext;
      theNeighFEDlowThresh = theNextFEDlowThresh;
      theNeighFEDhighThresh = theNextFEDhighThresh;
    }

    if (isAValidDigi()) {
      selectedSignal.push_back(SiStripDigi(strip, adc));
    }
  }
}

suppress() [2/4]

void SiStripFedZeroSuppression::suppress	(	const std::vector< SiStripDigi > &	in,
		std::vector< SiStripDigi > &	selectedSignal,
		uint32_t	detId
	)

Definition at line 21 of file SiStripFedZeroSuppression.cc.

References recoMuon::in, sistrip::extrainfo::noise_, and suppress.

                                                         {
  suppress(in, selectedSignal, detID, *noise_, *threshold_);
}

suppress() [3/4]

void SiStripFedZeroSuppression::suppress	(	const edm::DetSet< SiStripRawDigi > &	in,
		edm::DetSet< SiStripDigi > &	out
	)

Definition at line 129 of file SiStripFedZeroSuppression.cc.

References gpuClustering::adc, SiStripThreshold::Data::getHth(), SiStripThreshold::Data::getLth(), recoMuon::in, edm::isDebugEnabled(), LogTrace, sistrip::extrainfo::noise_, MillePedeFileConverter_cfg::out, digitizers_cfi::strip, particleFlowZeroSuppressionECAL_cff::thresholds, and HLT_2023v12_cff::truncate.

                                                                                                         {
  const uint32_t detID = out.id;
  SiStripNoises::Range detNoiseRange = noise_->getRange(detID);
  SiStripThreshold::Range detThRange = threshold_->getRange(detID);
#ifdef DEBUG_SiStripZeroSuppression_
  if (edm::isDebugEnabled())
    LogTrace("SiStripZeroSuppression")
        << "[SiStripFedZeroSuppression::suppress] Zero suppression on edm::DetSet<SiStripRawDigi>: detID " << detID
        << " size = " << in.data.size();
#endif
  edm::DetSet<SiStripRawDigi>::const_iterator in_iter = in.data.begin();
  for (; in_iter != in.data.end(); in_iter++) {
    const uint32_t strip = (uint32_t)(in_iter - in.data.begin());

#ifdef DEBUG_SiStripZeroSuppression_
    if (edm::isDebugEnabled())
      LogTrace("SiStripZeroSuppression") << "[SiStripFedZeroSuppression::suppress] detID= " << detID
                                         << " strip= " << strip << "  adc= " << in_iter->adc();
#endif
    adc = in_iter->adc();

    SiStripThreshold::Data thresholds = threshold_->getData(strip, detThRange);
    theFEDlowThresh = static_cast<int16_t>(thresholds.getLth() * noise_->getNoiseFast(strip, detNoiseRange) + 0.5);
    theFEDhighThresh = static_cast<int16_t>(thresholds.getHth() * noise_->getNoiseFast(strip, detNoiseRange) + 0.5);

    adcPrev = -9999;
    adcNext = -9999;
    /*
      If a strip is the last one on the chip
      set its next neighbor's thresholds to infinity
      because the FED does not merge clusters across
      chip boundaries right now
    */
    if (strip % 128 == 127) {
      adcNext = 0;
      theNextFEDlowThresh = 9999;
      theNextFEDhighThresh = 9999;
    } else {
      adcNext = (in_iter + 1)->adc();
      SiStripThreshold::Data thresholds_1 = threshold_->getData(strip + 1, detThRange);
      theNextFEDlowThresh =
          static_cast<int16_t>(thresholds_1.getLth() * noise_->getNoiseFast(strip + 1, detNoiseRange) + 0.5);
      theNextFEDhighThresh =
          static_cast<int16_t>(thresholds_1.getHth() * noise_->getNoiseFast(strip + 1, detNoiseRange) + 0.5);
    }
    /*
      Similarily, for the first strip 
      on a chip
    */
    if (strip % 128 == 0) {
      adcPrev = 0;
      thePrevFEDlowThresh = 9999;
      thePrevFEDhighThresh = 9999;
    } else {
      adcPrev = (in_iter - 1)->adc();
      SiStripThreshold::Data thresholds_1 = threshold_->getData(strip - 1, detThRange);
      thePrevFEDlowThresh =
          static_cast<int16_t>(thresholds_1.getLth() * noise_->getNoiseFast(strip - 1, detNoiseRange) + 0.5);
      thePrevFEDhighThresh =
          static_cast<int16_t>(thresholds_1.getHth() * noise_->getNoiseFast(strip - 1, detNoiseRange) + 0.5);
    }
    if (adcNext < adcPrev) {
      adcMaxNeigh = adcPrev;
      theNeighFEDlowThresh = thePrevFEDlowThresh;
      theNeighFEDhighThresh = thePrevFEDhighThresh;
    } else {
      adcMaxNeigh = adcNext;
      theNeighFEDlowThresh = theNextFEDlowThresh;
      theNeighFEDhighThresh = theNextFEDhighThresh;
    }

    //Find adc values for next neighbouring strips
    adcPrev2 = -9999;
    adcNext2 = -9999;
    thePrev2FEDlowThresh = 1;
    theNext2FEDlowThresh = 1;
    if (strip % 128 >= 126) {
      adcNext2 = 0;
      theNext2FEDlowThresh = 9999;
    } else if (strip % 128 < 126) {
      adcNext2 = (in_iter + 2)->adc();
      theNext2FEDlowThresh = static_cast<int16_t>(
          threshold_->getData(strip + 2, detThRange).getLth() * noise_->getNoiseFast(strip + 2, detNoiseRange) + 0.5);
    }
    if (strip % 128 <= 1) {
      adcPrev2 = 0;
      thePrev2FEDlowThresh = 9999;
    } else if (strip % 128 > 1) {
      adcPrev2 = (in_iter - 2)->adc();
      thePrev2FEDlowThresh = static_cast<int16_t>(
          threshold_->getData(strip - 2, detThRange).getLth() * noise_->getNoiseFast(strip - 2, detNoiseRange) + 0.5);
    }
    //GB 23/6/08: truncation should be done at the very beginning
    if (isAValidDigi())
      out.data.push_back(SiStripDigi(strip, truncate(in_iter->adc())));
  }
}

suppress() [4/4]

void SiStripFedZeroSuppression::suppress	(	const std::vector< int16_t > &	in,
		uint16_t	firstAPV,
		edm::DetSet< SiStripDigi > &	out
	)

Definition at line 254 of file SiStripFedZeroSuppression.cc.

References gpuClustering::adc, recoMuon::in, edm::isDebugEnabled(), LogTrace, MillePedeFileConverter_cfg::out, findQualityFiles::size, digitizers_cfi::strip, and HLT_2023v12_cff::truncate.

                                                                      {
  const uint32_t detID = out.id;
  size_t size = in.size();
#ifdef DEBUG_SiStripZeroSuppression_
  if (edm::isDebugEnabled())
    LogTrace("SiStripZeroSuppression")
        << "[SiStripFedZeroSuppression::suppress] Zero suppression on std::vector<int16_t>: detID " << detID
        << " size = " << in.size();
#endif

  fillThresholds_(detID, size + firstAPV * 128);  // want to decouple this from the other cost

  std::vector<int16_t>::const_iterator in_iter = in.begin();
  uint16_t strip = firstAPV * 128;
  for (; strip < size + firstAPV * 128; ++strip, ++in_iter) {
    size_t strip_mod_128 = strip & 127;
#ifdef DEBUG_SiStripZeroSuppression_
    if (edm::isDebugEnabled())
      LogTrace("SiStripZeroSuppression") << "[SiStripFedZeroSuppression::suppress]  detID= " << detID
                                         << " strip= " << strip << "  adc= " << *in_iter;
#endif
    adc = *in_iter;

    theFEDlowThresh = lowThr_[strip];
    theFEDhighThresh = highThr_[strip];

    //Find adc values for neighbouring strips

    /*
      If a strip is the last one on the chip
      set its next neighbor's thresholds to infinity
      because the FED does not merge clusters across
      chip boundaries right now
    */

    //adcPrev = -9999;  // useless, they are set
    //adcNext = -9999;  // in the next lines in any case
    if (strip_mod_128 == 127) {
      adcNext = 0;
      theNextFEDlowThresh = 9999;
      theNextFEDhighThresh = 9999;
    } else {
      adcNext = *(in_iter + 1);
      theNextFEDlowThresh = lowThr_[strip + 1];
      theNextFEDhighThresh = highThr_[strip + 1];
    }

    /*
      Similarily, for the first strip 
      on a chip
    */
    if (strip_mod_128 == 0) {
      adcPrev = 0;
      thePrevFEDlowThresh = 9999;
      thePrevFEDhighThresh = 9999;
    } else {
      adcPrev = *(in_iter - 1);
      thePrevFEDlowThresh = lowThr_[strip - 1];
      thePrevFEDhighThresh = highThr_[strip - 1];
    }

    if (adcNext < adcPrev) {
      adcMaxNeigh = adcPrev;
      theNeighFEDlowThresh = thePrevFEDlowThresh;
      theNeighFEDhighThresh = thePrevFEDhighThresh;
    } else {
      adcMaxNeigh = adcNext;
      theNeighFEDlowThresh = theNextFEDlowThresh;
      theNeighFEDhighThresh = theNextFEDhighThresh;
    }

    //Find adc values for next neighbouring strips
    //adcPrev2 = -9999;           //
    //adcNext2 = -9999;           // useless to set them here
    //thePrev2FEDlowThresh  = 1;  // they are overwritten always in the next 8 lines
    //theNext2FEDlowThresh  = 1;  //
    if (strip_mod_128 >= 126) {
      adcNext2 = 0;
      theNext2FEDlowThresh = 9999;
      //} else if ( strip_mod_128 < 126 ) { // if it's not >= then is <, no need to "if" again
    } else {
      adcNext2 = *(in_iter + 2);
      theNext2FEDlowThresh = lowThr_[strip + 2];
    }
    if (strip_mod_128 <= 1) {
      adcPrev2 = 0;
      thePrev2FEDlowThresh = 9999;
      //} else if ( strip_mod_128 > 1 ) { // same as above
    } else {
      adcPrev2 = *(in_iter - 2);
      thePrev2FEDlowThresh = lowThr_[strip - 2];
      ;
    }

    if (isAValidDigi()) {
#ifdef DEBUG_SiStripZeroSuppression_
      if (edm::isDebugEnabled())
        LogTrace("SiStripZeroSuppression")
            << "[SiStripFedZeroSuppression::suppress] DetId " << out.id << " strip " << strip << " adc " << *in_iter
            << " digiCollection size " << out.data.size();
#endif
      //GB 23/6/08: truncation should be done at the very beginning
      out.push_back(SiStripDigi(strip, (*in_iter < 0 ? 0 : truncate(*in_iter))));
    }
  }
}

truncate()

uint16_t SiStripFedZeroSuppression::truncate ( int16_t  adc ) const

inline

Definition at line 46 of file SiStripFedZeroSuppression.h.

References adc, doTruncate, and doTruncate10bits.

                                       {
    if (adc > 253 && doTruncate && !doTruncate10bits)
      return ((adc == 1023) ? 255 : 254);
    return adc;
  };

Friends And Related Function Documentation

SiStripRawProcessingFactory

friend class SiStripRawProcessingFactory

friend

Definition at line 21 of file SiStripFedZeroSuppression.h.

Member Data Documentation

adc

int16_t SiStripFedZeroSuppression::adc

private

Definition at line 68 of file SiStripFedZeroSuppression.h.

Referenced by truncate().

adcMaxNeigh

int16_t SiStripFedZeroSuppression::adcMaxNeigh

private

Definition at line 71 of file SiStripFedZeroSuppression.h.

adcNext

int16_t SiStripFedZeroSuppression::adcNext

private

Definition at line 70 of file SiStripFedZeroSuppression.h.

adcNext2

int16_t SiStripFedZeroSuppression::adcNext2

private

Definition at line 73 of file SiStripFedZeroSuppression.h.

adcPrev

int16_t SiStripFedZeroSuppression::adcPrev

private

Definition at line 69 of file SiStripFedZeroSuppression.h.

adcPrev2

int16_t SiStripFedZeroSuppression::adcPrev2

private

Definition at line 72 of file SiStripFedZeroSuppression.h.

doTruncate

bool SiStripFedZeroSuppression::doTruncate

private

Definition at line 63 of file SiStripFedZeroSuppression.h.

Referenced by truncate().

doTruncate10bits

bool SiStripFedZeroSuppression::doTruncate10bits

private

Definition at line 64 of file SiStripFedZeroSuppression.h.

Referenced by truncate().

highThr_

std::vector<int16_t> SiStripFedZeroSuppression::highThr_

private

Definition at line 87 of file SiStripFedZeroSuppression.h.

highThrSN_

std::vector<float> SiStripFedZeroSuppression::highThrSN_

private

Definition at line 88 of file SiStripFedZeroSuppression.h.

lowThr_

std::vector<int16_t> SiStripFedZeroSuppression::lowThr_

private

Definition at line 87 of file SiStripFedZeroSuppression.h.

lowThrSN_

std::vector<float> SiStripFedZeroSuppression::lowThrSN_

private

Definition at line 88 of file SiStripFedZeroSuppression.h.

noise_

const SiStripNoises* SiStripFedZeroSuppression::noise_

private

Definition at line 55 of file SiStripFedZeroSuppression.h.

noises_

std::vector<float> SiStripFedZeroSuppression::noises_

private

Definition at line 89 of file SiStripFedZeroSuppression.h.

noiseToken_

edm::ESGetToken<SiStripNoises, SiStripNoisesRcd> SiStripFedZeroSuppression::noiseToken_

private

Definition at line 50 of file SiStripFedZeroSuppression.h.

Referenced by SiStripFedZeroSuppression().

noiseWatcher_

edm::ESWatcher<SiStripNoisesRcd> SiStripFedZeroSuppression::noiseWatcher_

private

Definition at line 57 of file SiStripFedZeroSuppression.h.

theFEDalgorithm

uint16_t SiStripFedZeroSuppression::theFEDalgorithm

private

Definition at line 60 of file SiStripFedZeroSuppression.h.

theFEDhighThresh

int16_t SiStripFedZeroSuppression::theFEDhighThresh

private

Definition at line 66 of file SiStripFedZeroSuppression.h.

theFEDlowThresh

int16_t SiStripFedZeroSuppression::theFEDlowThresh

private

Definition at line 65 of file SiStripFedZeroSuppression.h.

theNeighFEDhighThresh

int16_t SiStripFedZeroSuppression::theNeighFEDhighThresh

private

Definition at line 81 of file SiStripFedZeroSuppression.h.

theNeighFEDlowThresh

int16_t SiStripFedZeroSuppression::theNeighFEDlowThresh

private

Definition at line 80 of file SiStripFedZeroSuppression.h.

theNext2FEDlowThresh

int16_t SiStripFedZeroSuppression::theNext2FEDlowThresh

private

Definition at line 84 of file SiStripFedZeroSuppression.h.

theNextFEDhighThresh

int16_t SiStripFedZeroSuppression::theNextFEDhighThresh

private

Definition at line 78 of file SiStripFedZeroSuppression.h.

theNextFEDlowThresh

int16_t SiStripFedZeroSuppression::theNextFEDlowThresh

private

Definition at line 77 of file SiStripFedZeroSuppression.h.

thePrev2FEDlowThresh

int16_t SiStripFedZeroSuppression::thePrev2FEDlowThresh

private

Definition at line 83 of file SiStripFedZeroSuppression.h.

thePrevFEDhighThresh

int16_t SiStripFedZeroSuppression::thePrevFEDhighThresh

private

Definition at line 76 of file SiStripFedZeroSuppression.h.

thePrevFEDlowThresh

int16_t SiStripFedZeroSuppression::thePrevFEDlowThresh

private

Definition at line 75 of file SiStripFedZeroSuppression.h.

threshold_

const SiStripThreshold* SiStripFedZeroSuppression::threshold_

private

Definition at line 56 of file SiStripFedZeroSuppression.h.

thresholdToken_

edm::ESGetToken<SiStripThreshold, SiStripThresholdRcd> SiStripFedZeroSuppression::thresholdToken_

private

Definition at line 54 of file SiStripFedZeroSuppression.h.

thresholdWatcher_

edm::ESWatcher<SiStripThresholdRcd> SiStripFedZeroSuppression::thresholdWatcher_

private

Definition at line 58 of file SiStripFedZeroSuppression.h.