Constrained fit

Gaussian fit to histogram data where the prefactor is frozen based on the histogram normalization.

Gaussian
χ²: 101.81 / 82 dof
amplitude: 10101 (frozen)
mu: -0.0018±0.0035
sigma: 0.9995±0.0025 (min=0.0)

<matplotlib.legend.Legend object at 0x7fb7d210b4d0>

import numpy as np

from aptapy.hist import Histogram1d
from aptapy.models import Gaussian
from aptapy.plotting import plt

hist = Histogram1d(np.linspace(-5., 5., 100), label="Random data", xlabel="z")
hist.fill(np.random.default_rng().normal(size=100000))
hist.plot()

model = Gaussian()
# Fix the amplitude. This is a generally useful technique, as you should
# never fit the normalization of a histogram.
model.amplitude.freeze(hist.area() / model.integral(-5., 5.))
model.fit(hist)
print(model)
model.plot(fit_output=True)

plt.legend()