Inheritance diagram for PhoAnyPFIsoWithEAAndExpoScalingCut:

Public Member Functions
CandidateType	candidateType () const override final

void	getEventContent (const edm::EventBase &) override final

result_type	operator() (const reco::PhotonPtr &) const override final

	PhoAnyPFIsoWithEAAndExpoScalingCut (const edm::ParameterSet &c)

void	setConsumes (edm::ConsumesCollector &) override final

double	value (const reco::CandidatePtr &cand) const override final

Public Member Functions inherited from CutApplicatorWithEventContentBase
	CutApplicatorWithEventContentBase ()

	CutApplicatorWithEventContentBase (const edm::ParameterSet &c)

	CutApplicatorWithEventContentBase (const CutApplicatorWithEventContentBase &)=delete

CutApplicatorWithEventContentBase &	operator= (const CutApplicatorWithEventContentBase &)=delete

virtual	~CutApplicatorWithEventContentBase ()
	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

virtual const std::string &	name () const

virtual 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 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

virtual	~CutApplicatorBase ()
	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 char	anyPFIsoWithEA_ [] = "anyPFIsoWithEA"

static char	rhoString_ [] = "rho"

Additional Inherited Members
Public Types inherited from CutApplicatorBase
enum	CandidateType { NONE, ELECTRON, MUON, PHOTON, TAU, PATELECTRON, PATMUON, PATPHOTON, PATTAU }

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 PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

Constructor & Destructor Documentation

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

Definition at line 50 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

References anyPFIsoWithEA_, CutApplicatorWithEventContentBase::contentTags_, edm::ParameterSet::getParameter(), rhoString_, and HiL1Corrector_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

inlinefinaloverridevirtual

Reimplemented from CutApplicatorBase.

Definition at line 17 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

References CutApplicatorBase::PHOTON.

                                                      { 
     return PHOTON; 
   }

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

finaloverridevirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 80 of file PhoAnyPFIsoWithEAAndExpoScalingCut.cc.

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

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

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

finaloverridevirtual

Reimplemented from CutApplicatorBase.

Definition at line 87 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_, JetChargeProducer_cfi::exp, EffectiveAreas::getEffectiveArea(), edm::Ptr< T >::id(), edm::ValueMap< T >::idSize(), edm::HandleBase::isValid(), edm::Ptr< T >::key(), hpstanc_transforms::max, EnergyCorrector::pt, reco::LeafCandidate::pt(), rho, AlCaHLTBitMon_QueryRunRegistry::string, reco::Photon::superCluster(), pat::PATObject< ObjectType >::userFloat(), and value().

Referenced by getEventContent().

                                            {  
 
   // 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 )

finaloverridevirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 71 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

finaloverridevirtual

Implements candidate_functions::CandidateCut.

Definition at line 134 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(), hpstanc_transforms::max, EnergyCorrector::pt, reco::LeafCandidate::pt(), rho, AlCaHLTBitMon_QueryRunRegistry::string, reco::Photon::superCluster(), and pat::PATObject< ObjectType >::userFloat().

Referenced by operator()().

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