L1ThrExtraInfo.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#include "TrigConfData/L1ThrExtraInfo.h"
 
#include <stdexcept>
#include <cmath>
 
 
using namespace std;
 
std::unique_ptr<TrigConf::L1ThrExtraInfoBase>
TrigConf::L1ThrExtraInfo::createExtraInfo(const std::string & thrTypeName, const boost::property_tree::ptree & data) {
   std::unique_ptr<TrigConf::L1ThrExtraInfoBase> extraInfo(nullptr);
 
   if( thrTypeName == "EM" )
      return std::make_unique<L1ThrExtraInfo_EMTAULegacy>(thrTypeName, data);
 
   if( thrTypeName == "TAU" )
      return std::make_unique<L1ThrExtraInfo_EMTAULegacy>(thrTypeName, data);      
 
   if( thrTypeName == "JET" )
      return std::make_unique<L1ThrExtraInfo_JETLegacy>(thrTypeName, data);
 
   if( thrTypeName == "XS" )
      return std::make_unique<L1ThrExtraInfo_XSLegacy>(thrTypeName, data);      
 
   if( thrTypeName == "MU" )
      return std::make_unique<L1ThrExtraInfo_MU>(thrTypeName, data);
 
   if( thrTypeName == "eEM" )
      return std::make_unique<L1ThrExtraInfo_eEM>(thrTypeName, data);
 
   if( thrTypeName == "jEM" )
      return std::make_unique<L1ThrExtraInfo_jEM>(thrTypeName, data);
 
   if( thrTypeName == "eTAU" )
      return std::make_unique<L1ThrExtraInfo_eTAU>(thrTypeName, data);
 
   if( thrTypeName == "jTAU" )
      return std::make_unique<L1ThrExtraInfo_jTAU>(thrTypeName, data);
 
   if( thrTypeName == "cTAU" )
      return std::make_unique<L1ThrExtraInfo_cTAU>(thrTypeName, data);
 
   if( thrTypeName == "jJ" )
      return std::make_unique<L1ThrExtraInfo_jJ>(thrTypeName, data);      
 
   if( thrTypeName == "jLJ" )
      return std::make_unique<L1ThrExtraInfo_jLJ>(thrTypeName, data);
 
   if( thrTypeName == "gJ" )
      return std::make_unique<L1ThrExtraInfo_gJ>(thrTypeName, data);
 
   if( thrTypeName == "gLJ" )
      return std::make_unique<L1ThrExtraInfo_gLJ>(thrTypeName, data);
 
   if( thrTypeName == "jXE" )
      return std::make_unique<L1ThrExtraInfo_jXE>(thrTypeName, data);
 
   if( thrTypeName == "jTE" )
      return std::make_unique<L1ThrExtraInfo_jTE>(thrTypeName, data);
 
   if( thrTypeName == "gXE" )
      return std::make_unique<L1ThrExtraInfo_gXE>(thrTypeName, data);
 
   if( thrTypeName == "gTE" )
      return std::make_unique<L1ThrExtraInfo_gTE>(thrTypeName, data);
   
   if( thrTypeName == "cXE" )
      return std::make_unique<L1ThrExtraInfo_cXE>(thrTypeName, data);
   
   // if no special extra information is supplied for the threshold type return base class
   return std::make_unique<L1ThrExtraInfoBase>(thrTypeName, data);
}
 
void
TrigConf::L1ThrExtraInfo::clear()
{
   m_thrExtraInfo.clear();
}
 
 
std::weak_ptr<TrigConf::L1ThrExtraInfoBase>
TrigConf::L1ThrExtraInfo::addExtraInfo(const std::string & thrTypeName, const boost::property_tree::ptree & data) {
   try {
      if( auto extraInfo = L1ThrExtraInfo::createExtraInfo( thrTypeName, data) ) {
         auto success = m_thrExtraInfo.emplace(thrTypeName, std::shared_ptr<TrigConf::L1ThrExtraInfoBase>(std::move(extraInfo)));
         return std::weak_ptr<TrigConf::L1ThrExtraInfoBase>( success.first->second );
      }
   }
   catch(std::exception & ex) {
      std::cerr << "L1ThrExtraInfo::addExtraInfo: exception occured when building extra info for " << thrTypeName << std::endl;
      throw;
   }
   return std::weak_ptr<TrigConf::L1ThrExtraInfoBase>( m_emptyInfo );
}
 
bool
TrigConf::L1ThrExtraInfo::hasInfo(const std::string & typeName) const
{
   return ( m_thrExtraInfo.find(typeName) != m_thrExtraInfo.end() );
}
 
const TrigConf::L1ThrExtraInfo_EMTAULegacy &
TrigConf::L1ThrExtraInfo::EM() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_EMTAULegacy&>( * m_thrExtraInfo.at("EM") );
}
 
const TrigConf::L1ThrExtraInfo_EMTAULegacy &
TrigConf::L1ThrExtraInfo::TAU() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_EMTAULegacy&>( * m_thrExtraInfo.at("TAU") );
}
 
const TrigConf::L1ThrExtraInfo_XSLegacy &
TrigConf::L1ThrExtraInfo::XS() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_XSLegacy&>( * m_thrExtraInfo.at("XS") );
}
 
const TrigConf::L1ThrExtraInfo_JETLegacy &
TrigConf::L1ThrExtraInfo::JET() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_JETLegacy&>( * m_thrExtraInfo.at("JET") );
}
 
const TrigConf::L1ThrExtraInfo_eEM &
TrigConf::L1ThrExtraInfo::eEM() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_eEM&>( * m_thrExtraInfo.at("eEM") );
}
 
const TrigConf::L1ThrExtraInfo_jEM &
TrigConf::L1ThrExtraInfo::jEM() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_jEM&>( * m_thrExtraInfo.at("jEM") );
}
 
const TrigConf::L1ThrExtraInfo_eTAU &
TrigConf::L1ThrExtraInfo::eTAU() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_eTAU&>( * m_thrExtraInfo.at("eTAU") );
}
 
const TrigConf::L1ThrExtraInfo_jTAU &
TrigConf::L1ThrExtraInfo::jTAU() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_jTAU&>( * m_thrExtraInfo.at("jTAU") );
}
 
const TrigConf::L1ThrExtraInfo_cTAU &
TrigConf::L1ThrExtraInfo::cTAU() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_cTAU&>( * m_thrExtraInfo.at("cTAU") );
}
 
const TrigConf::L1ThrExtraInfo_jJ &
TrigConf::L1ThrExtraInfo::jJ() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_jJ&>( * m_thrExtraInfo.at("jJ") );
}
 
const TrigConf::L1ThrExtraInfo_jLJ &
TrigConf::L1ThrExtraInfo::jLJ() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_jLJ&>( * m_thrExtraInfo.at("jLJ") );
}
 
const TrigConf::L1ThrExtraInfo_gJ &
TrigConf::L1ThrExtraInfo::gJ() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_gJ&>( * m_thrExtraInfo.at("gJ") );
}
 
const TrigConf::L1ThrExtraInfo_gLJ &
TrigConf::L1ThrExtraInfo::gLJ() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_gLJ&>( * m_thrExtraInfo.at("gLJ") );
}
 
const TrigConf::L1ThrExtraInfo_jXE &
TrigConf::L1ThrExtraInfo::jXE() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_jXE&>( * m_thrExtraInfo.at("jXE") );
}
 
const TrigConf::L1ThrExtraInfo_jTE &
TrigConf::L1ThrExtraInfo::jTE() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_jTE&>( * m_thrExtraInfo.at("jTE") );
}
 
const TrigConf::L1ThrExtraInfo_gXE &
TrigConf::L1ThrExtraInfo::gXE() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_gXE&>( * m_thrExtraInfo.at("gXE") );
}
 
const TrigConf::L1ThrExtraInfo_gTE &
TrigConf::L1ThrExtraInfo::gTE() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_gTE&>( * m_thrExtraInfo.at("gTE") );
}
 
const TrigConf::L1ThrExtraInfo_cXE &
TrigConf::L1ThrExtraInfo::cXE() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_cXE&>( * m_thrExtraInfo.at("cXE") );
}
 
const TrigConf::L1ThrExtraInfo_MU &
TrigConf::L1ThrExtraInfo::MU() const {
   return dynamic_cast<const TrigConf::L1ThrExtraInfo_MU&>( * m_thrExtraInfo.at("MU") );
}
 
const TrigConf::L1ThrExtraInfoBase &
TrigConf::L1ThrExtraInfo::thrExtraInfo(const std::string & thrTypeName) const
{
   try {
      return * m_thrExtraInfo.at(thrTypeName);
   }
   catch(std::exception & ex) {
      std::cerr << "Threshold type " << thrTypeName << " does not have extra info defined" << endl;
      throw;
   }
}

const TrigConf::IsolationLegacy &
TrigConf::L1ThrExtraInfo_EMTAULegacy::isolation(const std::string & thrType, size_t bit) const
{
   if(bit<1 or bit>5) {
      throw std::out_of_range("When accessing the legacy L1Calo EM or TAU isolation bit must be between 1 and 5, but bit=" 
                              + std::to_string(bit) + " was requested");
   }
   try {
      return m_isolation.at(thrType)[bit-1];
   }
   catch(std::exception & ex) {
      std::cerr << "Threshold type " << name() << " does not have isolation parameters for type " << thrType << endl;
      throw;
   }
}
 
void
TrigConf::L1ThrExtraInfo_EMTAULegacy::load()
{
   for( auto & x : m_extraInfo ) {
      if( x.first == "ptMinToTopo" ) {
         m_ptMinToTopoMeV = std::lround( 1000 * x.second.getValue<float>() );
      } else if( x.first == "isolation" ) {
         for( auto & y : x.second.data() ) {
            auto & isoV = m_isolation[y.first] = std::vector<IsolationLegacy>(5);
            for(auto & c : y.second.get_child("Parametrization") ) {
               auto iso = IsolationLegacy(c.second);
               isoV[iso.isobit()-1] = iso;
            }
         }
      }
   }
}
 
 
void
TrigConf::L1ThrExtraInfo_XSLegacy::load()
{
   if( hasExtraInfo("significance") ) {
      auto & sig = m_extraInfo["significance"];
      m_xeMin = sig.getAttribute<unsigned int>("xeMin");
      m_xeMax = sig.getAttribute<unsigned int>("xeMax");
      m_teSqrtMin = sig.getAttribute<unsigned int>("teSqrtMin");
      m_teSqrtMax = sig.getAttribute<unsigned int>("teSqrtMax");
      m_xsSigmaScale = sig.getAttribute<unsigned int>("xsSigmaScale");
      m_xsSigmaOffset = sig.getAttribute<unsigned int>("xsSigmaOffset");
   }
}
 
 

void
TrigConf::L1ThrExtraInfo_JETLegacy::load()
{
   for( auto & x : m_extraInfo ) {
      if( x.first == "ptMinToTopoLargeWindow" ) {
         m_ptMinToTopoLargeWindowMeV = std::lround( 1000 * x.second.getValue<float>() );
      } else if( x.first == "ptMinToTopoSmallWindow" ) {
         m_ptMinToTopoSmallWindowMeV = std::lround( 1000 * x.second.getValue<float>() );
      }
   }
}
 
 
/***********************************
 *
 * Extra info for new thresholds
 *
 ***********************************/
 
/*******
 * eEM
 *******/
TrigConf::L1ThrExtraInfo_eEM::WorkingPoints_eEM::WorkingPoints_eEM( const boost::property_tree::ptree & pt) {
   m_isDefined = true;
   m_reta_d  = pt.get_optional<float>("reta").get_value_or(0);
   m_wstot_d = pt.get_optional<float>("wstot").get_value_or(0);
   m_rhad_d  = pt.get_optional<float>("rhad").get_value_or(0);
   m_reta_fw  = pt.get_optional<int>("reta_fw").get_value_or(0);
   m_wstot_fw = pt.get_optional<int>("wstot_fw").get_value_or(0);
   m_rhad_fw  = pt.get_optional<int>("rhad_fw").get_value_or(0);
}
 
std::ostream &
TrigConf::operator<<(std::ostream & os, const TrigConf::L1ThrExtraInfo_eEM::WorkingPoints_eEM & iso) {
   os << "reta_fw=" << iso.reta_fw() << ", wstot_fw=" << iso.wstot_fw() << ", rhad_fw=" << iso.rhad_fw();
   return os;
}
 
void
TrigConf::L1ThrExtraInfo_eEM::load()
{
   for( auto & x : m_extraInfo ) {
      if( x.first == "maxEt" ){
         m_maxEt = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinToTopo" ) {
         m_ptMinToTopoMeV = lround(1000 * x.second.getValue<float>());
      } else if( x.first == "workingPoints" ) {
         for( auto & y : x.second.data() ) {
            auto wp = Selection::stringToWP(y.first);
            auto & iso = m_isolation.emplace(wp, string("eEM_WP_" + y.first)).first->second;
            for(auto & c : y.second ) {
               int etamin = c.second.get_optional<int>("etamin").get_value_or(-49);
               int etamax = c.second.get_optional<int>("etamax").get_value_or(49);
               unsigned int priority = c.second.get_optional<unsigned int>("priority").get_value_or(0);
               iso.addRangeValue(WorkingPoints_eEM(c.second), etamin, etamax, priority, /*symmetric=*/ false);
            }
         }
      } else if (x.first == "algoVersion") {
          m_algoVersion = x.second.getValue<unsigned int>();
      }
   }
}
 
/*******
 * jEM
 *******/
TrigConf::L1ThrExtraInfo_jEM::WorkingPoints_jEM::WorkingPoints_jEM( const boost::property_tree::ptree & pt ) {
   m_isDefined = true;
   m_iso_d  = pt.get_optional<float>("iso").get_value_or(0);
   m_frac_d = pt.get_optional<float>("frac").get_value_or(0);
   m_frac2_d  = pt.get_optional<float>("frac2").get_value_or(0);
   m_iso_fw  = pt.get_optional<int>("iso_fw").get_value_or(0);
   m_frac_fw = pt.get_optional<int>("frac_fw").get_value_or(0);
   m_frac2_fw  = pt.get_optional<int>("frac2_fw").get_value_or(0);
}
 
std::ostream &
TrigConf::operator<<(std::ostream & os, const TrigConf::L1ThrExtraInfo_jEM::WorkingPoints_jEM & iso) {
   os << "iso_fw=" << iso.iso_fw() << ", frac_fw=" << iso.frac_fw() << ", frac2_fw=" << iso.frac2_fw();
   return os;
}
 
void
TrigConf::L1ThrExtraInfo_jEM::load()
{
   for( auto & x : m_extraInfo ) {
      if( x.first == "maxEt" ){
         m_maxEt = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinToTopo1" ) {
         m_ptMinToTopoMeV1 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinToTopo2" ){
         m_ptMinToTopoMeV2 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinToTopo3" ){
         m_ptMinToTopoMeV3 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinxTOB1" ){
         m_ptMinxTOBMeV1 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinxTOB2" ){
         m_ptMinxTOBMeV2 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinxTOB3" ){
         m_ptMinxTOBMeV3 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "workingPoints" ) {
         for( auto & y : x.second.data() ) {
            auto wp = Selection::stringToWP(y.first);
            auto & iso = m_isolation.emplace(wp, string("jEM_WP_" + y.first)).first->second;
            for(auto & c : y.second ) {
               int etamin = c.second.get_optional<int>("etamin").get_value_or(-49);
               int etamax = c.second.get_optional<int>("etamax").get_value_or(49);
               unsigned int priority = c.second.get_optional<unsigned int>("priority").get_value_or(0);
               iso.addRangeValue(WorkingPoints_jEM(c.second), etamin, etamax, priority, /*symmetric=*/ false);
            }
         }
      }
   }
}
 
/*******
 * eTAU
 *******/
TrigConf::L1ThrExtraInfo_eTAU::WorkingPoints_eTAU::WorkingPoints_eTAU( const boost::property_tree::ptree & pt ) {
   m_isDefined = true;
   m_rCore_d   = pt.get_optional<float>("rCore").get_value_or(0);
   m_rHad_d    = pt.get_optional<float>("rHad").get_value_or(0);
   m_rCore_fw  = pt.get_optional<float>("rCore_fw").get_value_or(0);
   m_rHad_fw   = pt.get_optional<float>("rHad_fw").get_value_or(0);
}
 
std::ostream &
TrigConf::operator<<(std::ostream & os, const TrigConf::L1ThrExtraInfo_eTAU::WorkingPoints_eTAU & iso) {
   os << "rCore_fw=" << iso.rCore_fw() << ", rHad_fw=" << iso.rHad_fw() ;
   return os;
}
 
void
TrigConf::L1ThrExtraInfo_eTAU::load()
{
   for( auto & x : m_extraInfo ) {
      if( x.first == "maxEt" ){
         m_maxEt = 1000*x.second.getValue<unsigned int>(); // Original in units of GeV
      } else if( x.first == "minIsoEt" ){
         m_minIsoEt = lround(1000 * x.second.getValue<float>()); // Original in units of GeV
      } else if( x.first == "ptMinToTopo" ) {
         m_ptMinToTopoMeV = lround(1000 * x.second.getValue<float>());
      } else if( x.first == "workingPoints" ) {
         for( auto & y : x.second.data() ) {
            auto wp = TrigConf::Selection::stringToWP(y.first);
            auto & iso = m_isolation.emplace(wp, string("eTAU_WP_" + y.first)).first->second;
            for(auto & c : y.second ) {
               int etamin = c.second.get_optional<int>("etamin").get_value_or(-49);
               int etamax = c.second.get_optional<int>("etamax").get_value_or(49);
               unsigned int priority = c.second.get_optional<unsigned int>("priority").get_value_or(0);
               iso.addRangeValue(WorkingPoints_eTAU(c.second), etamin, etamax, priority, /*symmetric=*/ false);
            }
         }
      } else if (x.first == "algoVersion") {
          m_algoVersion = x.second.getValue<unsigned int>();
      }
   }
}
 
/*******
 * jTAU
 *******/
TrigConf::L1ThrExtraInfo_jTAU::WorkingPoints_jTAU::WorkingPoints_jTAU( const boost::property_tree::ptree & pt ) {
   m_isDefined = true;
   m_isolation_d    = pt.get_optional<float>("isolation").get_value_or(0);
   m_isolation_fw   = pt.get_optional<float>("isolation_fw").get_value_or(0);
}
 
std::ostream &
TrigConf::operator<<(std::ostream & os, const TrigConf::L1ThrExtraInfo_jTAU::WorkingPoints_jTAU & iso) {
   os << "isolation_fw=" << iso.isolation_fw() ;
   return os;
}
 
void
TrigConf::L1ThrExtraInfo_jTAU::load()
{
   for( auto & x : m_extraInfo ) {
      if( x.first == "maxEt" ){
         m_maxEt = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinToTopo1" ) {
         m_ptMinToTopoMeV1 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinToTopo2" ){
         m_ptMinToTopoMeV2 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinToTopo3" ){
         m_ptMinToTopoMeV3 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinxTOB1" ){
         m_ptMinxTOBMeV1 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinxTOB2" ){
         m_ptMinxTOBMeV2 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinxTOB3" ){
         m_ptMinxTOBMeV3 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "workingPoints" ) {
         for( auto & y : x.second.data() ) {
            auto wp = TrigConf::Selection::stringToWP(y.first);
            auto & iso = m_isolation.emplace(wp, string("jTAU_WP_" + y.first)).first->second;
            for(auto & c : y.second ) {
               int etamin = c.second.get_optional<int>("etamin").get_value_or(-49);
               int etamax = c.second.get_optional<int>("etamax").get_value_or(49);
               unsigned int priority = c.second.get_optional<unsigned int>("priority").get_value_or(0);
               iso.addRangeValue(WorkingPoints_jTAU(c.second), etamin, etamax, priority, /*symmetric=*/ false);
            }
         }
      }
   }
}
 
/*******
 * cTAU
 *******/
TrigConf::L1ThrExtraInfo_cTAU::WorkingPoints_cTAU::WorkingPoints_cTAU( const boost::property_tree::ptree & pt ) {
   m_isDefined = true;
   m_isolation_d    = pt.get_optional<float>("isolation").get_value_or(0);
   m_isolation_fw   = pt.get_optional<unsigned int>("isolation_fw").get_value_or(0);
   m_isolation_jTAUCoreScale_d    = pt.get_optional<float>("isolation_jTAUCoreScale").get_value_or(0);
   m_isolation_jTAUCoreScale_fw   = pt.get_optional<unsigned int>("isolation_jTAUCoreScale_fw").get_value_or(0);
   m_eTAU_rCoreMin_WP_d    = pt.get_optional<float>("eTAU_rCoreMin").get_value_or(0);
   m_eTAU_rCoreMin_WP_fw   = pt.get_optional<unsigned int>("eTAU_rCoreMin_WP_fw").get_value_or(0);
   m_eTAU_rHadMin_WP_d    = pt.get_optional<float>("eTAU_rHadMin").get_value_or(0);
   m_eTAU_rHadMin_WP_fw   = pt.get_optional<unsigned int>("eTAU_rHadMin_WP_fw").get_value_or(0);
}
 
std::ostream &
TrigConf::operator<<(std::ostream & os, const TrigConf::L1ThrExtraInfo_cTAU::WorkingPoints_cTAU & iso) {
   os << "isolation_fw=" << iso.isolation_fw() << ", isolation_jTAUCoreScale_fw=" << iso.isolation_jTAUCoreScale_fw();
   os << ", eTAU_rCoreMin_WP_fw=" << iso.eTAU_rCoreMin_WP_fw() << ", eTAU_rHadMin_WP_fw=" << iso.eTAU_rHadMin_WP_fw();
   return os;
}
 
void
TrigConf::L1ThrExtraInfo_cTAU::load()
{
   for( auto & x : m_extraInfo ) {
      if( x.first == "workingPoints" ) {
         for( auto & y : x.second.data() ) {
            auto wp = TrigConf::Selection::stringToWP(y.first);
            auto & iso = m_isolation.emplace(wp, string("cTAU_WP_" + y.first)).first->second;
            for(auto & c : y.second ) {
               int etamin = c.second.get_optional<int>("etamin").get_value_or(-49);
               int etamax = c.second.get_optional<int>("etamax").get_value_or(49);
               unsigned int priority = c.second.get_optional<unsigned int>("priority").get_value_or(0);
               iso.addRangeValue(WorkingPoints_cTAU(c.second), etamin, etamax, priority, /*symmetric=*/ false);
            }
         }
      }
   }
}
 
 
/*******
 * jJ
 *******/
void
TrigConf::L1ThrExtraInfo_jJ::load()
{
   for( auto & x : m_extraInfo ) {
      if( x.first == "ptMinToTopo1" ) {
         m_ptMinToTopoMeV1 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinToTopo2" ){
         m_ptMinToTopoMeV2 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinToTopo3" ){
         m_ptMinToTopoMeV3 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinxTOB1" ){
         m_ptMinxTOBMeV1 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinxTOB2" ){
         m_ptMinxTOBMeV2 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinxTOB3" ){
         m_ptMinxTOBMeV3 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "seedThreshold1" ){
         m_seedThresholdMeV1 = int(1000*x.second.getValue<float>());
      } else if( x.first == "seedThreshold2" ){
         m_seedThresholdMeV2 = int(1000*x.second.getValue<float>());
      } else if( x.first == "seedThreshold3" ){
         m_seedThresholdMeV3 = int(1000*x.second.getValue<float>());
      }
   }
}
 
/*******
 * jLJ
 *******/
void
TrigConf::L1ThrExtraInfo_jLJ::load()
{
   for( auto & x : m_extraInfo ) {
      if( x.first == "ptMinToTopo1" ) {
         m_ptMinToTopoMeV1 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinToTopo2" ){
         m_ptMinToTopoMeV2 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinToTopo3" ){
         m_ptMinToTopoMeV3 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinxTOB1" ){
         m_ptMinxTOBMeV1 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinxTOB2" ){
         m_ptMinxTOBMeV2 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinxTOB3" ){
         m_ptMinxTOBMeV3 = 1000*x.second.getValue<unsigned int>();
      } 
   }
}
 
/*******
 * gJ
 *******/
void
TrigConf::L1ThrExtraInfo_gJ::load()
{
   for( auto & x : m_extraInfo ) {
      if( x.first == "ptMinToTopo1" ) {
         m_ptMinToTopoMeV1 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinToTopo2" ){
         m_ptMinToTopoMeV2 = 1000*x.second.getValue<unsigned int>();
      }
   }
}
 
/*******
 * gLJ
 *******/
void
TrigConf::L1ThrExtraInfo_gLJ::load()
{
   for( auto & x : m_extraInfo ) {
      if( x.first == "ptMinToTopo1" ) {
         m_ptMinToTopoMeV1 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "ptMinToTopo2" ){
         m_ptMinToTopoMeV2 = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "seedThrA" ){
         m_seedThrMeVA = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "seedThrB" ){
         m_seedThrMeVB = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "seedThrC" ){
         m_seedThrMeVC = 1000*x.second.getValue<unsigned int>();
      } else if( x.first == "rhoTowerMinA" ){
         float rhoTower_tmp = 1000*x.second.getValue<float>();
         if( (int)rhoTower_tmp != rhoTower_tmp)
             throw std::runtime_error("gLJ: rhoTower param " + std::to_string(rhoTower_tmp/1000.) + " cannot be converted in MeV" ); 
         m_rhoTowerMinMeVA = (int)rhoTower_tmp;        
      } else if( x.first == "rhoTowerMinB" ){
         float rhoTower_tmp = 1000*x.second.getValue<float>();
         if( (int)rhoTower_tmp != rhoTower_tmp)
             throw std::runtime_error("gLJ: rhoTower param " + std::to_string(rhoTower_tmp/1000.) + " cannot be converted in MeV" );  
         m_rhoTowerMinMeVB = (int)rhoTower_tmp; 
      } else if( x.first == "rhoTowerMinC" ){
         float rhoTower_tmp = 1000*x.second.getValue<float>();
         if( (int)rhoTower_tmp != rhoTower_tmp)
             throw std::runtime_error("gLJ: rhoTower param " + std::to_string(rhoTower_tmp/1000.) + " cannot be converted in MeV" );  
         m_rhoTowerMinMeVC = (int)rhoTower_tmp; 
      } else if( x.first == "rhoTowerMaxA" ){
         float rhoTower_tmp = 1000*x.second.getValue<float>();
         if( (int)rhoTower_tmp != rhoTower_tmp)
             throw std::runtime_error("gLJ: rhoTower param " + std::to_string(rhoTower_tmp/1000.) + " cannot be converted in MeV" );  
         m_rhoTowerMaxMeVA = (int)rhoTower_tmp; 
      } else if( x.first == "rhoTowerMaxB" ){
         float rhoTower_tmp = 1000*x.second.getValue<float>();
         if( (int)rhoTower_tmp != rhoTower_tmp)
             throw std::runtime_error("gLJ: rhoTower param " + std::to_string(rhoTower_tmp/1000.) + " cannot be converted in MeV" );
         m_rhoTowerMaxMeVB = (int)rhoTower_tmp;
      } else if( x.first == "rhoTowerMaxC" ){
         float rhoTower_tmp = 1000*x.second.getValue<float>();
         if( (int)rhoTower_tmp != rhoTower_tmp)
             throw std::runtime_error("gLJ: rhoTower param " + std::to_string(rhoTower_tmp/1000.) + " cannot be converted in MeV" );
         m_rhoTowerMaxMeVC = (int)rhoTower_tmp;
      }
   }  
}
 
/*******
 * jXE
 *******/
void
TrigConf::L1ThrExtraInfo_jXE::load()
{
}
 
/*******
 * jTE
 *******/
void
TrigConf::L1ThrExtraInfo_jTE::load()
{
   for( auto & x : m_extraInfo ) {
      if( x.first == "etaBoundary1" ) {
         m_etaBoundary1 = x.second.getValue<unsigned int>();
      } else if( x.first == "etaBoundary1_fw" ) {
         m_etaBoundary1_fw = x.second.getValue<unsigned int>();
      } else if( x.first == "etaBoundary2" ) {
         m_etaBoundary2 = x.second.getValue<unsigned int>();
      } else if( x.first == "etaBoundary2_fw" ) {
         m_etaBoundary2_fw = x.second.getValue<unsigned int>();
      } else if( x.first == "etaBoundary3" ) {
         m_etaBoundary3 = x.second.getValue<unsigned int>();
      } else if( x.first == "etaBoundary3_fw" ) {
         m_etaBoundary3_fw = x.second.getValue<unsigned int>();
      } 
   }
}
 
/*******
 * gXE
 *******/
void
TrigConf::L1ThrExtraInfo_gXE::load()
{
   for( auto & x : m_extraInfo ) {
      if( x.first == "seedThrA" ){
         m_seedThrMeVA = 200*x.second.getValue<unsigned int>();
      } else if( x.first == "seedThrB" ){
         m_seedThrMeVB = 200*x.second.getValue<unsigned int>();
      } else if( x.first == "seedThrC" ){
         m_seedThrMeVC = 200*x.second.getValue<unsigned int>();
      } else if( x.first == "XERHO_sigmaPosA" ){
         m_XERHO_sigmaPosA = x.second.getValue<unsigned int>();
      } else if( x.first == "XERHO_sigmaPosB" ){
         m_XERHO_sigmaPosB = x.second.getValue<unsigned int>();
      } else if( x.first == "XERHO_sigmaPosC" ){
         m_XERHO_sigmaPosC = x.second.getValue<unsigned int>();
      } else if( x.first == "XERHO_sigmaNegA" ){
         m_XERHO_sigmaNegA = x.second.getValue<unsigned int>();
      } else if( x.first == "XERHO_sigmaNegB" ){
         m_XERHO_sigmaNegB = x.second.getValue<unsigned int>();
      } else if( x.first == "XERHO_sigmaNegC" ){
         m_XERHO_sigmaNegC = x.second.getValue<unsigned int>();
      } else if( x.first == "XEJWOJ_a_A" ){
         m_XEJWOJ_a_A = x.second.getValue<unsigned int>();
      } else if( x.first == "XEJWOJ_a_B" ){
         m_XEJWOJ_a_B = x.second.getValue<unsigned int>();
      } else if( x.first == "XEJWOJ_a_C" ){
         m_XEJWOJ_a_C = x.second.getValue<unsigned int>();
      } else if( x.first == "XEJWOJ_b_A" ){
         m_XEJWOJ_b_A = x.second.getValue<unsigned int>();
      } else if( x.first == "XEJWOJ_b_B" ){
         m_XEJWOJ_b_B = x.second.getValue<unsigned int>();
      } else if( x.first == "XEJWOJ_b_C" ){
         m_XEJWOJ_b_C = x.second.getValue<unsigned int>();
      } else if( x.first == "XEJWOJ_c_A" ){
         m_XEJWOJ_c_A = x.second.getValue<unsigned int>();
      } else if( x.first == "XEJWOJ_c_B" ){
         m_XEJWOJ_c_B = x.second.getValue<unsigned int>();
      } else if( x.first == "XEJWOJ_c_C" ){
         m_XEJWOJ_c_C = x.second.getValue<unsigned int>();
      }
   }
}
 
/*******
 * gTE
 *******/
void
TrigConf::L1ThrExtraInfo_gTE::load()
{
}
 
/*******
 * cXE
 *******/
void
TrigConf::L1ThrExtraInfo_cXE::load()
{
   for( auto & x : m_extraInfo ) {
      if( x.first == "jXeWeight" ){
         m_jXeWeight = x.second.getValue<float>();
      } else if( x.first == "gXeWeight" ){
         m_gXeWeight = x.second.getValue<float>();
      }
   }
}
 
/*******
 * MU
 *******/
unsigned int
TrigConf::L1ThrExtraInfo_MU::rpcIdxForPt(unsigned int pt) const
{
   try {
      return m_rpcPtMap.at(pt);
   }
   catch(std::exception & ex) {
      std::cerr << "No RPC index defined for pt " << pt << endl;
      throw;
   }   
}
 
 
unsigned int
TrigConf::L1ThrExtraInfo_MU::tgcIdxForPt(unsigned int pt) const
{
   try {
      return m_tgcPtMap.at(pt);
   }
   catch(std::exception & ex) {
      std::cerr << "No TGC index defined for pt " << pt << endl;
      throw;
   }
}
 
unsigned int
TrigConf::L1ThrExtraInfo_MU::ptForRpcIdx(unsigned int idx) const
{
 
   for(auto & x : m_rpcPtMap){
      if(x.second==idx) return x.first;
   } 
   throw std::runtime_error("index "+std::to_string(idx)+" not found for RPC roads"); 
 
}
 
unsigned int
TrigConf::L1ThrExtraInfo_MU::ptForTgcIdx(unsigned int idx) const
{
 
   for(auto & x : m_tgcPtMap){
      if(x.second==idx) return x.first;
   }
   throw std::runtime_error("index "+std::to_string(idx)+" not found for TGC roads");
 
}
 
unsigned int
TrigConf::L1ThrExtraInfo_MU::tgcIdxForRpcIdx(unsigned int rpcIdx) const
{
 
   int ptValue = ptForRpcIdx(rpcIdx);
   return tgcIdxForPt(ptValue);
}
 
std::vector<unsigned int>
TrigConf::L1ThrExtraInfo_MU::knownRpcPtValues() const
{
   std::vector<unsigned int> ptValues;   
   for( auto & x : m_rpcPtMap ) {
      ptValues.emplace_back(x.first);
   }
   return ptValues;
}
 
std::vector<unsigned int>
TrigConf::L1ThrExtraInfo_MU::knownTgcPtValues() const
{
   std::vector<unsigned int> ptValues;   
   for( auto & x : m_tgcPtMap ) {
      ptValues.emplace_back(x.first);
   }
   return ptValues;
}
 
 
std::vector<std::string>
TrigConf::L1ThrExtraInfo_MU::exclusionListNames() const
{
   std::vector<std::string> listNames;   
   for( auto & x : m_roiExclusionLists ) {
      listNames.emplace_back(x.first);
   }
   return listNames;
}
 
 
const std::map<std::string, std::vector<unsigned int> > &
TrigConf::L1ThrExtraInfo_MU::exclusionList(const std::string & listName) const
{
   try {
      return m_roiExclusionLists.at(listName);
   }
   catch(std::exception & ex) {
      std::cerr << "No exclusion list '" << listName << "' defined in MU threshold exlusionLists" << endl;
      throw;
   }
}
 
void
TrigConf::L1ThrExtraInfo_MU::load()
{
   {
      DataStructure ds = m_extraInfo["roads"].getObject("rpc");
      for( const auto & x : ds.data() ) {
        m_rpcPtMap.emplace( static_cast<unsigned int>(std::stoul(x.first)),
                            static_cast<unsigned int>(std::stoul(x.second.data())));
     }
   }
   {
      DataStructure ds = m_extraInfo["roads"].getObject("tgc");
      for( auto & x : ds.data() ) {
         m_tgcPtMap.emplace( static_cast<unsigned int>(std::stoul(x.first)),
                             static_cast<unsigned int>(std::stoul(x.second.data())));
      }
   }
   for( auto & x : m_extraInfo["exclusionLists"].data() ) {
      const std::string & listName = x.first;
      std::map<std::string, std::vector<unsigned int>> roisBySector;
      for( auto & list : x.second ) {
         const std::string & sectorName = list.second.get_child("sectorName").get_value<std::string>();
         std::vector<unsigned int> rois;
         for( auto & roi : list.second.get_child("rois") ) {
            rois.push_back( static_cast<unsigned int>(std::stoul( roi.second.data() )) );
         }
         roisBySector.emplace(sectorName, std::move(rois));
      }
      m_roiExclusionLists.emplace(listName, std::move(roisBySector));
   }
}