w = 0.1 * rho # calculate bin width (in force)
force_mean = alpha - gamma_e * rho
theory = Histogram()
- theory.bin_edges = arange(start=0, stop=3*force_mean, step=w)
+ theory.bin_edges = arange(start=0, stop=max(force_mean,0)+3*rho, step=w)
theory.bin_centers = theory.bin_edges[:-1] + w/2
theory.counts = w*num_domains*N*probability_distribution(
theory.bin_centers, [alpha, rho])