From b9080e560941b348bce044310ecf6ac8b63dfa0a Mon Sep 17 00:00:00 2001 From: gobbo_a Date: Fri, 5 May 2017 14:58:32 +0200 Subject: [PATCH] Script execution --- script/test/WireScanTest2.py | 20 +++++++++++++------- 1 file changed, 13 insertions(+), 7 deletions(-) diff --git a/script/test/WireScanTest2.py b/script/test/WireScanTest2.py index ae4aa56..019bdc8 100644 --- a/script/test/WireScanTest2.py +++ b/script/test/WireScanTest2.py @@ -1,3 +1,5 @@ +run("Diagnostics/sig_process_wrapper") + root = "2017/05/04/20170504_182411_WireScan.h5|/" wire = "x" blms = ["",""] @@ -6,10 +8,11 @@ cycles = 3 calc_values = [] -for i in range(len(blms)): +for i in range(len(blms)): path_back = root + "background/blm" + str(i+1) bg = get_attributes(path_back)["Mean"] - for cycle in range (cycles): + samples = [[], [], [], []] + for cycle in range (cycles): pos = load_data(root + wire+"_" + ("%04d" % (cycle+1)) + "/w_pos") path = root + wire+"_" + ("%04d" % (cycle+1)) + "/blm" + str(i+1) data = load_data(path) @@ -18,19 +21,22 @@ for i in range(len(blms)): 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 (rms_com is None) else rms_com) + #set_attribute(path, "RMS COM", 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]) + samples[0].append(rms_com);samples[1].append(rms);samples[2].append(com);samples[3].append(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)) - - - + stats = [] + for sample in samples: + stats.apppend((mean(sample), stdev(sample))) + print "-------------" + print i, stats + \ No newline at end of file