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)