#include <RecoEgamma/EgammaHLTProducers/interface/EgammaHLTElectronTrackIsolationProducers.h>

Inheritance diagram for EgammaHLTElectronTrackIsolationProducers:

Public Member Functions
	EgammaHLTElectronTrackIsolationProducers (const edm::ParameterSet &)

void	produce (edm::StreamID sid, edm::Event &, const edm::EventSetup &) const override

	~EgammaHLTElectronTrackIsolationProducers () override

Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Attributes
const edm::EDGetTokenT< reco::BeamSpot >	beamSpotProducer_

const double	egTrkIsoConeSize_

const double	egTrkIsoPtMin_

const double	egTrkIsoRSpan_

const double	egTrkIsoStripBarrel_

const double	egTrkIsoStripEndcap_

const double	egTrkIsoVetoConeSizeBarrel_

const double	egTrkIsoVetoConeSizeEndcap_

const double	egTrkIsoZSpan_

const edm::EDGetTokenT< reco::ElectronCollection >	electronProducer_

const edm::EDGetTokenT< reco::RecoEcalCandidateCollection >	recoEcalCandidateProducer_

const edm::EDGetTokenT< reco::TrackCollection >	trackProducer_

const bool	useGsfTrack_

const bool	useSCRefs_

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Author: Monica Vazquez Acosta (CERN)

Id: EgammaHLTElectronTrackIsolationProducers.cc,v 1.9 2011/12/19 11:17:28 sani Exp

Definition at line 37 of file EgammaHLTElectronTrackIsolationProducers.h.

Constructor & Destructor Documentation

EgammaHLTElectronTrackIsolationProducers::EgammaHLTElectronTrackIsolationProducers ( const edm::ParameterSet & config )

explicit

Definition at line 19 of file EgammaHLTElectronTrackIsolationProducers.cc.

References useSCRefs_.

                                                                                                                :
   electronProducer_           (consumes<reco::ElectronCollection>(config.getParameter<edm::InputTag>("electronProducer"))),
   trackProducer_              (consumes<reco::TrackCollection>(config.getParameter<edm::InputTag>("trackProducer"))),
   recoEcalCandidateProducer_  (consumes<reco::RecoEcalCandidateCollection>(config.getParameter<edm::InputTag>("recoEcalCandidateProducer"))),
   beamSpotProducer_           (consumes<reco::BeamSpot>(config.getParameter<edm::InputTag>("beamSpotProducer"))),
   useGsfTrack_                (config.getParameter<bool>("useGsfTrack")),
   useSCRefs_                  (config.getParameter<bool>("useSCRefs")),
   egTrkIsoPtMin_              (config.getParameter<double>("egTrkIsoPtMin")),
   egTrkIsoConeSize_           (config.getParameter<double>("egTrkIsoConeSize")),
   egTrkIsoZSpan_              (config.getParameter<double>("egTrkIsoZSpan")),
   egTrkIsoRSpan_              (config.getParameter<double>("egTrkIsoRSpan")),
   egTrkIsoVetoConeSizeBarrel_ (config.getParameter<double>("egTrkIsoVetoConeSizeBarrel")),
   egTrkIsoVetoConeSizeEndcap_ (config.getParameter<double>("egTrkIsoVetoConeSizeEndcap")),
   egTrkIsoStripBarrel_        (config.getParameter<double>("egTrkIsoStripBarrel")),
   egTrkIsoStripEndcap_        (config.getParameter<double>("egTrkIsoStripEndcap")) {
 
   //register your products
   if(useSCRefs_) 
     produces < reco::RecoEcalCandidateIsolationMap >();
   else 
     produces < reco::ElectronIsolationMap >();
 }

EgammaHLTElectronTrackIsolationProducers::~EgammaHLTElectronTrackIsolationProducers ( )

override

Definition at line 42 of file EgammaHLTElectronTrackIsolationProducers.cc.

43 {}

Member Function Documentation

void EgammaHLTElectronTrackIsolationProducers::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 45 of file EgammaHLTElectronTrackIsolationProducers.cc.

References edm::ConfigurationDescriptions::add(), and edm::ParameterSetDescription::add().

                                                                                                           {
   edm::ParameterSetDescription desc;  
   desc.add<edm::InputTag>("electronProducer", edm::InputTag("hltEleAnyWP80PixelMatchElectronsL1Seeded"));
   desc.add<edm::InputTag>("trackProducer", edm::InputTag("hltL1SeededEgammaRegionalCTFFinalFitWithMaterial"));
   desc.add<edm::InputTag>("recoEcalCandidateProducer", edm::InputTag()); 
   desc.add<edm::InputTag>("beamSpotProducer", edm::InputTag("hltOnlineBeamSpot"));
   desc.add<double>("egTrkIsoPtMin", 1.0);
   desc.add<double>("egTrkIsoConeSize", 0.3);
   desc.add<double>("egTrkIsoZSpan", 0.15);
   desc.add<double>("egTrkIsoRSpan", 999999.0);
   desc.add<double>("egTrkIsoVetoConeSizeBarrel", 0.03);
   desc.add<double>("egTrkIsoVetoConeSizeEndcap", 0.03);
   desc.add<double>("egTrkIsoStripBarrel", 0.03);
   desc.add<double>("egTrkIsoStripEndcap", 0.03);
   desc.add<bool>("useGsfTrack", false);
   desc.add<bool>("useSCRefs", false);
 
   descriptions.add("hltEgammaHLTElectronTrackIsolationProducers", desc);  
 }

void EgammaHLTElectronTrackIsolationProducers::produce	(	edm::StreamID	sid,
		edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)		const

override

Definition at line 64 of file EgammaHLTElectronTrackIsolationProducers.cc.

References beamSpotProducer_, egTrkIsoConeSize_, egTrkIsoPtMin_, egTrkIsoRSpan_, egTrkIsoStripBarrel_, egTrkIsoStripEndcap_, egTrkIsoVetoConeSizeBarrel_, egTrkIsoVetoConeSizeEndcap_, egTrkIsoZSpan_, electronProducer_, edm::Event::getByToken(), ElectronTkIsolation::getIso(), edm::AssociationMap< Tag >::insert(), reco::BeamSpot::position(), edm::Handle< T >::product(), edm::Event::put(), recoEcalCandidateProducer_, findElectronsInSiStrips_cfi::trackCollection, trackProducer_, useGsfTrack_, and useSCRefs_.

                                                                                                                            {
   
   edm::Handle<reco::ElectronCollection> electronHandle;
   iEvent.getByToken(electronProducer_,electronHandle);
 
  // Get the general tracks
   edm::Handle<reco::TrackCollection> trackHandle;
   iEvent.getByToken(trackProducer_, trackHandle);
   const reco::TrackCollection* trackCollection = trackHandle.product();
 
   reco::ElectronIsolationMap eleMap(electronHandle);
 
   edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
   iEvent.getByToken(beamSpotProducer_,recoBeamSpotHandle);
   
   const reco::BeamSpot::Point& beamSpotPosition = recoBeamSpotHandle->position(); 
 
   ElectronTkIsolation isoAlgo(egTrkIsoConeSize_,egTrkIsoVetoConeSizeBarrel_,egTrkIsoVetoConeSizeEndcap_,egTrkIsoStripBarrel_,egTrkIsoStripEndcap_,egTrkIsoPtMin_,egTrkIsoZSpan_,egTrkIsoRSpan_,trackCollection,beamSpotPosition);
   
   if(useSCRefs_){
 
     edm::Handle<reco::RecoEcalCandidateCollection> recoEcalCandHandle;
     iEvent.getByToken(recoEcalCandidateProducer_,recoEcalCandHandle);
     reco::RecoEcalCandidateIsolationMap recoEcalCandMap(recoEcalCandHandle);
 
     for(reco::RecoEcalCandidateCollection::const_iterator iRecoEcalCand = recoEcalCandHandle->begin(); iRecoEcalCand != recoEcalCandHandle->end(); iRecoEcalCand++){
       
       reco::RecoEcalCandidateRef recoEcalCandRef(recoEcalCandHandle,iRecoEcalCand-recoEcalCandHandle->begin());
       
       reco::ElectronRef eleRef;
       for(reco::ElectronCollection::const_iterator eleIt = electronHandle->begin(); eleIt != electronHandle->end(); eleIt++){
     if(eleIt->superCluster()==recoEcalCandRef->superCluster()){
       eleRef = reco::ElectronRef(electronHandle,eleIt - electronHandle->begin());
       break;
       }
       }
       float isol=999999;
       if(eleRef.isNonnull()){
     const reco::Track* eleTrk = useGsfTrack_ ? &*eleRef->gsfTrack() : &*eleRef->track();
     isol = isoAlgo.getIso(eleTrk).second;
       }
       recoEcalCandMap.insert(recoEcalCandRef,isol);
     }//end reco ecal candidate ref
 
     iEvent.put(std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandMap));
 
   }else{ //we are going to loop over electron instead
     for(reco::ElectronCollection::const_iterator iElectron = electronHandle->begin(); iElectron != electronHandle->end(); iElectron++){
       reco::ElectronRef eleRef(reco::ElectronRef(electronHandle,iElectron - electronHandle->begin()));
       const reco::Track* eleTrk = useGsfTrack_ ? &*eleRef->gsfTrack() : &*eleRef->track();
       float isol = isoAlgo.getIso(eleTrk).second;
       eleMap.insert(eleRef, isol);
     }
 
     iEvent.put(std::make_unique<reco::ElectronIsolationMap>(eleMap));
   }
 }