zebra/pyzebra
3
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Update fit2.py

Browse Source 
corrected value range for numerical fit, added export_fit parameter, which should in future decide if exported value is fitted or integrated one.
This commit is contained in:
JakHolzer 2020-09-14 13:43:56 +02:00 committed by Ivan Usov
parent 2788c77fa1
commit b157c2a3ae
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
pyzebra/fit2.py
View File

@ -122,8 +122,8 @@ def fitccl(
num_int_area = simps(y[numfit_min:numfit_max], x[numfit_min:numfit_max]) num_int_area = simps(y[numfit_min:numfit_max], x[numfit_min:numfit_max])
num_int_bacground = integrate.quad( num_int_bacground = integrate.quad(
background, background,
numfit_min, x[numfit_min],
numfit_max, x[numfit_max],
args=(result.params["slope"].value, result.params["intercept"].value), args=(result.params["slope"].value, result.params["intercept"].value),
) )
@ -131,6 +131,7 @@ def fitccl(
for pars in result.params: for pars in result.params:
d[str(pars)] = (result.params[str(pars)].value, result.params[str(pars)].vary) d[str(pars)] = (result.params[str(pars)].value, result.params[str(pars)].vary)
d["export_fit"] = False
d["int_area"] = num_int_area d["int_area"] = num_int_area
d["int_background"] = num_int_bacground d["int_background"] = num_int_bacground
d["full_report"] = result.fit_report() d["full_report"] = result.fit_report()