OMTFConfigMaker.cc

Go to the documentation of this file.

#include <iostream>

#include "L1Trigger/L1TMuonOverlap/interface/OMTFConfigMaker.h"
#include "L1Trigger/L1TMuonOverlap/interface/GoldenPattern.h"
#include "L1Trigger/L1TMuonOverlap/interface/XMLConfigReader.h"
#include "L1Trigger/L1TMuonOverlap/interface/OMTFinput.h"
#include "L1Trigger/L1TMuonOverlap/interface/OMTFResult.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

OMTFConfigMaker::OMTFConfigMaker(OMTFConfiguration * omtfConfig) : myOmtfConfig(omtfConfig) { 

  std::vector<int> refPhi1D(myOmtfConfig->nLogicRegions(),2*myOmtfConfig->nPhiBins());
  minRefPhi2D.assign(myOmtfConfig->nRefLayers(),refPhi1D);

  refPhi1D = std::vector<int>(myOmtfConfig->nLogicRegions(),-2*myOmtfConfig->nPhiBins());
  maxRefPhi2D.assign(myOmtfConfig->nRefLayers(),refPhi1D);
}
OMTFConfigMaker::~OMTFConfigMaker(){ }
void OMTFConfigMaker::fillCounts(unsigned int iProcessor,
                 const OMTFinput & aInput){

}
void OMTFConfigMaker::fillPhiMaps(unsigned int iProcessor,
                  const OMTFinput & aInput){

  for(unsigned int iRefLayer=0;iRefLayer<myOmtfConfig->nRefLayers();++iRefLayer){
    const OMTFinput::vector1D & refLayerHits = aInput.getLayerData(myOmtfConfig->getRefToLogicNumber()[iRefLayer]); 
    if(refLayerHits.empty()) continue;
    for(unsigned int iInput=0;iInput<myOmtfConfig->nInputs();++iInput){ 
      int phiRef = refLayerHits[iInput];
      unsigned int iRegion = myOmtfConfig->getRegionNumberFromMap(iInput,iRefLayer,phiRef);       
      if(phiRef>=(int)myOmtfConfig->nPhiBins()) continue;      
      if(phiRef<minRefPhi2D[iRefLayer][iRegion]) minRefPhi2D[iRefLayer][iRegion] = phiRef;      
      if(phiRef>maxRefPhi2D[iRefLayer][iRegion]) maxRefPhi2D[iRefLayer][iRegion] = phiRef;      
    }
  }
}
void OMTFConfigMaker::printPhiMap(std::ostream & out){
  
  out<<"min Phi in each logicRegion (X) in each ref Layer (Y): "<<std::endl; 
  for(unsigned int iRefLayer=0;iRefLayer<myOmtfConfig->nRefLayers();++iRefLayer){
    for(unsigned int iLogicRegion=0;iLogicRegion<myOmtfConfig->nLogicRegions();++iLogicRegion){
      out<<"          "<<minRefPhi2D[iRefLayer][iLogicRegion]<<"\t";
    }
    out<<std::endl;
  }
  out<<std::endl; 

  out<<"max Phi in each logicRegion (X) in each ref Layer (Y): "<<std::endl;
  for(unsigned int iRefLayer=0;iRefLayer<myOmtfConfig->nRefLayers();++iRefLayer){
    for(unsigned int iLogicRegion=0;iLogicRegion<myOmtfConfig->nLogicRegions();++iLogicRegion){
      out<<"          "<<maxRefPhi2D[iRefLayer][iLogicRegion]<<"\t";
    }
    out<<std::endl;
  }
  out<<std::endl; 

}
void OMTFConfigMaker::makeConnetionsMap(unsigned int iProcessor,
                    const OMTFinput & aInput){

  fillPhiMaps(iProcessor,aInput);
  
  for(unsigned int iRefLayer=0;iRefLayer<myOmtfConfig->nRefLayers();++iRefLayer){
    const OMTFinput::vector1D & refLayerHits = aInput.getLayerData(myOmtfConfig->getRefToLogicNumber()[iRefLayer]);    
    if(refLayerHits.empty()) continue;
    for(unsigned int iInput=0;iInput<refLayerHits.size();++iInput){
      int phiRef = refLayerHits[iInput];
      unsigned int iRegion = myOmtfConfig->getRegionNumberFromMap(iInput,iRefLayer,phiRef);     
      if(iRegion>=myOmtfConfig->nLogicRegions()) continue;      
      fillInputRange(iProcessor,iRegion,aInput);
      fillInputRange(iProcessor,iRegion,iRefLayer,iInput);
      if(iInput%2==0) fillInputRange(iProcessor,iRegion,iRefLayer,iInput+1);
    }      
  }
} 
void OMTFConfigMaker::fillInputRange(unsigned int iProcessor,
                   unsigned int iRegion,
                   const OMTFinput & aInput){

  for(unsigned int iLogicLayer=0;iLogicLayer<myOmtfConfig->nLayers();++iLogicLayer){
    for(unsigned int iInput=0;iInput<14;++iInput){
      bool isHit = aInput.getLayerData(iLogicLayer)[iInput]<(int)myOmtfConfig->nPhiBins();
      myOmtfConfig->getMeasurements4D()[iProcessor][iRegion][iLogicLayer][iInput]+=isHit;
    }
  }
}
void OMTFConfigMaker::fillInputRange(unsigned int iProcessor,
                   unsigned int iRegion,
                   unsigned int iRefLayer,
                   unsigned int iInput){

  ++myOmtfConfig->getMeasurements4Dref()[iProcessor][iRegion][iRefLayer][iInput]; 

}
void OMTFConfigMaker::printConnections(std::ostream & out,
                       unsigned int iProcessor,
                       unsigned int iRegion){

  out<<"iProcessor: "<<iProcessor
     <<" iRegion: "<<iRegion
     <<std::endl;

  out<<"Ref hits"<<std::endl;
  for(unsigned int iLogicLayer=0;iLogicLayer<myOmtfConfig->nLayers();++iLogicLayer){
    out<<"Logic layer: "<<iLogicLayer<<" Hits: ";
    for(unsigned int iInput=0;iInput<myOmtfConfig->nInputs();++iInput){
      out<<myOmtfConfig->getMeasurements4Dref()[iProcessor][iRegion][iLogicLayer][iInput]<<"\t";
    }
    out<<std::endl;
  }
  /*
  out<<"Measurement hits"<<std::endl;
  for(unsigned int iLogicLayer=0;iLogicLayer<myOmtfConfig->nLayers();++iLogicLayer){
    out<<"Logic layer: "<<iLogicLayer<<" Hits: ";
    for(unsigned int iInput=0;iInput<14;++iInput){
      out<<myOmtfConfig->getMeasurements4D()[iProcessor][iRegion][iLogicLayer][iInput]<<"\t";
    }
    out<<std::endl;
  }
  */
  out<<std::endl;
}