OMTFinput.cc

Go to the documentation of this file.

#include <cassert>
#include <iostream>
#include <cmath>
 
#include "L1Trigger/L1TMuonOverlap/interface/OMTFinput.h"
#include "L1Trigger/L1TMuonOverlap/interface/OMTFConfiguration.h"
#include "L1Trigger/L1TMuonOverlap/interface/XMLConfigReader.h"
 
OMTFinput::OMTFinput(const OMTFConfiguration *omtfConfig) {
  myOmtfConfig = omtfConfig;
  clear();
}
const OMTFinput::vector1D &OMTFinput::getLayerData(unsigned int iLayer, bool giveEta) const {
  assert(iLayer < measurementsPhi.size());
 
  if (giveEta)
    return measurementsEta[iLayer];
  return measurementsPhi[iLayer];
}
std::bitset<128> OMTFinput::getRefHits(unsigned int iProcessor) const {
  std::bitset<128> refHits;
 
  unsigned int iRefHit = 0;
  for (auto iRefHitDef : myOmtfConfig->getRefHitsDefs()[iProcessor]) {
    int iPhi = getLayerData(myOmtfConfig->getRefToLogicNumber()[iRefHitDef.iRefLayer])[iRefHitDef.iInput];
    int iEta = getLayerData(myOmtfConfig->getRefToLogicNumber()[iRefHitDef.iRefLayer], true)[iRefHitDef.iInput];
    if (iPhi < (int)myOmtfConfig->nPhiBins()) {
      refHits.set(iRefHit, iRefHitDef.fitsRange(iPhi));
      refHitsEta[iRefHit] = iEta;
    }
    iRefHit++;
  }
 
  return refHits;
}
bool OMTFinput::addLayerHit(unsigned int iLayer, unsigned int iInput, int iPhi, int iEta, bool allowOverwrite) {
  bool overwrite = false;
  assert(iLayer < myOmtfConfig->nLayers());
  assert(iInput < 14);
 
  if (iPhi >= (int)myOmtfConfig->nPhiBins())
    return true;
 
  if (allowOverwrite && measurementsPhi[iLayer][iInput] == iPhi && measurementsEta[iLayer][iInput] == iEta)
    return true;
 
  if (measurementsPhi[iLayer][iInput] != (int)myOmtfConfig->nPhiBins())
    ++iInput;
  if (measurementsPhi[iLayer][iInput] != (int)myOmtfConfig->nPhiBins())
    overwrite = true;
 
  if (iInput > 13)
    return true;
 
  measurementsPhi[iLayer][iInput] = iPhi;
  measurementsEta[iLayer][iInput] = iEta;
 
  return overwrite;
}
void OMTFinput::readData(XMLConfigReader *aReader, unsigned int iEvent, unsigned int iProcessor) {
  measurementsPhi = aReader->readEvent(iEvent, iProcessor);
  measurementsEta = aReader->readEvent(iEvent, iProcessor, true);
}
void OMTFinput::mergeData(const OMTFinput *aInput) {
  for (unsigned int iLayer = 0; iLayer < myOmtfConfig->nLayers(); ++iLayer) {
    const OMTFinput::vector1D &aPhiVec = aInput->getLayerData(iLayer, false);
    const OMTFinput::vector1D &aEtaVec = aInput->getLayerData(iLayer, true);
    if (aPhiVec.empty())
      continue;
 
    OMTFinput::vector1D layerData = getLayerData(iLayer, false);
    for (unsigned int iInput = 0; iInput < 14; ++iInput) {
      addLayerHit(iLayer, iInput, aPhiVec[iInput], aEtaVec[iInput]);
    }
  }
}
void OMTFinput::clear() {
  vector1D aLayer1D(14, myOmtfConfig->nPhiBins());
  measurementsPhi.assign(myOmtfConfig->nLayers(), aLayer1D);
  measurementsEta.assign(myOmtfConfig->nLayers(), aLayer1D);
  refHitsEta.assign(128, myOmtfConfig->nPhiBins());
}
void OMTFinput::shiftMyPhi(int phiShift) {
  int lowScaleEnd = std::pow(2, myOmtfConfig->nPhiBits() - 1);
  int highScaleEnd = lowScaleEnd - 1;
 
  for (unsigned int iLogicLayer = 0; iLogicLayer < measurementsPhi.size(); ++iLogicLayer) {
    for (unsigned int iHit = 0; iHit < measurementsPhi[iLogicLayer].size(); ++iHit) {
      if (!myOmtfConfig->getBendingLayers().count(iLogicLayer) &&
          measurementsPhi[iLogicLayer][iHit] < (int)myOmtfConfig->nPhiBins()) {
        if (measurementsPhi[iLogicLayer][iHit] < 0)
          measurementsPhi[iLogicLayer][iHit] += myOmtfConfig->nPhiBins();
        measurementsPhi[iLogicLayer][iHit] -= phiShift;
        if (measurementsPhi[iLogicLayer][iHit] < 0)
          measurementsPhi[iLogicLayer][iHit] += myOmtfConfig->nPhiBins();
        measurementsPhi[iLogicLayer][iHit] += -lowScaleEnd;
        if (measurementsPhi[iLogicLayer][iHit] < -lowScaleEnd || measurementsPhi[iLogicLayer][iHit] > highScaleEnd)
          measurementsPhi[iLogicLayer][iHit] = (int)myOmtfConfig->nPhiBins();
      }
    }
  }
}
std::ostream &operator<<(std::ostream &out, const OMTFinput &aInput) {
  for (unsigned int iLogicLayer = 0; iLogicLayer < aInput.measurementsPhi.size(); ++iLogicLayer) {
    out << "Logic layer: " << iLogicLayer << " Hits: ";
    for (unsigned int iHit = 0; iHit < aInput.measurementsPhi[iLogicLayer].size(); ++iHit) {
      out << aInput.measurementsPhi[iLogicLayer][iHit] << "\t";
    }
    out << std::endl;
  }
  return out;
}