GsfEleMinEcalEtCut.cc

Go to the documentation of this file.

#include "PhysicsTools/SelectorUtils/interface/CutApplicatorBase.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
 
class GsfEleMinEcalEtCut : public CutApplicatorBase {
public:
  GsfEleMinEcalEtCut(const edm::ParameterSet& c);
 
  result_type operator()(const reco::GsfElectronPtr&) const final;
 
  double value(const reco::CandidatePtr& cand) const final;
 
  CandidateType candidateType() const final { return ELECTRON; }
 
private:
  const double minEt_;
};
 
DEFINE_EDM_PLUGIN(CutApplicatorFactory, GsfEleMinEcalEtCut, "GsfEleMinEcalEtCut");
 
GsfEleMinEcalEtCut::GsfEleMinEcalEtCut(const edm::ParameterSet& c)
    : CutApplicatorBase(c), minEt_(c.getParameter<double>("minEt")) {}
 
CutApplicatorBase::result_type GsfEleMinEcalEtCut::operator()(const reco::GsfElectronPtr& cand) const {
  const reco::GsfElectron& ele = *cand;
  const float sinTheta = ele.p() != 0 ? ele.pt() / ele.p() : 0.;
  const float et = ele.ecalEnergy() * sinTheta;
  return et > minEt_;
}
 
double GsfEleMinEcalEtCut::value(const reco::CandidatePtr& cand) const {
  const reco::GsfElectronPtr elePtr(cand);
  const reco::GsfElectron& ele = *elePtr;
  const float sinTheta = ele.p() != 0 ? ele.pt() / ele.p() : 0.;
  const float et = ele.ecalEnergy() * sinTheta;
  return et;
}