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