MuonTrackCut.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 MuonTrackCut : public CutApplicatorBase {
public:
  MuonTrackCut(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:
  // inner track selection cuts
  int minTrackerLayersWithMeasurement_, minPixelLayersWithMeasurement_, minNumberOfValidPixelHits_;
  double minValidFraction_;
  reco::Track::TrackQuality trackQuality_;
  int minNumberOfValidMuonHits_;
 
  bool doInnerTrack_, doGlobalTrack_;
};
DEFINE_EDM_PLUGIN(CutApplicatorFactory, MuonTrackCut, "MuonTrackCut");
 
// Define constructors and initialization routines
MuonTrackCut::MuonTrackCut(const edm::ParameterSet& c)
    : CutApplicatorBase(c),
      minTrackerLayersWithMeasurement_(-1),
      minPixelLayersWithMeasurement_(-1),
      minNumberOfValidPixelHits_(-1),
      minValidFraction_(-1),
      trackQuality_(reco::Track::undefQuality),
      minNumberOfValidMuonHits_(-1),
      doInnerTrack_(false),
      doGlobalTrack_(false) {
  if (c.existsAs<edm::ParameterSet>("innerTrack")) {
    doInnerTrack_ = true;
 
    const edm::ParameterSet cc = c.getParameter<edm::ParameterSet>("innerTrack");
    if (cc.exists("minTrackerLayersWithMeasurement"))
      minTrackerLayersWithMeasurement_ = cc.getParameter<int>("minTrackerLayersWithMeasurement");
    if (cc.exists("minPixelLayersWithMeasurement"))
      minPixelLayersWithMeasurement_ = cc.getParameter<int>("minPixelLayersWithMeasurement");
    if (cc.exists("minNumberOfValidPixelHits"))
      minNumberOfValidPixelHits_ = cc.getParameter<int>("minNumberOfValidPixelHits");
    if (cc.exists("minValidFraction"))
      minValidFraction_ = cc.getParameter<double>("minValidFraction");
    const std::string trackQualityStr = cc.exists("trackQuality") ? cc.getParameter<std::string>("trackQuality") : "";
    trackQuality_ = reco::Track::qualityByName(trackQualityStr);
  }
  if (c.existsAs<edm::ParameterSet>("globalTrack")) {
    doGlobalTrack_ = true;
 
    const edm::ParameterSet cc = c.getParameter<edm::ParameterSet>("globalTrack");
    //if ( cc.exists("maxNormalizedChi2") ) maxNormalizedChi2_ = cc.getParameter<double>("maxNormalizedChi2");
    if (cc.exists("minNumberOfValidMuonHits"))
      minNumberOfValidMuonHits_ = cc.getParameter<int>("minNumberOfValidMuonHits");
  }
}
 
// Functors for evaluation
CutApplicatorBase::result_type MuonTrackCut::operator()(const reco::MuonPtr& muon) const {
  if (doInnerTrack_) {
    const reco::TrackRef t = muon->innerTrack();
    if (t.isNull())
      return false;
    const auto& h = t->hitPattern();
    if (trackQuality_ != reco::Track::undefQuality and !t->quality(trackQuality_))
      return false;
    if (h.trackerLayersWithMeasurement() < minTrackerLayersWithMeasurement_)
      return false;
    if (h.pixelLayersWithMeasurement() < minPixelLayersWithMeasurement_)
      return false;
    if (h.numberOfValidPixelHits() < minNumberOfValidPixelHits_)
      return false;
    if (t->validFraction() <= minValidFraction_)
      return false;
  }
  if (doGlobalTrack_) {
    const reco::TrackRef t = muon->globalTrack();
    if (t.isNull())
      return false;
    const auto& h = t->hitPattern();
    if (h.numberOfValidMuonHits() < minNumberOfValidMuonHits_)
      return false;
    // if ( t->normalizedChi2() > maxNormalizedChi2_ ) return false; Not used for
  }
 
  return true;
}
 
double MuonTrackCut::value(const reco::CandidatePtr& cand) const {
  const reco::MuonPtr muon(cand);
  if (doInnerTrack_) {
    const reco::TrackRef t = muon->innerTrack();
    if (t.isNull())
      return 0;
    const auto& h = t->hitPattern();
    if (trackQuality_ != reco::Track::undefQuality and !t->quality(trackQuality_))
      return t->quality(trackQuality_);
    if (h.trackerLayersWithMeasurement() < minTrackerLayersWithMeasurement_)
      return h.trackerLayersWithMeasurement();
    if (h.pixelLayersWithMeasurement() < minPixelLayersWithMeasurement_)
      return h.pixelLayersWithMeasurement();
    if (h.numberOfValidPixelHits() < minNumberOfValidPixelHits_)
      return h.numberOfValidPixelHits();
    if (t->validFraction() <= minValidFraction_)
      return t->validFraction();
 
    return t->validFraction();
  }
  if (doGlobalTrack_) {
    const reco::TrackRef t = muon->globalTrack();
    if (t.isNull())
      return 0;
    const auto& h = t->hitPattern();
    if (h.numberOfValidMuonHits() < minNumberOfValidMuonHits_)
      return h.numberOfValidMuonHits();
    // if ( t->normalizedChi2() > maxNormalizedChi2_ ) return false; Not used for
 
    return h.numberOfValidMuonHits();
  }
 
  return 0;
}