Convergent Bayesian global fits of 4D composite Higgs models