%pylab inline
import ae
Populating the interactive namespace from numpy and matplotlib
f = ae.open("files/M5.wfs")
print f.size
642798592
f.plot()
xlabel("Time [s]")
ylabel("Amplitude [V]")
grid()
100% 2.03s
events = f.get_events(0.02)
print events.size
100% 3.86s
6161
ae.hist(events.maxima)
xlabel("Event maximum [V]")
ylabel("Number of events")
xlim(0.015,3)
ylim(0.8,2000)
grid(True)
figure(figsize=(10,5))
for i,l in enumerate("durations energies maxima rise_times counts".split()):
subplot(2,3,i+1)
title(l.replace("_"," "))
ae.hist(getattr(events,l))
tight_layout()
figure(figsize=(10,7))
order = argsort(events.maxima)
for ax in events[order[::-1]].plot( ax=(subplot(3,3,x+1) for x in xrange(9)) ):
pass
#ax.set_xticks([])
#ax.set_yticks([])
ax.set_ylim(-3,3)
tight_layout()