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 ;
};

}

#endif