Gambit::Spectrum Class Reference

"Standard Model" (low-energy) plus high-energy model container class More...

#include <spectrum.hpp>

Collaboration diagram for Gambit::Spectrum:

Public Types
typedef std::vector< YAML::sdd >	mc_info
	Typedefs for making it easier to manipulate mass cut and mass ratio cut info. More...
typedef std::vector< YAML::ssdd >	mr_info

Public Member Functions
void	RunBothToScale (double scale)
	Linked running Only possible with non-const object. More...
void	check_mass_cuts ()
	Check the that the spectrum satisifies any mass cuts requested from the yaml file. More...
SLHAstruct	getSLHAea (int) const
	SLHAea object getter First constructs an SLHAea object from the SMINPUTS object, then adds the info from the LE subspectrum (if possible), followed by the HE subspectrum (if possible). More...
void	writeSLHAfile (int, const str &) const
	Output spectrum contents as an SLHA file, using getSLHAea. More...
void	drop_SLHAs_if_requested (const safe_ptr< Options > &, const str &)
	Helper function to drop SLHA files. More...
const std::map< int, int > &	PDG_translator () const
	PDG code translation map, for special cases where an SLHA file has been read in and the PDG codes changed. More...
	Spectrum ()
	Default constructor. More...
	Spectrum (const SubSpectrum &le, const SubSpectrum &he, const SMInputs &smi, const std::map< str, safe_ptr< const double > > *input_Param, const mc_info &, const mr_info &)
	Construct new object, cloning the SubSpectrum objects supplied and taking possession of them. More...
	Spectrum (SubSpectrum *const le, SubSpectrum *const he, const SMInputs &smi, const std::map< str, safe_ptr< const double > > *input_Param, const mc_info &, const mr_info &)
	Construct new object, wrapping existing SubSpectrum objects Make sure the original objects don't get deleted before this wrapper does! More...
	Spectrum (const SubSpectrum &he, const SMInputs &smi, const std::map< str, safe_ptr< const double > > *input_Param, const mc_info &, const mr_info &)
	Construct new object, automatically creating an SMSimpleSpec as the LE subspectrum, and cloning the HE SubSpectrum object supplied and taking possession of it. More...
	Spectrum (const Spectrum &other)
	Copy constructor, clones SubSpectrum objects. More...
Spectrum &	operator= (const Spectrum &other)
	Copy-assignment Using "copy-and-swap" idiom. More...
	Spectrum (Spectrum &&other)
	Move constructor. More...
SubSpectrum &	get_LE ()
	Standard getters Return references to internal data members. More...
SubSpectrum &	get_HE ()
SMInputs &	get_SMInputs ()
const SubSpectrum &	get_LE () const
const SubSpectrum &	get_HE () const
const SMInputs &	get_SMInputs () const
std::unique_ptr< SubSpectrum >	clone_LE () const
	Clone getters Note: If you want to clone the whole Spectrum object, just use copy constructor, not these. More...
std::unique_ptr< SubSpectrum >	clone_HE () const
bool	has (const Par::Tags partype, const std::string &mass) const
	Pole mass getters/checkers "Shortcut" getters/checkers to access pole masses in hosted SubSpectrum objects. More...
double	get (const Par::Tags partype, const std::string &mass) const
bool	has (const Par::Tags partype, const std::string &mass, const int index) const
double	get (const Par::Tags partype, const std::string &mass, const int index) const
bool	has (const Par::Tags partype, const std::string &mass, const int index1, const int index2) const
double	get (const Par::Tags partype, const std::string &mass, const int index1, const int index2) const
bool	has (const Par::Tags partype, const int pdg_code, const int context) const
	PDB getter/checker overloads. More...
double	get (const Par::Tags partype, const int pdg_code, const int context) const
bool	has (const Par::Tags partype, const std::pair< int, int > pdgpr) const
double	get (const Par::Tags partype, const std::pair< int, int > pdgpr) const
bool	has (const Par::Tags partype, const std::pair< str, int > shortpr) const
double	get (const Par::Tags partype, const std::pair< str, int > shortpr) const
double	safeget (const Par::Tags partype, const std::string &mass) const
	Getters which first check the sanity of the thing they are returning. More...
double	safeget (const Par::Tags partype, const std::string &mass, const int index) const
double	safeget (const Par::Tags partype, const int pdg_code, const int context) const
double	safeget (const Par::Tags partype, const std::pair< int, int > pdgpr) const
double	safeget (const Par::Tags partype, const std::pair< str, int > shortpr) const

Static Public Member Functions
static double	Wolf2V_ud (double, double, double, double)
	CKM Wolfenstein (lambda, A, rhobar, etabar) –> V_qq standard parameterisation convertors. More...
static double	Wolf2V_us (double, double, double, double)
	CKM Wolfenstein –> V_us standard parameterisation convertor. More...
static std::complex< double >	Wolf2V_ub (double, double, double, double)
	CKM Wolfenstein –> V_ub standard parameterisation convertor. More...
static std::complex< double >	Wolf2V_cd (double, double, double, double)
	CKM Wolfenstein –> V_cd standard parameterisation convertor. More...
static std::complex< double >	Wolf2V_cs (double, double, double, double)
	CKM Wolfenstein –> V_cs standard parameterisation convertor. More...
static double	Wolf2V_cb (double, double, double, double)
	CKM Wolfenstein –> V_cb standard parameterisation convertor. More...
static std::complex< double >	Wolf2V_td (double, double, double, double)
	CKM Wolfenstein –> V_td standard parameterisation convertor. More...
static std::complex< double >	Wolf2V_ts (double, double, double, double)
	CKM Wolfenstein –> V_ts standard parameterisation convertor. More...
static double	Wolf2V_tb (double, double, double, double)
	CKM Wolfenstein –> V_tb standard parameterisation convertor. More...

Private Member Functions
void	check_init () const
	Check if object has been fully initialised. More...

Static Private Member Functions
static std::complex< double >	rhoplusieta (double, double, double, double)
	Calculate Wolfenstein rho+i*eta from rhobar and etabar. More...

Private Attributes
std::unique_ptr< SubSpectrum >	LE_new
	Variables. More...
std::unique_ptr< SubSpectrum >	HE_new
SubSpectrum *	LE
SubSpectrum *	HE
SMInputs	SMINPUTS
const std::map< str, safe_ptr< const double > > *	input_Param
const mc_info *	mass_cuts
const mr_info *	mass_ratio_cuts
bool	initialised

Friends
void	swap (Spectrum &first, Spectrum &second)
	Friend function: swap resources of two Spectrum objects. More...

Detailed Description

"Standard Model" (low-energy) plus high-energy model container class

Definition at line 110 of file spectrum.hpp.

Member Typedef Documentation

mc_info

typedef std::vector<YAML::sdd> Gambit::Spectrum::mc_info

Typedefs for making it easier to manipulate mass cut and mass ratio cut info.

Definition at line 119 of file spectrum.hpp.

mr_info

typedef std::vector<YAML::ssdd> Gambit::Spectrum::mr_info

Definition at line 120 of file spectrum.hpp.

Constructor & Destructor Documentation

Spectrum() [1/6]

Gambit::Spectrum::Spectrum

(

)

Default constructor.

Constructors/destructors.

Constructors/Destructors Need custom copy and move constructors plus copy-assignment operator in order to manage the unique_ptrs properly.

Default constructor

Definition at line 74 of file spectrum.cpp.

: input_Param(NULL), mass_cuts(NULL), mass_ratio_cuts(NULL), initialised(false) {}

Spectrum() [2/6]

Gambit::Spectrum::Spectrum	(	const SubSpectrum &	le,
		const SubSpectrum &	he,
		const SMInputs &	smi,
		const std::map< str, safe_ptr< const double > > *	input_Param,
		const mc_info &	mci,
		const mr_info &	mri
	)

Construct new object, cloning the SubSpectrum objects supplied and taking possession of them.

Definition at line 77 of file spectrum.cpp.

References check_mass_cuts().

     : LE_new(le.clone())
     , HE_new(he.clone())
     , LE(LE_new.get())
     , HE(HE_new.get())
     , SMINPUTS(smi)
     , input_Param(params)
     , mass_cuts(&mci)
     , mass_ratio_cuts(&mri)
     , initialised(true)
   { check_mass_cuts(); }

Here is the call graph for this function:

Spectrum() [3/6]

Gambit::Spectrum::Spectrum	(	SubSpectrum *const	le,
		SubSpectrum *const	he,
		const SMInputs &	smi,
		const std::map< str, safe_ptr< const double > > *	input_Param,
		const mc_info &	mci,
		const mr_info &	mri
	)

Construct new object, wrapping existing SubSpectrum objects Make sure the original objects don't get deleted before this wrapper does!

Definition at line 105 of file spectrum.cpp.

References check_mass_cuts().

     : LE(le)
     , HE(he)
     , SMINPUTS(smi)
     , input_Param(params)
     , mass_cuts(&mci)
     , mass_ratio_cuts(&mri)
     , initialised(true)
   { check_mass_cuts(); }

Here is the call graph for this function:

Spectrum() [4/6]

Gambit::Spectrum::Spectrum	(	const SubSpectrum &	he,
		const SMInputs &	smi,
		const std::map< str, safe_ptr< const double > > *	input_Param,
		const mc_info &	mci,
		const mr_info &	mri
	)

Construct new object, automatically creating an SMSimpleSpec as the LE subspectrum, and cloning the HE SubSpectrum object supplied and taking possession of it.

(won't make a version of this taking a pointer, since this is an "advanced" task, let people use the full contructor to do it.)

Definition at line 91 of file spectrum.cpp.

References check_mass_cuts().

     : LE_new(SMSimpleSpec(smi).clone())
     , HE_new(he.clone())
     , LE(LE_new.get())
     , HE(HE_new.get())
     , SMINPUTS(smi)
     , input_Param(params)
     , mass_cuts(&mci)
     , mass_ratio_cuts(&mri)
     , initialised(true)
   { check_mass_cuts(); }

Here is the call graph for this function:

Spectrum() [5/6]

Gambit::Spectrum::Spectrum

(

const Spectrum &

other

)

Copy constructor, clones SubSpectrum objects.

Make a non-const copy in order to use e.g. RunBothToScale function.

Definition at line 118 of file spectrum.cpp.

     : LE_new(other.clone_LE())
     , HE_new(other.clone_HE())
     , LE(LE_new.get())
     , HE(HE_new.get())
     , SMINPUTS(other.SMINPUTS)
     , input_Param(other.input_Param)
     , mass_cuts(other.mass_cuts)
     , mass_ratio_cuts(other.mass_ratio_cuts)
     , initialised(other.initialised)
   {}

Spectrum() [6/6]

Gambit::Spectrum::Spectrum

(

Spectrum &&

other

)

Move constructor.

Definition at line 140 of file spectrum.cpp.

References swap.

   {
      swap(*this, other);
   }

Member Function Documentation

check_init()

void Gambit::Spectrum::check_init ( ) const

private

Check if object has been fully initialised.

Spectrum class member function definitions.

Check if object has been fully initialised

Definition at line 50 of file spectrum.cpp.

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

Referenced by clone_HE(), clone_LE(), get_HE(), get_LE(), and get_SMInputs().

   {
     if(not initialised) utils_error().raise(LOCAL_INFO,"Access or deepcopy of empty Spectrum object attempted!");
   }

Here is the call graph for this function:

Here is the caller graph for this function:

check_mass_cuts()

void Gambit::Spectrum::check_mass_cuts

(

)

Check the that the spectrum satisifies any mass cuts requested from the yaml file.

Definition at line 168 of file spectrum.cpp.

References has(), Gambit::invalid_point(), Gambit::is_abs(), LOCAL_INFO, generate_raster_scan_settings::m1, generate_raster_scan_settings::m2, mass_cuts, mass_ratio_cuts, Gambit::Par::Pole_Mass, Gambit::invalid_point_exception::raise(), and Gambit::utils_error().

Referenced by Spectrum().

   {
     if (mass_cuts != NULL and not mass_cuts->empty())
     {
       for (auto it = mass_cuts->begin(); it != mass_cuts->end(); ++it)
       {
         str p = it->first;
         bool absolute_value = is_abs(p);
         const double& low = it->second.first;
         const double& high = it->second.second;
         #ifdef SPECTRUM_DEBUG
           cout << "Applying mass cut " << low << " GeV < " << (absolute_value ? "|mass("+p+")|" : "mass("+p+")") << " < " << high << " GeV" << endl;
         #endif
         if (not has(Par::Pole_Mass, p)) utils_error().raise(LOCAL_INFO, "Cannot cut on mass of unrecognised particle: " + p);
         double m = get(Par::Pole_Mass, p);
         if (absolute_value) m = std::abs(m);
         #ifdef SPECTRUM_DEBUG
           cout << "Actual value: " << m << endl;
         #endif
         if (m < low or m > high) invalid_point().raise(p + " failed requested mass cut.");
       }
     }
     if (mass_ratio_cuts != NULL and not mass_ratio_cuts->empty())
     {
       for (auto it = mass_ratio_cuts->begin(); it != mass_ratio_cuts->end(); ++it)
       {
         str p1 = it->first.first;
         str p2 = it->first.second;
         bool absolute_value1 = is_abs(p1);
         bool absolute_value2 = is_abs(p2);
         const double& low = it->second.first;
         const double& high = it->second.second;
         #ifdef SPECTRUM_DEBUG
           cout << "Applying mass ratio cut " << low << " < "
                << (absolute_value1 ? "|mass("+p1+")|" : "mass("+p1+")") << " / "
                << (absolute_value2 ? "|mass("+p2+")|" : "mass("+p2+")")
                << " < " << high << endl;
         #endif
         if (not has(Par::Pole_Mass, p1)) utils_error().raise(LOCAL_INFO, "Cannot cut on ratio with mass of unrecognised particle: " + p1);
         if (not has(Par::Pole_Mass, p2)) utils_error().raise(LOCAL_INFO, "Cannot cut on ratio with mass of unrecognised particle: " + p2);
         double m1 = get(Par::Pole_Mass, p1);
         double m2 = get(Par::Pole_Mass, p2);
         if (absolute_value1) m1 = std::abs(m1);
         if (absolute_value2) m2 = std::abs(m2);
         double mratio = m1/m2;
         #ifdef SPECTRUM_DEBUG
           cout << "Actual value: " << mratio << endl;
         #endif
         if (mratio < low or mratio > high) invalid_point().raise(p1 + "/" + p2 +" failed requested mass ratio cut.");
       }
     }
     }

Here is the call graph for this function:

Here is the caller graph for this function:

clone_HE()

std::unique_ptr< SubSpectrum > Gambit::Spectrum::clone_HE

(

)

const

Definition at line 235 of file spectrum.cpp.

References check_init(), Gambit::SubSpectrum::clone(), and HE.

Referenced by Gambit::SpecBit::check_perturb(), Gambit::SpecBit::check_perturb_MDM(), Gambit::SpecBit::check_perturb_to_min_lambda(), and Gambit::SpecBit::find_min_lambda_Helper().

{check_init(); return HE->clone();}

Here is the call graph for this function:

Here is the caller graph for this function:

clone_LE()

std::unique_ptr< SubSpectrum > Gambit::Spectrum::clone_LE

(

)

const

Clone getters Note: If you want to clone the whole Spectrum object, just use copy constructor, not these.

Clone SubSpectrum getters To clone whole object, just use copy constructor.

Definition at line 234 of file spectrum.cpp.

References check_init(), Gambit::SubSpectrum::clone(), and LE.

Referenced by Gambit::SpecBit::Spectrum_test().

{check_init(); return LE->clone();}

Here is the call graph for this function:

Here is the caller graph for this function:

drop_SLHAs_if_requested()

void Gambit::Spectrum::drop_SLHAs_if_requested	(	const safe_ptr< Options > &	runOptions,
		const str &	default_name
	)

Helper function to drop SLHA files.

Definition at line 452 of file spectrum.cpp.

References writeSLHAfile().

Referenced by Gambit::SpecBit::get_MSSM_spectrum_SPheno(), Gambit::PrecisionBit::make_MSSM_precision_spectrum_4H_W(), Gambit::PrecisionBit::make_MSSM_precision_spectrum_H(), Gambit::PrecisionBit::make_MSSM_precision_spectrum_H_W(), Gambit::PrecisionBit::make_MSSM_precision_spectrum_none(), and Gambit::PrecisionBit::make_MSSM_precision_spectrum_W().

   {
      if (runOptions->getValueOrDef<bool>(false, "drop_SLHA_file"))
      {
         // Spit out the full spectrum as SLHA1 and SLHA2 files.
         str prefix   = runOptions->getValueOrDef<str>("", "SLHA_output_prefix");
         str filename = runOptions->getValueOrDef<str>(default_name, "SLHA_output_filename");
         writeSLHAfile(1,prefix+filename+".slha1");
         writeSLHAfile(2,prefix+filename+".slha2");
      }
   }

Here is the call graph for this function:

Here is the caller graph for this function:

get() [1/6]

double Gambit::Spectrum::get	(	const Par::Tags	partype,
		const std::string &	mass
	)		const

Definition at line 249 of file spectrum.cpp.

References Gambit::SubSpectrum::get(), Gambit::SubSpectrum::has(), HE, LE, LOCAL_INFO, and Gambit::utils_error().

Referenced by Gambit::ColliderBit::ALEPH_Selectron_Conservative_LLike(), Gambit::ColliderBit::ALEPH_Smuon_Conservative_LLike(), Gambit::ColliderBit::ALEPH_Stau_Conservative_LLike(), Gambit::ColliderBit::calc_susy_spectrum_scan_guide(), Gambit::SpecBit::check_EW_stability_ScalarSingletDM_Z3(), Gambit::DarkBit::DarkMatter_ID_MSSM(), Gambit::DarkBit::DarkMatterConj_ID_MSSM(), Gambit::DarkBit::DD_couplings_DiracSingletDM_Z2(), Gambit::DarkBit::DD_couplings_MajoranaSingletDM_Z2(), Gambit::DarkBit::DD_couplings_VectorSingletDM_Z2(), Gambit::DecayBit::DiracSingletDM_Higgs_decays(), Gambit::SpecBit::find_min_lambda_Helper(), Gambit::DarkBit::get_ScalarSingletDM_DD_couplings(), Gambit::SpecBit::has_neutralino_LSP(), Gambit::ColliderBit::L3_Chargino_All_Channels_Conservative_LLike(), Gambit::ColliderBit::L3_Chargino_Leptonic_Conservative_LLike(), Gambit::ColliderBit::L3_Gravitino_LLike(), Gambit::ColliderBit::L3_Neutralino_All_Channels_Conservative_LLike(), Gambit::ColliderBit::L3_Neutralino_Leptonic_Conservative_LLike(), Gambit::ColliderBit::L3_Selectron_Conservative_LLike(), Gambit::ColliderBit::L3_Smuon_Conservative_LLike(), Gambit::ColliderBit::L3_Stau_Conservative_LLike(), Gambit::DecayBit::MajoranaSingletDM_Higgs_decays(), Gambit::DecayBit::MSSM_inv_Higgs_BF(), Gambit::SpecBit::MSSMspectrum_test(), Gambit::ColliderBit::OPAL_Chargino_Hadronic_Conservative_LLike(), Gambit::ColliderBit::OPAL_Chargino_Leptonic_Conservative_LLike(), Gambit::ColliderBit::OPAL_Degenerate_Chargino_Conservative_LLike(), Gambit::ColliderBit::OPAL_Neutralino_Hadronic_Conservative_LLike(), Gambit::DecayBit::stau_1_decays_smallsplit(), Gambit::DarkBit::TH_ProcessCatalog_DiracSingletDM_Z2(), Gambit::DarkBit::TH_ProcessCatalog_MajoranaSingletDM_Z2(), Gambit::DarkBit::TH_ProcessCatalog_ScalarSingletDM_Z2(), Gambit::DarkBit::TH_ProcessCatalog_ScalarSingletDM_Z3(), Gambit::DarkBit::TH_ProcessCatalog_VectorSingletDM_Z2(), Gambit::DecayBit::VectorSingletDM_Higgs_decays(), and Gambit::DecayBit::Z_gamma_chi_0_MSSM_tree().

   {
     double result(-1);
     if( HE->has(partype,mass) )
     { result = HE->get(partype,mass); }
     else if( LE->has(partype,mass) )
     { result = LE->get(partype,mass); }
     else
     {
        std::ostringstream errmsg;
        errmsg << "Error retrieving particle spectrum data!" << std::endl;
        errmsg << "No pole mass with string reference '"<<mass<<"' could be found in either LE or HE SubSpectrum!" <<std::endl;
        utils_error().raise(LOCAL_INFO,errmsg.str());
     }
     // In c++11 we could add the [[noreturn]] attribute utils_error.raise()
     // to suppress the compiler warning about not returning anything (and enable
     // extra optimisations), however it isn't implemented in gcc until version
     // 4.8 (and we decided to support earlier versions).
     return result;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

get() [2/6]

double Gambit::Spectrum::get	(	const Par::Tags	partype,
		const std::string &	mass,
		const int	index
	)		const

Definition at line 275 of file spectrum.cpp.

References Gambit::SubSpectrum::get(), Gambit::SubSpectrum::has(), HE, LE, LOCAL_INFO, and Gambit::utils_error().

   {
     double result(-1);
     if( HE->has(partype,mass,index) )
     { result = HE->get(partype,mass,index); }
     else if( LE->has(partype,mass,index) )
     { result = LE->get(partype,mass,index); }
     else
     {
        std::ostringstream errmsg;
        errmsg << "Error retrieving particle spectrum data!" << std::endl;
        errmsg << "No pole mass with string reference '"<<mass<<"' and index '"<<index<<"' could be found in either LE or HE SubSpectrum!" <<std::endl;
        utils_error().raise(LOCAL_INFO,errmsg.str());
     }
     // [[noreturn]]
     return result;
   }

Here is the call graph for this function:

get() [3/6]

double Gambit::Spectrum::get	(	const Par::Tags	partype,
		const std::string &	mass,
		const int	index1,
		const int	index2
	)		const

Definition at line 298 of file spectrum.cpp.

References Gambit::SubSpectrum::get(), Gambit::SubSpectrum::has(), HE, LE, LOCAL_INFO, and Gambit::utils_error().

   {
     double result(-1);
     if( HE->has(partype,mass,index1,index2) )
     { result = HE->get(partype,mass,index1,index2); }
     else if( LE->has(partype,mass,index1,index2) )
     { result = LE->get(partype,mass,index1,index2); }
     else
     {
       std::ostringstream errmsg;
       errmsg << "Error retrieving particle spectrum data!" << std::endl;
       errmsg << "No pole mixing with string reference '"<<mass<<"' and indices '"<<index1<<"','"<<index2<<"' could be found in either LE or HE SubSepctrum!" << std::endl;
       utils_error().raise(LOCAL_INFO,errmsg.str());
     }
     // [[noreturn]]
     return result;
   }

Here is the call graph for this function:

get() [4/6]

double Gambit::Spectrum::get	(	const Par::Tags	partype,
		const int	pdg_code,
		const int	context
	)		const

Definition at line 326 of file spectrum.cpp.

   {
      return get( partype, std::make_pair(pdg_code,context) );
   }

get() [5/6]

double Gambit::Spectrum::get	(	const Par::Tags	partype,
		const std::pair< int, int >	pdgpr
	)		const

Definition at line 348 of file spectrum.cpp.

References Gambit::Models::partmap::long_name(), Gambit::Models::ParticleDB(), and Gambit::Models::partmap::short_name_pair().

   {
      /* If there is a short name, then retrieve that plus the index */
      if( Models::ParticleDB().has_short_name(pdgpr) )
      {
        return get( partype, Models::ParticleDB().short_name_pair(pdgpr) );
      }
      else /* Use the long name with no index instead */
      {
        return get( partype, Models::ParticleDB().long_name(pdgpr) );
      }
   }

Here is the call graph for this function:

get() [6/6]

double Gambit::Spectrum::get	(	const Par::Tags	partype,
		const std::pair< str, int >	shortpr
	)		const

Definition at line 370 of file spectrum.cpp.

   {
      return get( partype, shortpr.first, shortpr.second);
   }

get_HE() [1/2]

SubSpectrum & Gambit::Spectrum::get_HE

(

)

Definition at line 225 of file spectrum.cpp.

References check_init(), and HE.

Referenced by Gambit::ColliderBit::ALEPH_Selectron_Conservative_LLike(), Gambit::ColliderBit::ALEPH_Smuon_Conservative_LLike(), Gambit::ColliderBit::ALEPH_Stau_Conservative_LLike(), Gambit::DarkBit::DD_couplings_DiracSingletDM_Z2(), Gambit::DarkBit::DD_couplings_MajoranaSingletDM_Z2(), Gambit::DarkBit::DD_couplings_VectorSingletDM_Z2(), Gambit::DecayBit::DiracSingletDM_Higgs_decays(), Gambit::SpecBit::exampleRead(), Gambit::SpecBit::get_MSSM_spectrum_as_map(), Gambit::SpecBit::get_MSSM_spectrum_from_SLHAfile(), Gambit::DarkBit::get_ScalarSingletDM_DD_couplings(), Gambit::ColliderBit::get_sigma_ee_chi00(), Gambit::ColliderBit::get_sigma_ee_chipm(), Gambit::ColliderBit::get_sigma_ee_ll(), Gambit::SpecBit::get_unimproved_MSSM_spectrum_as_map(), Gambit::ColliderBit::L3_Selectron_Conservative_LLike(), Gambit::ColliderBit::L3_Smuon_Conservative_LLike(), Gambit::ColliderBit::L3_Stau_Conservative_LLike(), Gambit::DecayBit::MajoranaSingletDM_Higgs_decays(), Gambit::PrecisionBit::make_MSSM_precision_spectrum_4H_W(), Gambit::PrecisionBit::make_MSSM_precision_spectrum_H(), Gambit::PrecisionBit::make_MSSM_precision_spectrum_H_W(), Gambit::PrecisionBit::make_MSSM_precision_spectrum_W(), Gambit::SpecBit::MSSMspectrum_test(), Gambit::ColliderBit::OPAL_Chargino_Leptonic_Conservative_LLike(), Gambit::SpecBit::prepare_pass_MSSM_spectrum_to_vevacious(), Gambit::DarkBit::RD_spectrum_MSSM(), Gambit::DecayBit::stau_1_decays_smallsplit(), Gambit::DarkBit::TH_ProcessCatalog_DiracSingletDM_Z2(), Gambit::DarkBit::TH_ProcessCatalog_DS5_MSSM(), Gambit::DarkBit::TH_ProcessCatalog_DS_MSSM(), Gambit::DarkBit::TH_ProcessCatalog_MajoranaSingletDM_Z2(), Gambit::DarkBit::TH_ProcessCatalog_ScalarSingletDM_Z2(), Gambit::DarkBit::TH_ProcessCatalog_ScalarSingletDM_Z3(), Gambit::DarkBit::TH_ProcessCatalog_VectorSingletDM_Z2(), and Gambit::DecayBit::VectorSingletDM_Higgs_decays().

{check_init(); return *HE;}

Here is the call graph for this function:

Here is the caller graph for this function:

get_HE() [2/2]

const SubSpectrum & Gambit::Spectrum::get_HE

(

)

const

Definition at line 229 of file spectrum.cpp.

References check_init(), and HE.

{check_init(); return *HE;}

Here is the call graph for this function:

get_LE() [1/2]

SubSpectrum & Gambit::Spectrum::get_LE

(

)

Standard getters Return references to internal data members.

Standard SubSpectrum getters Return references to internal data members. Make sure original Spectrum object doesn't get destroyed before you finish using these or you will cause a segfault.

Make sure original Spectrum object doesn't get destroyed before you finish using these or you will cause a segfault.

Definition at line 224 of file spectrum.cpp.

References check_init(), and LE.

Referenced by Gambit::SpecBit::exampleRead(), Gambit::SpecBit::get_MSSM_spectrum_as_map(), Gambit::SpecBit::get_SM_SubSpectrum_from_MSSM_Spectrum(), Gambit::SpecBit::get_unimproved_MSSM_spectrum_as_map(), Gambit::PrecisionBit::make_MSSM_precision_spectrum_4H_W(), Gambit::PrecisionBit::make_MSSM_precision_spectrum_H_W(), Gambit::PrecisionBit::make_MSSM_precision_spectrum_W(), Gambit::SpecBit::MSSMspectrum_test(), Gambit::DarkBit::RD_spectrum_MSSM(), Gambit::SpecBit::Spectrum_test(), Gambit::DarkBit::TH_ProcessCatalog_DiracSingletDM_Z2(), Gambit::DarkBit::TH_ProcessCatalog_DS5_MSSM(), Gambit::DarkBit::TH_ProcessCatalog_DS_MSSM(), Gambit::DarkBit::TH_ProcessCatalog_MajoranaSingletDM_Z2(), Gambit::DarkBit::TH_ProcessCatalog_ScalarSingletDM_Z2(), Gambit::DarkBit::TH_ProcessCatalog_ScalarSingletDM_Z3(), and Gambit::DarkBit::TH_ProcessCatalog_VectorSingletDM_Z2().

{check_init(); return *LE;}

Here is the call graph for this function:

Here is the caller graph for this function:

get_LE() [2/2]

const SubSpectrum & Gambit::Spectrum::get_LE

(

)

const

Definition at line 228 of file spectrum.cpp.

References check_init(), and LE.

{check_init(); return *LE;}

Here is the call graph for this function:

get_SMInputs() [1/2]

SMInputs & Gambit::Spectrum::get_SMInputs

(

)

Definition at line 226 of file spectrum.cpp.

References check_init(), and SMINPUTS.

Referenced by Gambit::SpecBit::Spectrum_test(), Gambit::DarkBit::TH_ProcessCatalog_DiracSingletDM_Z2(), Gambit::DarkBit::TH_ProcessCatalog_DS5_MSSM(), Gambit::DarkBit::TH_ProcessCatalog_DS_MSSM(), Gambit::DarkBit::TH_ProcessCatalog_MajoranaSingletDM_Z2(), Gambit::DarkBit::TH_ProcessCatalog_ScalarSingletDM_Z2(), Gambit::DarkBit::TH_ProcessCatalog_ScalarSingletDM_Z3(), and Gambit::DarkBit::TH_ProcessCatalog_VectorSingletDM_Z2().

{check_init(); return SMINPUTS;}

Here is the call graph for this function:

Here is the caller graph for this function:

get_SMInputs() [2/2]

const SMInputs & Gambit::Spectrum::get_SMInputs

(

)

const

Definition at line 230 of file spectrum.cpp.

References check_init(), and SMINPUTS.

{check_init(); return SMINPUTS;}

Here is the call graph for this function:

getSLHAea()

SLHAstruct Gambit::Spectrum::getSLHAea

(

int

slha_version

)

const

SLHAea object getter First constructs an SLHAea object from the SMINPUTS object, then adds the info from the LE subspectrum (if possible), followed by the HE subspectrum (if possible).

Any duplicate entries are overwritten at each step, so HE takes precendence over LE, and LE takes precedence over SMINPUTS.

Definition at line 432 of file spectrum.cpp.

References Gambit::SubSpectrum::add_to_SLHAea(), Gambit::SMInputs::getSLHAea(), HE, LE, and SMINPUTS.

Referenced by Gambit::SpecBit::exampleRead(), Gambit::SpecBit::MSSMspectrum_test(), Gambit::SpecBit::SHD_HiggsMass(), Gambit::SpecBit::SPheno_MSSM_test(), and writeSLHAfile().

   {
      SLHAstruct slha(SMINPUTS.getSLHAea());
      LE->add_to_SLHAea(slha_version, slha);
      HE->add_to_SLHAea(slha_version, slha);
      return slha;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

has() [1/6]

bool Gambit::Spectrum::has	(	const Par::Tags	partype,
		const std::string &	mass
	)		const

Pole mass getters/checkers "Shortcut" getters/checkers to access pole masses in hosted SubSpectrum objects.

Pole mass getters "Shortcut" getters to access pole masses in hosted SubSpectrum objects. HE object given higher priority; if no match found, LE object will be checked. If still no match, error is thrown.

HE object given higher priority; if no match found, LE object will be checked. If still no match, error is thrown. TODO: These currently work for anything! Need to restrict them to only allow access to pole masses and their estimated uncertainties Also need to change error messages etc, plus the PDG overloads

Definition at line 244 of file spectrum.cpp.

References Gambit::SubSpectrum::has(), HE, and LE.

Referenced by check_mass_cuts(), and has().

   {
     return (HE->has(partype,mass) or LE->has(partype,mass));
   }

Here is the call graph for this function:

Here is the caller graph for this function:

has() [2/6]

bool Gambit::Spectrum::has	(	const Par::Tags	partype,
		const std::string &	mass,
		const int	index
	)		const

Definition at line 270 of file spectrum.cpp.

References Gambit::SubSpectrum::has(), HE, and LE.

   {
     return (HE->has(partype,mass,index) or LE->has(partype,mass,index));
   }

Here is the call graph for this function:

has() [3/6]

bool Gambit::Spectrum::has	(	const Par::Tags	partype,
		const std::string &	mass,
		const int	index1,
		const int	index2
	)		const

Definition at line 293 of file spectrum.cpp.

References Gambit::SubSpectrum::has(), HE, and LE.

   {
     return (HE->has(partype,mass,index1,index2) or LE->has(partype,mass,index1,index2));
   }

Here is the call graph for this function:

has() [4/6]

bool Gambit::Spectrum::has	(	const Par::Tags	partype,
		const int	pdg_code,
		const int	context
	)		const

PDB getter/checker overloads.

Definition at line 319 of file spectrum.cpp.

References has().

   {
      return has( partype, std::make_pair(pdg_code,context) );
   }

Here is the call graph for this function:

has() [5/6]

bool Gambit::Spectrum::has	(	const Par::Tags	partype,
		const std::pair< int, int >	pdgpr
	)		const

Definition at line 333 of file spectrum.cpp.

References has(), and Gambit::Models::ParticleDB().

   {
      /* If there is a short name, then retrieve that plus the index */
      if( Models::ParticleDB().has_short_name(pdgpr) )
      {
        return has( partype, Models::ParticleDB().short_name_pair(pdgpr) );
      }
      else /* Use the long name with no index instead */
      {
        return has( partype, Models::ParticleDB().long_name(pdgpr) );
      }
   }

Here is the call graph for this function:

has() [6/6]

bool Gambit::Spectrum::has	(	const Par::Tags	partype,
		const std::pair< str, int >	shortpr
	)		const

Definition at line 363 of file spectrum.cpp.

References has().

   {
      return has( partype, shortpr.first, shortpr.second);
   }

Here is the call graph for this function:

operator=()

Spectrum & Gambit::Spectrum::operator=

(

const Spectrum &

other

)

Copy-assignment Using "copy-and-swap" idiom.

Definition at line 132 of file spectrum.cpp.

References swap.

   {
      Spectrum temp(other);
      swap(*this, temp);
      return *this;
   }

PDG_translator()

const std::map< int, int > & Gambit::Spectrum::PDG_translator

(

)

const

PDG code translation map, for special cases where an SLHA file has been read in and the PDG codes changed.

Definition at line 465 of file spectrum.cpp.

References HE, and Gambit::SubSpectrum::PDG_translator().

   {
     return HE->PDG_translator();
   }

Here is the call graph for this function:

rhoplusieta()

std::complex< double > Gambit::Spectrum::rhoplusieta ( double  lambda, double  A, double  rhobar, double  etabar  )

staticprivate

Calculate Wolfenstein rho+i*eta from rhobar and etabar.

Helper function to calculate Wolfenstein rho+i*eta from rhobar and etabar.

Definition at line 473 of file spectrum.cpp.

References Gambit::Scanner::pow().

Referenced by Wolf2V_cb(), Wolf2V_cd(), Wolf2V_cs(), Wolf2V_tb(), Wolf2V_td(), Wolf2V_ts(), Wolf2V_ub(), Wolf2V_ud(), and Wolf2V_us().

   {
     std::complex<double> x(rhobar, etabar);
     double y = pow(A*lambda*lambda,2);
     return sqrt((1.0-y)/(1.0-lambda*lambda))*x/(1.0-x*y);
   }

Here is the call graph for this function:

Here is the caller graph for this function:

RunBothToScale()

void Gambit::Spectrum::RunBothToScale

(

double

scale

)

Linked running Only possible with non-const object.

Overloads for PDG types These just convert the types and then call the properly defined functions.

Linked running Only possible with non-const object

Definition at line 153 of file spectrum.cpp.

References HE, LE, and Gambit::SubSpectrum::RunToScale().

Referenced by Gambit::SpecBit::Spectrum_test().

   {
     LE->RunToScale(scale);
     HE->RunToScale(scale);
   }

Here is the call graph for this function:

Here is the caller graph for this function:

safeget() [1/5]

double Gambit::Spectrum::safeget	(	const Par::Tags	partype,
		const std::string &	mass
	)		const

Getters which first check the sanity of the thing they are returning.

Definition at line 380 of file spectrum.cpp.

References LOCAL_INFO, and Gambit::utils_error().

Referenced by Gambit::ColliderBit::get_sigma_ee_chi00(), Gambit::ColliderBit::get_sigma_ee_chipm(), Gambit::ColliderBit::get_sigma_ee_ll(), and Gambit::DecayBit::stau_1_decays_smallsplit().

   {
      double result = get(partype, mass);
      if (Utils::isnan(result))
         utils_error().raise(LOCAL_INFO,"SubSpectrum parameter is nan!!");
      return result;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

safeget() [2/5]

double Gambit::Spectrum::safeget	(	const Par::Tags	partype,
		const std::string &	mass,
		const int	index
	)		const

Definition at line 389 of file spectrum.cpp.

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

   {
      double result = get(partype, mass, index);
      if (Utils::isnan(result))
         utils_error().raise(LOCAL_INFO,"SubSpectrum parameter is nan!!");
      return result;
   }

Here is the call graph for this function:

safeget() [3/5]

double Gambit::Spectrum::safeget	(	const Par::Tags	partype,
		const int	pdg_code,
		const int	context
	)		const

Definition at line 398 of file spectrum.cpp.

References LOCAL_INFO, and Gambit::utils_error().

   {
      double result = get(partype, pdg_code, context);
      if (Utils::isnan(result))
         utils_error().raise(LOCAL_INFO,"SubSpectrum parameter is nan!!");
      return result;
   }

Here is the call graph for this function:

safeget() [4/5]

double Gambit::Spectrum::safeget	(	const Par::Tags	partype,
		const std::pair< int, int >	pdgpr
	)		const

Definition at line 407 of file spectrum.cpp.

References LOCAL_INFO, and Gambit::utils_error().

   {
      double result = get(partype, pdgpr);
      if (Utils::isnan(result))
         utils_error().raise(LOCAL_INFO,"SubSpectrum parameter is nan!!");
      return result;
   }

Here is the call graph for this function:

safeget() [5/5]

double Gambit::Spectrum::safeget	(	const Par::Tags	partype,
		const std::pair< str, int >	shortpr
	)		const

Definition at line 416 of file spectrum.cpp.

References LOCAL_INFO, and Gambit::utils_error().

   {
      double result = get(partype, shortpr);
      if (Utils::isnan(result))
         utils_error().raise(LOCAL_INFO,"SubSpectrum parameter is nan!!");
      return result;
   }

Here is the call graph for this function:

Wolf2V_cb()

double Gambit::Spectrum::Wolf2V_cb ( double  l, double  A, double  rhobar, double  etabar  )

static

CKM Wolfenstein –> V_cb standard parameterisation convertor.

Definition at line 519 of file spectrum.cpp.

References Gambit::Scanner::pow(), and rhoplusieta().

   {
     return A*l*l * (1.0 - 0.5*A*A*pow(l,6)*std::norm(rhoplusieta(l,A,rhobar,etabar)));
   }

Here is the call graph for this function:

Wolf2V_cd()

std::complex< double > Gambit::Spectrum::Wolf2V_cd ( double  l, double  A, double  rhobar, double  etabar  )

static

CKM Wolfenstein –> V_cd standard parameterisation convertor.

Definition at line 502 of file spectrum.cpp.

References Gambit::Scanner::pow(), and rhoplusieta().

   {
     std::complex<double> x(rhoplusieta(l,A,rhobar,etabar));
     return 0.5*pow(A*l,2)*(pow(l,3)*(1.0-2.0*x) + pow(l,5)*x) - l;
   }

Here is the call graph for this function:

Wolf2V_cs()

std::complex< double > Gambit::Spectrum::Wolf2V_cs ( double  l, double  A, double  rhobar, double  etabar  )

static

CKM Wolfenstein –> V_cs standard parameterisation convertor.

Definition at line 509 of file spectrum.cpp.

References Gambit::Scanner::pow(), and rhoplusieta().

   {
     double l2 = l*l;
     double fA2 = 4.0*A*A;
     return 1.0 - 0.5*l2 - 0.125*l2*l2*(1.0+fA2)
            - 0.0625*pow(l2,3)*(1.0-fA2+4.0*fA2*rhoplusieta(l,A,rhobar,etabar))
            - 0.0078125*pow(l2,4)*(5.0-fA2*(2.0+4.0*fA2));
   }

Here is the call graph for this function:

Wolf2V_tb()

double Gambit::Spectrum::Wolf2V_tb ( double  l, double  A, double  rhobar, double  etabar  )

static

CKM Wolfenstein –> V_tb standard parameterisation convertor.

Definition at line 539 of file spectrum.cpp.

References Gambit::Scanner::pow(), and rhoplusieta().

   {
     double norm = std::norm(rhoplusieta(l,A,rhobar,etabar));
     double l4 = pow(l,4);
     return 1.0 - 0.5*A*A*l4 * (1.0 + l*l*norm + 0.25*A*A*l4);
   }

Here is the call graph for this function:

Wolf2V_td()

std::complex< double > Gambit::Spectrum::Wolf2V_td ( double  l, double  A, double  rhobar, double  etabar  )

static

CKM Wolfenstein –> V_td standard parameterisation convertor.

Definition at line 525 of file spectrum.cpp.

References Gambit::Scanner::pow(), and rhoplusieta().

   {
     std::complex<double> x(rhoplusieta(l,A,rhobar,etabar));
     return A*l*l * (l*(1.0-x) + 0.5*pow(l,3)*x + 0.125*pow(l,5)*(1.0+4.0*A*A)*x);
   }

Here is the call graph for this function:

Wolf2V_ts()

std::complex< double > Gambit::Spectrum::Wolf2V_ts ( double  l, double  A, double  rhobar, double  etabar  )

static

CKM Wolfenstein –> V_ts standard parameterisation convertor.

Definition at line 532 of file spectrum.cpp.

References Gambit::Scanner::pow(), and rhoplusieta().

   {
     std::complex<double> x(rhoplusieta(l,A,rhobar,etabar));
     return A*l*l * (0.5*pow(l,2)*(1.0-2.0*x) + 0.125*pow(l,4) + 0.0625*pow(l,6)*(1.0+8.0*A*A*x) - 1.0);
   }

Here is the call graph for this function:

Wolf2V_ub()

std::complex< double > Gambit::Spectrum::Wolf2V_ub ( double  l, double  A, double  rhobar, double  etabar  )

static

CKM Wolfenstein –> V_ub standard parameterisation convertor.

Definition at line 496 of file spectrum.cpp.

References Gambit::Scanner::pow(), and rhoplusieta().

   {
     return A*pow(l,3)*std::conj(rhoplusieta(l,A,rhobar,etabar));
   }

Here is the call graph for this function:

Wolf2V_ud()

double Gambit::Spectrum::Wolf2V_ud ( double  l, double  A, double  rhobar, double  etabar  )

static

CKM Wolfenstein (lambda, A, rhobar, etabar) –> V_qq standard parameterisation convertors.

CKM Wolfenstein –> V_ud standard parameterisation convertor.

Definition at line 481 of file spectrum.cpp.

References Gambit::Scanner::pow(), and rhoplusieta().

   {
     double norm = std::norm(rhoplusieta(l,A,rhobar,etabar));
     return 1.0 - 0.5*pow(l,2) - 0.125*pow(l,4) - 0.0625*pow(l,6)*(1.0+8.0*A*A*norm)
            - 0.0078125*pow(l,8)*(5.0-32.0*A*A*norm);
   }

Here is the call graph for this function:

Wolf2V_us()

double Gambit::Spectrum::Wolf2V_us ( double  l, double  A, double  rhobar, double  etabar  )

static

CKM Wolfenstein –> V_us standard parameterisation convertor.

Definition at line 489 of file spectrum.cpp.

References Gambit::Scanner::pow(), and rhoplusieta().

   {
     double norm = std::norm(rhoplusieta(l,A,rhobar,etabar));
     return l - 0.5*A*A*pow(l,7)*norm;
   }

Here is the call graph for this function:

writeSLHAfile()

void Gambit::Spectrum::writeSLHAfile	(	int	slha_version,
		const str &	filename
	)		const

Output spectrum contents as an SLHA file, using getSLHAea.

Definition at line 441 of file spectrum.cpp.

References getSLHAea().

Referenced by drop_SLHAs_if_requested().

   {
      Utils::FileLock mylock(filename);
      mylock.get_lock();
      std::ofstream ofs(filename);
      ofs << getSLHAea(slha_version);
      ofs.close();
      mylock.release_lock();
   }

Here is the call graph for this function:

Here is the caller graph for this function:

Friends And Related Function Documentation

swap

void swap ( Spectrum &  first, Spectrum &  second  )

friend

Friend function: swap resources of two Spectrum objects.

Definition at line 57 of file spectrum.cpp.

Referenced by operator=(), Spectrum(), and Gambit::swap().

   {
       using std::swap; // enable ADL
       swap(first.LE, second.LE);
       swap(first.HE, second.HE);
       swap(first.LE_new, second.LE_new);
       swap(first.HE_new, second.HE_new);
       swap(first.SMINPUTS, second.SMINPUTS);
       swap(first.input_Param, second.input_Param);
       swap(first.mass_cuts, second.mass_cuts);
       swap(first.mass_ratio_cuts, second.mass_ratio_cuts);
       swap(first.initialised, second.initialised);
   }

Member Data Documentation

HE

SubSpectrum* Gambit::Spectrum::HE

private

Definition at line 130 of file spectrum.hpp.

Referenced by clone_HE(), get(), get_HE(), getSLHAea(), has(), PDG_translator(), RunBothToScale(), and Gambit::swap().

HE_new

std::unique_ptr<SubSpectrum> Gambit::Spectrum::HE_new

private

Definition at line 128 of file spectrum.hpp.

Referenced by Gambit::swap().

initialised

bool Gambit::Spectrum::initialised

private

Definition at line 135 of file spectrum.hpp.

Referenced by check_init(), and Gambit::swap().

input_Param

const std::map<str, safe_ptr<const double> >* Gambit::Spectrum::input_Param

private

Definition at line 132 of file spectrum.hpp.

Referenced by Gambit::swap().

LE

SubSpectrum* Gambit::Spectrum::LE

private

Definition at line 129 of file spectrum.hpp.

Referenced by clone_LE(), get(), get_LE(), getSLHAea(), has(), RunBothToScale(), and Gambit::swap().

LE_new

std::unique_ptr<SubSpectrum> Gambit::Spectrum::LE_new

private

Variables.

Definition at line 127 of file spectrum.hpp.

Referenced by Gambit::swap().

mass_cuts

const mc_info* Gambit::Spectrum::mass_cuts

private

Definition at line 133 of file spectrum.hpp.

Referenced by check_mass_cuts(), and Gambit::swap().

mass_ratio_cuts

const mr_info* Gambit::Spectrum::mass_ratio_cuts

private

Definition at line 134 of file spectrum.hpp.

Referenced by check_mass_cuts(), and Gambit::swap().

SMINPUTS

SMInputs Gambit::Spectrum::SMINPUTS

private

Definition at line 131 of file spectrum.hpp.

Referenced by get_SMInputs(), getSLHAea(), and Gambit::swap().

---

The documentation for this class was generated from the following files:

• Elements/include/gambit/Elements/spectrum.hpp
• Elements/src/spectrum.cpp