GsfEleDeltaBetaIsoCut Class Reference

Inheritance diagram for GsfEleDeltaBetaIsoCut:

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

Private Attributes
const float	_barrelCutOff
edm::Handle< edm::ValueMap< float > >	_chad_iso
const float	_deltaBetaConstant
const float	_isoCutEBHighPt
const float	_isoCutEBLowPt
const float	_isoCutEEHighPt
const float	_isoCutEELowPt
edm::Handle< edm::ValueMap< float > >	_nhad_iso
edm::Handle< edm::ValueMap< float > >	_ph_iso
const float	_ptCutOff
edm::Handle< edm::ValueMap< float > >	_PUchad_iso
bool	_relativeIso

Additional Inherited Members
Public Types inherited from CutApplicatorBase
enum	CandidateType {   NONE, ELECTRON, MUON, PHOTON,   TAU, PATELECTRON, PATMUON, PATPHOTON,   PATTAU }
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 8 of file GsfEleDeltaBetaIsoCut.cc.

Constructor & Destructor Documentation

GsfEleDeltaBetaIsoCut::GsfEleDeltaBetaIsoCut

(

const edm::ParameterSet &

c

)

Definition at line 34 of file GsfEleDeltaBetaIsoCut.cc.

References CutApplicatorWithEventContentBase::contentTags_, and edm::ParameterSet::getParameter().

                                                                     :
  CutApplicatorWithEventContentBase(c),
  _isoCutEBLowPt(c.getParameter<double>("isoCutEBLowPt")),
  _isoCutEBHighPt(c.getParameter<double>("isoCutEBHighPt")),
  _isoCutEELowPt(c.getParameter<double>("isoCutEELowPt")),
  _isoCutEEHighPt(c.getParameter<double>("isoCutEEHighPt")),
  _deltaBetaConstant(c.getParameter<double>("deltaBetaConstant")),
  _ptCutOff(c.getParameter<double>("ptCutOff")),
  _barrelCutOff(c.getParameter<double>("barrelCutOff")),
  _relativeIso(c.getParameter<bool>("isRelativeIso")) {  
  edm::InputTag chadtag = c.getParameter<edm::InputTag>("chargedHadronIsoDR03Src");
  edm::InputTag nhadtag = c.getParameter<edm::InputTag>("neutralHadronIsoDR03Src");
  edm::InputTag phtag = c.getParameter<edm::InputTag>("photonIsoDR03Src");
  edm::InputTag PUchadtag = c.getParameter<edm::InputTag>("puChargedHadronIsoDR03Src");
  contentTags_.emplace("chad",chadtag);
  contentTags_.emplace("nhad",nhadtag);
  contentTags_.emplace("pho",phtag);
  contentTags_.emplace("puchad",PUchadtag);
}

Member Function Documentation

CandidateType GsfEleDeltaBetaIsoCut::candidateType ( ) const

inlinefinaloverridevirtual

Reimplemented from CutApplicatorBase.

Definition at line 19 of file GsfEleDeltaBetaIsoCut.cc.

References CutApplicatorBase::ELECTRON.

                                                     { 
    return ELECTRON; 
  }

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

finaloverridevirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 69 of file GsfEleDeltaBetaIsoCut.cc.

References _chad_iso, _nhad_iso, _ph_iso, _PUchad_iso, CutApplicatorWithEventContentBase::contentTags_, edm::EventBase::getByLabel(), and operator()().

                                                                  {  
  ev.getByLabel(contentTags_["chad"],_chad_iso);
  ev.getByLabel(contentTags_["nhad"],_nhad_iso);
  ev.getByLabel(contentTags_["pho"],_ph_iso);
  ev.getByLabel(contentTags_["puchad"],_PUchad_iso);
}

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

finaloverridevirtual

Reimplemented from CutApplicatorBase.

Definition at line 78 of file GsfEleDeltaBetaIsoCut.cc.

References _barrelCutOff, _chad_iso, _deltaBetaConstant, _isoCutEBHighPt, _isoCutEBLowPt, _isoCutEEHighPt, _isoCutEELowPt, _nhad_iso, _ph_iso, _ptCutOff, _PUchad_iso, _relativeIso, funct::abs(), edm::ValueMap< T >::begin(), edm::ValueMap< T >::contains(), f, edm::Ptr< T >::id(), edm::ValueMap< T >::idSize(), edm::HandleBase::isValid(), edm::Ptr< T >::key(), hpstanc_transforms::max, reco::GsfElectron::p4(), reco::GsfElectron::pfIsolationVariables(), reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt, reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt, reco::GsfElectron::PflowIsolationVariables::sumPhotonEt, reco::GsfElectron::PflowIsolationVariables::sumPUPt, and reco::GsfElectron::superCluster().

Referenced by getEventContent().

                                                { 
  const float isoCut = 
    ( cand->p4().pt() < _ptCutOff ? 
      ( std::abs(cand->superCluster()->position().eta()) < _barrelCutOff ?
    _isoCutEBLowPt : _isoCutEELowPt ) :
      ( std::abs(cand->superCluster()->position().eta()) < _barrelCutOff ?
    _isoCutEBHighPt : _isoCutEEHighPt ) );
  const reco::GsfElectron::PflowIsolationVariables& pfIso =
    cand->pfIsolationVariables();

  float chad_val   = 0.0;
  float nhad_val   = 0.0;
  float pho_val    = 0.0;
  float puchad_val = 0.0;

  if( _chad_iso.isValid()   && _chad_iso->contains( cand.id() )   &&
      _nhad_iso.isValid()   && _nhad_iso->contains( cand.id() )   && 
      _ph_iso.isValid()     && _ph_iso->contains( cand.id() )     && 
      _PUchad_iso.isValid() && _PUchad_iso->contains( cand.id() )    ) {
    chad_val   = (*_chad_iso)[cand];
    nhad_val   = (*_nhad_iso)[cand];
    pho_val    = (*_ph_iso)[cand];
    puchad_val = (*_PUchad_iso)[cand];
  } else if ( _chad_iso.isValid()   && _chad_iso->idSize()   == 1 &&
              _nhad_iso.isValid()   && _nhad_iso->idSize()   == 1 &&
              _ph_iso.isValid()     && _ph_iso->idSize()     == 1 &&
              _PUchad_iso.isValid() && _PUchad_iso->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)
    chad_val   = _chad_iso->begin()[cand.key()];
    nhad_val   = _nhad_iso->begin()[cand.key()];
    pho_val    = _ph_iso->begin()[cand.key()];
    puchad_val = _PUchad_iso->begin()[cand.key()];
  } else if ( _chad_iso.isValid()   && _nhad_iso.isValid()   && 
              _ph_iso.isValid()     && _PUchad_iso.isValid()    ){ // throw an exception
    chad_val   = (*_chad_iso)[cand];
    nhad_val   = (*_nhad_iso)[cand];
    pho_val    = (*_ph_iso)[cand];
    puchad_val = (*_PUchad_iso)[cand];
  }
  
  const float chad   = _chad_iso.isValid()   ? chad_val   : pfIso.sumChargedHadronPt;
  const float nhad   = _nhad_iso.isValid()   ? nhad_val   : pfIso.sumNeutralHadronEt;
  const float pho    = _ph_iso.isValid()     ? pho_val    : pfIso.sumPhotonEt;
  const float puchad = _PUchad_iso.isValid() ? puchad_val : pfIso.sumPUPt;
  float iso = chad + std::max(0.0f, nhad + pho - _deltaBetaConstant*puchad);
  if( _relativeIso ) iso /= cand->p4().pt();
  return iso < isoCut;
}

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

finaloverridevirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 54 of file GsfEleDeltaBetaIsoCut.cc.

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

                                                                {
  auto chad = 
    cc.consumes<edm::ValueMap<float> >(contentTags_["chad"]);
  contentTokens_.emplace("chad",chad);
  auto nhad = 
    cc.consumes<edm::ValueMap<float> >(contentTags_["nhad"]);
  contentTokens_.emplace("nhad",nhad);
  auto pho = 
    cc.consumes<edm::ValueMap<float> >(contentTags_["pho"]);
  contentTokens_.emplace("pho",pho);
  auto puchad = 
    cc.consumes<edm::ValueMap<float> >(contentTags_["puchad"]);
  contentTokens_.emplace("puchad",puchad);
}

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

finaloverridevirtual

Implements candidate_functions::CandidateCut.

Definition at line 129 of file GsfEleDeltaBetaIsoCut.cc.

References _chad_iso, _deltaBetaConstant, _nhad_iso, _ph_iso, _PUchad_iso, _relativeIso, edm::ValueMap< T >::begin(), edm::ValueMap< T >::contains(), f, edm::Ptr< T >::id(), edm::ValueMap< T >::idSize(), edm::HandleBase::isValid(), edm::Ptr< T >::key(), hpstanc_transforms::max, reco::GsfElectron::pfIsolationVariables(), reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt, reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt, reco::GsfElectron::PflowIsolationVariables::sumPhotonEt, and reco::GsfElectron::PflowIsolationVariables::sumPUPt.

                                                                      {
  edm::Ptr<reco::GsfElectron> ele(cand);
  const reco::GsfElectron::PflowIsolationVariables& pfIso =
    ele->pfIsolationVariables();
  float chad_val   = 0.0;
  float nhad_val   = 0.0;
  float pho_val    = 0.0;
  float puchad_val = 0.0;

  if( _chad_iso.isValid()   && _chad_iso->contains( cand.id() )   &&
      _nhad_iso.isValid()   && _nhad_iso->contains( cand.id() )   &&
      _ph_iso.isValid()     && _ph_iso->contains( cand.id() )     &&
      _PUchad_iso.isValid() && _PUchad_iso->contains( cand.id() )    ) {
    chad_val   = (*_chad_iso)[cand];
    nhad_val   = (*_nhad_iso)[cand];
    pho_val    = (*_ph_iso)[cand];
    puchad_val = (*_PUchad_iso)[cand];
  } else if ( _chad_iso.isValid()   && _chad_iso->idSize()   == 1 &&
              _nhad_iso.isValid()   && _nhad_iso->idSize()   == 1 &&
              _ph_iso.isValid()     && _ph_iso->idSize()     == 1 &&
              _PUchad_iso.isValid() && _PUchad_iso->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)
    chad_val   = _chad_iso->begin()[cand.key()];
    nhad_val   = _nhad_iso->begin()[cand.key()];
    pho_val    = _ph_iso->begin()[cand.key()];
    puchad_val = _PUchad_iso->begin()[cand.key()];
  } else if ( _chad_iso.isValid()   && _nhad_iso.isValid()   &&
              _ph_iso.isValid()     && _PUchad_iso.isValid()    ){ // throw an exception
    chad_val   = (*_chad_iso)[cand];
    nhad_val   = (*_nhad_iso)[cand];
    pho_val    = (*_ph_iso)[cand];
    puchad_val = (*_PUchad_iso)[cand];
  }
  
  const float chad   = _chad_iso.isValid()   ? chad_val    : pfIso.sumChargedHadronPt;
  const float nhad   = _nhad_iso.isValid()   ? nhad_val    : pfIso.sumNeutralHadronEt;
  const float pho    = _ph_iso.isValid()     ? pho_val     : pfIso.sumPhotonEt;
  const float puchad = _PUchad_iso.isValid() ? puchad_val : pfIso.sumPUPt;
  float iso = chad + std::max(0.0f, nhad + pho - _deltaBetaConstant*puchad);
  if( _relativeIso ) iso /= cand->p4().pt();
  return iso;
}

Member Data Documentation

const float GsfEleDeltaBetaIsoCut::_barrelCutOff

private

Definition at line 25 of file GsfEleDeltaBetaIsoCut.cc.

Referenced by operator()().

edm::Handle<edm::ValueMap<float> > GsfEleDeltaBetaIsoCut::_chad_iso

private

Definition at line 27 of file GsfEleDeltaBetaIsoCut.cc.

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

const float GsfEleDeltaBetaIsoCut::_deltaBetaConstant

private

Definition at line 25 of file GsfEleDeltaBetaIsoCut.cc.

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

const float GsfEleDeltaBetaIsoCut::_isoCutEBHighPt

private

Definition at line 24 of file GsfEleDeltaBetaIsoCut.cc.

Referenced by operator()().

const float GsfEleDeltaBetaIsoCut::_isoCutEBLowPt

private

Definition at line 24 of file GsfEleDeltaBetaIsoCut.cc.

Referenced by operator()().

const float GsfEleDeltaBetaIsoCut::_isoCutEEHighPt

private

Definition at line 24 of file GsfEleDeltaBetaIsoCut.cc.

Referenced by operator()().

const float GsfEleDeltaBetaIsoCut::_isoCutEELowPt

private

Definition at line 24 of file GsfEleDeltaBetaIsoCut.cc.

Referenced by operator()().

edm::Handle<edm::ValueMap<float> > GsfEleDeltaBetaIsoCut::_nhad_iso

private

Definition at line 27 of file GsfEleDeltaBetaIsoCut.cc.

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

edm::Handle<edm::ValueMap<float> > GsfEleDeltaBetaIsoCut::_ph_iso

private

Definition at line 27 of file GsfEleDeltaBetaIsoCut.cc.

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

const float GsfEleDeltaBetaIsoCut::_ptCutOff

private

Definition at line 25 of file GsfEleDeltaBetaIsoCut.cc.

Referenced by operator()().

edm::Handle<edm::ValueMap<float> > GsfEleDeltaBetaIsoCut::_PUchad_iso

private

Definition at line 27 of file GsfEleDeltaBetaIsoCut.cc.

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

bool GsfEleDeltaBetaIsoCut::_relativeIso

private

Definition at line 26 of file GsfEleDeltaBetaIsoCut.cc.

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