TrackerDetToDTCELinkCablingMap.cc

Go to the documentation of this file.

#include "CondFormats/SiPhase2TrackerObjects/interface/TrackerDetToDTCELinkCablingMap.h"
#include "FWCore/Utilities/interface/Exception.h"

#include <utility>
#include <algorithm>
#include <iostream>

TrackerDetToDTCELinkCablingMap::TrackerDetToDTCELinkCablingMap() {}

TrackerDetToDTCELinkCablingMap::~TrackerDetToDTCELinkCablingMap() {}

std::unordered_map<DTCELinkId, uint32_t>::const_iterator TrackerDetToDTCELinkCablingMap::dtcELinkIdToDetId(
    DTCELinkId const& key) const {
  if (cablingMapDTCELinkIdToDetId_.find(key) == cablingMapDTCELinkIdToDetId_.end()) {
    throw cms::Exception(
        "TrackerDetToDTCELinkCablingMap has been asked to return a DetId associated to a DTCELinkId, but the latter is "
        "unknown to the map. ")
        << " (DTC, GBT, Elink) numbers = (" << key.dtc_id() << "," << key.gbtlink_id() << "," << key.elink_id() << ")"
        << std::endl;
  }

  return cablingMapDTCELinkIdToDetId_.find(key);
}

std::pair<std::unordered_multimap<uint32_t, DTCELinkId>::const_iterator,
          std::unordered_multimap<uint32_t, DTCELinkId>::const_iterator>
TrackerDetToDTCELinkCablingMap::detIdToDTCELinkId(uint32_t const key) const {
  auto const DTCELinkId_itpair = cablingMapDetIdToDTCELinkId_.equal_range(key);

  if (DTCELinkId_itpair.first == cablingMapDetIdToDTCELinkId_.end()) {
    throw cms::Exception(
        "TrackerDetToDTCELinkCablingMap has been asked to return a DTCELinkId associated to a DetId, but the latter is "
        "unknown to the map. ")
        << " DetId = " << key << std::endl;
  }

  return DTCELinkId_itpair;
}

bool TrackerDetToDTCELinkCablingMap::knowsDTCELinkId(DTCELinkId const& key) const {
  return cablingMapDTCELinkIdToDetId_.find(key) != cablingMapDTCELinkIdToDetId_.end();
}

bool TrackerDetToDTCELinkCablingMap::knowsDetId(uint32_t key) const {
  return cablingMapDetIdToDTCELinkId_.find(key) != cablingMapDetIdToDTCELinkId_.end();
}

std::vector<DTCELinkId> TrackerDetToDTCELinkCablingMap::getKnownDTCELinkIds() const {
  std::vector<DTCELinkId> knownDTCELinkIds(cablingMapDTCELinkIdToDetId_.size());

  // Unzip the map into a vector of DTCELinkId, discarding the DetIds
  std::transform(cablingMapDTCELinkIdToDetId_.begin(),
                 cablingMapDTCELinkIdToDetId_.end(),
                 knownDTCELinkIds.begin(),
                 [=](auto pair) { return pair.first; });

  return knownDTCELinkIds;
}

std::vector<uint32_t> TrackerDetToDTCELinkCablingMap::getKnownDetIds() const {
  std::vector<uint32_t> knownDetId;

  // To get the list of unique DetIds we need to iterate over the various equal_ranges
  // in the map associated to each unique key, and count them only once.

  for (auto allpairs_it = cablingMapDetIdToDTCELinkId_.begin(), allpairs_end = cablingMapDetIdToDTCELinkId_.end();
       allpairs_it != allpairs_end;) {
    // ***Store the first instance of the key***
    knownDetId.push_back(uint32_t(allpairs_it->first));

    // *** Skip to the end of the equal range ***
    // The following is just explicative, the bottom expression is equivalent
    //auto const current_key             = allpairs_it->first;
    //auto const current_key_equal_range = cablingMapDetIdToDTCELinkId_.equal_range(current_key);
    //auto const current_key_range_end   = current_key_equal_range.second;
    auto const current_key_range_end = cablingMapDetIdToDTCELinkId_.equal_range(allpairs_it->first).second;

    while (allpairs_it != current_key_range_end)
      ++allpairs_it;
  }

  return knownDetId;
}

void TrackerDetToDTCELinkCablingMap::insert(DTCELinkId const& dtcELinkId, uint32_t const detId) {
  cablingMapDTCELinkIdToDetId_.insert(std::make_pair(DTCELinkId(dtcELinkId), uint32_t(detId)));
  cablingMapDetIdToDTCELinkId_.insert(std::make_pair(uint32_t(detId), DTCELinkId(dtcELinkId)));
}

void TrackerDetToDTCELinkCablingMap::clear() {
  cablingMapDTCELinkIdToDetId_.clear();
  cablingMapDetIdToDTCELinkId_.clear();
}