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 overwrite = false;
  assert(iLayer<myOmtfConfig->nLayers());
  assert(iInput<14);

  if(iPhi>=(int)myOmtfConfig->nPhiBins()) return true;
//  unsigned int origInput = iInput;

  if(measurementsPhi[iLayer][iInput]!=(int)myOmtfConfig->nPhiBins()) ++iInput;
  if(measurementsPhi[iLayer][iInput]!=(int)myOmtfConfig->nPhiBins()) overwrite = true;

//  std::cout <<" addLayerHit :"<<iLayer<<" input orig :"<<inputorig<<" now:  "<< iInput <<" phi: "<<iPhi<< std::endl;
  
  if(iInput>13) return false;
  
  measurementsPhi[iLayer][iInput] = iPhi;
  measurementsEta[iLayer][iInput] = iEta;

//  return (origInput==iInput);                   
  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.size()) 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;
}