Module flavio.physics.higgs.signalstrength

Functions for Higgs signal strengths.

Functions

def higgs_signalstrength(wc_obj, par, name_prod, name_dec)

Expand source code

def higgs_signalstrength(wc_obj, par, name_prod, name_dec):
    scale = flavio.config['renormalization scale']['hdecays']
    C = wc_obj.get_wcxf(sector='all', scale=scale, par=par,
                        eft='SMEFT', basis='Warsaw')
    f_prod = getattr(production, name_prod)
    f_dec = getattr(decay, name_dec)
    return f_prod(C) * f_dec(C) / width.Gamma_h(par,  C)

def make_obs_higgs(name_prod, name_dec)

Expand source code

def make_obs_higgs(name_prod, name_dec):
    d_dec = decay_modes[name_dec]
    d_prod = prod_modes[name_prod]
    process_tex = r"h \to {}".format(d_dec['tex'])
    process_taxonomy = r'Process :: Higgs production and decay :: $' + d_dec['tex_class'] + r'$ :: $' + process_tex + r"$"
    obs_name = "mu_{}(h->{})".format(d_prod['str'], d_dec['str'])
    obs = flavio.classes.Observable(obs_name)
    obs.set_description(r"Signal strength of $" + process_tex + r"$ from " + d_prod['desc'])
    obs.tex = r"$\mu_{" + d_prod['tex'] + r"}(" + process_tex + r")$"
    obs.add_taxonomy(process_taxonomy)

    def obs_fct(wc_obj, par):
        return higgs_signalstrength(wc_obj, par, name_prod, name_dec)

    flavio.classes.Prediction(obs_name, obs_fct)