GEMDigiMatcher.h

Go to the documentation of this file.

#ifndef Validation_MuonGEMDigis_GEMDigiMatcher_h
#define Validation_MuonGEMDigis_GEMDigiMatcher_h

#include "Validation/MuonHits/interface/GEMSimHitMatcher.h"
#include "Geometry/GEMGeometry/interface/GEMGeometry.h"
#include "SimDataFormats/GEMDigiSimLink/interface/GEMDigiSimLink.h"
#include "DataFormats/GEMDigi/interface/GEMDigiCollection.h"
#include "DataFormats/GEMDigi/interface/GEMPadDigiCollection.h"
#include "DataFormats/GEMDigi/interface/GEMPadDigiClusterCollection.h"
#include "DataFormats/GEMDigi/interface/GEMCoPadDigiCollection.h"
#include "DataFormats/Common/interface/DetSetVector.h"

#include <vector>
#include <map>
#include <set>
#include <unordered_set>

typedef std::vector<GEMDigiSimLink> GEMDigiSimLinkContainer;
typedef std::vector<GEMDigi> GEMDigiContainer;
typedef std::vector<GEMPadDigi> GEMPadDigiContainer;
typedef std::vector<GEMPadDigiCluster> GEMPadDigiClusterContainer;
typedef std::vector<GEMCoPadDigi> GEMCoPadDigiContainer;

class GEMDigiMatcher {
public:
  // constructor
  GEMDigiMatcher(edm::ParameterSet const& iPS, edm::ConsumesCollector&& iC);

  // destructor
  ~GEMDigiMatcher() {}

  // initialize the event
  void init(const edm::Event& e, const edm::EventSetup& eventSetup);

  // do the matching
  void match(const SimTrack& t, const SimVertex& v);

  std::set<unsigned int> detIdsSimLink(int gem_type = MuonHitHelper::GEM_ALL) const;

  // partition GEM detIds with digis
  std::set<unsigned int> detIdsDigi(int gem_type = MuonHitHelper::GEM_ALL) const;
  std::set<unsigned int> detIdsPad(int gem_type = MuonHitHelper::GEM_ALL) const;
  std::set<unsigned int> detIdsCluster(int gem_type = MuonHitHelper::GEM_ALL) const;

  // chamber detIds with digis
  std::set<unsigned int> chamberIdsDigi(int gem_type = MuonHitHelper::GEM_ALL) const;
  std::set<unsigned int> chamberIdsPad(int gem_type = MuonHitHelper::GEM_ALL) const;
  std::set<unsigned int> chamberIdsCluster(int gem_type = MuonHitHelper::GEM_ALL) const;

  // superchamber detIds with digis
  std::set<unsigned int> superChamberIdsDigi(int gem_type = MuonHitHelper::GEM_ALL) const;
  std::set<unsigned int> superChamberIdsPad(int gem_type = MuonHitHelper::GEM_ALL) const;
  std::set<unsigned int> superChamberIdsCluster(int gem_type = MuonHitHelper::GEM_ALL) const;
  std::set<unsigned int> superChamberIdsCoPad(int gem_type = MuonHitHelper::GEM_ALL) const;

  const GEMDigiSimLinkContainer& simLinksInDetId(unsigned int) const;

  // GEM digis from a particular partition, chamber or superchamber
  const GEMDigiContainer& digisInDetId(unsigned int) const;
  const GEMDigiContainer& digisInChamber(unsigned int) const;
  const GEMDigiContainer& digisInSuperChamber(unsigned int) const;

  // GEM pads from a particular partition, chamber or superchamber
  const GEMPadDigiContainer& padsInDetId(unsigned int) const;
  const GEMPadDigiContainer& padsInChamber(unsigned int) const;
  const GEMPadDigiContainer& padsInSuperChamber(unsigned int) const;

  // GEM clusters from a particular partition, chamber or superchamber
  const GEMPadDigiClusterContainer& clustersInDetId(unsigned int) const;
  const GEMPadDigiClusterContainer& clustersInChamber(unsigned int) const;
  const GEMPadDigiClusterContainer& clustersInSuperChamber(unsigned int) const;

  // GEM co-pads from a particular partition or superchamber
  const GEMCoPadDigiContainer& coPadsInSuperChamber(unsigned int) const;

  const std::map<unsigned int, GEMPadDigiContainer> allPads() const { return detid_to_pads_; }

  // #layers with digis from this simtrack
  int nLayersWithDigisInSuperChamber(unsigned int) const;
  int nLayersWithPadsInSuperChamber(unsigned int) const;
  int nLayersWithClustersInSuperChamber(unsigned int) const;

  int nPads() const;

  int nCoPads() const;

  std::set<int> stripNumbersInDetId(unsigned int) const;
  std::set<int> padNumbersInDetId(unsigned int) const;

  int extrapolateHsfromGEMPad(unsigned int, int) const;
  int extrapolateHsfromGEMStrip(unsigned int, int) const;

  // what unique partitions numbers with digis from this simtrack?
  std::set<int> partitionNumbers() const;
  std::set<int> partitionNumbersWithCoPads() const;

  GlobalPoint getGlobalPointDigi(unsigned int rawId, const GEMDigi& d) const;
  GlobalPoint getGlobalPointPad(unsigned int rawId, const GEMPadDigi& tp) const;

  // get the GEMSimHitMatcher
  std::shared_ptr<GEMSimHitMatcher> muonSimHitMatcher() { return muonSimHitMatcher_; }

private:
  // match digi simlink to a simtrack (based on track Id and particle type)
  void matchDigisSLToSimTrack(const edm::DetSetVector<GEMDigiSimLink>&);
  // match digi to a simtrack (based on BX and strip number)
  void matchDigisToSimTrack(const GEMDigiCollection&);
  // match pad to a simtrack (based on BX and pad number)
  void matchPadsToSimTrack(const GEMPadDigiCollection&);
  // match pad to a simtrack (based on BX and pad number)
  void matchPadsInDetId(const GEMDetId& simhitId, const GEMPadDigiCollection&);
  // match cluster to a simtrack (based on BX and pad number)
  void matchClustersToSimTrack(const GEMPadDigiClusterCollection&);
  // match copad to a simtrack (based on previously matching pads)
  void matchCoPadsToSimTrack(const GEMCoPadDigiCollection&);

  void clear();

  edm::EDGetTokenT<edm::DetSetVector<GEMDigiSimLink>> gemSimLinkToken_;
  edm::EDGetTokenT<GEMDigiCollection> gemDigiToken_;
  edm::EDGetTokenT<GEMPadDigiCollection> gemPadToken_;
  edm::EDGetTokenT<GEMPadDigiClusterCollection> gemClusterToken_;
  edm::EDGetTokenT<GEMCoPadDigiCollection> gemCoPadToken_;

  edm::Handle<edm::DetSetVector<GEMDigiSimLink>> gemDigisSLH_;
  edm::Handle<GEMDigiCollection> gemDigisH_;
  edm::Handle<GEMPadDigiCollection> gemPadsH_;
  edm::Handle<GEMPadDigiClusterCollection> gemClustersH_;
  edm::Handle<GEMCoPadDigiCollection> gemCoPadsH_;

  std::shared_ptr<GEMSimHitMatcher> muonSimHitMatcher_;

  edm::ESGetToken<GEMGeometry, MuonGeometryRecord> geomToken_;
  const GEMGeometry* gemGeometry_;

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

  bool simMuOnly_;
  bool discardEleHits_;

  int minBXDigi_, maxBXDigi_;
  int minBXPad_, maxBXPad_;
  int minBXCluster_, maxBXCluster_;
  int minBXCoPad_, maxBXCoPad_;

  int matchDeltaStrip_;

  bool verboseSimLink_;
  bool verboseDigi_;
  bool verbosePad_;
  bool verboseCluster_;
  bool verboseCoPad_;

  bool matchToSimLink_;

  // map of raw GEMDetId to containers
  std::map<unsigned int, GEMDigiSimLinkContainer> detid_to_simLinks_;

  std::map<unsigned int, GEMDigiContainer> detid_to_digis_;
  std::map<unsigned int, GEMDigiContainer> chamber_to_digis_;
  std::map<unsigned int, GEMDigiContainer> superchamber_to_digis_;

  std::map<unsigned int, GEMPadDigiContainer> detid_to_pads_;
  std::map<unsigned int, GEMPadDigiContainer> chamber_to_pads_;
  std::map<unsigned int, GEMPadDigiContainer> superchamber_to_pads_;

  std::map<unsigned int, GEMPadDigiClusterContainer> detid_to_clusters_;
  std::map<unsigned int, GEMPadDigiClusterContainer> chamber_to_clusters_;
  std::map<unsigned int, GEMPadDigiClusterContainer> superchamber_to_clusters_;

  std::map<unsigned int, GEMCoPadDigiContainer> superchamber_to_copads_;

  GEMDigiSimLinkContainer no_gem_simLinks_;
  GEMDigiContainer no_gem_digis_;
  GEMPadDigiContainer no_gem_pads_;
  GEMPadDigiClusterContainer no_gem_clusters_;
  GEMCoPadDigiContainer no_gem_copads_;
};

template <class T>
std::set<unsigned int> GEMDigiMatcher::selectDetIds(const T& digis, int gem_type) const {
  std::set<unsigned int> result;
  for (const auto& p : digis) {
    const auto& id = p.first;
    if (gem_type > 0) {
      GEMDetId detId(id);
      if (MuonHitHelper::toGEMType(detId.station(), detId.ring()) != gem_type)
        continue;
    }
    result.insert(p.first);
  }
  return result;
}

#endif