EgammaElectronTkIsolationProducer.cc

Go to the documentation of this file.

//*****************************************************************************
// File:      EgammaElectronTkIsolationProducer.cc
// ----------------------------------------------------------------------------
// OrigAuth:  Matthias Mozer
// Institute: IIHE-VUB
//=============================================================================
//*****************************************************************************
 
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/global/EDProducer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "RecoEgamma/EgammaIsolationAlgos/interface/ElectronTkIsolation.h"
 
class EgammaElectronTkIsolationProducer : public edm::global::EDProducer<> {
public:
  explicit EgammaElectronTkIsolationProducer(const edm::ParameterSet&);
 
  void produce(edm::StreamID, edm::Event&, const edm::EventSetup&) const override;
 
private:
  const edm::EDGetTokenT<reco::GsfElectronCollection> electronProducer_;
  const edm::EDGetTokenT<reco::TrackCollection> trackProducer_;
  const edm::EDGetTokenT<reco::BeamSpot> beamspotProducer_;
 
  const double ptMin_;
  const double intRadiusBarrel_;
  const double intRadiusEndcap_;
  const double stripBarrel_;
  const double stripEndcap_;
  const double extRadius_;
  const double maxVtxDist_;
  const double drb_;
};
 
#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(EgammaElectronTkIsolationProducer);
 
EgammaElectronTkIsolationProducer::EgammaElectronTkIsolationProducer(const edm::ParameterSet& config)
    : electronProducer_{consumes(config.getParameter<edm::InputTag>("electronProducer"))},
      trackProducer_{consumes(config.getParameter<edm::InputTag>("trackProducer"))},
      beamspotProducer_{consumes(config.getParameter<edm::InputTag>("BeamspotProducer"))},
 
      ptMin_{config.getParameter<double>("ptMin")},
      intRadiusBarrel_{config.getParameter<double>("intRadiusBarrel")},
      intRadiusEndcap_{config.getParameter<double>("intRadiusEndcap")},
      stripBarrel_{config.getParameter<double>("stripBarrel")},
      stripEndcap_{config.getParameter<double>("stripEndcap")},
      extRadius_{config.getParameter<double>("extRadius")},
      maxVtxDist_{config.getParameter<double>("maxVtxDist")},
      drb_{config.getParameter<double>("maxVtxDistXY")}
 
{
  produces<edm::ValueMap<double>>();
}
 
void EgammaElectronTkIsolationProducer::produce(edm::StreamID, edm::Event& iEvent, const edm::EventSetup&) const {
  // Get the  filtered objects
  auto electronHandle = iEvent.getHandle(electronProducer_);
 
  //prepare product
  auto isoMap = std::make_unique<edm::ValueMap<double>>();
  edm::ValueMap<double>::Filler filler(*isoMap);
  std::vector<double> retV(electronHandle->size(), 0);
 
  ElectronTkIsolation myTkIsolation(extRadius_,
                                    intRadiusBarrel_,
                                    intRadiusEndcap_,
                                    stripBarrel_,
                                    stripEndcap_,
                                    ptMin_,
                                    maxVtxDist_,
                                    drb_,
                                    &iEvent.get(trackProducer_),
                                    iEvent.get(beamspotProducer_).position());
 
  for (unsigned int i = 0; i < electronHandle->size(); ++i) {
    double isoValue = myTkIsolation.getPtTracks(&(electronHandle->at(i)));
    retV[i] = isoValue;
  }
 
  //fill and insert valuemap
  filler.insert(electronHandle, retV.begin(), retV.end());
  filler.fill();
  iEvent.put(std::move(isoMap));
}