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 < (int)m_lut_CalMuEta.find(lutName)->second.size()) {
     value = m_lut_CalMuEta.find(lutName)->second.at(element);
   } else {
     edm::LogError("GlobalScales") << "Warning: Element Requested " << element << " too large for CalMuEta LUT (" << lutName <<") size = " << m_lut_CalMuEta.find(lutName)->second.size() << std::endl;
   }
   return value;
}

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

   if(element < (int)m_lut_CalMuPhi.find(lutName)->second.size()) {
     value = m_lut_CalMuPhi.find(lutName)->second.at(element);
   } else {
     edm::LogError("GlobalScales") << "Warning: Element Requested " << element << " too large for CalMuPhi LUT (" << lutName <<") size = " << m_lut_CalMuPhi.find(lutName)->second.size() << std::endl;
   }
   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 < (int)m_lut_DeltaEta.find(lutName)->second.size()) {
     value = m_lut_DeltaEta.find(lutName)->second.at(element);
   } else {
     edm::LogError("GlobalScales") << "Warning: Element Requested " << element << " too large for DeltaEta LUT (" << lutName <<") size = " << m_lut_DeltaEta.find(lutName)->second.size() << std::endl;
   }
   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 < (int)m_lut_DeltaPhi.find(lutName)->second.size()) {
     value = m_lut_DeltaPhi.find(lutName)->second.at(element);
   } else {
     edm::LogError("GlobalScales") << "Warning: Element Requested " << element << " too large for DeltaPhi LUT (" << lutName <<") size = " << m_lut_DeltaPhi.find(lutName)->second.size() << std::endl;
   }
   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 < (int)m_lut_Pt.find(lutName)->second.size()) {
     value = m_lut_Pt.find(lutName)->second.at(element);
   } else {
     edm::LogError("GlobalScales") << "Warning: Element Requested " << element << " too large for Pt LUT (" << lutName <<") size = " << m_lut_Pt.find(lutName)->second.size() << std::endl;
   }
   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 < (int)m_lut_Cosh.find(lutName)->second.size()) {
     value = m_lut_Cosh.find(lutName)->second.at(element);
   } else {
     edm::LogError("GlobalScales") << "Warning: Element Requested " << element << " too large for Cosh LUT (" << lutName <<") size = " << m_lut_Cosh.find(lutName)->second.size() << std::endl;
   }
   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 < (int)m_lut_Cos.find(lutName)->second.size()) {
     value = m_lut_Cos.find(lutName)->second.at(element);
   } else {
     edm::LogError("GlobalScales") << "Warning: Element Requested " << element << " too large for Cos LUT (" << lutName <<") size = " << m_lut_Cos.find(lutName)->second.size() << std::endl;
   }
   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 < (int)m_lut_Cos.find(lutName)->second.size()) {
     value = m_lut_Cos.find(lutName)->second.at(element);
   } else {
     edm::LogError("GlobalScales") << "Warning: Element Requested " << element << " too large for Cos LUT (" << lutName <<") size = " << m_lut_Cos.find(lutName)->second.size() << std::endl;
   }
   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 < (int)m_lut_Sin.find(lutName)->second.size()) {
     value = m_lut_Sin.find(lutName)->second.at(element);
   } else {
     edm::LogError("GlobalScales") << "Warning: Element Requested " << element << " too large for Sin LUT (" << lutName <<") size = " << m_lut_Sin.find(lutName)->second.size() << std::endl;
   }
   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()
           << "\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()
           << "\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;

}