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(const edm::ParameterSet & theConfig, OMTFConfiguration * omtf_config) : m_omtf_config(omtf_config) { 

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

  refPhi1D = std::vector<int>(OMTFConfiguration::instance()->nLogicRegions,-2*OMTFConfiguration::instance()->nPhiBins);
  maxRefPhi2D.assign(OMTFConfiguration::instance()->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<OMTFConfiguration::instance()->nRefLayers;++iRefLayer){
    const OMTFinput::vector1D & refLayerHits = aInput.getLayerData(OMTFConfiguration::instance()->refToLogicNumber[iRefLayer]); 
    if(!refLayerHits.size()) continue;
    for(unsigned int iInput=0;iInput<OMTFConfiguration::instance()->nInputs;++iInput){  
      int phiRef = refLayerHits[iInput];
      unsigned int iRegion = OMTFConfiguration::instance()->getRegionNumberFromMap(iInput,iRefLayer,phiRef);       
      if(phiRef>=(int)OMTFConfiguration::instance()->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<OMTFConfiguration::instance()->nRefLayers;++iRefLayer){
    for(unsigned int iLogicRegion=0;iLogicRegion<OMTFConfiguration::instance()->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<OMTFConfiguration::instance()->nRefLayers;++iRefLayer){
    for(unsigned int iLogicRegion=0;iLogicRegion<OMTFConfiguration::instance()->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<OMTFConfiguration::instance()->nRefLayers;++iRefLayer){
    const OMTFinput::vector1D & refLayerHits = aInput.getLayerData(OMTFConfiguration::instance()->refToLogicNumber[iRefLayer]);    
    if(!refLayerHits.size()) continue;
    for(unsigned int iInput=0;iInput<refLayerHits.size();++iInput){
      int phiRef = refLayerHits[iInput];
      unsigned int iRegion = OMTFConfiguration::instance()->getRegionNumberFromMap(iInput,iRefLayer,phiRef);     
      if(iRegion>=OMTFConfiguration::instance()->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<OMTFConfiguration::instance()->nLayers;++iLogicLayer){

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

  ++m_omtf_config->measurements4Dref[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<OMTFConfiguration::instance()->nLayers;++iLogicLayer){
    out<<"Logic layer: "<<iLogicLayer<<" Hits: ";
    for(unsigned int iInput=0;iInput<14;++iInput){
      out<<OMTFConfiguration::instance()->measurements4Dref[iProcessor][iRegion][iLogicLayer][iInput]<<"\t";
    }
    out<<std::endl;
  }
  /*
  out<<"Measurement hits"<<std::endl;
  for(unsigned int iLogicLayer=0;iLogicLayer<OMTFConfiguration::instance()->nLayers;++iLogicLayer){
    out<<"Logic layer: "<<iLogicLayer<<" Hits: ";
    for(unsigned int iInput=0;iInput<14;++iInput){
      out<<OMTFConfiguration::instance()->measurements4D[iProcessor][iRegion][iLogicLayer][iInput]<<"\t";
    }
    out<<std::endl;
  }
  */
  out<<std::endl;
}