CMGMuonCleanerBySegmentsAlgo.h

Go to the documentation of this file.

#ifndef PhysicsTools_Heppy_CMGMuonCleanerBySegmentsAlgo_h
#define PhysicsTools_Heppy_CMGMuonCleanerBySegmentsAlgo_h
 
#include "DataFormats/PatCandidates/interface/Muon.h"
#include "CommonTools/Utils/interface/StringCutObjectSelector.h"
#include <vector>
 
namespace heppy {
 
  class CMGMuonCleanerBySegmentsAlgo {
  public:
    CMGMuonCleanerBySegmentsAlgo(double sharedFraction = 0.499,
                                 const std::string &preselection = "track.isNonnull",
                                 const std::string &passthrough = "isGlobalMuon && numberOfMatches >= 2")
        : sharedFraction_(sharedFraction), preselection_(preselection, true), passthrough_(passthrough, true) {}
    ~CMGMuonCleanerBySegmentsAlgo();
 
    std::vector<bool> clean(const std::vector<pat::Muon> &muons) const;
 
  private:
    double sharedFraction_;
 
    StringCutObjectSelector<pat::Muon> preselection_;
    StringCutObjectSelector<pat::Muon> passthrough_;
 
    bool isSameMuon(const pat::Muon &mu1, const pat::Muon &mu2) const;
    bool isBetterMuon(const pat::Muon &mu1, bool mu1PF, const pat::Muon &mu2, bool mu2PF) const;
  };
 
}  // namespace heppy
 
#endif