CSCStubMatcher.h

Go to the documentation of this file.

#ifndef Validation_MuonCSCDigis_CSCStubMatcher_h
#define Validation_MuonCSCDigis_CSCStubMatcher_h

#include "Validation/MuonCSCDigis/interface/CSCDigiMatcher.h"
#include "Validation/MuonGEMDigis/interface/GEMDigiMatcher.h"

#include "DataFormats/CSCDigi/interface/CSCConstants.h"
#include "DataFormats/CSCDigi/interface/CSCALCTDigiCollection.h"
#include "DataFormats/CSCDigi/interface/CSCCLCTDigiCollection.h"
#include "DataFormats/CSCDigi/interface/CSCCorrelatedLCTDigiCollection.h"

typedef std::vector<CSCALCTDigi> CSCALCTDigiContainer;
typedef std::vector<CSCCLCTDigi> CSCCLCTDigiContainer;
typedef std::vector<CSCCorrelatedLCTDigi> CSCCorrelatedLCTDigiContainer;

class CSCStubMatcher {
public:
  CSCStubMatcher(edm::ParameterSet const& iPS, edm::ConsumesCollector&& iC);

  ~CSCStubMatcher() {}

  void init(const edm::Event& e, const edm::EventSetup& eventSetup);

  void match(const SimTrack& t, const SimVertex& v);

  std::set<unsigned int> chamberIdsAllCLCT(int csc_type = MuonHitHelper::CSC_ALL) const;
  std::set<unsigned int> chamberIdsAllALCT(int csc_type = MuonHitHelper::CSC_ALL) const;
  std::set<unsigned int> chamberIdsAllLCT(int csc_type = MuonHitHelper::CSC_ALL) const;
  std::set<unsigned int> chamberIdsAllMPLCT(int csc_type = MuonHitHelper::CSC_ALL) const;

  std::set<unsigned int> chamberIdsCLCT(int csc_type = MuonHitHelper::CSC_ALL) const;
  std::set<unsigned int> chamberIdsALCT(int csc_type = MuonHitHelper::CSC_ALL) const;
  std::set<unsigned int> chamberIdsLCT(int csc_type = MuonHitHelper::CSC_ALL) const;
  std::set<unsigned int> chamberIdsMPLCT(int csc_type = MuonHitHelper::CSC_ALL) const;

  const CSCCLCTDigiContainer& allCLCTsInChamber(unsigned int) const;
  const CSCALCTDigiContainer& allALCTsInChamber(unsigned int) const;
  const CSCCorrelatedLCTDigiContainer& allLCTsInChamber(unsigned int) const;
  const CSCCorrelatedLCTDigiContainer& allMPLCTsInChamber(unsigned int) const;

  const CSCCLCTDigiContainer& clctsInChamber(unsigned int) const;
  const CSCALCTDigiContainer& alctsInChamber(unsigned int) const;
  const CSCCorrelatedLCTDigiContainer& lctsInChamber(unsigned int) const;
  const CSCCorrelatedLCTDigiContainer& mplctsInChamber(unsigned int) const;

  std::map<unsigned int, CSCCLCTDigiContainer> clcts() const { return chamber_to_clcts_; }
  std::map<unsigned int, CSCALCTDigiContainer> alcts() const { return chamber_to_alcts_; }
  std::map<unsigned int, CSCCorrelatedLCTDigiContainer> lcts() const { return chamber_to_lcts_; }
  std::map<unsigned int, CSCCorrelatedLCTDigiContainer> mplcts() const { return chamber_to_mplcts_; }

  CSCCLCTDigi bestClctInChamber(unsigned int) const;
  CSCALCTDigi bestAlctInChamber(unsigned int) const;
  CSCCorrelatedLCTDigi bestLctInChamber(unsigned int) const;
  CSCCorrelatedLCTDigi bestMplctInChamber(unsigned int) const;

  //z position of  certain layer
  float zpositionOfLayer(unsigned int detid, int layer) const;

  int nChambersWithCLCT(int min_quality = 0) const;
  int nChambersWithALCT(int min_quality = 0) const;
  int nChambersWithLCT(int min_quality = 0) const;
  int nChambersWithMPLCT(int min_quality = 0) const;

  bool lctInChamber(const CSCDetId& id, const CSCCorrelatedLCTDigi& lct) const;

  // get the position of an LCT in global coordinates
  GlobalPoint getGlobalPosition(unsigned int rawId, const CSCCorrelatedLCTDigi& lct) const;

  std::shared_ptr<CSCDigiMatcher> cscDigiMatcher() { return cscDigiMatcher_; }
  std::shared_ptr<GEMDigiMatcher> gemDigiMatcher() { return gemDigiMatcher_; }

private:
  void matchCLCTsToSimTrack(const CSCCLCTDigiCollection&);
  void matchALCTsToSimTrack(const CSCALCTDigiCollection&);
  void matchLCTsToSimTrack(const CSCCorrelatedLCTDigiCollection&);
  void matchMPLCTsToSimTrack(const CSCCorrelatedLCTDigiCollection&);

  void clear();

  void addGhostLCTs(const CSCCorrelatedLCTDigi& lct11,
                    const CSCCorrelatedLCTDigi& lct22,
                    CSCCorrelatedLCTDigiContainer& lctcontainer) const;

  edm::InputTag clctInputTag_;
  edm::InputTag alctInputTag_;
  edm::InputTag lctInputTag_;
  edm::InputTag mplctInputTag_;

  edm::EDGetTokenT<CSCCLCTDigiCollection> clctToken_;
  edm::EDGetTokenT<CSCALCTDigiCollection> alctToken_;
  edm::EDGetTokenT<CSCCorrelatedLCTDigiCollection> lctToken_;
  edm::EDGetTokenT<CSCCorrelatedLCTDigiCollection> mplctToken_;

  edm::Handle<CSCCLCTDigiCollection> clctsH_;
  edm::Handle<CSCALCTDigiCollection> alctsH_;
  edm::Handle<CSCCorrelatedLCTDigiCollection> lctsH_;
  edm::Handle<CSCCorrelatedLCTDigiCollection> mplctsH_;

  std::shared_ptr<CSCDigiMatcher> cscDigiMatcher_;
  std::shared_ptr<GEMDigiMatcher> gemDigiMatcher_;

  edm::ESGetToken<CSCGeometry, MuonGeometryRecord> geomToken_;
  const CSCGeometry* cscGeometry_;

  // all stubs (not necessarily matching) in crossed chambers with digis
  std::map<unsigned int, CSCCLCTDigiContainer> chamber_to_clcts_all_;
  std::map<unsigned int, CSCALCTDigiContainer> chamber_to_alcts_all_;
  std::map<unsigned int, CSCCorrelatedLCTDigiContainer> chamber_to_lcts_all_;
  std::map<unsigned int, CSCCorrelatedLCTDigiContainer> chamber_to_mplcts_all_;

  // all matching stubs in crossed chambers with digis
  std::map<unsigned int, CSCCLCTDigiContainer> chamber_to_clcts_;
  std::map<unsigned int, CSCALCTDigiContainer> chamber_to_alcts_;
  std::map<unsigned int, CSCCorrelatedLCTDigiContainer> chamber_to_lcts_;
  std::map<unsigned int, CSCCorrelatedLCTDigiContainer> chamber_to_mplcts_;

  template <class D>
  std::set<unsigned int> selectDetIds(D&, int) const;

  bool addGhostLCTs_;
  bool useGEMs_;

  bool matchTypeTightLCT_;

  int minNHitsChamber_;
  int minNHitsChamberALCT_;
  int minNHitsChamberCLCT_;
  int minNHitsChamberLCT_;
  int minNHitsChamberMPLCT_;

  bool verboseALCT_;
  bool verboseCLCT_;
  bool verboseLCT_;
  bool verboseMPLCT_;

  int minBXCLCT_, maxBXCLCT_;
  int minBXALCT_, maxBXALCT_;
  int minBXLCT_, maxBXLCT_;
  int minBXMPLCT_, maxBXMPLCT_;

  CSCCLCTDigiContainer no_clcts_;
  CSCALCTDigiContainer no_alcts_;
  CSCCorrelatedLCTDigiContainer no_lcts_;
  CSCCorrelatedLCTDigiContainer no_mplcts_;
};

template <class D>
std::set<unsigned int> CSCStubMatcher::selectDetIds(D& digis, int csc_type) const {
  std::set<unsigned int> result;
  for (auto& p : digis) {
    auto id = p.first;
    if (csc_type > 0) {
      CSCDetId detId(id);
      if (MuonHitHelper::toCSCType(detId.station(), detId.ring()) != csc_type)
        continue;
    }
    result.insert(p.first);
  }
  return result;
}

#endif