L1MuBMTQualPatternLut.cc

Go to the documentation of this file.

//-------------------------------------------------
//
//   Class: L1MuDTQualPatternLut
//
//   Description: Look-up tables for eta matching unit (EMU)
//                stores lists of qualified patterns and
//                coarse eta values 
//
//
//   $Date: 2007/03/30 07:48:02 $
//   $Revision: 1.1 $
//
//   Author :
//   N. Neumeister            CERN EP
//   J. Troconiz              UAM Madrid
//
//--------------------------------------------------

//-----------------------
// This Class's Header --
//-----------------------

#include "L1Trigger/L1TMuonBarrel/interface/L1MuBMTQualPatternLut.h"

//---------------
// C++ Headers --
//---------------

#include <iostream>
#include <iomanip>
#include <string>

//-------------------------------
// Collaborating Class Headers --
//-------------------------------

#include "FWCore/ParameterSet/interface/FileInPath.h"
#include "CondFormats/L1TObjects/interface/L1TriggerLutFile.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

using namespace std;

// --------------------------------
//       class L1MuDTQualPatternLut
//---------------------------------

//----------------
// Constructors --
//----------------

L1MuBMTQualPatternLut::L1MuBMTQualPatternLut() {

  //  if ( load() != 0 ) {
  //    cout << "Can not open files to load eta matching look-up tables for DTTrackFinder!" << endl;
  //  }

  //  if ( L1MuDTTFConfig::Debug(6) ) print();

}


//--------------
// Destructor --
//--------------

L1MuBMTQualPatternLut::~L1MuBMTQualPatternLut() {

}


//--------------
// Operations --
//--------------

//
// reset look-up tables
//
void L1MuBMTQualPatternLut::reset() {

  m_lut.clear();

}


//
// load look-up tables for EMU
//
int L1MuBMTQualPatternLut::load() {

  // get directory name
  string defaultPath ="L1Trigger/"; //"L1TriggerConfig/DTTrackFinder/parameters/";
    string eau_dir = "L1TMuon/data/bmtf_luts/LUTs_Ass/";//"L1TriggerData/DTTrackFinder/Eau/";
  string emu_str = "";

  // loop over all sector processors
  for ( int sp = 0; sp < 6; sp++ ) { 

    emu_str="QualPatternList_SP"+std::to_string(sp+1);
  
    // assemble file name
    edm::FileInPath lut_f = edm::FileInPath(string(defaultPath + eau_dir + emu_str + ".lut"));
    string emu_file = lut_f.fullPath();

    // open file
    L1TriggerLutFile file(emu_file);
    if ( file.open() != 0 ) return -1;
    //    if ( L1MuDTTFConfig::Debug(1) ) cout << "Reading file : " 
    //                                         << file.getName() << endl; 

    // ignore comment lines 
    int skip2=getIgnoredLines(file);
    file.ignoreLines(skip2);

    // read file
    while ( file.good() ) {
  
      int id  = file.readInteger();
      if ( !file.good() ) break;
      int eta = file.readInteger();
      if ( !file.good() ) break;
      int num = file.readInteger();
      if ( !file.good() ) break;

      vector<short> patternlist;
      patternlist.reserve(num);

      for ( int i = 0; i < num; i++ ) {
        int pattern = file.readInteger();
        patternlist.push_back(pattern);
      }
      
      m_lut[make_pair(sp+1,id)] = make_pair(eta,patternlist);

      if ( !file.good() ) { file.close(); break; }
    
    }

    file.close();
  }  
    
  return 0;

}


//
// print look-up tables for EMU
//
void L1MuBMTQualPatternLut::print() const {

  cout << endl;
  cout << "L1 barrel Track Finder Qual Pattern look-up tables :" << endl;
  cout << "====================================================" << endl;
  cout << endl;

  int spold = 0;
  
  LUT::const_iterator iter = m_lut.begin();
  while ( iter != m_lut.end() ) {
    int sp = (*iter).first.first;
    if ( sp != spold ) {
      cout << endl;
      cout << "Qualified Patterns for Sector Processor " << setw(1) << sp << " :" << endl;
      cout << "===========================================" << endl;
      cout << endl;
      spold = sp;
    }
    cout << setw(2) << (*iter).first.second << "  "
         << setw(3) << (*iter).second.first << "  "
         << setw(5) << (*iter).second.second.size() << " : ";
    const vector<short>& patternlist = (*iter).second.second;
    vector<short>::const_iterator it;
    for ( it = patternlist.begin(); it != patternlist.end(); it++ ) {
      cout << setw(5) << (*it) << " ";
    }
    cout << endl;
    iter++;
  }
 
  cout << endl;

}


//
// get coarse eta value for a given sector processor [1-6] and address [1-22] 
//
int L1MuBMTQualPatternLut::getCoarseEta(int sp, int adr) const {

  LUT::const_iterator it = m_lut.find(make_pair(sp,adr));
  if ( it == m_lut.end() ) {
     edm::LogError ("L1MuBMTQualPatternLut") << "Error: L1MuBMTQualPatternLut: no coarse eta found for address " << adr << endl;
    return 0;
  }
  return (*it).second.first;

}


//
// get list of qualified patterns for a given sector processor [1-6] and address [1-22]
//
const vector<short>& L1MuBMTQualPatternLut::getQualifiedPatterns(int sp, int adr) const {

  LUT::const_iterator it = m_lut.find(make_pair(sp,adr));
  if ( it == m_lut.end() ) {
    edm::LogError ("L1MuBMTQualPatternLut") << "Error: L1MuBMTQualPatternLut: no pattern list found for address " << adr << endl;
  }
  return (*it).second.second;

}

int L1MuBMTQualPatternLut::getIgnoredLines(L1TriggerLutFile file) const{
   if ( file.open() != 0 ) return -1;
   int skip=0;
   while ( file.good() ) {

     string str=file.readString();
     if (str.find("#")==0) skip+=1;
           //cout<<"here "<<str<<" found "<<str.find("#")<<endl;
           if ( !file.good() ) { file.close(); break; }
     }
     file.close();

     // skip aditional lines of comments between "---".
     skip += 2;

     return skip;
}