Inheritance diagram for PhoAnyPFIsoWithEACut:

Public Member Functions
CandidateType	candidateType () const final

void	getEventContent (const edm::EventBase &) final

result_type	operator() (const reco::PhotonPtr &) const final

	PhoAnyPFIsoWithEACut (const edm::ParameterSet &c)

void	setConsumes (edm::ConsumesCollector &) final

double	value (const reco::CandidatePtr &cand) const final

Public Member Functions inherited from CutApplicatorWithEventContentBase
	CutApplicatorWithEventContentBase ()

	CutApplicatorWithEventContentBase (const CutApplicatorWithEventContentBase &)=delete

	CutApplicatorWithEventContentBase (const edm::ParameterSet &c)

CutApplicatorWithEventContentBase &	operator= (const CutApplicatorWithEventContentBase &)=delete

	~CutApplicatorWithEventContentBase () override
	Destructor. More...

Public Member Functions inherited from CutApplicatorBase
virtual result_type	asCandidate (const argument_type &) const

	CutApplicatorBase ()

	CutApplicatorBase (const CutApplicatorBase &)=delete

	CutApplicatorBase (const edm::ParameterSet &c)

const std::string &	name () const override

result_type	operator() (const argument_type &) const final

virtual result_type	operator() (const pat::ElectronPtr &) const

virtual result_type	operator() (const pat::MuonPtr &) const

virtual result_type	operator() (const reco::GsfElectronPtr &) const

virtual result_type	operator() (const reco::MuonPtr &) const

CutApplicatorBase &	operator= (const CutApplicatorBase &)=delete

	~CutApplicatorBase () override
	Destructor. More...

Public Member Functions inherited from candidate_functions::CandidateCut
	CandidateCut ()

virtual	~CandidateCut ()

Private Attributes
edm::Handle< edm::ValueMap< float > >	_anyPFIsoMap

float	_barrelCutOff

float	_C1_EB

float	_C1_EE

float	_C2_EB

float	_C2_EE

EffectiveAreas	_effectiveAreas

edm::Handle< double >	_rhoHandle

bool	_useRelativeIso

Static Private Attributes
constexpr static char	anyPFIsoWithEA_ [] = "anyPFIsoWithEA"

constexpr static char	rhoString_ [] = "rho"

Additional Inherited Members
Public Types inherited from CutApplicatorBase
enum	CandidateType { NONE, ELECTRON, MUON, PHOTON, TAU, PATELECTRON, PATMUON, PATPHOTON, PATTAU }

Public Types inherited from candidate_functions::CandidateCut
using	argument_type = reco::CandidatePtr

using	result_type = bool

Protected Attributes inherited from CutApplicatorWithEventContentBase
std::unordered_map< std::string, edm::InputTag >	contentTags_

std::unordered_map< std::string, edm::EDGetToken >	contentTokens_

Detailed Description

Definition at line 5 of file PhoAnyPFIsoWithEACut.cc.

Constructor & Destructor Documentation

◆ PhoAnyPFIsoWithEACut()

PhoAnyPFIsoWithEACut::PhoAnyPFIsoWithEACut ( const edm::ParameterSet & c )

Definition at line 43 of file PhoAnyPFIsoWithEACut.cc.

     : CutApplicatorWithEventContentBase(c),
       _C1_EB(c.getParameter<double>("C1_EB")),
       _C2_EB(c.getParameter<double>("C2_EB")),
       _C1_EE(c.getParameter<double>("C1_EE")),
       _C2_EE(c.getParameter<double>("C2_EE")),
       _barrelCutOff(c.getParameter<double>("barrelCutOff")),
       _useRelativeIso(c.getParameter<bool>("useRelativeIso")),
       _effectiveAreas((c.getParameter<edm::FileInPath>("effAreasConfigFile")).fullPath()) {
   edm::InputTag maptag = c.getParameter<edm::InputTag>("anyPFIsoMap");
   contentTags_.emplace(anyPFIsoWithEA_, maptag);
  
   edm::InputTag rhoTag = c.getParameter<edm::InputTag>("rho");
   contentTags_.emplace(rhoString_, rhoTag);
 }

References anyPFIsoWithEA_, HltBtagPostValidation_cff::c, CutApplicatorWithEventContentBase::contentTags_, rhoString_, and HLT_2018_cff::rhoTag.

Member Function Documentation

◆ candidateType()

CandidateType PhoAnyPFIsoWithEACut::candidateType ( ) const

inlinefinalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 16 of file PhoAnyPFIsoWithEACut.cc.

16 { return PHOTON; }

References CutApplicatorBase::PHOTON.

◆ getEventContent()

void PhoAnyPFIsoWithEACut::getEventContent ( const edm::EventBase & ev )

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 67 of file PhoAnyPFIsoWithEACut.cc.

                                                                  {
   ev.getByLabel(contentTags_[anyPFIsoWithEA_], _anyPFIsoMap);
   ev.getByLabel(contentTags_[rhoString_], _rhoHandle);
 }

References _anyPFIsoMap, _rhoHandle, anyPFIsoWithEA_, CutApplicatorWithEventContentBase::contentTags_, ev, and rhoString_.

◆ operator()()

CutApplicatorBase::result_type PhoAnyPFIsoWithEACut::operator() ( const reco::PhotonPtr & cand ) const

finalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 72 of file PhoAnyPFIsoWithEACut.cc.

                                                                                              {
   // in case we are by-value
   const std::string& inst_name = contentTags_.find(anyPFIsoWithEA_)->second.instance();
   edm::Ptr<pat::Photon> pat(cand);
   float anyisoval = -1.0;
   if (_anyPFIsoMap.isValid() && _anyPFIsoMap->contains(cand.id())) {
     anyisoval = (*_anyPFIsoMap)[cand];
   } else if (_anyPFIsoMap.isValid() && _anyPFIsoMap->idSize() == 1 && cand.id() == edm::ProductID()) {
     // in case we have spoofed a ptr
     //note this must be a 1:1 valuemap (only one product input)
     anyisoval = _anyPFIsoMap->begin()[cand.key()];
   } else if (_anyPFIsoMap.isValid()) {  // throw an exception
     anyisoval = (*_anyPFIsoMap)[cand];
   }
  
   // Figure out the cut value
   // The value is generally pt-dependent: C1 + pt * C2
   const double absEta = std::abs(cand->superCluster()->eta());
   const float anyPFIsoWithEACutValue =
       (absEta < _barrelCutOff ? _C1_EB + cand->pt() * _C2_EB : _C1_EE + cand->pt() * _C2_EE);
  
   // Retrieve the variable value for this particle
   float anyPFIso = _anyPFIsoMap.isValid() ? anyisoval : pat->userFloat(inst_name);
  
   // Apply pile-up correction
   const double eA = _effectiveAreas.getEffectiveArea(absEta);
   const double rho = _rhoHandle.isValid() ? *_rhoHandle : 0;
   const float anyPFIsoWithEA = std::max(0.0, anyPFIso - rho * eA);
  
   // Apply the cut and return the result
   // Scale by pT if the relative isolation is requested but avoid division by 0
   return anyPFIsoWithEA < anyPFIsoWithEACutValue * (_useRelativeIso ? cand->pt() : 1.);
 }

References _anyPFIsoMap, _barrelCutOff, _C1_EB, _C1_EE, _C2_EB, _C2_EE, _effectiveAreas, _rhoHandle, _useRelativeIso, funct::abs(), anyPFIsoWithEA_, edm::ValueMap< T >::begin(), edm::ValueMap< T >::contains(), CutApplicatorWithEventContentBase::contentTags_, EffectiveAreas::getEffectiveArea(), edm::ValueMap< T >::idSize(), edm::HandleBase::isValid(), SiStripPI::max, rho, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ setConsumes()

void PhoAnyPFIsoWithEACut::setConsumes ( edm::ConsumesCollector & cc )

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 59 of file PhoAnyPFIsoWithEACut.cc.

                                                                {
   auto anyPFIsoWithEA = cc.consumes<edm::ValueMap<float> >(contentTags_[anyPFIsoWithEA_]);
   contentTokens_.emplace(anyPFIsoWithEA_, anyPFIsoWithEA);
  
   auto rho = cc.consumes<double>(contentTags_[rhoString_]);
   contentTokens_.emplace(rhoString_, rho);
 }

References anyPFIsoWithEA_, CutApplicatorWithEventContentBase::contentTags_, CutApplicatorWithEventContentBase::contentTokens_, rho, and rhoString_.

◆ value()

double PhoAnyPFIsoWithEACut::value ( const reco::CandidatePtr & cand ) const

finalvirtual

Implements candidate_functions::CandidateCut.

Definition at line 106 of file PhoAnyPFIsoWithEACut.cc.

                                                                      {
   reco::PhotonPtr pho(cand);
  
   // in case we are by-value
   const std::string& inst_name = contentTags_.find(anyPFIsoWithEA_)->second.instance();
   edm::Ptr<pat::Photon> pat(cand);
   float anyisoval = -1.0;
   if (_anyPFIsoMap.isValid() && _anyPFIsoMap->contains(cand.id())) {
     anyisoval = (*_anyPFIsoMap)[cand];
   } else if (_anyPFIsoMap.isValid() && _anyPFIsoMap->idSize() == 1 && cand.id() == edm::ProductID()) {
     // in case we have spoofed a ptr
     //note this must be a 1:1 valuemap (only one product input)
     anyisoval = _anyPFIsoMap->begin()[cand.key()];
   } else if (_anyPFIsoMap.isValid()) {  // throw an exception
     anyisoval = (*_anyPFIsoMap)[cand];
   }
  
   // Figure out the cut value
   // The value is generally pt-dependent: C1 + pt * C2
   double absEta = std::abs(pho->superCluster()->eta());
  
   // Retrieve the variable value for this particle
   float anyPFIso = _anyPFIsoMap.isValid() ? anyisoval : pat->userFloat(inst_name);
  
   // Apply pile-up correction
   double eA = _effectiveAreas.getEffectiveArea(absEta);
   double rho = *_rhoHandle;
   float anyPFIsoWithEA = std::max(0.0, anyPFIso - rho * eA);
  
   // Divide by pT if the relative isolation is requested
   if (_useRelativeIso)
     anyPFIsoWithEA /= pho->pt();
  
   // Apply the cut and return the result
   return anyPFIsoWithEA;
 }

References _anyPFIsoMap, _effectiveAreas, _rhoHandle, _useRelativeIso, funct::abs(), anyPFIsoWithEA_, edm::ValueMap< T >::begin(), edm::ValueMap< T >::contains(), CutApplicatorWithEventContentBase::contentTags_, EffectiveAreas::getEffectiveArea(), edm::ValueMap< T >::idSize(), edm::HandleBase::isValid(), SiStripPI::max, rho, and AlCaHLTBitMon_QueryRunRegistry::string.

Member Data Documentation

◆ _anyPFIsoMap

edm::Handle<edm::ValueMap<float> > PhoAnyPFIsoWithEACut::_anyPFIsoMap

private

Definition at line 30 of file PhoAnyPFIsoWithEACut.cc.

Referenced by getEventContent(), operator()(), and value().

◆ _barrelCutOff

float PhoAnyPFIsoWithEACut::_barrelCutOff

private

Definition at line 25 of file PhoAnyPFIsoWithEACut.cc.

Referenced by operator()().

◆ _C1_EB

float PhoAnyPFIsoWithEACut::_C1_EB

private

Definition at line 20 of file PhoAnyPFIsoWithEACut.cc.

Referenced by operator()().

◆ _C1_EE

float PhoAnyPFIsoWithEACut::_C1_EE

private

Definition at line 22 of file PhoAnyPFIsoWithEACut.cc.

Referenced by operator()().

◆ _C2_EB

float PhoAnyPFIsoWithEACut::_C2_EB

private

Definition at line 21 of file PhoAnyPFIsoWithEACut.cc.

Referenced by operator()().

◆ _C2_EE

float PhoAnyPFIsoWithEACut::_C2_EE

private

Definition at line 23 of file PhoAnyPFIsoWithEACut.cc.

Referenced by operator()().

◆ _effectiveAreas

EffectiveAreas PhoAnyPFIsoWithEACut::_effectiveAreas

private

Definition at line 28 of file PhoAnyPFIsoWithEACut.cc.

Referenced by operator()(), and value().

◆ _rhoHandle

edm::Handle<double> PhoAnyPFIsoWithEACut::_rhoHandle

private

Definition at line 32 of file PhoAnyPFIsoWithEACut.cc.

Referenced by getEventContent(), operator()(), and value().

◆ _useRelativeIso

bool PhoAnyPFIsoWithEACut::_useRelativeIso

private

Definition at line 26 of file PhoAnyPFIsoWithEACut.cc.

Referenced by operator()(), and value().

◆ anyPFIsoWithEA_

constexpr char PhoAnyPFIsoWithEACut::anyPFIsoWithEA_ = "anyPFIsoWithEA"

staticconstexprprivate

Definition at line 34 of file PhoAnyPFIsoWithEACut.cc.

Referenced by getEventContent(), operator()(), PhoAnyPFIsoWithEACut(), setConsumes(), and value().

◆ rhoString_

constexpr char PhoAnyPFIsoWithEACut::rhoString_ = "rho"

staticconstexprprivate

Definition at line 35 of file PhoAnyPFIsoWithEACut.cc.

Referenced by getEventContent(), PhoAnyPFIsoWithEACut(), and setConsumes().

Public Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ PhoAnyPFIsoWithEACut()

Member Function Documentation

◆ candidateType()

◆ getEventContent()

◆ operator()()

◆ setConsumes()

◆ value()

Member Data Documentation

◆ _anyPFIsoMap

◆ _barrelCutOff

◆ _C1_EB

◆ _C1_EE

◆ _C2_EB

◆ _C2_EE

◆ _effectiveAreas

◆ _rhoHandle

◆ _useRelativeIso

◆ anyPFIsoWithEA_

◆ rhoString_