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)