#include <IsolatedPFCandidateSelectorDefinition.h>

Inheritance diagram for pf2pat::IsolatedPFCandidateSelectorDefinition:

Public Types
typedef edm::ValueMap< double >	IsoMap

Public Types inherited from pf2pat::PFCandidateSelectorDefinition
typedef reco::PFCandidateCollection	collection

typedef boost::transform_iterator < Pointer, container::const_iterator >	const_iterator

typedef std::vector < reco::PFCandidate >	container

typedef edm::Handle< collection >	HandleToCollection

Public Member Functions
	deltaBetaFactor_ (cfg.getParameter< double >("deltaBetaFactor"))

	deltaBetaIsolationValueMapToken_ (iC.mayConsume< IsoMap >(cfg.getParameter< edm::InputTag >("deltaBetaIsolationValueMap")))

	doDeltaBetaCorrection_ (cfg.getParameter< bool >("doDeltaBetaCorrection"))

	IsolatedPFCandidateSelectorDefinition (const edm::ParameterSet &cfg, edm::ConsumesCollector &&iC)

	isolationCut_ (cfg.getParameter< double >("isolationCut"))

	isolationValueMapNeutralTokens_ (edm::vector_transform(cfg.getParameter< std::vector< edm::InputTag > >("isolationValueMapsNeutral"), [&](edm::InputTag const &tag){return iC.consumes< IsoMap >(tag);}))

	isRelative_ (cfg.getParameter< bool >("isRelative"))

void	select (const HandleToCollection &hc, const edm::Event &e, const edm::EventSetup &s)

Public Member Functions inherited from pf2pat::PFCandidateSelectorDefinition
const_iterator	begin () const

const_iterator	end () const

	PFCandidateSelectorDefinition ()

const container &	selected () const

size_t	size () const

Private Attributes
double	deltaBetaFactor_

edm::EDGetTokenT< IsoMap >	deltaBetaIsolationValueMapToken_

bool	doDeltaBetaCorrection_

double	isolationCut_

std::vector< edm::EDGetTokenT < IsoMap > >	isolationValueMapChargedTokens_

std::vector< edm::EDGetTokenT < IsoMap > >	isolationValueMapNeutralTokens_

bool	isRelative_

Additional Inherited Members
Protected Attributes inherited from pf2pat::PFCandidateSelectorDefinition
container	selected_

Detailed Description

Definition at line 16 of file IsolatedPFCandidateSelectorDefinition.h.

Member Typedef Documentation

typedef edm::ValueMap<double> pf2pat::IsolatedPFCandidateSelectorDefinition::IsoMap

Definition at line 18 of file IsolatedPFCandidateSelectorDefinition.h.

Constructor & Destructor Documentation

pf2pat::IsolatedPFCandidateSelectorDefinition::IsolatedPFCandidateSelectorDefinition	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &&	iC
	)

inline

Definition at line 20 of file IsolatedPFCandidateSelectorDefinition.h.

References GlobalPosition_Frontier_DevDB_cff::tag.

         : isolationValueMapChargedTokens_(
               edm::vector_transform(cfg.getParameter<std::vector<edm::InputTag> >("isolationValueMapsCharged"),
                                     [&](edm::InputTag const& tag) { return iC.consumes<IsoMap>(tag); })),

Member Function Documentation

pf2pat::IsolatedPFCandidateSelectorDefinition::deltaBetaFactor_ ( cfg.getParameter< double > "deltaBetaFactor" )

pf2pat::IsolatedPFCandidateSelectorDefinition::deltaBetaIsolationValueMapToken_ ( iC.mayConsume< IsoMap > cfg.getParameter< edm::InputTag >("deltaBetaIsolationValueMap") )

pf2pat::IsolatedPFCandidateSelectorDefinition::doDeltaBetaCorrection_ ( cfg.getParameter< bool > "doDeltaBetaCorrection" )

pf2pat::IsolatedPFCandidateSelectorDefinition::isolationCut_ ( cfg.getParameter< double > "isolationCut" )

inline

Definition at line 32 of file IsolatedPFCandidateSelectorDefinition.h.

32 {}

pf2pat::IsolatedPFCandidateSelectorDefinition::isolationValueMapNeutralTokens_ ( edm:: vector_transformcfg.getParameter< std::vector< edm::InputTag > >("isolationValueMapsNeutral"),[&](edm::InputTag const &tag){return iC.consumes< IsoMap >(tag);} )

pf2pat::IsolatedPFCandidateSelectorDefinition::isRelative_ ( cfg.getParameter< bool > "isRelative" )

void pf2pat::IsolatedPFCandidateSelectorDefinition::select	(	const HandleToCollection &	hc,
		const edm::Event &	e,
		const edm::EventSetup &	s
	)

inline

Definition at line 34 of file IsolatedPFCandidateSelectorDefinition.h.

References edm::Event::getByToken(), objects.IsoTrackAnalyzer::isoSum, submitPVResolutionJobs::key, and hgcalPerformanceValidation::val.

                                                                                          {
       selected_.clear();
 
       // read all charged isolation value maps
       std::vector<edm::Handle<IsoMap> > isoMapsCharged(isolationValueMapChargedTokens_.size());
       for (unsigned iMap = 0; iMap < isolationValueMapChargedTokens_.size(); ++iMap) {
         e.getByToken(isolationValueMapChargedTokens_[iMap], isoMapsCharged[iMap]);
       }
 
       // read all neutral isolation value maps
       std::vector<edm::Handle<IsoMap> > isoMapsNeutral(isolationValueMapNeutralTokens_.size());
       for (unsigned iMap = 0; iMap < isolationValueMapNeutralTokens_.size(); ++iMap) {
         e.getByToken(isolationValueMapNeutralTokens_[iMap], isoMapsNeutral[iMap]);
       }
 
       edm::Handle<IsoMap> dBetaH;
       if (doDeltaBetaCorrection_) {
         e.getByToken(deltaBetaIsolationValueMapToken_, dBetaH);
       }
 
       unsigned key = 0;
       for (collection::const_iterator pfc = hc->begin(); pfc != hc->end(); ++pfc, ++key) {
         reco::PFCandidateRef candidate(hc, key);
 
         bool passed = true;
         double isoSumCharged = 0.0;
         double isoSumNeutral = 0.0;
 
         for (unsigned iMap = 0; iMap < isoMapsCharged.size(); ++iMap) {
           const IsoMap& isoMap = *(isoMapsCharged[iMap]);
           double val = isoMap[candidate];
           isoSumCharged += val;
         }
 
         for (unsigned iMap = 0; iMap < isoMapsNeutral.size(); ++iMap) {
           const IsoMap& isoMap = *(isoMapsNeutral[iMap]);
           double val = isoMap[candidate];
           isoSumNeutral += val;
         }
 
         if (doDeltaBetaCorrection_) {
           const IsoMap& isoMap = *dBetaH;
           double dBetaVal = isoMap[candidate];
           double dBetaCorIsoSumNeutral = isoSumNeutral + deltaBetaFactor_ * dBetaVal;
           isoSumNeutral = dBetaCorIsoSumNeutral > 0 ? dBetaCorIsoSumNeutral : 0;  //follow muon POG definition in 2012
         }
 
         double isoSum = isoSumCharged + isoSumNeutral;
 
         if (isRelative_) {
           isoSum /= candidate->pt();
         }
 
         if (isoSum > isolationCut_) {
           passed = false;
         }
 
         if (passed) {
           // passed all cuts, selected
           selected_.push_back(reco::PFCandidate(*pfc));
           reco::PFCandidatePtr ptrToMother(hc, key);
           selected_.back().setSourceCandidatePtr(ptrToMother);
         }
       }
     }