PhoSingleTowerHadOverEmCut.cc

Go to the documentation of this file.

#include "PhysicsTools/SelectorUtils/interface/CutApplicatorBase.h"
#include "DataFormats/EgammaCandidates/interface/Photon.h"
 
class PhoSingleTowerHadOverEmCut : public CutApplicatorBase {
public:
  PhoSingleTowerHadOverEmCut(const edm::ParameterSet& c)
      : CutApplicatorBase(c),
        _hadronicOverEMCutValueEB(c.getParameter<double>("hadronicOverEMCutValueEB")),
        _hadronicOverEMCutValueEE(c.getParameter<double>("hadronicOverEMCutValueEE")),
        _barrelCutOff(c.getParameter<double>("barrelCutOff")) {}
 
  result_type operator()(const reco::PhotonPtr&) const final;
 
  double value(const reco::CandidatePtr& cand) const final;
 
  CandidateType candidateType() const final { return PHOTON; }
 
private:
  const float _hadronicOverEMCutValueEB, _hadronicOverEMCutValueEE, _barrelCutOff;
};
 
DEFINE_EDM_PLUGIN(CutApplicatorFactory, PhoSingleTowerHadOverEmCut, "PhoSingleTowerHadOverEmCut");
 
CutApplicatorBase::result_type PhoSingleTowerHadOverEmCut::operator()(const reco::PhotonPtr& cand) const {
  const float hadronicOverEMCutValue =
      (std::abs(cand->superCluster()->eta()) < _barrelCutOff ? _hadronicOverEMCutValueEB : _hadronicOverEMCutValueEE);
 
  return cand->hadTowOverEm() < hadronicOverEMCutValue;
}
 
double PhoSingleTowerHadOverEmCut::value(const reco::CandidatePtr& cand) const {
  reco::PhotonPtr pho(cand);
  return pho->hadTowOverEm();
}