pshell_config/sf-op
0
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Script execution

Browse Source
This commit is contained in:
gobbo_a
2017-05-05 14:58:32 +02:00
parent 38811dd315
commit b9080e5609
1 changed files with 13 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
script/test/WireScanTest2.py
View File

@@ -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