MuonCosmicCompatibilityFiller.h

Go to the documentation of this file.

#ifndef MuonIdentification_MuonCosmicCompatibilityFiller_h
#define MuonIdentification_MuonCosmicCompatibilityFiller_h

#include <vector>

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "DataFormats/Common/interface/ValueMap.h"

#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonCosmicCompatibility.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/Utilities/interface/ESGetToken.h"

#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "Geometry/Records/interface/GlobalTrackingGeometryRecord.h"
#include "Geometry/CommonDetUnit/interface/GlobalTrackingGeometry.h"

namespace edm {
  class ParameterSet;
  class Event;
  class EventSetup;
}  // namespace edm
class GlobalMuonRefitter;

class MuonCosmicCompatibilityFiller {
public:
  MuonCosmicCompatibilityFiller(const edm::ParameterSet&, edm::ConsumesCollector&);
  ~MuonCosmicCompatibilityFiller();

  reco::MuonCosmicCompatibility fillCompatibility(const reco::Muon& muon, edm::Event&, const edm::EventSetup&);

private:
  float muonTiming(const edm::Event& iEvent, const reco::Muon& muon, bool isLoose) const;

  unsigned int backToBack2LegCosmic(const edm::Event&, const reco::Muon&) const;

  unsigned int pvMatches(const edm::Event&, const reco::Muon&, bool) const;

  bool isOverlappingMuon(const edm::Event&, const edm::EventSetup& iSetup, const reco::Muon&) const;

  unsigned int nMuons(const edm::Event&) const;

  unsigned int eventActivity(const edm::Event&, const reco::Muon&) const;

  float combinedCosmicID(
      const edm::Event&, const edm::EventSetup& iSetup, const reco::Muon&, bool CheckMuonID, bool checkVertex) const;

  bool checkMuonID(const reco::Muon&) const;

  bool checkMuonSegments(const reco::Muon& muon) const;

private:
  std::vector<edm::InputTag> inputMuonCollections_;
  std::vector<edm::InputTag> inputTrackCollections_;
  edm::InputTag inputCosmicMuonCollection_;
  edm::InputTag inputVertexCollection_;

  std::vector<edm::EDGetTokenT<reco::MuonCollection> > muonTokens_;
  std::vector<edm::EDGetTokenT<reco::TrackCollection> > trackTokens_;
  edm::EDGetTokenT<reco::MuonCollection> cosmicToken_;
  edm::EDGetTokenT<reco::VertexCollection> vertexToken_;
  edm::ESGetToken<GlobalTrackingGeometry, GlobalTrackingGeometryRecord> geometryToken_;

  double maxdxyLoose_;
  double maxdzLoose_;
  double maxdxyTight_;
  double maxdzTight_;
  double maxdxyLooseMult_;
  double maxdzLooseMult_;
  double maxdxyTightMult_;
  double maxdzTightMult_;
  double largedxyMult_;
  double largedxy_;
  double hIpTrdxy_;
  double hIpTrvProb_;
  double minvProb_;
  double maxvertZ_;
  double maxvertRho_;
  unsigned int nTrackThreshold_;
  double offTimePosLoose_;
  double offTimeNegLoose_;
  double offTimePosTight_;
  double offTimeNegTight_;
  double offTimePosLooseMult_;
  double offTimeNegLooseMult_;
  double offTimePosTightMult_;
  double offTimeNegTightMult_;
  double corrTimePos_;
  double corrTimeNeg_;
  double deltaPt_;
  double angleThreshold_;
  int sharedHits_;
  double sharedFrac_;
  double ipThreshold_;
  int nChamberMatches_;
  double segmentComp_;
};
#endif