def draw_times(ax, wf, wt):
#wf -= wf[0:1000].mean()
edges = find_edges(wf, BASE, THR, CFR, DEADTIME, LEADINGEDGE)
# pairs of (amplitude,sampleNumber)
print 'MCP edges'
print edges
for (amp, ind) in edges:
x0 = wt[int(ind)]
xarr = (x0, x0)
yarr = (amp, -amp)
gg.drawLine(ax, xarr, yarr, s=10, linewidth=1, color='k')
def draw_times(ax, wf, wt):
#wf -= wf[0:1000].mean()
t0_sec = time()
edges = find_edges(wf, BASE, THR, CFR, DEADTIME, LEADINGEDGE)
print ' consumed time = %10.6f(sec)' % (time() - t0_sec)
# pairs of (amplitude,sampleNumber)
print ' nhits:', edges.shape[0]
#print ' edges:', edges
for (amp, ind) in edges:
x0 = wt[int(ind)]
xarr = (x0, x0)
yarr = (amp, -amp)
gg.drawLine(ax, xarr, yarr, s=10, linewidth=1, color='k')
def draw_times(ax, wf, wt):
#extrema = 0.1* wf_extremas(ax, wf, wt, rank=10)
#gg.drawLine(ax, wt, extrema, s=10, linewidth=1, color='k')
#return
#wf -= wf[0:1000].mean()
#edges = find_edges(wf, BASE, THR, CFR, DEADTIME, LEADINGEDGE)
edges = wf_extremas(ax, wf, wt, rank=10)
# pairs of (amplitude,sampleNumber)
print ' nhits:', edges.shape[0],
print ' edges:', edges
for (amp, ind) in edges:
x0 = wt[int(ind)]
xarr = (x0, x0)
yarr = (amp, -amp)
gg.drawLine(ax, xarr, yarr, s=10, linewidth=1, color='k')
amp, ind = edges[0]
#if amp<THRU : continue
if (ind < NBINSOFF) or (BNUM - ind < NBINSWL): continue
sp.w[i, :] += wfsmo[ind - NBINSOFF:ind + NBINSWL - NBINSOFF]
#ax[i].plot(wtsel, wfsel, gfmt[i], linewidth=lw)
ax[i].plot(wtsel, wfsmo, gfmt[0], linewidth=lw)
#ax[i].plot(wtsel, grad2, gfmt[0], linewidth=lw)
gg.drawLine(ax[i],
ax[i].get_xlim(), (THR, THR),
s=10,
linewidth=1,
color='k')
#wf2sel = wf2[ch[i],BBEG:BEND]
#wt2sel = wt2[ch[i],BBEG:BEND]
#ax[i].plot(wt2sel, wf2sel, gfmt[i], linewidth=lw)
#edges = draw_times(ax[i], wf2sel, wt2sel)
gr.draw_fig(fig)
gr.show(mode='non-hold')
for i in range(NAXES):
ax[i].clear()
ax[i].plot(wttot[:NBINSWL], sp.w[i, :], gfmt[0], linewidth=lw)
