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

Static Private Attributes
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 6 of file GsfEleCalPFClusterIsoCut.cc.

Member Enumeration Documentation

enum GsfEleCalPFClusterIsoCut::IsoType

Enumerator
UNDEF
ISO_ECAL
ISO_HCAL

Definition at line 9 of file GsfEleCalPFClusterIsoCut.cc.

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

Constructor & Destructor Documentation

GsfEleCalPFClusterIsoCut::GsfEleCalPFClusterIsoCut

(

const edm::ParameterSet &

c

)

Definition at line 46 of file GsfEleCalPFClusterIsoCut.cc.

References CutApplicatorWithEventContentBase::contentTags_, edm::ParameterSet::getParameter(), rhoString_, and regressionModifier_cfi::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 20 of file GsfEleCalPFClusterIsoCut.cc.

References CutApplicatorBase::ELECTRON.

                                            { 
    return ELECTRON; 
  }

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

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 70 of file GsfEleCalPFClusterIsoCut.cc.

References CutApplicatorWithEventContentBase::contentTags_, edm::EventBase::getByLabel(), operator()(), 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 77 of file GsfEleCalPFClusterIsoCut.cc.

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

Referenced by getEventContent().

                                                {  

  // 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 64 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 114 of file GsfEleCalPFClusterIsoCut.cc.

References funct::abs(), effectiveAreas_, Exception, f, objects.autophobj::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

const float GsfEleCalPFClusterIsoCut::barrelCutOff_

private

Definition at line 30 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

EffectiveAreas GsfEleCalPFClusterIsoCut::effectiveAreas_

private

Definition at line 33 of file GsfEleCalPFClusterIsoCut.cc.

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

const float GsfEleCalPFClusterIsoCut::isoCutEBHighPt_

private

Definition at line 26 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

const float GsfEleCalPFClusterIsoCut::isoCutEBLowPt_

private

Definition at line 26 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

const float GsfEleCalPFClusterIsoCut::isoCutEEHighPt_

private

Definition at line 26 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

const float GsfEleCalPFClusterIsoCut::isoCutEELowPt_

private

Definition at line 26 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

const int GsfEleCalPFClusterIsoCut::isoType_

private

Definition at line 28 of file GsfEleCalPFClusterIsoCut.cc.

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

bool GsfEleCalPFClusterIsoCut::isRelativeIso_

private

Definition at line 31 of file GsfEleCalPFClusterIsoCut.cc.

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

const float GsfEleCalPFClusterIsoCut::ptCutOff_

private

Definition at line 29 of file GsfEleCalPFClusterIsoCut.cc.

Referenced by operator()().

edm::Handle< double > GsfEleCalPFClusterIsoCut::rhoHandle_

private

Definition at line 35 of file GsfEleCalPFClusterIsoCut.cc.

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

char GsfEleCalPFClusterIsoCut::rhoString_ = "rho"

staticprivate

Definition at line 37 of file GsfEleCalPFClusterIsoCut.cc.

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