Histogram any kinematic distribution for the current process.
2→2 cos θ uses deterministic scipy.quad per bin.
2→N and pp processes use Monte Carlo, with
optional NLO running-coupling rescaling. Bin axis ranges are
automatically chosen for the selected observable.