--- /dev/null
+#!/usr/bin/env python
+
+from pylab import *
+
+import os
+from glob import glob
+from random import sample as rand_sample
+
+
+
+fig = figure(figsize=(5.0, 5.0))
+subplots_adjust(top=0.98, bottom=0.08, left=0.08, right=0.98)
+
+
+
+d = [False]
+for i in range(1, 5):
+ data = np.load('chip%02i/chip%02i-caldac.npz' % (i, i))['data']
+ N = len(data)
+
+ ax = subplot(4, 1, i)
+ n, bins, patches = hist(data, bins=range(128), align='left')
+
+ xlim((0, 128))
+ xticks(range(0, 129, 16))
+ ax.set_xticklabels('') # no labels by default
+
+ #print n.max(), N, 1.0*n.max()/N
+ ylim((0, 0.08*N))
+ ax.set_yticklabels('')
+ #ylabel('Bin count')
+ #title('48x2 channels x %i runs' % (len(infiles),))
+
+ text(0.04, 0.7,
+ 'Chip \#%02i, %i runs, N=%i' % (i, N/96, N),
+ transform=ax.transAxes)
+
+ax = subplot(4, 1, 4)
+xlabel('Offset tuning DAC code')
+ax.set_xticklabels(ax.get_xticks())
+
+fig.text(0.02, 0.52,
+ 'Relative bin counts',
+ rotation=90,
+ va='center')
+
+savefig('caldac-hist.pdf')
+
figure(figsize=(5.0, 3.5))
subplots_adjust(top=0.93, bottom=0.12, left=0.10, right=0.98)
data = cv.flatten()
- hist(data, bins=range(128), align='left')
- ylim((0, 7000))
- xlabel('OTA tuning DAC code')
+ n, bins, patches = hist(data, bins=range(128), align='left')
+ #ylim((0, 7000))
+ ylabel('Bin count')
+ xlim((0, 128))
+ xticks(range(0, 129, 16))
+ xlabel('Offset tuning DAC code')
title('48x2 channels x %i runs' % (len(infiles),))
savefig(os.path.basename(os.getcwd()) + '-caldac-hist.pdf')
+ np.savez(os.path.basename(os.getcwd()) + '-caldac.npz',
+ data=data)
figure(figsize=(5.0, 3.5))
data = co.flatten() * 1e-3
hist(data, bins=arange(-1.500, 1.500, 0.025), align='left')
xlabel('Residual $V_{oos}$ (V)')
- ylim((0, 7000))
+ #ylim((0, 7000))
title('48x2 channels x %i runs' % (len(infiles),))
savefig(os.path.basename(os.getcwd()) + '-calout-hist.pdf')
fig = figure(figsize=(6.0, 3.5))
subplots_adjust(top=0.97, bottom=0.15, left=0.15, right=0.98)
+zphase_outs = zeros((2, 7))
+
+
for name in glob('arb0-*5.0-*.npz'):
# ['f0', 'fs', 'vout', 'phase', 'Nperiods', 'fh', 'fsin', 'duration']
data = np.load(name)
else:
style = '+g'
- plot(wphase, vout, style)
+ zphase_outs[col, row] = vout[0] - vout.mean()
+
+ #plot(wphase, vout, style)
+ plot(wphase, vout-vout.mean(), style)
ylim([-0.8, 0.8])
xlim([0, 360])
horizontalalignment='center',
rotation='horizontal')
-savefig('harmonic-phase.pdf')
+#savefig('harmonic-phase.pdf')
+
+
+
+
+x = 5.0 * arange(1, 8)
+
+figure()
+subplot(211)
+plot(x, zphase_outs[0,:], 'o')
+plot([0, 70], [0, 0], '-k')
+ylim((-0.1, 0.8))
+xlim((4, 31))
+
+subplot(212)
+plot(x, zphase_outs[1,:], 'o')
+plot([0, 70], [0, 0], '-k')
+ylim((-0.1, 0.8))
+xlim((4, 31))
+savefig('tmp.pdf')