Main addition is the 1/2 sigma determination from the crossing point (using bisection from boost) Additionally some tweaks to plotting (forcing limits in the custom plot, vertical position for labels)
@zewolffs @brkortma