gambit is hosted by Hepforge, IPPP Durham
GAMBIT  v1.5.0-2191-ga4742ac
a Global And Modular Bsm Inference Tool
Functions
Gambit::DecayBit::MSSM_H Namespace Reference

Functions

double lambda (double x, double y, double z)
 
double gamma_h_chi_gL (std::array< double, 2 > m, double gL, double mh=125., double mw=80.385, double GF=1.1663787e-5)
 
double gamma_h_chi_0 (int i, int j, std::array< double, 4 > m, std::array< std::array< double, 4 >, 4 > Z, double alpha, double mh=125., double mw=80.385, double GF=1.1663787e-5, double sw2=0.22)
 
double gamma_h_chi_pm (int i, int j, std::array< double, 2 > m, std::array< std::array< double, 2 >, 2 > U, std::array< std::array< double, 2 >, 2 > V, double alpha, double mh=125., double mw=80.385, double GF=1.1663787e-5)
 
double gamma_h_chi (std::array< double, 2 > m_pm, std::array< double, 4 > m_0, std::array< std::array< double, 2 >, 2 > U, std::array< std::array< double, 2 >, 2 > V, std::array< std::array< double, 4 >, 4 > Z, double alpha, double mh=125., double mw=80.385, double GF=1.1663787e-5, double sw2=0.22)
 

Function Documentation

◆ gamma_h_chi()

double Gambit::DecayBit::MSSM_H::gamma_h_chi ( std::array< double, 2 >  m_pm,
std::array< double, 4 >  m_0,
std::array< std::array< double, 2 >, 2 >  U,
std::array< std::array< double, 2 >, 2 >  V,
std::array< std::array< double, 4 >, 4 >  Z,
double  alpha,
double  mh = 125.,
double  mw = 80.385,
double  GF = 1.1663787e-5,
double  sw2 = 0.22 
)

Lightest Higgs boson decay to neutralinos and charginos at tree-level in GeV

Warning
Tree-level formula
Returns
$\Gamma(h \to \chi\chi)$ in GeV
Parameters
m_pmChargino masses
m_0Neutralino mases with phases
UReal chargino mixing matrix
VReal chargino mixing matrix
ZReal neutralino mixing matrix
alpha$\alpha$, Higgs mixing angle
mhLightest Higgs mass, $m_h$
mwW-boson mass, $M_W$
GFFermi constant, $G_F$
sw2Weinberg angle, $\sin^2\theta_W$

Definition at line 174 of file MSSM_H.hpp.

References gamma_h_chi_0(), gamma_h_chi_pm(), mh, and mw.

Referenced by Gambit::DecayBit::MSSM_inv_Higgs_BF().

183  {
202  double gamma = 0.;
203 
204  for (int i = 0; i <= 3; i += 1) {
205  // Do not double count e.g. 12 and 21 - they are not distinct
206  for (int j = i; j <= 3; j += 1) {
207  gamma += MSSM_H::gamma_h_chi_0(i, j, m_0, Z, alpha, mh, mw, GF, sw2);
208  }
209  }
210 
211  for (int i = 0; i <= 1; i += 1) {
212  // Do count e.g. 12 and 21 - they are distinct
213  for (int j = 0; j <= 1; j += 1) {
214  gamma += MSSM_H::gamma_h_chi_pm(i, j, m_pm, U, V, alpha, mh, mw, GF);
215  }
216  }
217 
218  return gamma;
219  }
mw
Spectrum Spectrum Spectrum mw
Definition: PrecisionBit_rollcall.hpp:114
alpha
START_MODEL alpha
Definition: HaloModels_gNFW.hpp:30
Gambit::DecayBit::MSSM_H::gamma_h_chi_0
double gamma_h_chi_0(int i, int j, std::array< double, 4 > m, std::array< std::array< double, 4 >, 4 > Z, double alpha, double mh=125., double mw=80.385, double GF=1.1663787e-5, double sw2=0.22)
Definition: MSSM_H.hpp:90
mh
Spectrum Spectrum Spectrum Spectrum Spectrum Spectrum mh
Definition: PrecisionBit_rollcall.hpp:120
Gambit::DecayBit::MSSM_H::gamma_h_chi_pm
double gamma_h_chi_pm(int i, int j, std::array< double, 2 > m, std::array< std::array< double, 2 >, 2 > U, std::array< std::array< double, 2 >, 2 > V, double alpha, double mh=125., double mw=80.385, double GF=1.1663787e-5)
Definition: MSSM_H.hpp:136
Here is the call graph for this function:
Here is the caller graph for this function:

◆ gamma_h_chi_0()

double Gambit::DecayBit::MSSM_H::gamma_h_chi_0 ( int  i,
int  j,
std::array< double, 4 >  m,
std::array< std::array< double, 4 >, 4 >  Z,
double  alpha,
double  mh = 125.,
double  mw = 80.385,
double  GF = 1.1663787e-5,
double  sw2 = 0.22 
)

Lightest Higgs boson decay to neutralinos at tree-level in GeV

Warning
Arguments $i$ and $j$ are zero-based
Tree-level formula
Returns
$\Gamma(h \to \chi_i \chi_j)$ in GeV
Parameters
iNeutralino $\chi_i$ in final state
jNeutralino $\chi_j$ in final state
mNeutralino mases with phases
ZReal neutralino mixing matrix
alpha$\alpha$, Higgs mixing angle
mhLightest Higgs mass, $m_h$
mwW-boson mass, $M_W$
GFFermi constant, $G_F$
sw2Weinberg angle, $\sin^2\theta_W$

Definition at line 90 of file MSSM_H.hpp.

References gamma_h_chi_gL(), mh, and mw.

Referenced by gamma_h_chi(), and Gambit::DecayBit::MSSM_inv_Higgs_BF().

98  {
116  // Weinberg angle
117  const double tw = std::sqrt(sw2) / std::sqrt(1. - sw2);
118 
119  // Eq. 1.113
120  const double e2 = -std::sin(alpha);
121  const double d2 = -std::cos(alpha);
122 
123  // Eq. 1.112 (without \f$\sin\theta_W\f$ which cancels)
124  const double gL = 0.5 * (
125  (Z[j][1] - tw * Z[j][0]) * (e2 * Z[i][2] + d2 * Z[i][3]) +
126  (Z[i][1] - tw * Z[i][0]) * (e2 * Z[j][2] + d2 * Z[j][3]));
127 
128  std::array<double, 2> mf{{m[i], m[j]}};
129  double gamma = gamma_h_chi_gL(mf, gL, mh, mw, GF);
130  if (i == j) {
131  gamma *= 0.5;
132  }
133  return gamma;
134  }
mw
Spectrum Spectrum Spectrum mw
Definition: PrecisionBit_rollcall.hpp:114
alpha
START_MODEL alpha
Definition: HaloModels_gNFW.hpp:30
mh
Spectrum Spectrum Spectrum Spectrum Spectrum Spectrum mh
Definition: PrecisionBit_rollcall.hpp:120
Gambit::DecayBit::MSSM_H::gamma_h_chi_gL
double gamma_h_chi_gL(std::array< double, 2 > m, double gL, double mh=125., double mw=80.385, double GF=1.1663787e-5)
Definition: MSSM_H.hpp:53
Here is the call graph for this function:
Here is the caller graph for this function:

◆ gamma_h_chi_gL()

double Gambit::DecayBit::MSSM_H::gamma_h_chi_gL ( std::array< double, 2 >  m,
double  gL,
double  mh = 125.,
double  mw = 80.385,
double  GF = 1.1663787e-5 
)

Eq. 2.56 for charginos or neutralinos

Returns
$\Gamma(h \to \chi_i \chi_j)$ in GeV
Parameters
mMasses of final state particles
gLCoupling $g_L$ in Eq. 1.112
mhLightest Higgs mass, $m_h$
mwW-boson mass, $M_W$
GFFermi constant, $G_F$

Definition at line 53 of file MSSM_H.hpp.

References lambda(), mh, mw, and Gambit::Scanner::pow().

Referenced by gamma_h_chi_0(), and gamma_h_chi_pm().

57  {
67  const double eps_2 = 1.;
68  const double gR = gL * eps_2;
69 
70  // Phase-space
71 
72  if (std::fabs(m[0]) + std::fabs(m[1]) >= mh) {
73  return 0.;
74  }
75 
76  const double l = lambda(pow(m[0], 2), pow(m[1], 2), pow(mh, 2));
77  if (l <= 0.) {
78  return 0.;
79  }
80 
81  // Eq. 2.56 without common factor (without f$\sin\theta_W\f$ which cancels)
82  const double gamma_no_prefactor = std::sqrt(l) * ((
83  pow(gL, 2) + pow(gR, 2)) * (1. - (pow(m[0], 2) + pow(m[1], 2)) / pow(mh, 2))
84  -4. * gL * gR * m[0] * m[1] / pow(mh, 2));
85 
86  // Eq. 2.56 with common factor
87  return GF * pow(mw, 2) / (2. * M_SQRT2 * M_PI) * mh * gamma_no_prefactor;
88  }
mw
Spectrum Spectrum Spectrum mw
Definition: PrecisionBit_rollcall.hpp:114
Gambit::DecayBit::MSSM_H::lambda
double lambda(double x, double y, double z)
Definition: MSSM_H.hpp:38
Gambit::Scanner::pow
double pow(const double &a)
Outputs a^i.
Definition: scanner_utils.hpp:310
mh
Spectrum Spectrum Spectrum Spectrum Spectrum Spectrum mh
Definition: PrecisionBit_rollcall.hpp:120
Here is the call graph for this function:
Here is the caller graph for this function:

◆ gamma_h_chi_pm()

double Gambit::DecayBit::MSSM_H::gamma_h_chi_pm ( int  i,
int  j,
std::array< double, 2 >  m,
std::array< std::array< double, 2 >, 2 >  U,
std::array< std::array< double, 2 >, 2 >  V,
double  alpha,
double  mh = 125.,
double  mw = 80.385,
double  GF = 1.1663787e-5 
)

Lightest Higgs boson decay to charginos at tree-level in GeV

Warning
Arguments $i$ and $j$ are zero-based
Tree-level formula
Returns
$\Gamma(h \to \chi^-_i \chi^+_j)$ in GeV
Parameters
iNegative chargino $\chi^-_i$ in final state
jPositive chargino $\chi^+_j$ in final state
mChargino masses
UReal chargino mixing matrix
VReal chargino mixing matrix
alpha$\alpha$, Higgs mixing angle
mhLightest Higgs mass, $m_h$
mwW-boson mass, $M_W$
GFFermi constant, $G_F$

Definition at line 136 of file MSSM_H.hpp.

References gamma_h_chi_gL(), mh, and mw.

Referenced by gamma_h_chi().

144  {
162  // Eq. 1.113
163  const double e2 = -std::sin(alpha);
164  const double d2 = -std::cos(alpha);
165 
166  // Eq. 1.111 (without \f$\sin\theta_W\f$ which cancels)
167  const double gL = M_SQRT1_2 *
168  (e2 * V[j][0] * U[i][1] - d2 * V[j][1] * U[i][0]);
169 
170  std::array<double, 2> mf{{m[i], m[j]}};
171  return gamma_h_chi_gL(mf, gL, mh, mw, GF);
172  }
mw
Spectrum Spectrum Spectrum mw
Definition: PrecisionBit_rollcall.hpp:114
alpha
START_MODEL alpha
Definition: HaloModels_gNFW.hpp:30
mh
Spectrum Spectrum Spectrum Spectrum Spectrum Spectrum mh
Definition: PrecisionBit_rollcall.hpp:120
Gambit::DecayBit::MSSM_H::gamma_h_chi_gL
double gamma_h_chi_gL(std::array< double, 2 > m, double gL, double mh=125., double mw=80.385, double GF=1.1663787e-5)
Definition: MSSM_H.hpp:53
Here is the call graph for this function:
Here is the caller graph for this function:

◆ lambda()

double Gambit::DecayBit::MSSM_H::lambda ( double  x,
double  y,
double  z 
)

Phase-space function (similar to Kallen function)

\[ \lambda(x, y; z) = (1 - x / z - y / z)^2 - 4 x y / z^2 \]

Definition at line 38 of file MSSM_H.hpp.

References Gambit::Scanner::pow().

Referenced by Gambit::DecayBit::DiracSingletDM_Higgs_decays(), gamma_h_chi_gL(), Gambit::DecayBit::MajoranaSingletDM_Higgs_decays(), Gambit::SpecBit::run_lambda(), Gambit::DecayBit::ScalarSingletDM_Higgs_decays(), scanner_plugin(), Gambit::DecayBit::VectorSingletDM_Higgs_decays(), and Gambit::FlavBit::Vertices::VuL().

38  {
46  if (z == 0.) {
47  throw std::invalid_argument("z must be greater than zero");
48  }
49 
50  return pow(1. - x / z - y / z, 2) - 4. * x * y / pow(z, 2);
51  }
Gambit::Scanner::pow
double pow(const double &a)
Outputs a^i.
Definition: scanner_utils.hpp:310
Here is the call graph for this function:
Here is the caller graph for this function: