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 override 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 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 6 of file GsfEleEffAreaPFIsoCut.cc.

Constructor & Destructor Documentation

GsfEleEffAreaPFIsoCut::GsfEleEffAreaPFIsoCut

(

const edm::ParameterSet &

c

)

Definition at line 42 of file GsfEleEffAreaPFIsoCut.cc.

References CutApplicatorWithEventContentBase::contentTags_, edm::ParameterSet::getParameter(), rhoString_, and lowPtGsfElectrons_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 GsfEleEffAreaPFIsoCut::candidateType ( ) const

inlinefinalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 17 of file GsfEleEffAreaPFIsoCut.cc.

References CutApplicatorBase::ELECTRON.

                                            { 
    return ELECTRON; 
  }

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

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 65 of file GsfEleEffAreaPFIsoCut.cc.

References _rhoHandle, CutApplicatorWithEventContentBase::contentTags_, edm::EventBase::getByLabel(), operator()(), and rhoString_.

                                                                  {  

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

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

finalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 72 of file GsfEleEffAreaPFIsoCut.cc.

References _barrelCutOff, _effectiveAreas, _isoCutEBHighPt, _isoCutEBLowPt, _isoCutEEHighPt, _isoCutEELowPt, _isRelativeIso, _ptCutOff, _rhoHandle, funct::abs(), f, objects.autophobj::float, EffectiveAreas::getEffectiveArea(), edm::HandleBase::isValid(), SiStripPI::max, reco::GsfElectron::pfIsolationVariables(), reco::LeafCandidate::pt(), rho, reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt, reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt, reco::GsfElectron::PflowIsolationVariables::sumPhotonEt, and reco::GsfElectron::superCluster().

Referenced by getEventContent().

                                                {  

  // 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.);
}

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

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 59 of file GsfEleEffAreaPFIsoCut.cc.

References edm::ConsumesCollector::consumes(), CutApplicatorWithEventContentBase::contentTags_, CutApplicatorWithEventContentBase::contentTokens_, rho, and rhoString_.

                                                                {

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

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

finalvirtual

Implements candidate_functions::CandidateCut.

Definition at line 97 of file GsfEleEffAreaPFIsoCut.cc.

References _effectiveAreas, _isRelativeIso, _rhoHandle, funct::abs(), f, objects.autophobj::float, EffectiveAreas::getEffectiveArea(), 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

const float GsfEleEffAreaPFIsoCut::_barrelCutOff

private

Definition at line 26 of file GsfEleEffAreaPFIsoCut.cc.

Referenced by operator()().

EffectiveAreas GsfEleEffAreaPFIsoCut::_effectiveAreas

private

Definition at line 29 of file GsfEleEffAreaPFIsoCut.cc.

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

const float GsfEleEffAreaPFIsoCut::_isoCutEBHighPt

private

Definition at line 23 of file GsfEleEffAreaPFIsoCut.cc.

Referenced by operator()().

const float GsfEleEffAreaPFIsoCut::_isoCutEBLowPt

private

Definition at line 23 of file GsfEleEffAreaPFIsoCut.cc.

Referenced by operator()().

const float GsfEleEffAreaPFIsoCut::_isoCutEEHighPt

private

Definition at line 23 of file GsfEleEffAreaPFIsoCut.cc.

Referenced by operator()().

const float GsfEleEffAreaPFIsoCut::_isoCutEELowPt

private

Definition at line 23 of file GsfEleEffAreaPFIsoCut.cc.

Referenced by operator()().

bool GsfEleEffAreaPFIsoCut::_isRelativeIso

private

Definition at line 27 of file GsfEleEffAreaPFIsoCut.cc.

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

const float GsfEleEffAreaPFIsoCut::_ptCutOff

private

Definition at line 25 of file GsfEleEffAreaPFIsoCut.cc.

Referenced by operator()().

edm::Handle< double > GsfEleEffAreaPFIsoCut::_rhoHandle

private

Definition at line 31 of file GsfEleEffAreaPFIsoCut.cc.

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

char GsfEleEffAreaPFIsoCut::rhoString_ = "rho"

staticprivate

Definition at line 33 of file GsfEleEffAreaPFIsoCut.cc.

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