ESZeroSuppressionProducer.cc

Go to the documentation of this file.

#include "DataFormats/EcalDigi/interface/ESDataFrame.h"
#include "Geometry/Records/interface/IdealGeometryRecord.h"
#include "SimCalorimetry/EcalZeroSuppressionProducers/interface/ESZeroSuppressionProducer.h"
 
ESZeroSuppressionProducer::ESZeroSuppressionProducer(const edm::ParameterSet &ps) {
  digiProducer_ = ps.getParameter<std::string>("digiProducer");
  ESdigiCollection_ = ps.getParameter<std::string>("ESdigiCollection");
  ESZSdigiCollection_ = ps.getParameter<std::string>("ESZSdigiCollection");
 
  produces<ESDigiCollection>(ESZSdigiCollection_);
 
  ES_token = consumes<ESDigiCollection>(edm::InputTag(digiProducer_));
  ;
}
 
ESZeroSuppressionProducer::~ESZeroSuppressionProducer() {}
 
void ESZeroSuppressionProducer::produce(edm::Event &event, const edm::EventSetup &eventSetup) {
  eventSetup.get<ESThresholdsRcd>().get(esthresholds_);
  const ESThresholds *thresholds = esthresholds_.product();
 
  eventSetup.get<ESPedestalsRcd>().get(espeds_);
  const ESPedestals *pedestals = espeds_.product();
 
  float ts2Threshold = thresholds->getTS2Threshold();
 
  edm::Handle<ESDigiCollection> ESDigis;
 
  bool fullESDigis = true;
  event.getByToken(ES_token, ESDigis);
  if (!ESDigis.isValid()) {
    edm::LogError("ZeroSuppressionError") << "Error! can't get the product " << ESdigiCollection_.c_str();
    fullESDigis = false;
  }
 
  std::unique_ptr<ESDigiCollection> ESZSDigis(new ESDigiCollection());
 
  if (fullESDigis) {
    for (ESDigiCollection::const_iterator i(ESDigis->begin()); i != ESDigis->end(); ++i) {
      ESDataFrame dataframe = (*i);
 
      ESPedestals::const_iterator it_ped = pedestals->find(dataframe.id());
 
      if (dataframe.sample(1).adc() > (ts2Threshold + it_ped->getMean())) {
        // std::cout<<dataframe.sample(1).adc()<<"
        // "<<ts2Threshold+it_ped->getMean()<<std::endl;
        (*ESZSDigis).push_back(*i);
      }
    }
  }
 
  event.put(std::move(ESZSDigis), ESZSdigiCollection_);
}