Relative difference

Compare two 1D histograms using the relative difference [(h1-h2)/h2].

from plothist_utils import get_dummy_data

df = get_dummy_data()

name = "variable_1"

x1 = df[name][df["category"] == 2]
x2 = df[name][df["category"] == 3]

x_range = (min(*x1, *x2), max(*x1, *x2))

from plothist import make_hist

h1 = make_hist(x1, bins=50, range=x_range)
h2 = make_hist(x2, bins=50, range=x_range)

###
from plothist import plot_two_hist_comparison

fig, ax_main, ax_comparison = plot_two_hist_comparison(
    h1,
    h2,
    xlabel=name,
    ylabel="Entries",
    h1_label=r"$\mathbf{H\,\,1}$",
    h2_label=r"$\mathbf{H\,\,2}$",
    comparison="relative_difference",  # <--
)

fig.savefig("1d_comparison_relative_difference.svg", bbox_inches="tight")