ESRecHitProducer.cc

Go to the documentation of this file.

// ESRecHitProducer author : Chia-Ming, Kuo

#include "DataFormats/Common/interface/EDCollection.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHit.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/stream/EDProducer.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "RecoLocalCalo/EcalRecProducers/interface/ESRecHitWorkerFactory.h"

#include "ESRecHitWorker.h"

class ESRecHitProducer : public edm::stream::EDProducer<> {
public:
  explicit ESRecHitProducer(const edm::ParameterSet& ps);
  void produce(edm::Event& e, const edm::EventSetup& es) override;

private:
  const edm::EDGetTokenT<ESDigiCollection> digiToken_;
  const std::string rechitCollection_;  // secondary name to be given to collection of hits

  std::unique_ptr<ESRecHitWorkerBaseClass> worker_;
};

ESRecHitProducer::ESRecHitProducer(edm::ParameterSet const& ps)
    : digiToken_(consumes<ESDigiCollection>(ps.getParameter<edm::InputTag>("ESdigiCollection"))),
      rechitCollection_(ps.getParameter<std::string>("ESrechitCollection")),
      worker_{ESRecHitWorkerFactory::get()->create(ps.getParameter<std::string>("algo"), ps, consumesCollector())} {
  produces<ESRecHitCollection>(rechitCollection_);
}

void ESRecHitProducer::produce(edm::Event& e, const edm::EventSetup& es) {
  edm::Handle<ESDigiCollection> digiHandle;
  const ESDigiCollection* digi = nullptr;
  e.getByToken(digiToken_, digiHandle);

  digi = digiHandle.product();
  LogDebug("ESRecHitInfo") << "total # ESdigis: " << digi->size();

  // Create empty output
  auto rec = std::make_unique<ESRecHitCollection>();

  if (digi) {
    rec->reserve(digi->size());

    worker_->set(es);

    // run the algorithm
    for (ESDigiCollection::const_iterator i(digi->begin()); i != digi->end(); i++) {
      worker_->run(i, *rec);
    }
  }

  e.put(std::move(rec), rechitCollection_);
}

#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(ESRecHitProducer);