MuonMomQualityCut.cc

Go to the documentation of this file.

#include "PhysicsTools/SelectorUtils/interface/CutApplicatorBase.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonSelectors.h"

class MuonMomQualityCut : public CutApplicatorBase
{
public:
  MuonMomQualityCut(const edm::ParameterSet& c);

  result_type operator()(const reco::MuonPtr&) const final;
  CandidateType candidateType() const final { return MUON; }
  double value(const reco::CandidatePtr&) const final;

private:
  const double maxRelPtErr_;
};
DEFINE_EDM_PLUGIN(CutApplicatorFactory, MuonMomQualityCut, "MuonMomQualityCut");

// Define constructors and initialization routines
MuonMomQualityCut::MuonMomQualityCut(const edm::ParameterSet& c):
  CutApplicatorBase(c),
  maxRelPtErr_(c.getParameter<double>("maxRelPtErr"))
{
}

// Functors for evaluation
CutApplicatorBase::result_type MuonMomQualityCut::operator()(const reco::MuonPtr& cand) const
{
  const auto trackRef = cand->muonBestTrack();
  return trackRef.isNonnull() and trackRef->ptError() <= maxRelPtErr_*trackRef->pt();

  return true;
}

double MuonMomQualityCut::value(const reco::CandidatePtr& cand) const
{
  const reco::MuonPtr muon(cand);
  const auto trackRef = muon->muonBestTrack();
  if ( trackRef.isNull() or trackRef->pt() <= 0 ) return -1;

  return trackRef->ptError()/trackRef->pt();
}