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

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

(

const edm::ParameterSet &

c

)

Definition at line 40 of file GsfEleCalPFClusterIsoCut.cc.

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

inlinefinalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 18 of file GsfEleCalPFClusterIsoCut.cc.

References CutApplicatorBase::ELECTRON.

{ return ELECTRON; }

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

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 60 of file GsfEleCalPFClusterIsoCut.cc.

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

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

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, validate-o2o-wbm::f, edm::Ptr< T >::get(), 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.);
}

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

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 55 of file GsfEleCalPFClusterIsoCut.cc.

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

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

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, validate-o2o-wbm::f, edm::Ptr< T >::get(), EffectiveAreas::getEffectiveArea(), ISO_ECAL, ISO_HCAL, isoType_, isRelativeIso_, edm::HandleBase::isValid(), SiStripPI::max, rho, and rhoHandle_.

Referenced by Types.int32::__nonzero__(), Types.uint32::__nonzero__(), Types.int64::__nonzero__(), Types.uint64::__nonzero__(), Types.double::__nonzero__(), Types.bool::__nonzero__(), Types.string::__nonzero__(), Types.string::configValue(), Types.FileInPath::configValue(), Types.int32::insertInto(), Types.uint32::insertInto(), Types.int64::insertInto(), Types.uint64::insertInto(), Types.double::insertInto(), Types.bool::insertInto(), Types.string::insertInto(), Types.FileInPath::insertInto(), Types.vint32::insertInto(), Types.vuint32::insertInto(), Types.vint64::insertInto(), Types.vuint64::insertInto(), Types.vdouble::insertInto(), Types.vbool::insertInto(), and Types.vstring::insertInto().

                                                                         {
  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

const float GsfEleCalPFClusterIsoCut::barrelCutOff_

private

Definition at line 26 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

EffectiveAreas GsfEleCalPFClusterIsoCut::effectiveAreas_

private

Definition at line 29 of file GsfEleCalPFClusterIsoCut.cc.

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

const float GsfEleCalPFClusterIsoCut::isoCutEBHighPt_

private

Definition at line 22 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

const float GsfEleCalPFClusterIsoCut::isoCutEBLowPt_

private

Definition at line 22 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

const float GsfEleCalPFClusterIsoCut::isoCutEEHighPt_

private

Definition at line 22 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

const float GsfEleCalPFClusterIsoCut::isoCutEELowPt_

private

Definition at line 22 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

const int GsfEleCalPFClusterIsoCut::isoType_

private

Definition at line 24 of file GsfEleCalPFClusterIsoCut.cc.

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

bool GsfEleCalPFClusterIsoCut::isRelativeIso_

private

Definition at line 27 of file GsfEleCalPFClusterIsoCut.cc.

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

const float GsfEleCalPFClusterIsoCut::ptCutOff_

private

Definition at line 25 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

edm::Handle<double> GsfEleCalPFClusterIsoCut::rhoHandle_

private

Definition at line 31 of file GsfEleCalPFClusterIsoCut.cc.

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

constexpr char GsfEleCalPFClusterIsoCut::rhoString_ = "rho"

staticprivate

Definition at line 33 of file GsfEleCalPFClusterIsoCut.cc.

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