#include <EgammaElectronTkIsolationProducer.h>

Inheritance diagram for EgammaElectronTkIsolationProducer:

Public Member Functions
	EgammaElectronTkIsolationProducer (const edm::ParameterSet &)
virtual void	produce (edm::Event &, const edm::EventSetup &)
	~EgammaElectronTkIsolationProducer ()
Private Attributes
edm::InputTag	beamspotProducer_
edm::ParameterSet	conf_
double	drb_
edm::InputTag	electronProducer_
double	extRadius_
double	intRadiusBarrel_
double	intRadiusEndcap_
double	maxVtxDist_
double	ptMin_
double	stripBarrel_
double	stripEndcap_
edm::InputTag	trackProducer_

Detailed Description

Definition at line 20 of file EgammaElectronTkIsolationProducer.h.

Constructor & Destructor Documentation

EgammaElectronTkIsolationProducer::EgammaElectronTkIsolationProducer ( const edm::ParameterSet & config ) [explicit]

Definition at line 27 of file EgammaElectronTkIsolationProducer.cc.

References beamspotProducer_, conf_, drb_, electronProducer_, extRadius_, edm::ParameterSet::getParameter(), intRadiusBarrel_, intRadiusEndcap_, maxVtxDist_, ptMin_, stripBarrel_, stripEndcap_, and trackProducer_.

                                                                                                  : conf_(config)
{
  // use configuration file to setup input/output collection names
  electronProducer_               = conf_.getParameter<edm::InputTag>("electronProducer");
  
  trackProducer_           = conf_.getParameter<edm::InputTag>("trackProducer");
  beamspotProducer_        = conf_.getParameter<edm::InputTag>("BeamspotProducer");

  ptMin_                = conf_.getParameter<double>("ptMin");
  intRadiusBarrel_      = conf_.getParameter<double>("intRadiusBarrel");
  intRadiusEndcap_      = conf_.getParameter<double>("intRadiusEndcap");
  stripBarrel_          = conf_.getParameter<double>("stripBarrel");
  stripEndcap_          = conf_.getParameter<double>("stripEndcap");
  extRadius_            = conf_.getParameter<double>("extRadius");
  maxVtxDist_           = conf_.getParameter<double>("maxVtxDist");
  drb_                  = conf_.getParameter<double>("maxVtxDistXY");

  //register your products
  produces < edm::ValueMap<double> >();

}

EgammaElectronTkIsolationProducer::~EgammaElectronTkIsolationProducer ( )

Definition at line 49 of file EgammaElectronTkIsolationProducer.cc.

{}

Member Function Documentation

void EgammaElectronTkIsolationProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)		`[virtual]`

Implements edm::EDProducer.

Definition at line 53 of file EgammaElectronTkIsolationProducer.cc.

References beamspotProducer_, drb_, electronProducer_, extRadius_, edm::helper::Filler< Map >::fill(), edm::Event::getByLabel(), ElectronTkIsolation::getPtTracks(), i, edm::helper::Filler< Map >::insert(), intRadiusBarrel_, intRadiusEndcap_, maxVtxDist_, edm::Handle< T >::product(), ptMin_, edm::Event::put(), stripBarrel_, stripEndcap_, trackProducer_, and testEve_cfg::tracks.

{
  // Get the  filtered objects
  edm::Handle< reco::GsfElectronCollection> electronHandle;
  iEvent.getByLabel(electronProducer_,electronHandle);
  
  //get the tracks
  edm::Handle<reco::TrackCollection> tracks;
  iEvent.getByLabel(trackProducer_,tracks);
  const reco::TrackCollection* trackCollection = tracks.product();
  
  //prepare product
  std::auto_ptr<edm::ValueMap<double> > isoMap(new edm::ValueMap<double>());
  edm::ValueMap<double>::Filler filler(*isoMap);
  std::vector<double> retV(electronHandle->size(),0);
 
  edm::Handle<reco::BeamSpot> beamSpotH;
  iEvent.getByLabel(beamspotProducer_,beamSpotH);
  reco::TrackBase::Point beamspot = beamSpotH->position();
 
  ElectronTkIsolation myTkIsolation (extRadius_,intRadiusBarrel_,intRadiusEndcap_,stripBarrel_,stripEndcap_,ptMin_,maxVtxDist_,drb_,trackCollection,beamspot) ;
  
  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(isoMap);
}