CastorJetIDProducer.cc

Go to the documentation of this file.

#include "RecoJets/JetProducers/plugins/CastorJetIDProducer.h"
#include "DataFormats/JetReco/interface/CastorJetID.h"

#include <vector>

//
// constants, enums and typedefs
//


//
// static data member definitions
//

//
// constructors and destructor
//
CastorJetIDProducer::CastorJetIDProducer(const edm::ParameterSet& iConfig) :
  src_       ( iConfig.getParameter<edm::InputTag>("src") ),
  helper_    ( )
{
  produces< reco::CastorJetIDValueMap >();

  input_jet_token_ = consumes<edm::View<reco::BasicJet> >(src_);

}


CastorJetIDProducer::~CastorJetIDProducer()
{
}


//
// member functions
//

// ------------ method called to produce the data  ------------
void
CastorJetIDProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup)
{

  // get the input jets
  edm::Handle< edm::View<reco::BasicJet> > h_jets;
  iEvent.getByToken( input_jet_token_, h_jets );

  // allocate the jet--->jetid value map
  auto castorjetIdValueMap = std::make_unique<reco::CastorJetIDValueMap>();
  // instantiate the filler with the map
  reco::CastorJetIDValueMap::Filler filler(*castorjetIdValueMap);
  
  // allocate the vector of ids
  size_t njets = h_jets->size();
  std::vector<reco::CastorJetID>  ids (njets);
   
  // loop over the jets
  for ( edm::View<reco::BasicJet>::const_iterator jetsBegin = h_jets->begin(),
      jetsEnd = h_jets->end(),
      ijet = jetsBegin;
    ijet != jetsEnd; ++ijet ) {

    // get the id from each jet
    helper_.calculate( iEvent, *ijet );

    ids[ijet-jetsBegin].emEnergy               =  helper_.emEnergy();
    ids[ijet-jetsBegin].hadEnergy               =  helper_.hadEnergy();
    ids[ijet-jetsBegin].fem            =  helper_.fem();
    ids[ijet-jetsBegin].depth      =  helper_.depth();
    ids[ijet-jetsBegin].width      =  helper_.width();
    ids[ijet-jetsBegin].fhot      =  helper_.fhot();
    ids[ijet-jetsBegin].sigmaz      =  helper_.sigmaz();
    ids[ijet-jetsBegin].nTowers      =  helper_.nTowers(); 


  }
  
  // set up the map
  filler.insert( h_jets, ids.begin(), ids.end() );
 
  // fill the vals
  filler.fill();

  // write map to the event
  iEvent.put(std::move(castorjetIdValueMap));
}

//define this as a plug-in
DEFINE_FWK_MODULE(CastorJetIDProducer);