SLHASimpleSpec.cpp

Go to the documentation of this file.

//   GAMBIT: Global and Modular BSM Inference Tool
//   *********************************************
//

#include "gambit/Models/SimpleSpectra/SLHASimpleSpec.hpp"

#include <boost/preprocessor/tuple/to_seq.hpp>
#include <boost/preprocessor/seq/elem.hpp>
#include <boost/preprocessor/seq/for_each_product.hpp>

// Relies on "tmp_map" being used as the variable name for the temporary maps
// inside the fill_map functions.
#define addtomap_EL(r, PRODUCT)                                         \
{                                                                       \
   str key      = BOOST_PP_SEQ_ELEM(0,PRODUCT); /* string map key */    \
   tmp_map[key] = BOOST_PP_SEQ_ELEM(1,PRODUCT); /* function pointer */  \
}

#define addtomap(__KEYS,FPTR) BOOST_PP_SEQ_FOR_EACH_PRODUCT(addtomap_EL, (BOOST_PP_TUPLE_TO_SEQ(__KEYS))((FPTR)) )

using namespace SLHAea;

namespace Gambit
{


      SLHAeaModel::SLHAeaModel()
        : data()
      {}

      SLHAeaModel::SLHAeaModel(const SLHAea::Coll& input)
        : data(input)
      {
         // No idea what kind of model this is, so cannot set valid slha_version. Use "0" for "unknown" or "n/a"
         // If you write a model-specific derived class from this, be sure to add a sensible check for this and
         // overwrite what we set here in the base class.
         wrapped_slha_version = 0;
      }

      int SLHAeaModel::slha_version() const { return wrapped_slha_version; }

      const SLHAea::Coll& SLHAeaModel::get_slhaea() const { return data; }

      const std::map<int, int>& SLHAeaModel::PDG_translator() const
      {
        return PDG_translation_map;
      }


      double SLHAeaModel::getdata(const std::string& block, int index) const
      {
         double output = 0.0;
         try
         {
           output = to<double>(data.at(block).at(index).at(1));
         }
         catch (const std::out_of_range& e)
         {
           std::ostringstream errmsg;
           errmsg << "Error accessing data at index "<<index<<" of block "<<block<<". Please check that the SLHAea object was properly filled." << std::endl;
           errmsg  << "(Received out_of_range error from SLHAea class with message: " << e.what() << ")";
           utils_error().raise(LOCAL_INFO,errmsg.str());
         }
         return output;
      }

      double SLHAeaModel::getdata(const std::string& block, int i, int j) const
      {
         double output = 0.0;
         try
         {
           output = to<double>(data.at(block).at(i,j).at(2));
         }
         catch (const std::out_of_range& e)
         {
           std::ostringstream errmsg;
           errmsg << "Error accessing data at index "<<i<<","<<j<<" of block "<<block<<". Please check that the SLHAea object was properly filled." << std::endl;
           errmsg  << "(Received out_of_range error from SLHAea class with message: " << e.what() << ")";
           utils_error().raise(LOCAL_INFO,errmsg.str());
         }
         return output;
      }

      // Check if block and entry exist, 1 index
      bool SLHAeaModel::checkdata(const std::string& block, int i) const
      {
        if (data.find(block) == data.end()) return false;

        SLHAea::Block::key_type key(1);
        key[0] = std::to_string(i);
        if (data.at(block).find(key) != data.at(block).end())
          return true;

        return false;
      }




} // end Gambit namespace