GlobalScales.cc

Go to the documentation of this file.

#include "L1Trigger/L1TGlobal/interface/GlobalScales.h"

#include <iostream>
#include <fstream>
#include <iomanip>

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/MessageLogger/interface/MessageDrop.h"

// constructor
l1t::GlobalScales::GlobalScales() {}

// destructor
l1t::GlobalScales::~GlobalScales() {}

std::string l1t::GlobalScales::getScalesName() const { return m_ScaleSetName; }

void l1t::GlobalScales::setLUT_CalMuEta(const std::string& lutName, std::vector<long long> lut) {
  if (m_lut_CalMuEta.count(lutName) != 0) {
    LogTrace("GlobalScales") << "      LUT \"" << lutName << "\"already exists in the LUT map- not inserted!"
                             << std::endl;
    return;
  }

  // Insert this LUT into the Table
  m_lut_CalMuEta.insert(std::map<std::string, std::vector<long long>>::value_type(lutName, lut));

  return;
}

void l1t::GlobalScales::setLUT_CalMuPhi(const std::string& lutName, std::vector<long long> lut) {
  if (m_lut_CalMuPhi.count(lutName) != 0) {
    LogTrace("GlobalScales") << "      LUT \"" << lutName << "\"already exists in the LUT map- not inserted!"
                             << std::endl;
    return;
  }

  // Insert this LUT into the Table
  m_lut_CalMuPhi.insert(std::map<std::string, std::vector<long long>>::value_type(lutName, lut));

  return;
}

void l1t::GlobalScales::setLUT_DeltaEta(const std::string& lutName,
                                        std::vector<long long> lut,
                                        unsigned int precision) {
  if (m_lut_DeltaEta.count(lutName) != 0) {
    LogTrace("GlobalScales") << "      LUT \"" << lutName << "\"already exists in the LUT map- not inserted!"
                             << std::endl;
    return;
  }

  // Insert this LUT into the Table
  m_lut_DeltaEta.insert(std::map<std::string, std::vector<long long>>::value_type(lutName, lut));
  m_Prec_DeltaEta.insert(std::map<std::string, unsigned int>::value_type(lutName, precision));

  return;
}

void l1t::GlobalScales::setLUT_DeltaPhi(const std::string& lutName,
                                        std::vector<long long> lut,
                                        unsigned int precision) {
  if (m_lut_DeltaPhi.count(lutName) != 0) {
    LogTrace("GlobalScales") << "      LUT \"" << lutName << "\"already exists in the LUT map- not inserted!"
                             << std::endl;
    return;
  }

  // Insert this LUT into the Table
  m_lut_DeltaPhi.insert(std::map<std::string, std::vector<long long>>::value_type(lutName, lut));
  m_Prec_DeltaPhi.insert(std::map<std::string, unsigned int>::value_type(lutName, precision));

  return;
}

void l1t::GlobalScales::setLUT_Pt(const std::string& lutName, std::vector<long long> lut, unsigned int precision) {
  if (m_lut_Pt.count(lutName) != 0) {
    LogTrace("GlobalScales") << "      LUT \"" << lutName << "\"already exists in the LUT map- not inserted!"
                             << std::endl;
    return;
  }

  // Insert this LUT into the Table
  m_lut_Pt.insert(std::map<std::string, std::vector<long long>>::value_type(lutName, lut));
  m_Prec_Pt.insert(std::map<std::string, unsigned int>::value_type(lutName, precision));

  return;
}

void l1t::GlobalScales::setLUT_Cosh(const std::string& lutName, std::vector<long long> lut, unsigned int precision) {
  if (m_lut_Cosh.count(lutName) != 0) {
    LogTrace("GlobalScales") << "      LUT \"" << lutName << "\"already exists in the LUT map- not inserted!"
                             << std::endl;
    return;
  }

  // Insert this LUT into the Table
  m_lut_Cosh.insert(std::map<std::string, std::vector<long long>>::value_type(lutName, lut));
  m_Prec_Cosh.insert(std::map<std::string, unsigned int>::value_type(lutName, precision));

  return;
}

void l1t::GlobalScales::setLUT_Cos(const std::string& lutName, std::vector<long long> lut, unsigned int precision) {
  if (m_lut_Cos.count(lutName) != 0) {
    LogTrace("GlobalScales") << "      LUT \"" << lutName << "\"already exists in the LUT map- not inserted!"
                             << std::endl;
    return;
  }

  // Insert this LUT into the Table
  m_lut_Cos.insert(std::map<std::string, std::vector<long long>>::value_type(lutName, lut));
  m_Prec_Cos.insert(std::map<std::string, unsigned int>::value_type(lutName, precision));

  return;
}

void l1t::GlobalScales::setLUT_Sin(const std::string& lutName, std::vector<long long> lut, unsigned int precision) {
  if (m_lut_Sin.count(lutName) != 0) {
    LogTrace("GlobalScales") << "      LUT \"" << lutName << "\"already exists in the LUT map- not inserted!"
                             << std::endl;
    return;
  }

  // Insert this LUT into the Table
  m_lut_Sin.insert(std::map<std::string, std::vector<long long>>::value_type(lutName, lut));
  m_Prec_Sin.insert(std::map<std::string, unsigned int>::value_type(lutName, precision));

  return;
}

long long l1t::GlobalScales::getLUT_CalMuEta(const std::string& lutName, int element) const {
  long long value = 0;

  if (element < 0) {
    edm::LogError("GlobalScales") << "Error: Negative index, " << element << ", requested for CalMuEta LUT ( "
                                  << lutName << ")" << std::endl;
  } else if (element >= (int)m_lut_CalMuEta.find(lutName)->second.size()) {
    edm::LogError("GlobalScales") << "Error: Element Requested " << element << " too large for CalMuEta LUT ("
                                  << lutName << ") size = " << m_lut_CalMuEta.find(lutName)->second.size() << std::endl;
  } else {
    value = m_lut_CalMuEta.find(lutName)->second.at(element);
  }
  return value;
}

long long l1t::GlobalScales::getLUT_CalMuPhi(const std::string& lutName, int element) const {
  long long value = 0;

  if (element < 0) {
    edm::LogError("GlobalScales") << "Error: Negative index, " << element << ", requested for CalMuPhi LUT ( "
                                  << lutName << ")" << std::endl;
  } else if (element >= (int)m_lut_CalMuPhi.find(lutName)->second.size()) {
    edm::LogError("GlobalScales") << "Error: Element Requested " << element << " too large for CalMuPhi LUT ("
                                  << lutName << ") size = " << m_lut_CalMuPhi.find(lutName)->second.size() << std::endl;
  } else {
    value = m_lut_CalMuPhi.find(lutName)->second.at(element);
  }
  return value;
}

long long l1t::GlobalScales::getLUT_DeltaEta(std::string lutName, int element) const {
  long long value = 0;

  //first check whether this LUT exists
  if (m_lut_DeltaEta.find(lutName) == m_lut_DeltaEta.end()) {
    //does not exist. Check for oppoisite ordering
    std::size_t pos = lutName.find("-");
    std::string name = lutName.substr(pos + 1);
    name += "-";
    name += lutName.substr(0, pos);

    //check again
    if (m_lut_DeltaEta.find(name) == m_lut_DeltaEta.end()) {
      edm::LogError("GlobalScales") << "Warning: No LUT by name " << lutName << " or " << name << " exists! "
                                    << std::endl;
      return value;
    } else {
      lutName = name;
    }
  }

  if (element < 0) {
    edm::LogError("GlobalScales") << "Error: Negative index, " << element << ", requested for DeltaEta LUT ( "
                                  << lutName << ")" << std::endl;
  } else if (element >= (int)m_lut_DeltaEta.find(lutName)->second.size()) {
    edm::LogError("GlobalScales") << "Error: Element Requested " << element << " too large for DeltaEta LUT ("
                                  << lutName << ") size = " << m_lut_DeltaEta.find(lutName)->second.size() << std::endl;
  } else {
    value = m_lut_DeltaEta.find(lutName)->second.at(element);
  }
  return value;
}
unsigned int l1t::GlobalScales::getPrec_DeltaEta(const std::string& lutName) const {
  unsigned int value = 0;

  if (m_Prec_DeltaEta.find(lutName) != m_Prec_DeltaEta.end()) {
    value = m_Prec_DeltaEta.find(lutName)->second;
  } else {
    //does not exist. Check for oppoisite ordering
    std::size_t pos = lutName.find("-");
    std::string name = lutName.substr(pos + 1);
    name += "-";
    name += lutName.substr(0, pos);

    //check again
    if (m_Prec_DeltaEta.find(name) != m_Prec_DeltaEta.end()) {
      value = m_Prec_DeltaEta.find(name)->second;
    } else {
      edm::LogError("GlobalScales") << "Warning: LUT " << lutName << " for DeltaEta not found" << std::endl;
    }
  }
  return value;
}

long long l1t::GlobalScales::getLUT_DeltaPhi(std::string lutName, int element) const {
  long long value = 0;

  //first check whether this LUT exists
  if (m_lut_DeltaPhi.find(lutName) == m_lut_DeltaPhi.end()) {
    //does not exist. Check for oppoisite ordering
    std::size_t pos = lutName.find("-");
    std::string name = lutName.substr(pos + 1);
    name += "-";
    name += lutName.substr(0, pos);

    //check again
    if (m_lut_DeltaPhi.find(name) == m_lut_DeltaPhi.end()) {
      edm::LogError("GlobalScales") << "Warning: No LUT by name " << lutName << " or " << name << " exists! "
                                    << std::endl;
      return value;
    } else {
      lutName = name;
    }
  }

  if (element < 0) {
    edm::LogError("GlobalScales") << "Error: Negative index, " << element << ", requested for DeltaPhi LUT ( "
                                  << lutName << ")" << std::endl;
  } else if (element >= (int)m_lut_DeltaPhi.find(lutName)->second.size()) {
    edm::LogError("GlobalScales") << "Error: Element Requested " << element << " too large for DeltaPhi LUT ("
                                  << lutName << ") size = " << m_lut_DeltaPhi.find(lutName)->second.size() << std::endl;
  } else {
    value = m_lut_DeltaPhi.find(lutName)->second.at(element);
  }
  return value;
}
unsigned int l1t::GlobalScales::getPrec_DeltaPhi(const std::string& lutName) const {
  unsigned int value = 0;

  if (m_Prec_DeltaPhi.find(lutName) != m_Prec_DeltaPhi.end()) {
    value = m_Prec_DeltaPhi.find(lutName)->second;
  } else {
    //does not exist. Check for oppoisite ordering
    std::size_t pos = lutName.find("-");
    std::string name = lutName.substr(pos + 1);
    name += "-";
    name += lutName.substr(0, pos);

    //check again
    if (m_Prec_DeltaPhi.find(name) != m_Prec_DeltaPhi.end()) {
      value = m_Prec_DeltaPhi.find(name)->second;
    } else {
      edm::LogError("GlobalScales") << "Warning: LUT " << lutName << " for DeltaPhi not found" << std::endl;
    }
  }
  return value;
}

long long l1t::GlobalScales::getLUT_Pt(const std::string& lutName, int element) const {
  long long value = 0;

  if (element < 0) {
    edm::LogError("GlobalScales") << "Error: Negative index, " << element << ", requested for Pt LUT ( " << lutName
                                  << ")" << std::endl;
  } else if (element >= (int)m_lut_Pt.find(lutName)->second.size()) {
    edm::LogError("GlobalScales") << "Error: Element Requested " << element << " too large for Pt LUT (" << lutName
                                  << ") size = " << m_lut_Pt.find(lutName)->second.size() << std::endl;
  } else {
    value = m_lut_Pt.find(lutName)->second.at(element);
  }
  return value;
}
unsigned int l1t::GlobalScales::getPrec_Pt(const std::string& lutName) const {
  unsigned int value = 0;

  if (m_Prec_Pt.find(lutName) != m_Prec_Pt.end()) {
    value = m_Prec_Pt.find(lutName)->second;
  } else {
    edm::LogError("GlobalScales") << "Warning: LUT " << lutName << " for Pt not found" << std::endl;
  }
  return value;
}

long long l1t::GlobalScales::getLUT_DeltaEta_Cosh(std::string lutName, int element) const {
  long long value = 0;

  //first check whether this LUT exists
  if (m_lut_Cosh.find(lutName) == m_lut_Cosh.end()) {
    //does not exist. Check for oppoisite ordering
    std::size_t pos = lutName.find("-");
    std::string name = lutName.substr(pos + 1);
    name += "-";
    name += lutName.substr(0, pos);

    //check again
    if (m_lut_Cosh.find(name) == m_lut_Cosh.end()) {
      edm::LogError("GlobalScales") << "Warning: No LUT by name " << lutName << " or " << name << " exists! "
                                    << std::endl;
      return value;
    } else {
      lutName = name;
    }
  }

  if (element < 0) {
    edm::LogError("GlobalScales") << "Error: Negative index, " << element << ", requested for Cosh LUT ( " << lutName
                                  << ")" << std::endl;
  } else if (element >= (int)m_lut_Cosh.find(lutName)->second.size()) {
    edm::LogError("GlobalScales") << "Error: Element Requested " << element << " too large for Cosh LUT (" << lutName
                                  << ") size = " << m_lut_Cosh.find(lutName)->second.size() << std::endl;
  } else {
    value = m_lut_Cosh.find(lutName)->second.at(element);
  }
  return value;
}
unsigned int l1t::GlobalScales::getPrec_DeltaEta_Cosh(const std::string& lutName) const {
  unsigned int value = 0;

  if (m_Prec_Cosh.find(lutName) != m_Prec_Cosh.end()) {
    value = m_Prec_Cosh.find(lutName)->second;
  } else {
    //does not exist. Check for oppoisite ordering
    std::size_t pos = lutName.find("-");
    std::string name = lutName.substr(pos + 1);
    name += "-";
    name += lutName.substr(0, pos);

    //check again
    if (m_Prec_Cosh.find(name) != m_Prec_Cosh.end()) {
      value = m_Prec_Cosh.find(name)->second;
    } else {
      edm::LogError("GlobalScales") << "Warning: LUT " << lutName << " for Cosh not found" << std::endl;
    }
  }
  return value;
}

long long l1t::GlobalScales::getLUT_DeltaPhi_Cos(std::string lutName, int element) const {
  long long value = 0;

  //first check whether this LUT exists
  if (m_lut_Cos.find(lutName) == m_lut_Cos.end()) {
    //does not exist. Check for oppoisite ordering
    std::size_t pos = lutName.find("-");
    std::string name = lutName.substr(pos + 1);
    name += "-";
    name += lutName.substr(0, pos);

    //check again
    if (m_lut_Cos.find(name) == m_lut_Cos.end()) {
      edm::LogError("GlobalScales") << "Warning: No LUT by name " << lutName << " or " << name << " exists! "
                                    << std::endl;
      return value;
    } else {
      lutName = name;
    }
  }

  if (element < 0) {
    edm::LogError("GlobalScales") << "Error: Negative index, " << element << ", requested for Cos LUT ( " << lutName
                                  << ")" << std::endl;
  } else if (element >= (int)m_lut_DeltaPhi.find(lutName)->second.size()) {
    edm::LogError("GlobalScales") << "Error: Element Requested " << element << " too large for Cos LUT (" << lutName
                                  << ") size = " << m_lut_Cos.find(lutName)->second.size() << std::endl;
  } else {
    value = m_lut_Cos.find(lutName)->second.at(element);
  }
  return value;
}

long long l1t::GlobalScales::getLUT_Cos(const std::string& lutName, int element) const {
  long long value = 0;

  //first check whether this LUT exists
  if (m_lut_Cos.find(lutName) == m_lut_Cos.end()) {
    edm::LogError("GlobalScales") << "Warning: No LUT by name " << lutName << " exists! " << std::endl;
    return value;
  }

  if (element < 0) {
    edm::LogError("GlobalScales") << "Error: Negative index, " << element << ", requested for Cos LUT ( " << lutName
                                  << ")" << std::endl;
  } else if (element >= (int)m_lut_Cos.find(lutName)->second.size()) {
    edm::LogError("GlobalScales") << "Error: Element Requested " << element << " too large for Cos LUT (" << lutName
                                  << ") size = " << m_lut_Cos.find(lutName)->second.size() << std::endl;
  } else {
    value = m_lut_Cos.find(lutName)->second.at(element);
  }
  return value;
}

long long l1t::GlobalScales::getLUT_Sin(const std::string& lutName, int element) const {
  long long value = 0;

  //first check whether this LUT exists
  if (m_lut_Sin.find(lutName) == m_lut_Sin.end()) {
    edm::LogError("GlobalScales") << "Warning: No LUT by name " << lutName << " exists! " << std::endl;
    return value;
  }

  if (element < 0) {
    edm::LogError("GlobalScales") << "Error: Negative index, " << element << ", requested for Sin LUT ( " << lutName
                                  << ")" << std::endl;
  } else if (element >= (int)m_lut_Sin.find(lutName)->second.size()) {
    edm::LogError("GlobalScales") << "Error: Element Requested " << element << " too large for Sin LUT (" << lutName
                                  << ") size = " << m_lut_Sin.find(lutName)->second.size() << std::endl;
  } else {
    value = m_lut_Sin.find(lutName)->second.at(element);
  }
  return value;
}

unsigned int l1t::GlobalScales::getPrec_DeltaPhi_Cos(const std::string& lutName) const {
  unsigned int value = 0;

  if (m_Prec_Cos.find(lutName) != m_Prec_Cos.end()) {
    value = m_Prec_Cos.find(lutName)->second;
  } else {
    //does not exist. Check for oppoisite ordering
    std::size_t pos = lutName.find("-");
    std::string name = lutName.substr(pos + 1);
    name += "-";
    name += lutName.substr(0, pos);

    //check again
    if (m_Prec_Cos.find(name) != m_Prec_Cos.end()) {
      value = m_Prec_Cos.find(name)->second;
    } else {
      edm::LogError("GlobalScales") << "Warning: LUT " << lutName << " for Cos not found" << std::endl;
    }
  }
  return value;
}

unsigned int l1t::GlobalScales::getPrec_Cos(const std::string& lutName) const {
  unsigned int value = 0;

  if (m_Prec_Sin.find(lutName) != m_Prec_Sin.end()) {
    value = m_Prec_Sin.find(lutName)->second;
  } else {
    edm::LogError("GlobalScales") << "Warning: LUT " << lutName << " for Sin not found" << std::endl;
  }
  return value;
}

unsigned int l1t::GlobalScales::getPrec_Sin(const std::string& lutName) const {
  unsigned int value = 0;

  if (m_Prec_Sin.find(lutName) != m_Prec_Sin.end()) {
    value = m_Prec_Sin.find(lutName)->second;
  } else {
    edm::LogError("GlobalScales") << "Warning: LUT " << lutName << " for Sin not found" << std::endl;
  }
  return value;
}

void l1t::GlobalScales::dumpAllLUTs(std::ostream& myCout) const {
  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_CalMuEta.begin();
       itr != m_lut_CalMuEta.end();
       itr++) {
    dumpLUT(myCout, 1, itr->first);
  }

  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_CalMuPhi.begin();
       itr != m_lut_CalMuPhi.end();
       itr++) {
    dumpLUT(myCout, 2, itr->first);
  }

  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_DeltaEta.begin();
       itr != m_lut_DeltaEta.end();
       itr++) {
    dumpLUT(myCout, 3, itr->first);
  }

  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_DeltaPhi.begin();
       itr != m_lut_DeltaPhi.end();
       itr++) {
    dumpLUT(myCout, 4, itr->first);
  }

  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_Cosh.begin(); itr != m_lut_Cosh.end();
       itr++) {
    dumpLUT(myCout, 5, itr->first);
  }

  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_Cos.begin(); itr != m_lut_Cos.end();
       itr++) {
    dumpLUT(myCout, 6, itr->first);
  }

  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_Sin.begin(); itr != m_lut_Sin.end();
       itr++) {
    dumpLUT(myCout, 7, itr->first);
  }

  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_Pt.begin(); itr != m_lut_Pt.end();
       itr++) {
    dumpLUT(myCout, 8, itr->first);
  }
}

void l1t::GlobalScales::dumpLUT(std::ostream& myCout, int LUTtype, std::string name) const {
  std::vector<long long> dumpV;
  unsigned int prec = 0;
  std::string type = "Null";
  switch (LUTtype) {
    case 1: {
      dumpV = m_lut_CalMuEta.find(name)->second;
      type = "Cal-Mu Eta";
      break;
    }
    case 2: {
      dumpV = m_lut_CalMuPhi.find(name)->second;
      type = "Cal-Mu Phi";
      break;
    }
    case 3: {
      dumpV = m_lut_DeltaEta.find(name)->second;
      prec = m_Prec_DeltaEta.find(name)->second;
      type = "Delta Eta";
      break;
    }
    case 4: {
      dumpV = m_lut_DeltaPhi.find(name)->second;
      prec = m_Prec_DeltaPhi.find(name)->second;
      type = "Delta Phi";
      break;
    }
    case 5: {
      dumpV = m_lut_Cosh.find(name)->second;
      prec = m_Prec_Cosh.find(name)->second;
      type = "Cosh";
      break;
    }
    case 6: {
      dumpV = m_lut_Cos.find(name)->second;
      prec = m_Prec_Cos.find(name)->second;
      type = "Cos";
      break;
    }
    case 7: {
      dumpV = m_lut_Sin.find(name)->second;
      prec = m_Prec_Sin.find(name)->second;
      type = "Sin";
      break;
    }
    case 8: {
      dumpV = m_lut_Pt.find(name)->second;
      prec = m_Prec_Pt.find(name)->second;
      type = "Pt";
      break;
    }
  }

  myCout << "=========================================" << std::endl;
  myCout << "   LUT: " << type << "  Table: " << name << "   Size = " << dumpV.size() << " Precision " << prec
         << std::endl;
  myCout << "=========================================" << std::endl;
  for (int i = 0; i < (int)dumpV.size(); i++) {
    myCout << " Element " << i << " " << dumpV.at(i) << std::endl;
  }
}

void l1t::GlobalScales::print(std::ostream& myCout) const {
  myCout << "\n  *************  L1T Global Scales ************" << std::endl;

  myCout << "   Muon Scales: " << std::endl;
  printScale(m_muScales, myCout);

  myCout << "   EG Scales: " << std::endl;
  printScale(m_egScales, myCout);

  myCout << "   Tau Scales: " << std::endl;
  printScale(m_tauScales, myCout);

  myCout << "   Jet Scales: " << std::endl;
  printScale(m_jetScales, myCout);

  myCout << "   HTT Scales: " << std::endl;
  printScale(m_httScales, myCout);

  myCout << "   ETT Scales: " << std::endl;
  printScale(m_ettScales, myCout);

  myCout << "   HTM Scales: " << std::endl;
  printScale(m_htmScales, myCout);

  myCout << "   ETM Scales: " << std::endl;
  printScale(m_etmScales, myCout);

  myCout << std::endl;
  myCout << "   LUTs Stored: " << std::endl;
  myCout << " CalMuEta:";
  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_CalMuEta.begin();
       itr != m_lut_CalMuEta.end();
       itr++) {
    myCout << " " << itr->first;
  }
  myCout << std::endl;

  myCout << " CalMuPhi:";
  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_CalMuPhi.begin();
       itr != m_lut_CalMuPhi.end();
       itr++) {
    myCout << " " << itr->first;
  }
  myCout << std::endl;

  myCout << " DeltaEta:";
  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_DeltaEta.begin();
       itr != m_lut_DeltaEta.end();
       itr++) {
    myCout << " " << itr->first;
  }
  myCout << std::endl;

  myCout << " DeltaPhi:";
  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_DeltaPhi.begin();
       itr != m_lut_DeltaPhi.end();
       itr++) {
    myCout << " " << itr->first;
  }
  myCout << std::endl;

  myCout << " Cos:     ";
  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_Cos.begin(); itr != m_lut_Cos.end();
       itr++) {
    myCout << " " << itr->first;
  }
  myCout << std::endl;

  myCout << " Sin:     ";
  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_Sin.begin(); itr != m_lut_Sin.end();
       itr++) {
    myCout << " " << itr->first;
  }
  myCout << std::endl;

  myCout << " Cosh:    ";
  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_Cosh.begin(); itr != m_lut_Cosh.end();
       itr++) {
    myCout << " " << itr->first;
  }
  myCout << std::endl;

  myCout << " Pt:      ";
  for (std::map<std::string, std::vector<long long>>::const_iterator itr = m_lut_Pt.begin(); itr != m_lut_Pt.end();
       itr++) {
    myCout << " " << itr->first;
  }
  myCout << std::endl;
}
void l1t::GlobalScales::printScale(ScaleParameters scale, std::ostream& myCout) const {
  myCout << "    Pt Min   = " << std::setw(10) << scale.etMin << "      Pt Max   = " << std::setw(10) << scale.etMax
         << "      Pt Step  = " << std::setw(10) << scale.etStep << "      Number   = " << std::setw(10)
         << scale.etBins.size() << std::endl;
  myCout << "\n";
  for (int i = 0; i < int(scale.etBins.size()); i++) {
    std::pair<double, double> binEdges = scale.etBins.at(i);
    myCout << "        etBins[" << i << "]\trange:\t" << binEdges.first << " - " << binEdges.second << std::endl;
  }

  myCout << "\n    Phi Min  = " << std::setw(10) << scale.phiMin << "      Phi Max  = " << std::setw(10) << scale.phiMax
         << "      Phi Step = " << std::setw(10) << scale.phiStep << "      Number   = " << std::setw(10)
         << scale.phiBins.size() << std::endl;
  myCout << "\n";
  for (int i = 0; i < int(scale.phiBins.size()); i++) {
    std::pair<double, double> binEdges = scale.phiBins.at(i);
    myCout << "        phiBins[" << i << "]\trange:\t" << binEdges.first << " - " << binEdges.second << std::endl;
  }

  myCout << "\n    Eta Min  = " << std::setw(10) << scale.etaMin << "      Eta Max  = " << std::setw(10) << scale.etaMax
         << "      Eta Step = " << std::setw(10) << scale.etaStep << "      Number   = " << std::setw(10)
         << scale.etaBins.size() << std::endl;
  myCout << "\n";
  for (int i = 0; i < int(scale.etaBins.size()); i++) {
    std::pair<double, double> binEdges = scale.etaBins.at(i);
    myCout << "        etaBins[" << i << "]\trange:\t" << binEdges.first << " - " << binEdges.second << std::endl;
  }
}