SiStripFedZeroSuppression.h

Go to the documentation of this file.

#ifndef RECOLOCALTRACKER_SISTRIPZEROSUPPRESSION_SISTRIPZEROSUPPRESSOR_H
#define RECOLOCALTRACKER_SISTRIPZEROSUPPRESSION_SISTRIPZEROSUPPRESSOR_H

#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/SiStripDigi/interface/SiStripDigi.h"
#include "DataFormats/SiStripDigi/interface/SiStripRawDigi.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"

#include <vector>
class SiStripNoises;
class SiStripThreshold;

class SiStripFedZeroSuppression {
  friend class SiStripRawProcessingFactory;

public:
  SiStripFedZeroSuppression(uint16_t fedalgo, bool trunc = true, bool trunc10bits = false)
      : noise_cache_id(0),
        threshold_cache_id(0),
        theFEDalgorithm(fedalgo),
        doTruncate(trunc),
        doTruncate10bits(trunc10bits) {}
  ~SiStripFedZeroSuppression(){};
  void init(const edm::EventSetup& es);
  void suppress(const std::vector<SiStripDigi>& in,
                std::vector<SiStripDigi>& selectedSignal,
                uint32_t detId,
                edm::ESHandle<SiStripNoises>&,
                edm::ESHandle<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 {
    if (adc > 253 && doTruncate && !doTruncate10bits)
      return ((adc == 1023) ? 255 : 254);
    return adc;
  };

private:
  edm::ESHandle<SiStripNoises> noiseHandle;
  edm::ESHandle<SiStripThreshold> thresholdHandle;
  uint32_t noise_cache_id, threshold_cache_id;

  uint16_t theFEDalgorithm;
  bool isAValidDigi();

  bool doTruncate;
  bool doTruncate10bits;
  int16_t theFEDlowThresh;
  int16_t theFEDhighThresh;

  int16_t adc;
  int16_t adcPrev;
  int16_t adcNext;
  int16_t adcMaxNeigh;
  int16_t adcPrev2;
  int16_t adcNext2;

  int16_t thePrevFEDlowThresh;
  int16_t thePrevFEDhighThresh;
  int16_t theNextFEDlowThresh;
  int16_t theNextFEDhighThresh;

  int16_t theNeighFEDlowThresh;
  int16_t theNeighFEDhighThresh;

  int16_t thePrev2FEDlowThresh;
  int16_t theNext2FEDlowThresh;

  // working caches
  std::vector<int16_t> highThr_, lowThr_;    // thresholds in adc counts
  std::vector<float> highThrSN_, lowThrSN_;  // thresholds as S/N
  std::vector<float> noises_;

  void fillThresholds_(const uint32_t detID, size_t size);
};
#endif