PhoAnyPFIsoWithEAAndExpoScalingCut Class Reference

Inheritance diagram for PhoAnyPFIsoWithEAAndExpoScalingCut:

Public Member Functions
CandidateType	candidateType () const final
void	getEventContent (const edm::EventBase &) final
result_type	operator() (const reco::PhotonPtr &) const final
	PhoAnyPFIsoWithEAAndExpoScalingCut (const edm::ParameterSet &c)
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 reco::GsfElectronPtr &) const
virtual result_type	operator() (const pat::ElectronPtr &) const
virtual result_type	operator() (const reco::MuonPtr &) const
virtual result_type	operator() (const pat::MuonPtr &) const
CutApplicatorBase &	operator= (const CutApplicatorBase &)=delete
	~CutApplicatorBase () override
	Destructor. More...
Public Member Functions inherited from candidate_functions::CandidateCut
	CandidateCut ()
virtual	~CandidateCut ()

Private Attributes
edm::Handle< edm::ValueMap < float > >	_anyPFIsoMap
float	_barrelCutOff
float	_C1_EB
float	_C1_EE
float	_C2_EB
float	_C2_EE
float	_C3_EB
float	_C3_EE
EffectiveAreas	_effectiveAreas
edm::Handle< double >	_rhoHandle
bool	_useRelativeIso

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

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

Constructor & Destructor Documentation

PhoAnyPFIsoWithEAAndExpoScalingCut::PhoAnyPFIsoWithEAAndExpoScalingCut

(

const edm::ParameterSet &

c

)

Definition at line 45 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

References anyPFIsoWithEA_, CutApplicatorWithEventContentBase::contentTags_, edm::ParameterSet::getParameter(), rhoString_, and HLT_FULL_cff::rhoTag.

    : CutApplicatorWithEventContentBase(c),
      _C1_EB(c.getParameter<double>("C1_EB")),
      _C2_EB(c.getParameter<double>("C2_EB")),
      _C3_EB(c.getParameter<double>("C3_EB")),
      _C1_EE(c.getParameter<double>("C1_EE")),
      _C2_EE(c.getParameter<double>("C2_EE")),
      _C3_EE(c.getParameter<double>("C3_EE")),
      _barrelCutOff(c.getParameter<double>("barrelCutOff")),
      _useRelativeIso(c.getParameter<bool>("useRelativeIso")),
      _effectiveAreas((c.getParameter<edm::FileInPath>("effAreasConfigFile")).fullPath()) {
  edm::InputTag maptag = c.getParameter<edm::InputTag>("anyPFIsoMap");
  contentTags_.emplace(anyPFIsoWithEA_, maptag);

  edm::InputTag rhoTag = c.getParameter<edm::InputTag>("rho");
  contentTags_.emplace(rhoString_, rhoTag);
}

Member Function Documentation

CandidateType PhoAnyPFIsoWithEAAndExpoScalingCut::candidateType ( ) const

inlinefinalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 16 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

References CutApplicatorBase::PHOTON.

{ return PHOTON; }

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

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 71 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

References _anyPFIsoMap, _rhoHandle, anyPFIsoWithEA_, CutApplicatorWithEventContentBase::contentTags_, edm::EventBase::getByLabel(), and rhoString_.

                                                                               {
  ev.getByLabel(contentTags_[anyPFIsoWithEA_], _anyPFIsoMap);
  ev.getByLabel(contentTags_[rhoString_], _rhoHandle);
}

CutApplicatorBase::result_type PhoAnyPFIsoWithEAAndExpoScalingCut::operator() ( const reco::PhotonPtr &  cand ) const

finalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 76 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

References _anyPFIsoMap, _barrelCutOff, _C1_EB, _C1_EE, _C2_EB, _C2_EE, _C3_EB, _C3_EE, _effectiveAreas, _rhoHandle, _useRelativeIso, funct::abs(), anyPFIsoWithEA_, edm::ValueMap< T >::begin(), edm::ValueMap< T >::contains(), CutApplicatorWithEventContentBase::contentTags_, funct::exp(), EffectiveAreas::getEffectiveArea(), edm::Ptr< T >::id(), edm::ValueMap< T >::idSize(), edm::HandleBase::isValid(), edm::Ptr< T >::key(), SiStripPI::max, DiDispStaMuonMonitor_cfi::pt, rho, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                           {
  // in case we are by-value
  const std::string& inst_name = contentTags_.find(anyPFIsoWithEA_)->second.instance();
  edm::Ptr<pat::Photon> pat(cand);
  float anyisoval = -1.0;
  if (_anyPFIsoMap.isValid() && _anyPFIsoMap->contains(cand.id())) {
    anyisoval = (*_anyPFIsoMap)[cand];
  } else if (_anyPFIsoMap.isValid() && _anyPFIsoMap->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)
    anyisoval = _anyPFIsoMap->begin()[cand.key()];
  } else if (_anyPFIsoMap.isValid()) {  // throw an exception
    anyisoval = (*_anyPFIsoMap)[cand];
  }

  // Figure out the cut value
  // The value is generally pt-dependent: C1 + pt * C2
  const float pt = cand->pt();

  // In this version of the isolation cut we apply
  // exponential pt scaling to the barrel isolation cut,
  // and linear pt scaling to the endcap isolation cut.
  double absEta = std::abs(cand->superCluster()->eta());
  const float isolationCutValue =
      (absEta < _barrelCutOff ? _C1_EB + exp(pt * _C2_EB + _C3_EB) : _C1_EE + exp(pt * _C2_EE + _C3_EE));

  // Retrieve the variable value for this particle
  float anyPFIso = _anyPFIsoMap.isValid() ? anyisoval : pat->userFloat(inst_name);

  // Apply pile-up correction
  double eA = _effectiveAreas.getEffectiveArea(absEta);
  double rho = *_rhoHandle;
  float anyPFIsoWithEA = std::max(0.0, anyPFIso - rho * eA);

  // Divide by pT if the relative isolation is requested
  if (_useRelativeIso)
    anyPFIsoWithEA /= pt;

  // Apply the cut and return the result
  return anyPFIsoWithEA < isolationCutValue;
}

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

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 63 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

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

                                                                             {
  auto anyPFIsoWithEA = cc.consumes<edm::ValueMap<float> >(contentTags_[anyPFIsoWithEA_]);
  contentTokens_.emplace(anyPFIsoWithEA_, anyPFIsoWithEA);

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

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

finalvirtual

Implements candidate_functions::CandidateCut.

Definition at line 118 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

References _anyPFIsoMap, _effectiveAreas, _rhoHandle, _useRelativeIso, funct::abs(), anyPFIsoWithEA_, edm::ValueMap< T >::begin(), edm::ValueMap< T >::contains(), CutApplicatorWithEventContentBase::contentTags_, EffectiveAreas::getEffectiveArea(), edm::Ptr< T >::id(), edm::ValueMap< T >::idSize(), edm::HandleBase::isValid(), edm::Ptr< T >::key(), SiStripPI::max, DiDispStaMuonMonitor_cfi::pt, rho, and AlCaHLTBitMon_QueryRunRegistry::string.

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::PhotonPtr pho(cand);

  // in case we are by-value
  const std::string& inst_name = contentTags_.find(anyPFIsoWithEA_)->second.instance();
  edm::Ptr<pat::Photon> pat(cand);
  float anyisoval = -1.0;
  if (_anyPFIsoMap.isValid() && _anyPFIsoMap->contains(cand.id())) {
    anyisoval = (*_anyPFIsoMap)[cand];
  } else if (_anyPFIsoMap.isValid() && _anyPFIsoMap->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)
    anyisoval = _anyPFIsoMap->begin()[cand.key()];
  } else if (_anyPFIsoMap.isValid()) {  // throw an exception
    anyisoval = (*_anyPFIsoMap)[cand];
  }

  // Figure out the cut value
  // The value is generally pt-dependent: C1 + pt * C2
  const float pt = pho->pt();

  // In this version of the isolation cut we apply
  // exponential pt scaling to the barrel isolation cut,
  // and linear pt scaling to the endcap isolation cut.
  double absEta = std::abs(pho->superCluster()->eta());

  // Retrieve the variable value for this particle
  float anyPFIso = _anyPFIsoMap.isValid() ? anyisoval : pat->userFloat(inst_name);

  // Apply pile-up correction
  double eA = _effectiveAreas.getEffectiveArea(absEta);
  double rho = *_rhoHandle;
  float anyPFIsoWithEA = std::max(0.0, anyPFIso - rho * eA);

  // Divide by pT if the relative isolation is requested
  if (_useRelativeIso)
    anyPFIsoWithEA /= pt;

  // Apply the cut and return the result
  return anyPFIsoWithEA;
}

Member Data Documentation

edm::Handle<edm::ValueMap<float> > PhoAnyPFIsoWithEAAndExpoScalingCut::_anyPFIsoMap

private

Definition at line 32 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

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

float PhoAnyPFIsoWithEAAndExpoScalingCut::_barrelCutOff

private

Definition at line 27 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

Referenced by operator()().

float PhoAnyPFIsoWithEAAndExpoScalingCut::_C1_EB

private

Definition at line 20 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

Referenced by operator()().

float PhoAnyPFIsoWithEAAndExpoScalingCut::_C1_EE

private

Definition at line 23 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

Referenced by operator()().

float PhoAnyPFIsoWithEAAndExpoScalingCut::_C2_EB

private

Definition at line 21 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

Referenced by operator()().

float PhoAnyPFIsoWithEAAndExpoScalingCut::_C2_EE

private

Definition at line 24 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

Referenced by operator()().

float PhoAnyPFIsoWithEAAndExpoScalingCut::_C3_EB

private

Definition at line 22 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

Referenced by operator()().

float PhoAnyPFIsoWithEAAndExpoScalingCut::_C3_EE

private

Definition at line 25 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

Referenced by operator()().

EffectiveAreas PhoAnyPFIsoWithEAAndExpoScalingCut::_effectiveAreas

private

Definition at line 30 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

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

edm::Handle<double> PhoAnyPFIsoWithEAAndExpoScalingCut::_rhoHandle

private

Definition at line 34 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

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

bool PhoAnyPFIsoWithEAAndExpoScalingCut::_useRelativeIso

private

Definition at line 28 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

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

constexpr char PhoAnyPFIsoWithEAAndExpoScalingCut::anyPFIsoWithEA_ = "anyPFIsoWithEA"

staticprivate

Definition at line 36 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

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

constexpr char PhoAnyPFIsoWithEAAndExpoScalingCut::rhoString_ = "rho"

staticprivate

Definition at line 37 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

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