GsfEleDEtaInLinearCut.cc

Go to the documentation of this file.

#include "PhysicsTools/SelectorUtils/interface/CutApplicatorBase.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "RecoEgamma/ElectronIdentification/interface/EBEECutValues.h"

class GsfEleDEtaInLinearCut : public CutApplicatorBase {
public:
  GsfEleDEtaInLinearCut(const edm::ParameterSet& param) :
    CutApplicatorBase(param),
    slopeTerm_(param,"slopeTerm"),
    constTerm_(param,"constTerm"),
    minValue_(param,"minValue")
  {
  }
  
  result_type operator()(const reco::GsfElectronPtr&) const final;

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

  CandidateType candidateType() const final { 
    return ELECTRON; 
  }

private:
  const EBEECutValues slopeTerm_;
  const EBEECutValues constTerm_;
  const EBEECutValues minValue_;
};

DEFINE_EDM_PLUGIN(CutApplicatorFactory,
          GsfEleDEtaInLinearCut,
          "GsfEleDEtaInLinearCut");

CutApplicatorBase::result_type 
GsfEleDEtaInLinearCut::
operator()(const reco::GsfElectronPtr& cand) const
{  
  float et = cand->energy()!=0. ? cand->et()/cand->energy()*cand->caloEnergy() : 0.;
  double cutValue = std::max(constTerm_(cand)+slopeTerm_(cand)*et,minValue_(cand));
  return std::abs(cand->deltaEtaSuperClusterTrackAtVtx())<cutValue;
 
}

double GsfEleDEtaInLinearCut::value(const reco::CandidatePtr& cand) const {
  reco::GsfElectronPtr ele(cand);
  return std::abs(ele->deltaEtaSuperClusterTrackAtVtx());
}