From f5c405bee8047317dc1b059cd95b7eafbf3db279 Mon Sep 17 00:00:00 2001 From: JakHolzer <53743814+JakHolzer@users.noreply.github.com> Date: Fri, 18 Dec 2020 15:40:37 +0100 Subject: [PATCH] Numerically integrate whole area Added two lines to integrate full area on the request of Oksana and Romain, might be changed back after some more testing, ergo I did not delete the previous code and only change the final result regardless of what the code produces. I find this least invasive, but a bit messy. --- pyzebra/fit2.py | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/pyzebra/fit2.py b/pyzebra/fit2.py index 7ac0c8e..585d84b 100644 --- a/pyzebra/fit2.py +++ b/pyzebra/fit2.py @@ -190,7 +190,9 @@ def fitccl( print("Maximal integration value lower than minimal") else: pass - + + numfit_min = 0 + numfit_max = len(x) count_errors = create_uncertanities(y, y_err) # create error vector for numerical integration propagation num_int_area = simps(count_errors[numfit_min:numfit_max], x[numfit_min:numfit_max])