L1MuBMTEtaPatternLut.cc

Go to the documentation of this file.

//-------------------------------------------------
//
//   Class: L1MuDTEtaPatternLut
//
//   Description: Look-up table for eta track finder
//
//
//   $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/L1MuBMTEtaPatternLut.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 L1MuDTEtaPatternLut
//---------------------------------
 
//----------------
// Constructors --
//----------------
 
L1MuBMTEtaPatternLut::L1MuBMTEtaPatternLut() {
  //  if ( load() != 0 ) {
  //    cout << "Can not open files to load eta track finder look-up tables for DTTrackFinder!" << endl;
  //  }
 
  //  if ( L1MuDTTFConfig::Debug(6) ) print();
}
 
//--------------
// Destructor --
//--------------
 
L1MuBMTEtaPatternLut::~L1MuBMTEtaPatternLut() { m_lut.clear(); }
 
//--------------
// Operations --
//--------------
 
//
// reset look-up table
//
void L1MuBMTEtaPatternLut::reset() { m_lut.clear(); }
 
//
// load pattern look-up table for ETF
//
int L1MuBMTEtaPatternLut::load() {
  // get directory name
  string defaultPath = "L1Trigger/";
  string eau_dir = "L1TMuon/data/bmtf_luts/LUTs_Ass/";
 
  // assemble file name
  edm::FileInPath lut_f = edm::FileInPath(string(defaultPath + eau_dir + "ETFPatternList.lut"));
  string etf_file = lut_f.fullPath();
 
  // open file
  L1TriggerLutFile file(etf_file);
  if (file.open() != 0)
    return -1;
  //  if ( L1MuDTTFConfig::Debug(1) ) cout << "Reading file : "
  //                                       << file.getName() << endl;
 
  // ignore comment lines (always at the beginning)
  int skip2 = getIgnoredLines(file);
  file.ignoreLines(skip2);
 
  // read patterns
  while (file.good()) {
    int id = file.readInteger();
    if (!file.good())
      break;
    string pat = file.readString();
    if (!file.good())
      break;
    int qual = file.readInteger();
    if (!file.good())
      break;
    int eta = file.readInteger();
    if (!file.good())
      break;
    L1MuDTEtaPattern pattern(id, pat, eta, qual);
 
    m_lut[pattern.id()] = pattern;
 
    if (!file.good()) {
      file.close();
      break;
    }
  }
 
  file.close();
 
  return 0;
}
 
//
// print pattern look-up table for ETF
//
void L1MuBMTEtaPatternLut::print() const {
  cout << endl;
  cout << "L1 barrel Track Finder ETA Pattern look-up table :" << endl;
  cout << "==================================================" << endl;
  cout << endl;
 
  cout << "ETF Patterns : " << m_lut.size() << endl;
  cout << "======================" << endl;
  cout << endl;
 
  LUT::const_iterator iter = m_lut.begin();
  while (iter != m_lut.end()) {
    cout << (*iter).second << endl;
    iter++;
  }
 
  cout << endl;
}
 
//
// get pattern with a given ID
//
L1MuDTEtaPattern L1MuBMTEtaPatternLut::getPattern(int id) const {
  LUT::const_iterator it = m_lut.find(id);
  if (it == m_lut.end()) {
    edm::LogError("L1MuBMTEtaPatternLut: fine eta not found")
        << "Error: L1MuBMTEtaPatternLut: pattern not found : " << id << endl;
  }
  return (*it).second;
}
 
int L1MuBMTEtaPatternLut::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();
  return skip;
}