print "plotting"
common._import_pylab()
common._pylab.figure(config.BASE_FIGNUM+2)
- common._pylab.hold(False)
# plot time series
common._pylab.subplot(311)
+ common._pylab.hold(False)
common._pylab.plot(deflection_bits, 'r.')
common._pylab.title("free oscillation")
# plot histogram distribution and gaussian fit
common._pylab.subplot(312)
+ common._pylab.hold(False)
n, bins, patches = \
common._pylab.hist(deflection_bits, bins=30,
normed=1, align='center')
# plot FFTed data
axes = common._pylab.subplot(313)
+ common._pylab.hold(False)
common._pylab.semilogy(freq_axis, power, 'r.-')
fitdata = C / ((A**2-freq_axis**2)**2 + (B*freq_axis)**2)
xmin,xmax = axes.get_xbound()