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);
   std::auto_ptr<double> outputRho(new double(bge_.rho()));
   iEvent.put(outputRho);
}

DEFINE_FWK_MODULE(FixedGridRhoProducerFastjet);