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();
}