GsfEleEffAreaPFIsoCut Class Reference

Inheritance diagram for GsfEleEffAreaPFIsoCut:

Public Member Functions
CandidateType	candidateType () const final
void	getEventContent (const edm::EventBase &) final
	GsfEleEffAreaPFIsoCut (const edm::ParameterSet &c)
result_type	operator() (const reco::GsfElectronPtr &) const final
void	setConsumes (edm::ConsumesCollector &) final
double	value (const reco::CandidatePtr &cand) const final
Public Member Functions inherited from CutApplicatorWithEventContentBase
	CutApplicatorWithEventContentBase ()
	CutApplicatorWithEventContentBase (const edm::ParameterSet &c)
	CutApplicatorWithEventContentBase (const CutApplicatorWithEventContentBase &)=delete
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 edm::ParameterSet &c)
	CutApplicatorBase (const CutApplicatorBase &)=delete
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 reco::PhotonPtr &) const
virtual result_type	operator() (const pat::PhotonPtr &) const
virtual result_type	operator() (const reco::MuonPtr &) const
virtual result_type	operator() (const pat::MuonPtr &) const
virtual result_type	operator() (const reco::PFTauPtr &) const
virtual result_type	operator() (const pat::TauPtr &) const
CutApplicatorBase &	operator= (const CutApplicatorBase &)=delete
	~CutApplicatorBase () override
	Destructor. More...
Public Member Functions inherited from candidate_functions::CandidateCut
	CandidateCut ()
virtual	~CandidateCut ()

Private Attributes
const float	_barrelCutOff
EffectiveAreas	_effectiveAreas
const float	_isoCutEBHighPt
const float	_isoCutEBLowPt
const float	_isoCutEEHighPt
const float	_isoCutEELowPt
bool	_isRelativeIso
const float	_ptCutOff
edm::Handle< double >	_rhoHandle

Static Private Attributes
static constexpr 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 GsfEleEffAreaPFIsoCut.cc.

Constructor & Destructor Documentation

GsfEleEffAreaPFIsoCut()

GsfEleEffAreaPFIsoCut::GsfEleEffAreaPFIsoCut

(

const edm::ParameterSet &

c

)

Definition at line 37 of file GsfEleEffAreaPFIsoCut.cc.

References HltBtagPostValidation_cff::c, CutApplicatorWithEventContentBase::contentTags_, rhoString_, and HLT_2024v13_cff::rhoTag.

    : CutApplicatorWithEventContentBase(c),
      _isoCutEBLowPt(c.getParameter<double>("isoCutEBLowPt")),
      _isoCutEBHighPt(c.getParameter<double>("isoCutEBHighPt")),
      _isoCutEELowPt(c.getParameter<double>("isoCutEELowPt")),
      _isoCutEEHighPt(c.getParameter<double>("isoCutEEHighPt")),
      _ptCutOff(c.getParameter<double>("ptCutOff")),
      _barrelCutOff(c.getParameter<double>("barrelCutOff")),
      _isRelativeIso(c.getParameter<bool>("isRelativeIso")),
      _effectiveAreas((c.getParameter<edm::FileInPath>("effAreasConfigFile")).fullPath()) {
  edm::InputTag rhoTag = c.getParameter<edm::InputTag>("rho");
  contentTags_.emplace(rhoString_, rhoTag);
}

Member Function Documentation

candidateType()

CandidateType GsfEleEffAreaPFIsoCut::candidateType ( ) const

inlinefinalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 16 of file GsfEleEffAreaPFIsoCut.cc.

References CutApplicatorBase::ELECTRON.

{ return ELECTRON; }

getEventContent()

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

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 56 of file GsfEleEffAreaPFIsoCut.cc.

References _rhoHandle, CutApplicatorWithEventContentBase::contentTags_, makeMEIFBenchmarkPlots::ev, and rhoString_.

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

operator()()

CutApplicatorBase::result_type GsfEleEffAreaPFIsoCut::operator() ( const reco::GsfElectronPtr &  cand ) const

finalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 60 of file GsfEleEffAreaPFIsoCut.cc.

References _barrelCutOff, _effectiveAreas, _isoCutEBHighPt, _isoCutEBLowPt, _isoCutEEHighPt, _isoCutEELowPt, _isRelativeIso, _ptCutOff, _rhoHandle, funct::abs(), f, nano_mu_digi_cff::float, EffectiveAreas::getEffectiveArea(), genparticles_cff::iso, edm::HandleBase::isValid(), SiStripPI::max, rho, reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt, reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt, and reco::GsfElectron::PflowIsolationVariables::sumPhotonEt.

                                                                                                   {
  // Establish the cut value
  double absEta = std::abs(cand->superCluster()->eta());
  const float isoCut = (cand->pt() < _ptCutOff ? (absEta < _barrelCutOff ? _isoCutEBLowPt : _isoCutEELowPt)
                                               : (absEta < _barrelCutOff ? _isoCutEBHighPt : _isoCutEEHighPt));

  // Compute the combined isolation with effective area correction
  const reco::GsfElectron::PflowIsolationVariables& pfIso = cand->pfIsolationVariables();
  const float chad = pfIso.sumChargedHadronPt;
  const float nhad = pfIso.sumNeutralHadronEt;
  const float pho = pfIso.sumPhotonEt;
  const float eA = _effectiveAreas.getEffectiveArea(absEta);
  const float rho = _rhoHandle.isValid() ? (float)(*_rhoHandle) : 0;  // std::max likes float arguments
  const float iso = chad + std::max(0.0f, nhad + pho - rho * eA);

  // Apply the cut and return the result
  // Scale by pT if the relative isolation is requested but avoid division by 0
  return iso < isoCut * (_isRelativeIso ? cand->pt() : 1.);
}

setConsumes()

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

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 51 of file GsfEleEffAreaPFIsoCut.cc.

References gpuPixelDoublets::cc, CutApplicatorWithEventContentBase::contentTags_, CutApplicatorWithEventContentBase::contentTokens_, rho, and rhoString_.

                                                                {
  auto rho = cc.consumes<double>(contentTags_[rhoString_]);
  contentTokens_.emplace(rhoString_, rho);
}

value()

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

finalvirtual

Implements candidate_functions::CandidateCut.

Definition at line 80 of file GsfEleEffAreaPFIsoCut.cc.

References _effectiveAreas, _isRelativeIso, _rhoHandle, funct::abs(), f, nano_mu_digi_cff::float, EffectiveAreas::getEffectiveArea(), genparticles_cff::iso, SiStripPI::max, reco::GsfElectron::pfIsolationVariables(), reco::LeafCandidate::pt(), rho, reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt, and reco::GsfElectron::superCluster().

                                                                      {
  reco::GsfElectronPtr ele(cand);
  // Establish the cut value
  double absEta = std::abs(ele->superCluster()->eta());

  // Compute the combined isolation with effective area correction
  const reco::GsfElectron::PflowIsolationVariables& pfIso = ele->pfIsolationVariables();
  const float chad = pfIso.sumChargedHadronPt;
  const float nhad = pfIso.sumNeutralHadronEt;
  const float pho = pfIso.sumPhotonEt;
  float eA = _effectiveAreas.getEffectiveArea(absEta);
  float rho = (float)(*_rhoHandle);  // std::max likes float arguments
  float iso = chad + std::max(0.0f, nhad + pho - rho * eA);

  // Divide by pT if the relative isolation is requested
  if (_isRelativeIso)
    iso /= ele->pt();

  // Apply the cut and return the result
  return iso;
}

Member Data Documentation

_barrelCutOff

const float GsfEleEffAreaPFIsoCut::_barrelCutOff

private

Definition at line 23 of file GsfEleEffAreaPFIsoCut.cc.

Referenced by operator()().

_effectiveAreas

EffectiveAreas GsfEleEffAreaPFIsoCut::_effectiveAreas

private

Definition at line 26 of file GsfEleEffAreaPFIsoCut.cc.

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

_isoCutEBHighPt

const float GsfEleEffAreaPFIsoCut::_isoCutEBHighPt

private

Definition at line 20 of file GsfEleEffAreaPFIsoCut.cc.

Referenced by operator()().

_isoCutEBLowPt

const float GsfEleEffAreaPFIsoCut::_isoCutEBLowPt

private

Definition at line 20 of file GsfEleEffAreaPFIsoCut.cc.

Referenced by operator()().

_isoCutEEHighPt

const float GsfEleEffAreaPFIsoCut::_isoCutEEHighPt

private

Definition at line 20 of file GsfEleEffAreaPFIsoCut.cc.

Referenced by operator()().

_isoCutEELowPt

const float GsfEleEffAreaPFIsoCut::_isoCutEELowPt

private

Definition at line 20 of file GsfEleEffAreaPFIsoCut.cc.

Referenced by operator()().

_isRelativeIso

bool GsfEleEffAreaPFIsoCut::_isRelativeIso

private

Definition at line 24 of file GsfEleEffAreaPFIsoCut.cc.

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

_ptCutOff

const float GsfEleEffAreaPFIsoCut::_ptCutOff

private

Definition at line 22 of file GsfEleEffAreaPFIsoCut.cc.

Referenced by operator()().

_rhoHandle

edm::Handle<double> GsfEleEffAreaPFIsoCut::_rhoHandle

private

Definition at line 28 of file GsfEleEffAreaPFIsoCut.cc.

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

rhoString_

constexpr char GsfEleEffAreaPFIsoCut::rhoString_ = "rho"

staticprivate

Definition at line 30 of file GsfEleEffAreaPFIsoCut.cc.

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