OmtfLinkMappingRpc.cc

Go to the documentation of this file.

#include "EventFilter/L1TRawToDigi/interface/OmtfLinkMappingRpc.h"
 
#include <fstream>
 
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESTransientHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
#include "CondFormats/RPCObjects/interface/RPCReadOutMapping.h"
#include "CondFormats/RPCObjects/interface/RPCEMap.h"
#include "CondFormats/DataRecord/interface/RPCEMapRcd.h"
#include "CondFormats/RPCObjects/interface/DccSpec.h"
#include "CondFormats/RPCObjects/interface/LinkBoardElectronicIndex.h"
#include "CondFormats/RPCObjects/interface/LinkBoardPackedStrip.h"
#include "CondFormats/RPCObjects/interface/TriggerBoardSpec.h"
#include "CondFormats/RPCObjects/interface/LinkBoardSpec.h"
 
#include "CondFormats/DataRecord/interface/RPCOMTFLinkMapRcd.h"
#include "CondFormats/RPCObjects/interface/RPCAMCLinkMap.h"
 
namespace omtf {
 
  bool lessLinkBoardElectronicIndex::operator()(const LinkBoardElectronicIndex& o1,
                                                const LinkBoardElectronicIndex& o2) const {
    if (o1.dccId < o2.dccId)
      return true;
    else if (o1.dccId == o2.dccId && o1.dccInputChannelNum < o2.dccInputChannelNum)
      return true;
    else if (o1.dccId == o2.dccId && o1.dccInputChannelNum == o2.dccInputChannelNum &&
             o1.tbLinkInputNum < o2.tbLinkInputNum)
      return true;
    else
      return false;
  }
 
  MapEleIndex2LBIndex translateOmtf2Pact(const RpcLinkMap& omtfLink2Ele, const RPCReadOutMapping* pactCabling) {
    MapEleIndex2LBIndex omtf2rpc;
 
    std::vector<const DccSpec*> dccs = pactCabling->dccList();
    for (auto it1 : dccs) {
      const std::vector<TriggerBoardSpec>& rmbs = it1->triggerBoards();
      for (auto const& it2 : rmbs) {
        const std::vector<LinkConnSpec>& links = it2.linkConns();
        for (auto const& it3 : links) {
          const std::vector<LinkBoardSpec>& lbs = it3.linkBoards();
          for (std::vector<LinkBoardSpec>::const_iterator it4 = lbs.begin(); it4 != lbs.end(); ++it4) {
            std::string lbNameCH = it4->linkBoardName();
            std::string lbName = lbNameCH.substr(0, lbNameCH.size() - 4);
            std::vector<EleIndex> omtfEles = omtfLink2Ele.omtfEleIndex(lbName);
            //        if (!omtfEles.empty()) std::cout <<"  isOK ! " <<  it4->linkBoardName() <<" has: " << omtfEles.size() << " first: "<< omtfEles[0] << std::endl;
            LinkBoardElectronicIndex rpcEle = {
                it1->id(), it2.dccInputChannelNum(), it3.triggerBoardInputNumber(), it4->linkBoardNumInLink()};
            for (const auto& omtfEle : omtfEles)
              omtf2rpc[omtfEle] = rpcEle;
          }
        }
      }
    }
    LogTrace(" ") << " SIZE OF OMTF to RPC map TRANSLATION is: " << omtf2rpc.size() << std::endl;
    return omtf2rpc;
  }
 
  MapLBIndex2EleIndex translatePact2Omtf(const RpcLinkMap& omtfLink2Ele, const RPCReadOutMapping* pactCabling) {
    MapLBIndex2EleIndex pact2omtfs;
    MapEleIndex2LBIndex omtf2rpcs = translateOmtf2Pact(omtfLink2Ele, pactCabling);
    for (const auto& omtf2rpc : omtf2rpcs) {
      std::pair<EleIndex, EleIndex>& omtfs = pact2omtfs[omtf2rpc.second];
      if (omtfs.first.fed() == 0)
        omtfs.first = omtf2rpc.first;
      else if (omtfs.second.fed() == 0)
        omtfs.second = omtf2rpc.first;
      else
        edm::LogError(" translatePact2Omtf ") << " PROBLEM LinkBoardElectronicIndex already USED!!!! ";
    }
    return pact2omtfs;
  }
 
  void RpcLinkMap::init(const edm::EventSetup& es) {
    edm::ESHandle<RPCAMCLinkMap> amcMapping;
    es.get<RPCOMTFLinkMapRcd>().get(amcMapping);
    const RPCAMCLinkMap::map_type& amcMap = amcMapping->getMap();
 
    for (const auto& item : amcMap) {
      unsigned int fedId = item.first.getFED();
      unsigned int amcSlot = item.first.getAMCNumber();
      unsigned int link = item.first.getAMCInput();
      std::string lbName = item.second.getName();
 
      std::string processorNameStr = "OMTF";
      ;
      if (fedId == 1380)
        processorNameStr = "OMTFn";
      else
        processorNameStr = "OMTFp";
      processorNameStr += std::to_string(amcSlot / 2 + 1);
      std::string processorName = processorNameStr;
 
      std::map<unsigned int, std::string>& li2lb = link2lbName[processorName];
      std::map<std::string, unsigned int>& lb2li = lbName2link[processorName];
      li2lb[link] = lbName;
      lb2li[lbName] = link;
      EleIndex ele(processorName, link);
      lbName2OmtfIndex[lbName].push_back(ele);
    }
  }
 
  void RpcLinkMap::init(const std::string& fName) {
    std::ifstream inFile;
    inFile.open(fName);
    if (inFile) {
      LogTrace("") << " reading OmtfRpcLinksMap from: " << fName;
    } else {
      LogTrace("") << " Unable to open file " << fName;
 
      throw std::runtime_error("Unable to open OmtfRpcLinksMap file " + fName);
    }
 
    std::string line;
    while (std::getline(inFile, line)) {
      line.erase(0, line.find_first_not_of(" \t\r\n"));  //cut first character
      if (line.empty() || !line.compare(0, 2, "--"))
        continue;  // empty or comment line
      std::stringstream ss(line);
      std::string processorName, lbName;
      unsigned int link, dbId;
      if (ss >> processorName >> link >> lbName >> dbId) {
        std::map<unsigned int, std::string>& li2lb = link2lbName[processorName];
        std::map<std::string, unsigned int>& lb2li = lbName2link[processorName];
        li2lb[link] = lbName;
        lb2li[lbName] = link;
        EleIndex ele(processorName, link);
        lbName2OmtfIndex[lbName].push_back(ele);
      }
    }
    inFile.close();
  }
 
  std::vector<EleIndex> RpcLinkMap::omtfEleIndex(const std::string& lbName) const {
    const auto pos = lbName2OmtfIndex.find(lbName);
    return pos != lbName2OmtfIndex.end() ? pos->second : std::vector<EleIndex>();
  }
 
}  // namespace omtf