From 5a9065c64900f048dceded259860b192cde335bc Mon Sep 17 00:00:00 2001 From: Dmitry Kalinkin Date: Mon, 28 Oct 2024 03:02:11 -0400 Subject: [PATCH] tau_plots.py: maxfev=10000 --- benchmarks/insert_tau/analysis/tau_plots.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/benchmarks/insert_tau/analysis/tau_plots.py b/benchmarks/insert_tau/analysis/tau_plots.py index f9f5c86..77c2dd3 100644 --- a/benchmarks/insert_tau/analysis/tau_plots.py +++ b/benchmarks/insert_tau/analysis/tau_plots.py @@ -135,7 +135,7 @@ def gauss(x, A,mu, sigma): fnc=gauss sigma=np.sqrt(y[slc])+(y[slc]==0) - coeff, var_matrix = curve_fit(fnc, list(bc[slc]), list(y[slc]), p0=p0,sigma=list(sigma)) + coeff, var_matrix = curve_fit(fnc, list(bc[slc]), list(y[slc]), p0=p0, sigma=list(sigma), maxfev=10000) res.append(abs(coeff[2])) dres.append(np.sqrt(var_matrix[2][2])) emid.append((emin+emax)/2) @@ -147,7 +147,7 @@ def gauss(x, A,mu, sigma): fnc=lambda E,B:B/E p0=[1,] -coeff, var_matrix = curve_fit(fnc, emid, res, p0=p0,sigma=list(dres)) +coeff, var_matrix = curve_fit(fnc, emid, res, p0=p0, sigma=list(dres), maxfev=10000) xx=np.linspace(10, 100, 100) plt.plot(xx, 1000*fnc(xx, *coeff), label=f"fit: ${coeff[0]:.2f}/E$ mrad") plt.legend()