Script execution
This commit is contained in:
@@ -17,15 +17,15 @@ for i in range(len(blms)):
|
||||
[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 i, j, [com, rms]
|
||||
print [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 i, j, [off, amp, com, sigma]
|
||||
print [off, amp, com, sigma]
|
||||
|
||||
from mathutils import Gaussian
|
||||
g = Gaussian(amp, com, sigma)
|
||||
gauss = [g.value(i)+off for i in pos]
|
||||
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))
|
||||
|
||||
|
||||
Reference in New Issue
Block a user