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;
}