FixedGridRhoProducer.cc

Go to the documentation of this file.

#include "RecoJets/JetProducers/plugins/FixedGridRhoProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

using namespace std;

FixedGridRhoProducer::FixedGridRhoProducer(const edm::ParameterSet& iConfig) {
  pfCandidatesTag_ = iConfig.getParameter<edm::InputTag>("pfCandidatesTag");
  string etaRegion = iConfig.getParameter<string>("EtaRegion");
  if (etaRegion == "Central")
    myEtaRegion = FixedGridEnergyDensity::Central;
  else if (etaRegion == "Forward")
    myEtaRegion = FixedGridEnergyDensity::Forward;
  else if (etaRegion == "All")
    myEtaRegion = FixedGridEnergyDensity::All;
  else {
    edm::LogWarning("FixedGridRhoProducer")
        << "Wrong EtaRegion parameter: " << etaRegion << ". Using EtaRegion = Central";
    myEtaRegion = FixedGridEnergyDensity::Central;
  }
  produces<double>();

  input_pfcoll_token_ = consumes<reco::PFCandidateCollection>(pfCandidatesTag_);
}

FixedGridRhoProducer::~FixedGridRhoProducer() {}

void FixedGridRhoProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {
  edm::Handle<reco::PFCandidateCollection> pfColl;
  iEvent.getByToken(input_pfcoll_token_, pfColl);

  algo = new FixedGridEnergyDensity(pfColl.product());

  double result = algo->fixedGridRho(myEtaRegion);
  iEvent.put(std::make_unique<double>(result));

  delete algo;
}

DEFINE_FWK_MODULE(FixedGridRhoProducer);