GsfEleCalPFClusterIsoCut Class Reference

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 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_
const int	isoType_
bool	isRelativeIso_
const float	ptCutOff_
edm::Handle< double >	rhoHandle_

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

{ 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.

References DummyCfis::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")),
      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);
}

Member Function Documentation

candidateType()

CandidateType GsfEleCalPFClusterIsoCut::candidateType ( ) const

inlinefinalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 18 of file GsfEleCalPFClusterIsoCut.cc.

References CutApplicatorBase::ELECTRON.

{ return ELECTRON; }

getEventContent()

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

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 60 of file GsfEleCalPFClusterIsoCut.cc.

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

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

operator()()

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

finalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 64 of file GsfEleCalPFClusterIsoCut.cc.

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

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

setConsumes()

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

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 55 of file GsfEleCalPFClusterIsoCut.cc.

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

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

value()

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

finalvirtual

Implements candidate_functions::CandidateCut.

Definition at line 99 of file GsfEleCalPFClusterIsoCut.cc.

References funct::abs(), effectiveAreas_, Exception, f, nano_mu_digi_cff::float, edm::Ptr< T >::get(), EffectiveAreas::getEffectiveArea(), ISO_ECAL, ISO_HCAL, isoType_, isRelativeIso_, edm::HandleBase::isValid(), SiStripPI::max, reco::LeafCandidate::pt(), rho, rhoHandle_, and reco::GsfElectron::superCluster().

                                                                         {
  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;
}

Member Data Documentation

barrelCutOff_

const float GsfEleCalPFClusterIsoCut::barrelCutOff_

private

Definition at line 26 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

effectiveAreas_

EffectiveAreas GsfEleCalPFClusterIsoCut::effectiveAreas_

private

Definition at line 29 of file GsfEleCalPFClusterIsoCut.cc.

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

isoCutEBHighPt_

const float GsfEleCalPFClusterIsoCut::isoCutEBHighPt_

private

Definition at line 22 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

isoCutEBLowPt_

const float GsfEleCalPFClusterIsoCut::isoCutEBLowPt_

private

Definition at line 22 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

isoCutEEHighPt_

const float GsfEleCalPFClusterIsoCut::isoCutEEHighPt_

private

Definition at line 22 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

isoCutEELowPt_

const float GsfEleCalPFClusterIsoCut::isoCutEELowPt_

private

Definition at line 22 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

isoType_

const int GsfEleCalPFClusterIsoCut::isoType_

private

Definition at line 24 of file GsfEleCalPFClusterIsoCut.cc.

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

isRelativeIso_

bool GsfEleCalPFClusterIsoCut::isRelativeIso_

private

Definition at line 27 of file GsfEleCalPFClusterIsoCut.cc.

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

ptCutOff_

const float GsfEleCalPFClusterIsoCut::ptCutOff_

private

Definition at line 25 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

rhoHandle_

edm::Handle<double> GsfEleCalPFClusterIsoCut::rhoHandle_

private

Definition at line 31 of file GsfEleCalPFClusterIsoCut.cc.

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

rhoString_

constexpr char GsfEleCalPFClusterIsoCut::rhoString_ = "rho"

staticprivate

Definition at line 33 of file GsfEleCalPFClusterIsoCut.cc.

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