[submodule "test-pcb/waeda-sym"]
path = sch-pcb/waeda-sym
- url = ssh://dan@twain.unl.edu/space/git/waeda-sym.git
+ ;url = ssh://dan@twain.unl.edu/space/git/waeda-sym.git
+ url = git@git.whiteaudio.com:waeda-sym.git
[submodule "test-pcb/waeda-fp"]
path = sch-pcb/waeda-fp
- url = ssh://dan@twain.unl.edu/space/git/waeda-fp.git
+ ;url = ssh://dan@twain.unl.edu/space/git/waeda-fp.git
+ url = git@git.whiteaudio.com:waeda-fp.git
cal_outs = []
infiles = sorted(glob('*.npz'))
-#N_RUNS = min(len(infiles), 8)
-N_RUNS = 0
+N_RUNS = min(len(infiles), 8)
+#N_RUNS = 0
sample_runs = sorted(rand_sample(range(len(infiles)), N_RUNS))
-#N_PLOTS = 4
-N_PLOTS = 0
-#sample_channels = sorted(rand_sample(range(48), N_PLOTS))
-sample_channels = range(N_PLOTS)
-#N_COLS = 4
-N_COLS = 0
+N_PLOTS = 4
+#N_PLOTS = 0
+sample_channels = sorted(rand_sample(range(48), N_PLOTS))
+#sample_channels = range(N_PLOTS)
+N_COLS = 4
+#N_COLS = 0
nrun = 0
for infile in infiles:
#for infile in sorted(glob('*.npz')):
hist(data, bins=100, range=(-dm, dm), align='left')
xlim((-dm, dm))
#ylim((0, 50))
+ print i,k,data.std()
subplot(N_PLOTS, N_COLS, 1)
title('Offset DAC value')
interactive(False)
-if 1:
+if 0:
figure()
data = cv.flatten()
hist(data, bins=range(128), align='left')
(len(infiles),))
savefig('calout-hist.pdf')
-print sample_runs
-print [infiles[s] for s in sample_runs]
+
+if 1:
+ # stddev for each calibrated channel
+ means = []
+ stdds = []
+ c = 0
+ for k in range(co.shape[2]):
+ data = 1e3*(2.5/2/2**13)*co[:, 0, k]
+ m = data.mean()
+ s = data.std()
+ means.append(m)
+ stdds.append(s)
+ if abs(m) < 10.0 and abs(s) < 10.0:
+ c += 1
+ print '%02i: mean: %5.1f, std: %5.1f' % (k, m, s)
+
+ print 'good channels:', c
+#print sample_runs
+#print [infiles[s] for s in sample_runs]