spectrum_factories.hpp

Go to the documentation of this file.

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

#include <sstream>

#include "gambit/Elements/spectrum.hpp"
#include "gambit/Elements/sminputs.hpp"
#include "gambit/Models/SimpleSpectra/SMSimpleSpec.hpp"
#include "gambit/Models/SimpleSpectra/ScalarSingletDMSimpleSpec.hpp"

namespace Gambit
{
   template <typename HE, typename HEmod>
   Spectrum spectrum_from_SLHAea(HEmod hemod, SLHAstruct slhaea, const Spectrum::mc_info& mci, const Spectrum::mr_info& mri)
   {
     // Start by stripping out any DECAY and DCINFO blocks, as we're only
     // interested in holding spectrum information internally in the
     // spectrum object, not decay info.
     SLHAstruct::key_matches is_dcinfo("DCINFO");
     for (SLHAstruct::iterator block = slhaea.begin(); block != slhaea.end();)
     {
       bool delete_block = false;
       if (is_dcinfo(*block)) delete_block = true;
       else
       {
         auto block_def = block->find_block_def();
         if (block_def != block->end())
         {
           if (block_def->at(0) == "DECAY") delete_block = true;
         }
       }
       if (delete_block) slhaea.erase(block);
       else ++block;
     }

     // Create HE simple SubSpectrum object from the SLHAea object
     // (interacts with MSSM blocks in MSSM case)
     HE he(hemod);

     // Create SMInputs object from the SLHAea object
     SMInputs sminputs(slhaea);

     // Create SMSimpleSpec SubSpectrum object from the SLHAea object
     // (basically just interacts with SMINPUTS block)
     SMSimpleSpec sm(slhaea);

     // Create full Spectrum object from components above
     // Note subtlety! There are TWO constructors for the Spectrum object:
     // If pointers to SubSpectrum objects are passed, it is assumed that
     // these objects are managed EXTERNALLY! So if we were to do this:
     //   matched_spectra = Spectrum(&smskel,&mssmskel,sminputs);
     // then the SubSpectrum objects would end up DELETED at the end of
     // this scope, and we will get a segfault if we try to access them
     // later. INSTEAD, we should just pass the objects themselves, and
     // then they will be CLONED and the Spectrum object will take
     // possession of them:
     return Spectrum(sm,he,sminputs,NULL,mci,mri);
   }

   template <typename HE>
   Spectrum spectrum_from_SLHA(str slha, const Spectrum::mc_info& mci, const Spectrum::mr_info& mri)
   {
     // Read the SLHA file in to an SLHAea object
     SLHAstruct slhaea = read_SLHA(slha);
     // Create the final object from the SLHAea object
     return spectrum_from_SLHAea<HE, SLHAstruct>(slhaea, slhaea, mci, mri);
   }
}