Inheritance diagram for GsfEleCalPFClusterIsoCut:

Public Types
enum	IsoType { UNDEF = -1, ISO_ECAL = 0, ISO_HCAL = 1 }

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

Public Member Functions
CandidateType	candidateType () const final

void	getEventContent (const edm::EventBase &) final

	GsfEleCalPFClusterIsoCut (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 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 pat::PhotonPtr &) const

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

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

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

virtual result_type	operator() (const reco::PhotonPtr &) 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_

const int	isoType_

bool	isRelativeIso_

const float	ptCutOff_

edm::Handle< double >	rhoHandle_

Static Private Attributes
constexpr static char	rhoString_ [] = "rho"

Additional Inherited Members
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 GsfEleCalPFClusterIsoCut.cc.

Member Enumeration Documentation

◆ IsoType

enum GsfEleCalPFClusterIsoCut::IsoType

Enumerator
UNDEF
ISO_ECAL
ISO_HCAL

Definition at line 7 of file GsfEleCalPFClusterIsoCut.cc.

7 { UNDEF = -1, ISO_ECAL = 0, ISO_HCAL = 1 };

Constructor & Destructor Documentation

◆ GsfEleCalPFClusterIsoCut()

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

Definition at line 40 of file GsfEleCalPFClusterIsoCut.cc.

     : CutApplicatorWithEventContentBase(c),
       isoCutEBLowPt_(c.getParameter<double>("isoCutEBLowPt")),
       isoCutEBHighPt_(c.getParameter<double>("isoCutEBHighPt")),
       isoCutEELowPt_(c.getParameter<double>("isoCutEELowPt")),
       isoCutEEHighPt_(c.getParameter<double>("isoCutEEHighPt")),
       isoType_(c.getParameter<int>("isoType")),
       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);
 }

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

Member Function Documentation

◆ candidateType()

CandidateType GsfEleCalPFClusterIsoCut::candidateType ( ) const

inlinefinalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 18 of file GsfEleCalPFClusterIsoCut.cc.

18 { return ELECTRON; }

References CutApplicatorBase::ELECTRON.

◆ getEventContent()

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

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 60 of file GsfEleCalPFClusterIsoCut.cc.

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

References CutApplicatorWithEventContentBase::contentTags_, ev, rhoHandle_, and rhoString_.

◆ operator()()

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

finalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 64 of file GsfEleCalPFClusterIsoCut.cc.

                                                                                                       {
   // Establish the cut value
   double absEta = std::abs(cand->superCluster()->eta());
  
   const pat::Electron* elPat = dynamic_cast<const pat::Electron*>(cand.get());
   if (!elPat) {
     throw cms::Exception("ERROR: this VID selection is meant to be run on miniAOD/PAT only")
         << std::endl
         << "Change input format to PAT/miniAOD or contact Egamma experts" << std::endl
         << std::endl;
   }
  
   const float isoCut = (cand->pt() < ptCutOff_ ? (absEta < barrelCutOff_ ? isoCutEBLowPt_ : isoCutEELowPt_)
                                                : (absEta < barrelCutOff_ ? isoCutEBHighPt_ : isoCutEEHighPt_));
  
   const float eA = effectiveAreas_.getEffectiveArea(absEta);
   const float rho = rhoHandle_.isValid() ? (float)(*rhoHandle_) : 0;  // std::max likes float arguments
  
   float isoValue = -999;
   if (isoType_ == ISO_ECAL) {
     isoValue = elPat->ecalPFClusterIso();
   } else if (isoType_ == ISO_HCAL) {
     isoValue = elPat->hcalPFClusterIso();
   } else {
     throw cms::Exception("ERROR: unknown type requested for PF cluster isolation.")
         << std::endl
         << "Check VID configuration." << std::endl;
   }
   float isoValueCorr = std::max(0.0f, isoValue - rho * eA);
  
   // Apply the cut and return the result
   // Scale by pT if the relative isolation is requested but avoid division by 0
   return isoValueCorr < isoCut * (isRelativeIso_ ? cand->pt() : 1.);
 }

References funct::abs(), barrelCutOff_, effectiveAreas_, Exception, f, dqmMemoryStats::float, EffectiveAreas::getEffectiveArea(), ISO_ECAL, ISO_HCAL, zSelection_cfi::isoCut, isoCutEBHighPt_, isoCutEBLowPt_, isoCutEEHighPt_, isoCutEELowPt_, isoType_, isRelativeIso_, edm::HandleBase::isValid(), SiStripPI::max, ptCutOff_, rho, and rhoHandle_.

◆ setConsumes()

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

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 55 of file GsfEleCalPFClusterIsoCut.cc.

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

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

◆ value()

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

finalvirtual

Implements candidate_functions::CandidateCut.

Definition at line 99 of file GsfEleCalPFClusterIsoCut.cc.

                                                                          {
   reco::GsfElectronPtr ele(cand);
   // Establish the cut value
   double absEta = std::abs(ele->superCluster()->eta());
  
   const pat::Electron* elPat = dynamic_cast<const pat::Electron*>(ele.get());
   if (!elPat) {
     throw cms::Exception("ERROR: this VID selection is meant to be run on miniAOD/PAT only")
         << std::endl
         << "Change input format to PAT/miniAOD or contact Egamma experts" << std::endl
         << std::endl;
   }
  
   const float eA = effectiveAreas_.getEffectiveArea(absEta);
   const float rho = rhoHandle_.isValid() ? (float)(*rhoHandle_) : 0;  // std::max likes float arguments
  
   float isoValue = -999;
   if (isoType_ == ISO_ECAL) {
     isoValue = elPat->ecalPFClusterIso();
   } else if (isoType_ == ISO_HCAL) {
     isoValue = elPat->hcalPFClusterIso();
   } else {
     throw cms::Exception("ERROR: unknown type requested for PF cluster isolation.")
         << std::endl
         << "Check VID configuration." << std::endl;
   }
   float isoValueCorr = std::max(0.0f, isoValue - rho * eA);
  
   // Divide by pT if the relative isolation is requested
   if (isRelativeIso_)
     isoValueCorr /= ele->pt();
  
   // Apply the cut and return the result
   return isoValueCorr;
 }