HGCSiliconDetIdToROC.cc

Go to the documentation of this file.

#include "DataFormats/ForwardDetId/interface/HGCSiliconDetIdToROC.h"

HGCSiliconDetIdToROC::HGCSiliconDetIdToROC() {
  const int rocMax(6), cellRows(4), cellColumn(2);
  const int type[rocMax] = {0, 0, 1, 1, 2, 2};
  const int start1[rocMax][cellColumn] = {{3, 2}, {1, 0}, {1, 2}, {3, 4}, {7, 6}, {5, 4}};
  const int start2[rocMax] = {0, 0, 0, 1, 4, 4};

  // key = ROC #; value = vector of trigger cell U.V
  for (int roc = 0; roc < rocMax; ++roc) {
    std::vector<std::pair<int, int> > cells;
    if (type[roc] == 1) {
      for (int i1 = 0; i1 < cellColumn; ++i1) {
        int u = start1[roc][i1];
        int v = start2[roc] + u / 2;
        for (int i2 = 0; i2 < cellRows; ++i2) {
          cells.emplace_back(u, v);
          ++u;
        }
      }
    } else if (type[roc] == 2) {
      for (int i1 = 0; i1 < cellColumn; ++i1) {
        int v = start2[roc];
        int u = start1[roc][i1];
        for (int i2 = 0; i2 < cellRows; ++i2) {
          cells.emplace_back(u, v);
          ++v;
        }
      }
    } else {
      for (int i1 = 0; i1 < cellColumn; ++i1) {
        int u = start2[roc];
        int v = start1[roc][i1];
        for (int i2 = 0; i2 < cellRows; ++i2) {
          cells.emplace_back(u, v);
          ++u;
          ++v;
        }
      }
    }
    triggerIdFromROC_[roc + 1] = cells;
  }

  // key = trigger cell U,V; value = roc #
  for (auto const& itr : triggerIdFromROC_) {
    for (auto cell : itr.second)
      triggerIdToROC_[cell] = itr.first;
  }
}

int HGCSiliconDetIdToROC::getROCNumber(int triggerCellU, int triggerCellV, int type) const {
  auto itr = triggerIdToROC_.find(std::make_pair(triggerCellU, triggerCellV));
  int roc = ((itr == triggerIdToROC_.end()) ? -1 : ((type == 0) ? itr->second : (1 + itr->second) / 2));
  return roc;
}

std::vector<std::pair<int, int> > HGCSiliconDetIdToROC::getTriggerId(int roc, int type) const {
  std::vector<std::pair<int, int> > list;
  if (type == 0) {
    auto itr = triggerIdFromROC_.find(roc);
    return ((itr == triggerIdFromROC_.end()) ? list : itr->second);
  } else {
    for (int k = 0; k < 2; ++k) {
      int rocx = 2 * roc + k - 1;
      auto itr = triggerIdFromROC_.find(rocx);
      if (itr != triggerIdFromROC_.end()) {
        for (auto cell : itr->second)
          list.emplace_back(cell);
      }
    }
    return list;
  }
}

void HGCSiliconDetIdToROC::print() const {
  for (auto const& itr : triggerIdFromROC_) {
    std::cout << "ROC " << itr.first << " with " << (itr.second).size() << " trigger cells:";
    for (auto cell : itr.second)
      std::cout << " (" << cell.first << "," << cell.second << ")";
    std::cout << std::endl;
  }
  for (auto const& itr : triggerIdToROC_)
    std::cout << "Trigger cell (" << (itr.first).first << "," << (itr.first).second << ") : ROC " << itr.second
              << std::endl;
}