Moving from GnuplotBiDir fitting to scipy.optimize.leastsq fitting. Major improvement in fit quality on our problem case: $ cc_vib_analyze.py -m 300 -M 7000 -f50000 -p -G -d 20080919/20080919121041_zp_time_50000Hz