fit1, covar1 = curve_fit(gaus, centers, bins[0], p0=[0.1, 0.2])
err1 = np.sqrt(np.diag(covar1))
fit2, covar2 = curve_fit(gaus, centers, bins[0], p0=[1, 0.1])
err2 = np.sqrt(np.diag(covar2))