def plot_counts_vs_time(windowdata,
savepath,
window=1,
bins=700,
save_prefix='',
save_postfix=''):
gated = hht3.filter_counts_on_marker(windowdata, 2, [window - 2])
print 'length of window data = ', len(windowdata)
fig = plt.figure()
t0, t1 = hht3.get_click_times(windowdata)
h0y, h0x = np.histogram(t0, bins=bins, range=(0, bins - 1))
h1y, h1x = np.histogram(t1, bins=bins, range=(0, bins - 1))
g_t0, g_t1 = hht3.get_click_times(gated)
g_h0y, g_h0x = np.histogram(g_t0, bins=bins, range=(0, bins - 1))
g_h1y, g_h1x = np.histogram(g_t1, bins=bins, range=(0, bins - 1))
ax = plt.subplot(121)
ax.semilogy(g_h0x[:-1], g_h0y, color='g', label='ch0 gated')
ax.semilogy(h0x[:-1], h0y, label='ch0', color='b')
ax.set_xlabel('Bins')
ax.legend()
ax = plt.subplot(122)
ax.semilogy(h1x[:-1], h1y, label='ch1')
ax.semilogy(g_h1x[:-1], g_h1y, label='ch1 gated')
ax.set_xlabel('Bins')
ax.legend()
plt.suptitle('Clicks after optical pulse, window %d' % window)
fig.savefig(os.path.join(savepath,'%shistogram_window%s%s.pdf' % \
(save_prefix,window,save_postfix)))
np.savez(os.path.join(savepath, 'histogram_window%s.npz' % window),
window_number=window,
window_data=windowdata,
ch0=h0y,
ch0_xaxis=h0x[:-1],
ch0_gated=g_h0y,
ch0_gated_xaxis=g_h0x[:-1],
ch1=h1y,
ch1_xaxis=h1x[:-1],
ch1_gated=g_h1y,
ch1_gated_xaxis=g_h1x[:-1])
return True
def plot_counts_vs_time(windowdata, savepath, window=1, bins=700,
save_prefix='', save_postfix=''):
gated = hht3.filter_counts_on_marker(windowdata, 2, [window-2])
print 'length of window data = ',len(windowdata)
fig = plt.figure()
t0, t1 = hht3.get_click_times(windowdata)
h0y, h0x = np.histogram(t0, bins=bins, range=(0,bins-1))
h1y, h1x = np.histogram(t1, bins=bins, range=(0,bins-1))
g_t0, g_t1 = hht3.get_click_times(gated)
g_h0y, g_h0x = np.histogram(g_t0, bins=bins, range=(0,bins-1))
g_h1y, g_h1x = np.histogram(g_t1, bins=bins, range=(0,bins-1))
ax = plt.subplot(121)
ax.semilogy(g_h0x[:-1], g_h0y, color = 'g', label='ch0 gated')
ax.semilogy(h0x[:-1], h0y, label='ch0', color = 'b')
ax.set_xlabel('Bins')
ax.legend()
ax = plt.subplot(122)
ax.semilogy(h1x[:-1], h1y, label='ch1')
ax.semilogy(g_h1x[:-1], g_h1y, label='ch1 gated')
ax.set_xlabel('Bins')
ax.legend()
plt.suptitle('Clicks after optical pulse, window %d' % window)
fig.savefig(os.path.join(savepath,'%shistogram_window%s%s.pdf' % \
(save_prefix,window,save_postfix)))
np.savez(os.path.join(savepath, 'histogram_window%s.npz'%window),
window_number = window,
window_data = windowdata,
ch0 = h0y,
ch0_xaxis = h0x[:-1],
ch0_gated = g_h0y,
ch0_gated_xaxis = g_h0x[:-1],
ch1 = h1y,
ch1_xaxis = h1x[:-1],
ch1_gated = g_h1y,
ch1_gated_xaxis = g_h1x[:-1])
return True
def plot_counts_vs_time(data, rng=500.):
fig = plt.figure()
t0, t1 = hht3.get_click_times(data)
h0y, h0x = np.histogram(t0, bins=rng+1, range=(0.,rng))
h1y, h1x = np.histogram(t1, bins=rng+1, range=(0.,rng))
ax = plt.subplot(121)
ax.semilogy(h0x[:-1], h0y, 'o-', c='k', mfc='None', mec='k', label='ch0')
ax.set_xlabel('Bins')
ax.legend()
ax = plt.subplot(122)
ax.semilogy(h1x[:-1], h1y, 'o-', c='k', mfc='None', mec='k', label='ch1')
ax.set_xlabel('Bins')
ax.legend()
def plot_counts_vs_time(data, rng=500.0):
fig = plt.figure()
t0, t1 = hht3.get_click_times(data)
h0y, h0x = np.histogram(t0, bins=rng + 1, range=(0.0, rng))
h1y, h1x = np.histogram(t1, bins=rng + 1, range=(0.0, rng))
ax = plt.subplot(121)
ax.semilogy(h0x[:-1], h0y, "o-", c="k", mfc="None", mec="k", label="ch0")
ax.set_xlabel("Bins")
ax.legend()
ax = plt.subplot(122)
ax.semilogy(h1x[:-1], h1y, "o-", c="k", mfc="None", mec="k", label="ch1")
ax.set_xlabel("Bins")
ax.legend()
def plot_counts_vs_time(data, rng=500.):
fig = plt.figure()
t0, t1 = hht3.get_click_times(data)
h0y, h0x = np.histogram(t0, bins=rng+1, range=(0.,rng))
h1y, h1x = np.histogram(t1, bins=rng+1, range=(0.,rng))
ax = plt.subplot(121)
ax.semilogy(h0x[:-1], h0y, 'o-', c='k', mfc='None', mec='k', label='ch0')
ax.set_xlabel('Bins')
ax.legend()
ax = plt.subplot(122)
ax.semilogy(h1x[:-1], h1y, 'o-', c='k', mfc='None', mec='k', label='ch1')
ax.set_xlabel('Bins')
ax.legend()
