FixedGridRhoProducerFastjet.cc

Go to the documentation of this file.

#include "RecoJets/JetProducers/plugins/FixedGridRhoProducerFastjet.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Common/interface/View.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"

using namespace std;

FixedGridRhoProducerFastjet::FixedGridRhoProducerFastjet(const edm::ParameterSet& iConfig)
    : bge_(iConfig.getParameter<double>("maxRapidity"), iConfig.getParameter<double>("gridSpacing")) {
  pfCandidatesTag_ = iConfig.getParameter<edm::InputTag>("pfCandidatesTag");
  produces<double>();

  input_pfcoll_token_ = consumes<edm::View<reco::Candidate> >(pfCandidatesTag_);
}

FixedGridRhoProducerFastjet::~FixedGridRhoProducerFastjet() {}

void FixedGridRhoProducerFastjet::produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {
  edm::Handle<edm::View<reco::Candidate> > pfColl;
  iEvent.getByToken(input_pfcoll_token_, pfColl);
  std::vector<fastjet::PseudoJet> inputs;
  for (edm::View<reco::Candidate>::const_iterator ibegin = pfColl->begin(), iend = pfColl->end(), i = ibegin; i != iend;
       ++i) {
    inputs.push_back(fastjet::PseudoJet(i->px(), i->py(), i->pz(), i->energy()));
  }
  bge_.set_particles(inputs);
  iEvent.put(std::make_unique<double>(bge_.rho()));
}

DEFINE_FWK_MODULE(FixedGridRhoProducerFastjet);