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/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::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_;
 
  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