SiPixelFedCablingMap.cc

Go to the documentation of this file.

#include "CondFormats/SiPixelObjects/interface/SiPixelFedCablingMap.h"
#include "CondFormats/SiPixelObjects/interface/SiPixelFedCablingTree.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"

#include <vector>
#include <iostream>
#include <algorithm>
#include <iostream>

using namespace sipixelobjects;

void SiPixelFedCablingMap::initializeRocs() {
  // OLD Method
  //for (auto & v : theMap) v.second.initFrameConversion();
  // below is the new code, works for phase0 and phase1

  // Decide if it is phase0 or phase1 based on the first fed, 0-phase0, 1200-phase1
  unsigned int fedId = (theMap.begin())->first.fed;  // get the first fed

  // Specifically for CMSSW_9_0_X, we need to call a different version of the frame
  // conversion steered by the version name in the cabling map
  if (theVersion.find("CMSSW_9_0_X") != std::string::npos) {
    for (auto &v : theMap)
      v.second.initFrameConversionPhase1_CMSSW_9_0_X();  // works
    std::cout << "*** Found CMSSW_9_0_X specific cabling map\n";
    return;
  }

  if (fedId >= FEDNumbering::MINSiPixeluTCAFEDID) {  // phase1 >= 1200
    for (auto &v : theMap)
      v.second.initFrameConversionPhase1();  // works
  } else {                                   // phase0
    for (auto &v : theMap)
      v.second.initFrameConversion();  // works
  }

  // if(0) {  // for testing
  //   for (Map::iterator im = theMap.begin(); im != theMap.end(); im++) {
  //     unsigned int fedId = im->first.fed;
  //     unsigned int linkId = im->first.link;
  //     unsigned int rocId = im->first.roc;
  //     auto rawDetID = im->second.rawId();
  //     auto idInDetUnit = im->second.idInDetUnit();
  //     auto idInLink = im->second.idInLink();
  //     //auto v = *im;
  //     if(fedId>=1200) {
  //    //v.second.initFrameConversionPhase1(); //
  //    im->second.initFrameConversionPhase1(); //
  //     } else {
  //    im->second.initFrameConversion();
  //     }
  //   } //
  // } // end if(0)
}

bool SiPixelFedCablingMap::Key::operator<(const Key &other) const {
  if (fed < other.fed)
    return true;
  if (fed > other.fed)
    return false;

  if (link < other.link)
    return true;
  if (link > other.link)
    return false;

  if (roc < other.roc)
    return true;
  if (roc > other.roc)
    return false;

  return false;
}

SiPixelFedCablingMap::SiPixelFedCablingMap(const SiPixelFedCablingTree *cab) : theVersion(cab->version()) {
  // Never called
  std::vector<const PixelFEDCabling *> fedList = cab->fedList();
  for (std::vector<const PixelFEDCabling *>::const_iterator ifed = fedList.begin(); ifed != fedList.end(); ifed++) {
    unsigned int fed = (**ifed).id();
    unsigned int numLink = (**ifed).numberOfLinks();
    for (unsigned int link = 1; link <= numLink; link++) {
      const PixelFEDLink *pLink = (**ifed).link(link);
      if (pLink == nullptr)
        continue;
      //unsigned int linkId = pLink->id();
      //if (linkId != 0 && linkId!= link)
      //  std::cout << "PROBLEM WITH LINK NUMBER!!!!" << std::endl;
      unsigned int numberROC = pLink->numberOfROCs();

      for (unsigned int roc = 1; roc <= numberROC; roc++) {
        const PixelROC *pROC = pLink->roc(roc);
        if (pROC == nullptr)
          continue;
        //if (pROC->idInLink() != roc)
        //  std::cout << "PROBLEM WITH ROC NUMBER!!!!" << std::endl;
        Key key = {fed, link, roc};
        theMap[key] = (*pROC);
      }
    }
  }  // fed loop
}

std::unique_ptr<SiPixelFedCablingTree> SiPixelFedCablingMap::cablingTree() const {
  std::unique_ptr<SiPixelFedCablingTree> tree(new SiPixelFedCablingTree(theVersion));
  for (Map::const_iterator im = theMap.begin(); im != theMap.end(); im++) {
    const sipixelobjects::PixelROC &roc = im->second;
    unsigned int fedId = im->first.fed;
    unsigned int linkId = im->first.link;
    tree->addItem(fedId, linkId, roc);
  }
  return tree;
}

std::vector<unsigned int> SiPixelFedCablingMap::fedIds() const {
  std::vector<unsigned int> result;
  for (Map::const_iterator im = theMap.begin(); im != theMap.end(); im++) {
    unsigned int fedId = im->first.fed;
    if (find(result.begin(), result.end(), fedId) == result.end())
      result.push_back(fedId);
  }
  return result;
}

const sipixelobjects::PixelROC *SiPixelFedCablingMap::findItem(const sipixelobjects::CablingPathToDetUnit &path) const {
  const PixelROC *roc = nullptr;
  Key key = {path.fed, path.link, path.roc};
  Map::const_iterator inMap = theMap.find(key);
  if (inMap != theMap.end())
    roc = &(inMap->second);
  return roc;
}

std::unordered_map<uint32_t, unsigned int> SiPixelFedCablingMap::det2fedMap() const {
  std::unordered_map<uint32_t, unsigned int> result;
  for (auto im = theMap.begin(); im != theMap.end(); ++im) {
    result[im->second.rawId()] = im->first.fed;  // we know: a det is in only one fed!
  }
  return result;
}

std::map<uint32_t, std::vector<sipixelobjects::CablingPathToDetUnit> > SiPixelFedCablingMap::det2PathMap() const {
  std::map<uint32_t, std::vector<sipixelobjects::CablingPathToDetUnit> > result;
  for (auto im = theMap.begin(); im != theMap.end(); ++im) {
    CablingPathToDetUnit path = {im->first.fed, im->first.link, im->first.roc};
    result[im->second.rawId()].push_back(path);
  }
  return result;
}

std::vector<sipixelobjects::CablingPathToDetUnit> SiPixelFedCablingMap::pathToDetUnit(uint32_t rawDetId) const {
  std::vector<sipixelobjects::CablingPathToDetUnit> result;
  for (auto im = theMap.begin(); im != theMap.end(); ++im) {
    if (im->second.rawId() == rawDetId) {
      CablingPathToDetUnit path = {im->first.fed, im->first.link, im->first.roc};
      result.push_back(path);
    }
  }
  return result;
}

bool SiPixelFedCablingMap::pathToDetUnitHasDetUnit(uint32_t rawDetId, unsigned int fedId) const {
  auto end = theMap.end();
  for (auto im = theMap.lower_bound({fedId, 0, 0}); im != end and im->first.fed == fedId; ++im) {
    if (im->second.rawId() == rawDetId) {
      return true;
    }
  }
  return false;
}