1 files changed, 32 insertions, 0 deletions
diff --git a/calib_plots.py b/calib_plots.py
new file mode 100644
index 0000000..a905d91
--- /dev/null
+++ b/calib_plots.py
@@ -0,0 +1,32 @@
+
+def find_peaks(x,num=8):
+    d = x;
+    #figure();
+    #plot(d);
+
+    peaks = r_[0.:num];
+    i = 0;
+    mmm = max(d)
+    for p in range(0,int(num)):
+        while (d[i] < mmm*.2):
+            i += 1;
+        peaks[p] = 0.;
+        for n in range(-1,10):
+            peaks[p] += r_[i+n]*d[i+n]
+        peaks[p] = peaks[p]/sum(d[i-1:i+9])
+        i += 13
+    return peaks;
+
+
+def fit_peaks(p):
+    l=lambda a,x:a[0]*x+a[1]
+    print p; 
+    plotfit(fittodata_linear(r_[0.:size(p)],p,[25.,1.],err=.5*ones(size(p)),ret_all=True));
+
+def plot_mcacalib():
+    figure();
+    plot(fmca_calib[0:450]);
+    txt = text(0.22,.75,"delta-T/",
+            horizontalalignment='center',
+            verticalalignment='center',
+            transform = ax.transAxes)