L1MuDTQualPatternLut.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 "CondFormats/L1TObjects/interface/L1MuDTQualPatternLut.h"

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

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

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

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

using namespace std;

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

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

L1MuDTQualPatternLut::L1MuDTQualPatternLut() {

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

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

}


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

L1MuDTQualPatternLut::~L1MuDTQualPatternLut() {

  LUT::iterator iter = m_lut.begin();
  while ( iter != m_lut.end() ) {
    (*iter).second.second.clear(); 
    iter++;
  }

  m_lut.clear();

}


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

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

  m_lut.clear();

}


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

  // get directory name
  string defaultPath = "L1TriggerConfig/DTTrackFinder/parameters/";
  string eau_dir = "L1TriggerData/DTTrackFinder/Eau/";
  string emu_str = "";

  // loop over all sector processors
  for ( int sp = 0; sp < 6; sp++ ) { 
    switch ( sp ) {
      case 0  : { emu_str = "QualPatternList_SP1"; break; }
      case 1  : { emu_str = "QualPatternList_SP2"; break; }
      case 2  : { emu_str = "QualPatternList_SP3"; break; }
      case 3  : { emu_str = "QualPatternList_SP4"; break; }
      case 4  : { emu_str = "QualPatternList_SP5"; break; }
      case 5  : { emu_str = "QualPatternList_SP6"; break; }
    }  
  
    // 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 
    file.ignoreLines(14);

    // 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;
      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 L1MuDTQualPatternLut::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 L1MuDTQualPatternLut::getCoarseEta(int sp, int adr) const {

  LUT::const_iterator it = m_lut.find(make_pair(sp,adr));
  if ( it == m_lut.end() ) {
    cerr << "Error: L1MuDTQualPatternLut: 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>& L1MuDTQualPatternLut::getQualifiedPatterns(int sp, int adr) const {

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

}