def test_histogram_fits(histogram_file): """ Write to fits,read back, and check """ hist = Histogram(nbins=[5, 11], ranges=[[-2.5, 2.5], [-1, 1]]) hist.fill(np.array([[0, 0], [0, 0.5]])) hist.to_fits().writeto(histogram_file, overwrite=True) newhist = Histogram.from_fits(histogram_file) # check that the values are the same compare_histograms(hist, newhist)