namespace Gambit::DecayBit #

[No description available]

Namespaces #

Functions #

Functions Documentation #

function sqrt_lambda #

double sqrt_lambda(
    double a,
    double b,
    double c
)

Square root of the standard kinematic function lambda(a,b,c)

function BW #

std::complex< double > BW(
    double q2,
    double m2,
    double imag_term
)

Breit-Wigner pole (complex)

function check_width #

void check_width(
    const str & info,
    double & w,
    bool raise_invalid_pt_negative_width =false,
    bool raise_invalid_pt_large_width =false
)

Check if a width is negative or suspiciously large and raise an error.

function compute_SM_higgs_decays #

void compute_SM_higgs_decays(
    DecayTable::Entry & result,
    double mh
)

Populate SM Higgs decay channels for a higgs mass of m_h.

function set_FH_neutral_h_decay #

void set_FH_neutral_h_decay(
    DecayTable::Entry & result,
    int iH,
    const fh_Couplings_container & FH_input,
    const mass_es_pseudonyms & psn,
    bool invalidate,
    bool SM
)

Set neutral h decays computed by FeynHiggs.

function combine_decay_table_entries #

void combine_decay_table_entries(
    DecayTable::Entry & entry1,
    DecayTable::Entry & entry2,
    DecayTable::Entry & result
)

Function combining two Decaytable entries with different decay channels of the same particle into a single entry.

function W_plus_decays #

void W_plus_decays(
    DecayTable::Entry & result
)

SM decays: W+ Reference: 2017 PDG

function Z_decays #

void Z_decays(
    DecayTable::Entry & result
)

SM decays: Z Reference: 2019 PDG

function t_decays #

void t_decays(
    DecayTable::Entry & result
)

SM decays: t Reference: 2017 PDG

function mu_plus_decays #

void mu_plus_decays(
    DecayTable::Entry & result
)

SM decays: mu+ Reference: 2017 PDG

function tau_plus_decays #

void tau_plus_decays(
    DecayTable::Entry & result
)

SM decays: tau+ Reference: 2017 PDG

function pi_0_decays #

void pi_0_decays(
    DecayTable::Entry & result
)

SM decays: pi0 Reference: 2017 PDG

function pi_plus_decays #

void pi_plus_decays(
    DecayTable::Entry & result
)

SM decays: pi+ Reference: 2017 PDG

function eta_decays #

void eta_decays(
    DecayTable::Entry & result
)

SM decays: eta Reference: 2017 PDG

function rho_0_decays #

void rho_0_decays(
    DecayTable::Entry & result
)

SM decays: rho0 Reference: 2017 PDG

function rho_plus_decays #

void rho_plus_decays(
    DecayTable::Entry & result
)

SM decays: rho+ Reference: 2017 PDG

function omega_decays #

void omega_decays(
    DecayTable::Entry & result
)

SM decays: omega Reference: 2017 PDG

function rho1450_decays #

void rho1450_decays(
    DecayTable::Entry & result
)

SM decays: rho1450 Reference: 2017 PDG

function W_minus_decays #

void W_minus_decays(
    DecayTable::Entry & result
)

SM decays: conjugates

function tbar_decays #

void tbar_decays(
    DecayTable::Entry & result
)

function mu_minus_decays #

void mu_minus_decays(
    DecayTable::Entry & result
)

function tau_minus_decays #

void tau_minus_decays(
    DecayTable::Entry & result
)

function pi_minus_decays #

void pi_minus_decays(
    DecayTable::Entry & result
)

function rho_minus_decays #

void rho_minus_decays(
    DecayTable::Entry & result
)

function Ref_SM_Higgs_decays_table #

void Ref_SM_Higgs_decays_table(
    DecayTable::Entry & result
)

Reference SM Higgs decays.

These functions are given a different capability to regular decay functions, to allow other module functions to specifically depend on the SM values for reference, even when scanning another model. Reference SM Higgs decays from LHCHiggsXSWG: most SM-like Higgs

function Ref_SM_other_Higgs_decays_table #

void Ref_SM_other_Higgs_decays_table(
    DecayTable::Entry & result
)

Reference SM Higgs decays from LHCHiggsXSWG: least SM-like Higgs.

function Ref_SM_A0_decays_table #

void Ref_SM_A0_decays_table(
    DecayTable::Entry & result
)

Reference SM Higgs decays from LHCHiggsXSWG: A0.

function Ref_SM_Higgs_decays_FeynHiggs #

void Ref_SM_Higgs_decays_FeynHiggs(
    DecayTable::Entry & result
)

Reference SM Higgs decays from FeynHiggs: h0_1.

function Ref_SM_other_Higgs_decays_FH #

void Ref_SM_other_Higgs_decays_FH(
    DecayTable::Entry & result
)

Reference SM Higgs decays from FeynHiggs: h0_2.

function Ref_SM_A0_decays_FH #

void Ref_SM_A0_decays_FH(
    DecayTable::Entry & result
)

Reference SM Higgs decays from FeynHiggs: A0.

function H_minus_decays #

void H_minus_decays(
    DecayTable::Entry & result
)

MSSM decays: conjugates

function stopbar_1_decays #

void stopbar_1_decays(
    DecayTable::Entry & result
)

function stopbar_2_decays #

void stopbar_2_decays(
    DecayTable::Entry & result
)

function sbottombar_1_decays #

void sbottombar_1_decays(
    DecayTable::Entry & result
)

function sbottombar_2_decays #

void sbottombar_2_decays(
    DecayTable::Entry & result
)

function supbar_l_decays #

void supbar_l_decays(
    DecayTable::Entry & result
)

function supbar_r_decays #

void supbar_r_decays(
    DecayTable::Entry & result
)

function sdownbar_l_decays #

void sdownbar_l_decays(
    DecayTable::Entry & result
)

function sdownbar_r_decays #

void sdownbar_r_decays(
    DecayTable::Entry & result
)

function scharmbar_l_decays #

void scharmbar_l_decays(
    DecayTable::Entry & result
)

function scharmbar_r_decays #

void scharmbar_r_decays(
    DecayTable::Entry & result
)

function sstrangebar_l_decays #

void sstrangebar_l_decays(
    DecayTable::Entry & result
)

function sstrangebar_r_decays #

void sstrangebar_r_decays(
    DecayTable::Entry & result
)

function selectronbar_l_decays #

void selectronbar_l_decays(
    DecayTable::Entry & result
)

function selectronbar_r_decays #

void selectronbar_r_decays(
    DecayTable::Entry & result
)

function smuonbar_l_decays #

void smuonbar_l_decays(
    DecayTable::Entry & result
)

function smuonbar_r_decays #

void smuonbar_r_decays(
    DecayTable::Entry & result
)

function staubar_1_decays #

void staubar_1_decays(
    DecayTable::Entry & result
)

function staubar_2_decays #

void staubar_2_decays(
    DecayTable::Entry & result
)

function snubar_electronl_decays #

void snubar_electronl_decays(
    DecayTable::Entry & result
)

function snubar_muonl_decays #

void snubar_muonl_decays(
    DecayTable::Entry & result
)

function snubar_taul_decays #

void snubar_taul_decays(
    DecayTable::Entry & result
)

function chargino_minus_1_decays #

void chargino_minus_1_decays(
    DecayTable::Entry & result
)

function chargino_minus_2_decays #

void chargino_minus_2_decays(
    DecayTable::Entry & result
)

function CP_conjugate #

DecayTable::Entry CP_conjugate(
    const DecayTable::Entry & in
)

Construct a decay table entry for a particle from the entry for its antiparticle.

Construct a decay table entry for a particle, from the entry for its antiparticle.

function SM_Higgs_decays #

void SM_Higgs_decays(
    DecayTable::Entry & result
)

SM decays: Higgs.

function FeynHiggs_t_decays #

void FeynHiggs_t_decays(
    DecayTable::Entry & result
)

FeynHiggs MSSM decays: t Reference for total width: 2017 PDG

function FeynHiggs_MSSM_h0_1_decays #

void FeynHiggs_MSSM_h0_1_decays(
    DecayTable::Entry & result
)

FeynHiggs MSSM decays: h0_1.

function FeynHiggs_h0_2_decays #

void FeynHiggs_h0_2_decays(
    DecayTable::Entry & result
)

FeynHiggs MSSM decays: h0_2.

function FeynHiggs_A0_decays #

void FeynHiggs_A0_decays(
    DecayTable::Entry & result
)

FeynHiggs MSSM decays: A0.

function FeynHiggs_H_plus_decays #

void FeynHiggs_H_plus_decays(
    DecayTable::Entry & result
)

FeynHiggs MSSM decays: H+.

function MSSM_h0_1_decays #

void MSSM_h0_1_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: h0_1.

function h0_2_decays #

void h0_2_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: h0_2.

function A0_decays #

void A0_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: A0.

function H_plus_decays #

void H_plus_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: H_plus.

function gluino_decays #

void gluino_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: gluino.

function stop_1_decays #

void stop_1_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: stop_1.

function stop_2_decays #

void stop_2_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: stop_2.

function sbottom_1_decays #

void sbottom_1_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: sbottom_1.

function sbottom_2_decays #

void sbottom_2_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: sbottom_2.

function sup_l_decays #

void sup_l_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: sup_l.

function sup_r_decays #

void sup_r_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: sup_r.

function sdown_l_decays #

void sdown_l_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: sdown_l.

function sdown_r_decays #

void sdown_r_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: sdown_r.

function scharm_l_decays #

void scharm_l_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: scharm_l.

function scharm_r_decays #

void scharm_r_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: scharm_r.

function sstrange_l_decays #

void sstrange_l_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: sstrange_l.

function sstrange_r_decays #

void sstrange_r_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: sstrange_r.

function selectron_l_decays #

void selectron_l_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: selectron_l.

function selectron_r_decays #

void selectron_r_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: selectron_r.

function smuon_l_decays #

void smuon_l_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: smuon_l.

function smuon_r_decays #

void smuon_r_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: smuon_r.

function stau_1_decays #

void stau_1_decays(
    DecayTable::Entry & result
)

MSSM decays: stau_1 (Uses SUSY-HIT results or dedicated DecayBit calculation for small mass splittings)

function stau_1_decays_SH #

void stau_1_decays_SH(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: stau_1.

function stau_2_decays #

void stau_2_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: stau_2.

function snu_electronl_decays #

void snu_electronl_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: snu_electronl.

function snu_muonl_decays #

void snu_muonl_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: snu_muonl.

function snu_taul_decays #

void snu_taul_decays(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: snu_taul Note that SUSY-HIT calls ~nu_tau_L “snutau1” even though it has no RH (~)nus.

function chargino_plus_1_decays #

void chargino_plus_1_decays(
    DecayTable::Entry & result
)

MSSM decays: chargino_plus_1 without decays to gravitinos (Uses SUSY-HIT results or dedicated DecayBit calculation for small mass splittings) Aggregator function providing chargino_plus_1_decay_rates for MSSM models without gravitinos

function chargino_plus_2_decays #

void chargino_plus_2_decays(
    DecayTable::Entry & result
)

MSSM decays: chargino_plus_2 without decays to gravitinos (Uses SUSY-HIT results) Aggregator function providing chargino_plus_2_decay_rates for MSSM models without gravitinos

function chargino_plus_1_decays_all #

void chargino_plus_1_decays_all(
    DecayTable::Entry & result
)

MSSM decays: chargino_plus_1 incl. decays to gravitinos (Uses SUSY-HIT results and DecayBit calculations for small mass splittings and decays to gravitinos) Aggregator function providing chargino_plus_1_decay_rates for MSSM models with gravitinos

function chargino_plus_2_decays_all #

void chargino_plus_2_decays_all(
    DecayTable::Entry & result
)

MSSM decays: chargino_plus_2 incl. decays to gravitinos (Uses SUSY-HIT results and DecayBit calculations for decays to gravitinos) Aggregator function providing chargino_plus_2_decay_rates for MSSM models with gravitinos

function chargino_plus_1_decays_SH_or_smallsplit #

void chargino_plus_1_decays_SH_or_smallsplit(
    DecayTable::Entry & result
)

MSSM decays: chargino_plus_1 (Uses SUSY-HIT results or dedicated DecayBit calculation for small mass splittings)

function chargino_plus_1_decays_SH #

void chargino_plus_1_decays_SH(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: chargino_plus_1.

function chargino_plus_2_decays_SH #

void chargino_plus_2_decays_SH(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: chargino_plus_2.

function neutralino_1_decays #

void neutralino_1_decays(
    DecayTable::Entry & result
)

MSSM decays: neutralino_1 without decays to gravitinos (Uses SUSY-HIT results) Aggregator function providing neutralino_1_decay_rates for MSSM models without gravitinos

function neutralino_2_decays #

void neutralino_2_decays(
    DecayTable::Entry & result
)

MSSM decays: neutralino_2 without decays to gravitinos (Uses SUSY-HIT results) Aggregator function providing neutralino_2_decay_rates for MSSM models without gravitinos

function neutralino_3_decays #

void neutralino_3_decays(
    DecayTable::Entry & result
)

MSSM decays: neutralino_3 without decays to gravitinos (Uses SUSY-HIT results) Aggregator function providing neutralino_3_decay_rates for MSSM models without gravitinos

function neutralino_4_decays #

void neutralino_4_decays(
    DecayTable::Entry & result
)

MSSM decays: neutralino_4 without decays to gravitinos (Uses SUSY-HIT results) Aggregator function providing neutralino_4_decay_rates for MSSM models without gravitinos

function neutralino_1_decays_all #

void neutralino_1_decays_all(
    DecayTable::Entry & result
)

MSSM decays: neutralino_1 incl. decays to gravitinos (Uses SUSY-HIT results and dedicated DecayBit calculation for decays to gravitinos) Aggregator function providing neutralino_1_decay_rates for MSSM models with gravitinos

function neutralino_2_decays_all #

void neutralino_2_decays_all(
    DecayTable::Entry & result
)

MSSM decays: neutralino_2 incl. decays to gravitinos (Uses SUSY-HIT results and dedicated DecayBit calculation for decays to gravitinos) Aggregator function providing neutralino_2_decay_rates for MSSM models with gravitinos

function neutralino_3_decays_all #

void neutralino_3_decays_all(
    DecayTable::Entry & result
)

MSSM decays: neutralino_3 incl. decays to gravitinos (Uses SUSY-HIT results and dedicated DecayBit calculation for decays to gravitinos) Aggregator function providing neutralino_3_decay_rates for MSSM models with gravitinos

function neutralino_4_decays_all #

void neutralino_4_decays_all(
    DecayTable::Entry & result
)

MSSM decays: neutralino_4 incl. decays to gravitinos (Uses SUSY-HIT results and dedicated DecayBit calculation for decays to gravitinos) Aggregator function providing neutralino_4_decay_rates for MSSM models with gravitinos

function neutralino_1_decays_SH #

void neutralino_1_decays_SH(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: neutralino_1.

function neutralino_2_decays_SH #

void neutralino_2_decays_SH(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: neutralino_2.

function neutralino_3_decays_SH #

void neutralino_3_decays_SH(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: neutralino_3.

function neutralino_4_decays_SH #

void neutralino_4_decays_SH(
    DecayTable::Entry & result
)

SUSY-HIT MSSM decays: neutralino_4.

function chargino_plus_1_decays_smallsplit #

void chargino_plus_1_decays_smallsplit(
    DecayTable::Entry & result
)

MSSM decays: chargino decays for small chargino–neutralino mass splitting. Using results from hep-ph/9607421.

function stau_1_decays_smallsplit #

void stau_1_decays_smallsplit(
    DecayTable::Entry & result
)

MSSM decays: stau decays for small stau–neutralino mass splitting. Using results from T. Jittoh, J. Sato, T. Shimomura, M. Yamanaka, Phys. Rev. D73 (2006), hep-ph/0512197

function neutralino_decays_gravitino #

DecayTable::Entry neutralino_decays_gravitino(
    int n_Neu,
    const Spectrum & spec,
    const DecayTable::Entry & Z_decays
)

MSSM decays: neutralino decays to gravitinos. Using results from hep-ph/0404231v2 and 0908.3399. Includes 3-body decays via the Z boson to gravitinos and SM fermion pairs, for all of the latter. Using results from http://www-library.desy.de/preparch/desy/thesis/desy-thesis-09-016.pdf. Complete analytic formulas, not just the light-gravitino limit.

Neutralino three-body decays to gravitinos and fermion pairs, for all of the latter. Using results from http://www-library.desy.de/preparch/desy/thesis/desy-thesis-09-016.pdf. Complete analytic formulas, not just the light-gravitino limit.

function neutralino_1_decays_gravitino #

void neutralino_1_decays_gravitino(
    DecayTable::Entry & result
)

MSSM decays: neutralino_1 decays to gravitinos.

function neutralino_2_decays_gravitino #

void neutralino_2_decays_gravitino(
    DecayTable::Entry & result
)

MSSM decays: neutralino_2 decays to gravitinos.

function neutralino_3_decays_gravitino #

void neutralino_3_decays_gravitino(
    DecayTable::Entry & result
)

MSSM decays: neutralino_3 decays to gravitinos.

function neutralino_4_decays_gravitino #

void neutralino_4_decays_gravitino(
    DecayTable::Entry & result
)

MSSM decays: neutralino_4 decays to gravitinos.

function chargino_decays_gravitino #

DecayTable::Entry chargino_decays_gravitino(
    int n_Cha,
    const Spectrum & spec,
    const DecayTable::Entry & W_decays
)

MSSM decays: chargino decays to gravitinos. Using results from hep-ph/9605398. The decay width formulas are valid in the light-gravitino limit.

function chargino_plus_1_decays_gravitino #

void chargino_plus_1_decays_gravitino(
    DecayTable::Entry & result
)

MSSM decays: chargino_1 decay to gravitinos.

function chargino_plus_2_decays_gravitino #

void chargino_plus_2_decays_gravitino(
    DecayTable::Entry & result
)

MSSM decays: chargino_2 decay to gravitinos.

function ScalarSingletDM_Higgs_decays #

void ScalarSingletDM_Higgs_decays(
    DecayTable::Entry & result
)

Add the decay of Higgs to singlets for the ScalarSingletDM models.

function VectorSingletDM_Higgs_decays #

void VectorSingletDM_Higgs_decays(
    DecayTable::Entry & result
)

Add the decay of Higgs to vectors for the VectorSingletDM models (see arXiv:1512.06458v4)

function MajoranaSingletDM_Higgs_decays #

void MajoranaSingletDM_Higgs_decays(
    DecayTable::Entry & result
)

Add the decay of Higgs to Majorana fermions for the MajoranaSingletDM models (see arXiv:1512.06458v4)

function DiracSingletDM_Higgs_decays #

void DiracSingletDM_Higgs_decays(
    DecayTable::Entry & result
)

Add the decay of Higgs to Dirac fermions for the DiracSingletDM models (see arXiv:1512.06458v4)

function SubGeVDM_dark_photon_decays #

void SubGeVDM_dark_photon_decays(
    DecayTable::Entry & result
)

function get_dark_photon_visible_branching #

void get_dark_photon_visible_branching(
    double & result
)

function get_dark_photon_decay_length #

void get_dark_photon_decay_length(
    double & result
)

function CH_DMsimpVectorMedDiracDM_Y1_decays #

void CH_DMsimpVectorMedDiracDM_Y1_decays(
    DecayTable::Entry & result
)

Add the decay of Higgs to vectors for the VectorSingletDM models (see arXiv:1512.06458v4) ///////// Dirac DM Simplified Model /////////////////////.

function CH_DMsimpVectorMedMajoranaDM_Y1_decays #

void CH_DMsimpVectorMedMajoranaDM_Y1_decays(
    DecayTable::Entry & result
)

function CH_DMsimpVectorMedScalarDM_Y1_decays #

void CH_DMsimpVectorMedScalarDM_Y1_decays(
    DecayTable::Entry & result
)

function CH_DMsimpVectorMedVectorDM_Y1_decays #

void CH_DMsimpVectorMedVectorDM_Y1_decays(
    DecayTable::Entry & result
)

function all_decays #

void all_decays(
    DecayTable & decays
)

Collect all the DecayTable entries into an actual DecayTable.

Spit out the full decay table as SLHA1 and SLHA2 files.

Invalidate MSSM points that have a stable charged particle?

Check all particles for negative decay width

Spit out the full decay table as an SLHA file.

function all_decays_from_SLHA #

void all_decays_from_SLHA(
    DecayTable & decays
)

Read an SLHA2 file in and use it to create a GAMBIT DecayTable. Note that creating a DecayTable from an SLHA1 file is not possible at present.

function get_decaytable_as_map #

void get_decaytable_as_map(
    map_str_dbl & map
)

Convert the DecayTable to a format where we can print each individual channel’s BF.

Otherwise just print the specific, named channels

function get_mass_es_pseudonyms #

void get_mass_es_pseudonyms(
    mass_es_pseudonyms & result
)

Get MSSM mass eigenstate pseudonyms for the gauge eigenstates.

function check_first_sec_gen_mixing #

void check_first_sec_gen_mixing(
    int & result
)

SLHA1 mixing check flag This gets set non-zero if first or second generation left-right mixing exceeds the specified tolerance (which means that SLHA1 is an invalid format for this model point). 1 = u, 2 = d, 3 = c, 4 = s, 5 = e, 6 = mu

function get_Higgs_invWidth_chi2 #

daFunk::Funk get_Higgs_invWidth_chi2(
    std::string filename
)

function MSSM_inv_Higgs_BF #

void MSSM_inv_Higgs_BF(
    double & BF
)

Branching fraction for Higgs into pair of lightest neutralinos

Tree-level formulas

Assumes decoupling limit for Higgs mixing angle [ \alpha = \beta - \frac12 \pi ]

Only includes neutralinos, charginos and SM-like width in total width

BF(\textrm{BR}(h\to\chi^0_1\chi^0_1))

function ScalarSingletDM_inv_Higgs_BF #

void ScalarSingletDM_inv_Higgs_BF(
    double & BF
)

Branching fraction for Higgs into scalar singlet DM BF(\textrm{BR}(h\to S S))

function VectorSingletDM_inv_Higgs_BF #

void VectorSingletDM_inv_Higgs_BF(
    double & BF
)

Branching fraction for Higgs into vector singlet DM BF(\textrm{BR}(h\to V V))

function MajoranaSingletDM_inv_Higgs_BF #

void MajoranaSingletDM_inv_Higgs_BF(
    double & BF
)

Branching fraction for Higgs into Majorana singlet DM BF(\textrm{BR}(h\to X X))

function DiracSingletDM_inv_Higgs_BF #

void DiracSingletDM_inv_Higgs_BF(
    double & BF
)

Branching fraction for Higgs into Dirac singlet DM BF(\textrm{BR}(h\to F F))

function lnL_Higgs_invWidth_SMlike #

void lnL_Higgs_invWidth_SMlike(
    double & lnL
)

Log-likelihood for Higgs invisible branching ratio

We assume that the log likelihood as a function of the invisible branching fraction BF can be written in the form: - 2 log L = a * (BF - b)^2

We can explicitly test this assumption for CMS, where the chi2 is given in figure 12 of arXiv:2201.11585. Indeed, the assumed functional form gives a perfect fit for a_CMS = 339 and b_CMS = 0.089.

For ATLAS we extract the fit parameters from the two values stated in arXiv:2202.07953: BR < 0.145 at 95% CL and BR < 0.127 at 90% CL. These values correspond to a_ATLAS = 303 and b_ATLAS = 0.032.

The combined log likelihood gives BR < 0.14 at 95% CL. The best-fit value is BR = 0.06, which is preferred over BR = 0 with a significance of 1.2 sigma.

This typically assumes that the Higgs is otherwise SM-like, i.e., no changes to production cross sections or any other decays.

lnLLog-likelihood for Higgs invisible branching ratio

function lnL_Z_inv #

void lnL_Z_inv(
    double & lnL
)

Log-likelihood from LEP measurements of (Z)-boson invisible width

This is valid for SM, RHN models and for MSSM-like models with Z-boson invisible decays to neutralinos and neutrinos

lnLLog-likelihood

function Z_gamma_nu_2l #

void Z_gamma_nu_2l(
    triplet< double > & gamma
)

Calculate width of (Z) decays to neutrinos (with RHN correction if present), (\Gamma(Z\to\nu\nu)) at two-loop in GeV

This uses input (\alpha(M_Z)) - does not include input light-quark thresholds.

gamma(\Gamma(Z\to\chi\chi))

function Z_gamma_chi_0_MSSM_tree #

void Z_gamma_chi_0_MSSM_tree(
    triplet< double > & gamma
)

Calculate width of (Z) decays to the lightest neutralinos, (\Gamma(Z\to\chi^0_1\chi^0_1)) in GeV

Tree-level formula with 10% error

gamma(\Gamma(Z\to\chi^0_1\chi^0_1))

function RHN_W_to_l_decays #

void RHN_W_to_l_decays(
    std::vector< double > & result
)

function lnL_W_decays_chi2 #

void lnL_W_decays_chi2(
    double & result
)

Updated on 2024-07-18 at 13:53:31 +0000