GAMBIT native higgs coupling table class. More...

#include <higgs_couplings_table.hpp>

Collaboration diagram for Gambit::HiggsCouplingsTable:

Public Types

typedef const DecayTable::Entry *	h0_decay_array_type[max_neutral_higgses]
	Types to make returning decay arrays easier. More...

typedef const DecayTable::Entry *	hp_decay_array_type[max_charged_higgses]

Public Member Functions
	HiggsCouplingsTable ()
	Constructor. More...


template<typename T >
double	compute_effective_coupling (int index, const T &p1, const T &p2)
	Compute a neutral higgs effective coupling from the current two-body neutral higgs decays. More...

template<typename T >
double	compute_effective_coupling (str name, const T &p1, const T &p2)


void	set_neutral_decays_SM (int, const str &, const DecayTable::Entry &)
	Assign decay entries to the various table components. More...

void	set_neutral_decays (int, const str &, const DecayTable::Entry &)
	Assign decay entries to neutral higgses. More...

void	set_charged_decays (int, const str &, const DecayTable::Entry &)
	Assign decay entries to charged higgses. More...

void	set_t_decays (const DecayTable::Entry &)
	Assign decay entries to top. More...


const h0_decay_array_type &	get_neutral_decays_SM_array (int) const
	Retrieve decay sets. More...

const DecayTable::Entry &	get_neutral_decays_SM (int) const
	Retrieve SM decays of a specific neutral Higgs, by index. More...

const DecayTable::Entry &	get_neutral_decays_SM (const str &) const
	Retrieve SM decays of a specific neutral Higgs, by name. More...

const h0_decay_array_type &	get_neutral_decays_array (int) const
	Retrieve decays of all neutral higgses. More...

const DecayTable::Entry &	get_neutral_decays (int) const
	Retrieve decays of a specific neutral Higgs, by index. More...

const DecayTable::Entry &	get_neutral_decays (const str &) const
	Retrieve decays of a specific neutral Higgs, by name. More...

const hp_decay_array_type &	get_charged_decays_array (int) const
	Retrieve decays of all charged higgses. More...

const DecayTable::Entry &	get_charged_decays (int) const
	Retrieve decays of a specific charged Higgs, by index. More...

const DecayTable::Entry &	get_charged_decays (const str &) const
	Retrieve decays of a specific charged Higgs, by name. More...

const DecayTable::Entry &	get_t_decays () const
	Retrieve decays of the top quark. More...

Public Attributes
double	CP [max_neutral_higgses]
	CP of neutral higgses. More...

std::vector< str >	invisibles
	Particles that higgses can decay invisibly to. More...


double	C_WW2 [max_neutral_higgses]
	Effective couplings for neutral higgses. More...

double	C_ZZ2 [max_neutral_higgses]

double	C_tt2 [max_neutral_higgses]

double	C_bb2 [max_neutral_higgses]

double	C_cc2 [max_neutral_higgses]

double	C_tautau2 [max_neutral_higgses]

double	C_gaga2 [max_neutral_higgses]

double	C_gg2 [max_neutral_higgses]

double	C_mumu2 [max_neutral_higgses]

double	C_Zga2 [max_neutral_higgses]

double	C_ss2 [max_neutral_higgses]

double	C_hiZ2 [max_neutral_higgses][max_neutral_higgses]

Static Public Attributes

static const int	max_neutral_higgses = 3
	The maximal Higgs sector manageable by this class; update as needed. More...

static const int	max_charged_higgses = 1

Private Attributes
const DecayTable::Entry *	t_decays
	t decays (as t can decay to h) More...


const DecayTable::Entry *	neutral_decays_SM_array [max_neutral_higgses]
	SM neutral higgs decays. More...

std::map< str, const DecayTable::Entry & >	neutral_decays_SM_map


const DecayTable::Entry *	neutral_decays_array [max_neutral_higgses]
	BSM higgs decays. More...

const DecayTable::Entry *	charged_decays_array [max_charged_higgses]

std::map< str, const DecayTable::Entry & >	neutral_decays_map

std::map< str, const DecayTable::Entry & >	charged_decays_map

Detailed Description

GAMBIT native higgs coupling table class.

Definition at line 28 of file higgs_couplings_table.hpp.

Member Typedef Documentation

◆ h0_decay_array_type

typedef const DecayTable::Entry* Gambit::HiggsCouplingsTable::h0_decay_array_type[max_neutral_higgses]

Types to make returning decay arrays easier.

Definition at line 61 of file higgs_couplings_table.hpp.

◆ hp_decay_array_type

typedef const DecayTable::Entry* Gambit::HiggsCouplingsTable::hp_decay_array_type[max_charged_higgses]

Definition at line 62 of file higgs_couplings_table.hpp.

Constructor & Destructor Documentation

◆ HiggsCouplingsTable()

Gambit::HiggsCouplingsTable::HiggsCouplingsTable ( )

Constructor.

Constructor: initialise all effective couplings to 1.

Definition at line 24 of file higgs_couplings_table.cpp.

References C_bb2, C_cc2, C_gaga2, C_gg2, C_hiZ2, C_mumu2, C_ss2, C_tautau2, C_tt2, C_WW2, C_Zga2, C_ZZ2, and max_neutral_higgses.

   {
     for (int i = 0; i < max_neutral_higgses; i++)
     {
       C_WW2[i] = 1.0;
       C_ZZ2[i] = 1.0;
       C_tt2[i] = 1.0;
       C_bb2[i] = 1.0;
       C_cc2[i] = 1.0;
       C_tautau2[i] = 1.0;
       C_gaga2[i] = 1.0;
       C_gg2[i] = 1.0;
       C_mumu2[i] = 1.0;
       C_Zga2[i] = 1.0;
       C_ss2[i] = 1.0;
       for(int j = 0; j < max_neutral_higgses; j++) C_hiZ2[i][j] = 1.0;
     }
   }

Member Function Documentation

◆ compute_effective_coupling() [1/2]

template<typename T >

double Gambit::HiggsCouplingsTable::compute_effective_coupling	(	int	index,
		const T &	p1,
		const T &	p2
	)

inline

Compute a neutral higgs effective coupling from the current two-body neutral higgs decays.

Definition at line 93 of file higgs_couplings_table.hpp.

References Gambit::DecayTable::Entry::BF(), combine_hdf5::index, LOCAL_INFO, Gambit::utils_error(), and Gambit::DecayTable::Entry::width_in_GeV.

Referenced by Gambit::SpecBit::MSSM_higgs_couplings_FH(), and Gambit::SpecBit::MSSM_higgs_couplings_pwid().

       {
         if (index > max_neutral_higgses - 1) utils_error().raise(LOCAL_INFO, "Requested index beyond max_neutral_higgses.");
         // If channel is missing from either SM or BSM decays, return unity.
         if (!neutral_decays_SM_array[index]->has_channel(p1, p2) or !neutral_decays_array[index]->has_channel(p1, p2)) return 1.;
         double smwidth = neutral_decays_SM_array[index]->width_in_GeV;
         double smbf = neutral_decays_SM_array[index]->BF(p1, p2);
         if (smwidth <= 0. or smbf <= 0.) return 1.;
         double total_width_ratio = neutral_decays_array[index]->width_in_GeV / smwidth;
         double BF_ratio = neutral_decays_array[index]->BF(p1, p2) / smbf;
         return total_width_ratio * BF_ratio;
       }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ compute_effective_coupling() [2/2]

template<typename T >

double Gambit::HiggsCouplingsTable::compute_effective_coupling	(	str	name,
		const T &	p1,
		const T &	p2
	)

inline

Definition at line 106 of file higgs_couplings_table.hpp.

References get_charged_decays(), get_charged_decays_array(), get_neutral_decays(), get_neutral_decays_array(), get_neutral_decays_SM(), get_neutral_decays_SM_array(), get_t_decays(), LOCAL_INFO, set_charged_decays(), set_neutral_decays(), set_neutral_decays_SM(), set_t_decays(), and Gambit::utils_error().

       {
         if (neutral_decays_SM_map.find(name) == neutral_decays_SM_map.end() or
             neutral_decays_map.find(name) == neutral_decays_map.end())
          utils_error().raise(LOCAL_INFO, "Requested higgs not found.");
         // If channel is missing from either SM or BSM decays, return unity.
         if (!neutral_decays_SM_map.at(name).has_channel(p1, p2) or !neutral_decays_map.at(name).has_channel(p1, p2)) return 1.;
         double smwidth = neutral_decays_SM_map.at(name).width_in_GeV;
         double smbf = neutral_decays_SM_map.at(name).BF(p1, p2);
         if (smwidth <= 0. or smbf <= 0.) return 1.;
         double total_width_ratio = neutral_decays_map.at(name).width_in_GeV / smwidth;
         double BF_ratio = neutral_decays_map.at(name).BF(p1, p2) / smbf;
         return total_width_ratio * BF_ratio;
       }

Here is the call graph for this function:

◆ get_charged_decays() [1/2]

const DecayTable::Entry & Gambit::HiggsCouplingsTable::get_charged_decays ( int index ) const

Retrieve decays of a specific charged Higgs, by index.

Definition at line 120 of file higgs_couplings_table.cpp.

References charged_decays_array, combine_hdf5::index, LOCAL_INFO, max_charged_higgses, and Gambit::utils_error().

Referenced by compute_effective_coupling().

   {
     if (index > max_charged_higgses - 1) utils_error().raise(LOCAL_INFO, "Requested index beyond max_charged_higgses.");
     return *charged_decays_array[index];
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ get_charged_decays() [2/2]

const DecayTable::Entry & Gambit::HiggsCouplingsTable::get_charged_decays ( const str & name ) const

Retrieve decays of a specific charged Higgs, by name.

Definition at line 127 of file higgs_couplings_table.cpp.

References charged_decays_map, LOCAL_INFO, and Gambit::utils_error().

   {
     if (charged_decays_map.find(name) == charged_decays_map.end()) utils_error().raise(LOCAL_INFO, "Requested higgs not found.");
     return charged_decays_map.at(name);
   }

Here is the call graph for this function:

◆ get_charged_decays_array()

const HiggsCouplingsTable::hp_decay_array_type & Gambit::HiggsCouplingsTable::get_charged_decays_array ( int array_size ) const

Retrieve decays of all charged higgses.

Definition at line 113 of file higgs_couplings_table.cpp.

References charged_decays_array, LOCAL_INFO, max_charged_higgses, and Gambit::utils_error().

Referenced by compute_effective_coupling().

   {
     if (array_size > max_charged_higgses) utils_error().raise(LOCAL_INFO, "Requested array size beyond max_charged_higgses.");
     return charged_decays_array;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ get_neutral_decays() [1/2]

const DecayTable::Entry & Gambit::HiggsCouplingsTable::get_neutral_decays ( int index ) const

Retrieve decays of a specific neutral Higgs, by index.

Definition at line 99 of file higgs_couplings_table.cpp.

References combine_hdf5::index, LOCAL_INFO, max_neutral_higgses, neutral_decays_array, and Gambit::utils_error().

Referenced by compute_effective_coupling(), Gambit::SpecBit::MSSM_higgs_couplings_pwid(), and Gambit::ColliderBit::set_SMLikeHiggs_ModelParameters().

   {
     if (index > max_neutral_higgses - 1) utils_error().raise(LOCAL_INFO, "Requested index beyond max_neutral_higgses.");
     return *neutral_decays_array[index];
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ get_neutral_decays() [2/2]

const DecayTable::Entry & Gambit::HiggsCouplingsTable::get_neutral_decays ( const str & name ) const

Retrieve decays of a specific neutral Higgs, by name.

Definition at line 106 of file higgs_couplings_table.cpp.

References LOCAL_INFO, neutral_decays_map, and Gambit::utils_error().

   {
     if (neutral_decays_map.find(name) == neutral_decays_map.end()) utils_error().raise(LOCAL_INFO, "Requested higgs not found.");
     return neutral_decays_map.at(name);
   }

Here is the call graph for this function:

◆ get_neutral_decays_array()

const HiggsCouplingsTable::h0_decay_array_type & Gambit::HiggsCouplingsTable::get_neutral_decays_array ( int array_size ) const

Retrieve decays of all neutral higgses.

Definition at line 92 of file higgs_couplings_table.cpp.

References LOCAL_INFO, max_neutral_higgses, neutral_decays_array, and Gambit::utils_error().

Referenced by compute_effective_coupling().

   {
     if (array_size > max_neutral_higgses) utils_error().raise(LOCAL_INFO, "Requested arrray size beyond max_neutral_higgses.");
     return neutral_decays_array;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ get_neutral_decays_SM() [1/2]

const DecayTable::Entry & Gambit::HiggsCouplingsTable::get_neutral_decays_SM ( int index ) const

Retrieve SM decays of a specific neutral Higgs, by index.

Definition at line 78 of file higgs_couplings_table.cpp.

References combine_hdf5::index, LOCAL_INFO, max_neutral_higgses, neutral_decays_SM_array, and Gambit::utils_error().

Referenced by compute_effective_coupling().

   {
     if (index > max_neutral_higgses - 1) utils_error().raise(LOCAL_INFO, "Requested index beyond max_neutral_higgses.");
     return *neutral_decays_SM_array[index];
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ get_neutral_decays_SM() [2/2]

const DecayTable::Entry & Gambit::HiggsCouplingsTable::get_neutral_decays_SM ( const str & name ) const

Retrieve SM decays of a specific neutral Higgs, by name.

Definition at line 85 of file higgs_couplings_table.cpp.

References LOCAL_INFO, neutral_decays_SM_map, and Gambit::utils_error().

   {
     if (neutral_decays_SM_map.find(name) == neutral_decays_SM_map.end()) utils_error().raise(LOCAL_INFO, "Requested higgs not found.");
     return neutral_decays_SM_map.at(name);
   }

Here is the call graph for this function:

◆ get_neutral_decays_SM_array()

const HiggsCouplingsTable::h0_decay_array_type & Gambit::HiggsCouplingsTable::get_neutral_decays_SM_array ( int array_size ) const

Retrieve decay sets.

Retrieve SM decays of all neutral higgses.

Definition at line 71 of file higgs_couplings_table.cpp.

References LOCAL_INFO, max_neutral_higgses, neutral_decays_SM_array, and Gambit::utils_error().

Referenced by compute_effective_coupling().

   {
     if (array_size > max_neutral_higgses) utils_error().raise(LOCAL_INFO, "Requested array size beyond max_neutral_higgses.");
     return neutral_decays_SM_array;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ get_t_decays()

const DecayTable::Entry & Gambit::HiggsCouplingsTable::get_t_decays ( ) const

Retrieve decays of the top quark.

Definition at line 134 of file higgs_couplings_table.cpp.

References t_decays.

Referenced by compute_effective_coupling().

   {
     return *t_decays;
   }

Here is the caller graph for this function:

◆ set_charged_decays()

void Gambit::HiggsCouplingsTable::set_charged_decays	(	int	index,
		const str &	name,
		const DecayTable::Entry &	entry
	)

Assign decay entries to charged higgses.

Definition at line 60 of file higgs_couplings_table.cpp.

References charged_decays_array, charged_decays_map, combine_hdf5::index, LOCAL_INFO, max_charged_higgses, and Gambit::utils_error().

Referenced by compute_effective_coupling(), Gambit::SpecBit::MSSM_higgs_couplings_FH(), and Gambit::SpecBit::MSSM_higgs_couplings_pwid().

   {
     if (index > max_charged_higgses - 1) utils_error().raise(LOCAL_INFO, "Requested index beyond max_charged_higgses.");
     charged_decays_array[index] = &entry;
     charged_decays_map.insert(std::pair<str, const DecayTable::Entry&>(name,entry));
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ set_neutral_decays()

void Gambit::HiggsCouplingsTable::set_neutral_decays	(	int	index,
		const str &	name,
		const DecayTable::Entry &	entry
	)

Assign decay entries to neutral higgses.

Definition at line 52 of file higgs_couplings_table.cpp.

References combine_hdf5::index, LOCAL_INFO, max_neutral_higgses, neutral_decays_array, neutral_decays_map, and Gambit::utils_error().

Referenced by compute_effective_coupling(), Gambit::SpecBit::MSSM_higgs_couplings_FH(), Gambit::SpecBit::MSSM_higgs_couplings_pwid(), Gambit::SpecBit::ScalarSingletDM_higgs_couplings_pwid(), and Gambit::SpecBit::SM_higgs_couplings().

   {
     if (index > max_neutral_higgses - 1) utils_error().raise(LOCAL_INFO, "Requested index beyond max_neutral_higgses.");
     neutral_decays_array[index] = &entry;
     neutral_decays_map.insert(std::pair<str, const DecayTable::Entry&>(name,entry));
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ set_neutral_decays_SM()

void Gambit::HiggsCouplingsTable::set_neutral_decays_SM	(	int	index,
		const str &	name,
		const DecayTable::Entry &	entry
	)

Assign decay entries to the various table components.

Assign SM decay entries to neutral higgses.

Definition at line 44 of file higgs_couplings_table.cpp.

References combine_hdf5::index, LOCAL_INFO, max_neutral_higgses, neutral_decays_SM_array, neutral_decays_SM_map, and Gambit::utils_error().

Referenced by compute_effective_coupling(), Gambit::SpecBit::MSSM_higgs_couplings_FH(), Gambit::SpecBit::MSSM_higgs_couplings_pwid(), Gambit::SpecBit::ScalarSingletDM_higgs_couplings_pwid(), and Gambit::SpecBit::SM_higgs_couplings().

   {
     if (index > max_neutral_higgses - 1) utils_error().raise(LOCAL_INFO, "Requested index beyond max_neutral_higgses.");
     neutral_decays_SM_array[index] = &entry;
     neutral_decays_SM_map.insert(std::pair<str, const DecayTable::Entry&>(name,entry));
   }