flexibledecay_settings.cpp

Go to the documentation of this file.

// ====================================================================
// This file is part of FlexibleSUSY.
//
// FlexibleSUSY is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published
// by the Free Software Foundation, either version 3 of the License,
// or (at your option) any later version.
//
// FlexibleSUSY is distributed in the hope that it will be useful, but
// WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with FlexibleSUSY.  If not, see
// <http://www.gnu.org/licenses/>.
// ====================================================================
 
#include "flexibledecay_settings.hpp"
#include "error.hpp"
#include "string_format.hpp"
 
#include <cmath>
#include <iostream>
#include <string>
 
namespace flexiblesusy {
 
namespace {
const std::array<std::string, FlexibleDecay_settings::NUMBER_OF_OPTIONS> descriptions = {
   "calculate particle decays",
   "minimum BR to print",
   "include higher order corrections in decays",
   "use Thomson alpha(0) instead of alpha(m) in decays to γγ and γZ",
   "off-shell decays into VV pair",
   "print loop-induced Higgs couplings",
   "calculate effective couplings normalized to SM",
   "call HiggsTools",
   "call Lilith",
   "use pole Higgs mixings in decays"
};
 
bool is_integer(double value)
{
   double intpart;
   return std::modf(value, &intpart) == 0.0;
}
 
void assert_bool(double value, const char* quantity)
{
   if (value != 0.0 && value != 1.0) {
      throw SetupError(std::string(quantity) + " must either 0 or 1");
   }
}
 
void assert_integer(double value, const char* quantity)
{
   if (!is_integer(value)) {
      throw SetupError(std::string(quantity) + " must be an integer");
   }
}
 
void assert_ge(double value, double lower_bound, const char* quantity)
{
   if (value < lower_bound) {
      throw SetupError(std::string(quantity) +
                       " must be greater than or equal to " +
                       flexiblesusy::to_string(lower_bound));
   }
}
 
void assert_le(double value, double upper_bound, const char* quantity)
{
   if (value > upper_bound) {
      throw SetupError(std::string(quantity) +
                       " must be lower than or equal to " +
                       flexiblesusy::to_string(upper_bound));
   }
}
 
} // anonymous namespace
 
FlexibleDecay_settings::FlexibleDecay_settings()
{
   reset();
}
 
double FlexibleDecay_settings::get(Settings o) const
{
   return values.at(o);
}
 
FlexibleDecay_settings::Settings_t FlexibleDecay_settings::get() const
{
   Settings_t s(&values[0]);
   return s;
}
 
std::string FlexibleDecay_settings::get_description(Settings o) const
{
   return descriptions.at(o);
}
 
void FlexibleDecay_settings::set(Settings o, double value)
{
   switch (o) {
   case calculate_decays: // 0 [bool]
      assert_bool(value, descriptions.at(o).c_str());
      break;
   case min_br_to_print: // 1 [double >= 0 and <= 1]
      assert_ge(value, 0, descriptions.at(o).c_str());
      assert_le(value, 1, descriptions.at(o).c_str());
      break;
   case include_higher_order_corrections: // 2 [int >= 0 and <= 4]
      assert_integer(value, descriptions.at(o).c_str());
      assert_ge(value, 0, descriptions.at(o).c_str());
      assert_le(value, 4, descriptions.at(o).c_str());
      break;
   case use_Thomson_alpha_in_Phigamgam_and_PhigamZ: // 3 [bool]
      assert_bool(value, descriptions.at(o).c_str());
      break;
   case offshell_VV_decays: // 4 [int >= 0 and <= 2]
      assert_integer(value, descriptions.at(o).c_str());
      assert_ge(value, 0, descriptions.at(o).c_str());
      assert_le(value, 2, descriptions.at(o).c_str());
      break;
   case print_effc_block: // 5 [bool]
      assert_bool(value, descriptions.at(o).c_str());
      break;
   case calculate_normalized_effc: // 6 [bool]
      assert_bool(value, descriptions.at(o).c_str());
      break;
   case call_higgstools: // 7 [bool]
      assert_bool(value, descriptions.at(o).c_str());
      break;
   case call_lilith: // 8 [bool]
      assert_bool(value, descriptions.at(o).c_str());
      break;
   case use_pole_higgs_mixings: // 9 [bool]
      assert_bool(value, descriptions.at(o).c_str());
      break;
   default:
      break;
   }
 
   values.at(o) = value;
}
 
void FlexibleDecay_settings::set(const FlexibleDecay_settings::Settings_t& s)
{
   std::copy(s.data(), s.data() + s.size(), values.begin());
}
 
void FlexibleDecay_settings::reset()
{
   values[calculate_decays]                           = 1.0;
   values[min_br_to_print]                            = 1e-5;
   values[include_higher_order_corrections]           = 4.0;
   values[use_Thomson_alpha_in_Phigamgam_and_PhigamZ] = 1.0;
   values[offshell_VV_decays]                         = 2.0;
   values[print_effc_block]                           = 1.0;
   values[calculate_normalized_effc]                  = 0.0;
   values[call_higgstools]                            = 0.0;
   values[call_lilith]                                = 0.0;
   values[use_pole_higgs_mixings]                     = 1.0;
}
 
bool is_integer(double value)
{
   double intpart;
   return std::modf(value, &intpart) == 0.0;
}
 
}