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 12:57:55 +02:00
parent 29c8f9ec99
commit 38811dd315
1 changed files with 9 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files
script/test/WireScanTest2.py
+9 -5
View File
@@ -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))