From 38811dd315fa4d2c99a2171ec87477dc0e9fec05 Mon Sep 17 00:00:00 2001 From: gobbo_a Date: Fri, 5 May 2017 12:57:55 +0200 Subject: [PATCH] Script execution --- script/test/WireScanTest2.py | 14 +++++++++----- 1 file changed, 9 insertions(+), 5 deletions(-) diff --git a/script/test/WireScanTest2.py b/script/test/WireScanTest2.py index ecf52c7..ae4aa56 100644 --- a/script/test/WireScanTest2.py +++ b/script/test/WireScanTest2.py @@ -3,6 +3,8 @@ wire = "x" blms = ["",""] cycles = 3 +calc_values = [] + for i in range(len(blms)): path_back = root + "background/blm" + str(i+1) @@ -14,19 +16,21 @@ for i in range(len(blms)): sp = blm_remove_spikes(data) sig = sp if bg is None else [v-bg for v in sp] - [com, rms] = profile_rms_stats(pos, sig,noise_std=0, n_sigma=3.5) + rms_com, rms = profile_rms_stats(pos, sig,noise_std=0, n_sigma=3.5) #set_attribute(path, "RMS", float("nan") if (rms is None) else rms) - #set_attribute(path, "RMS centroid", float("nan") if (com is None) else com) - print [com, rms] - [off, amp, com, sigma] = profile_gauss_stats(pos, sig, off=None, amp=None, com=None, sigma=None) + #set_attribute(path, "RMS centroid", float("nan") if (rms_com is None) else rms_com) + print [rms_com, rms] + off, amp, com, sigma = profile_gauss_stats(pos, sig, off=None, amp=None, com=None, sigma=None) #set_attribute(path, "Gauss COM", float("nan") if (com is None) else com) #set_attribute(path, "Gauss Sigma", float("nan") if (sigma is None) else sigma) print [off, amp, com, sigma] + calc_values.append([rms_com, rms, com, sigma]) + from mathutils import Gaussian g = Gaussian(amp, com, sigma) gauss = [g.value(v)+off for v in pos] - plot([data, sp, sig, gauss], ["data", "sp", "signal", "gauss", ], xdata = pos, title="Fit blm" + str(i+1) + " - " + str(cycle+1)) +