reference_psd.py: add functions to apply a moving median to the spectrum and a polynomial fit to the log

far.py: clip the SNR of events to be less than 100 prior to putting them into the likelihood ratio bins

gstlal_plot_psd: only ever compute horizon up to 90% nyquist to avoid filter rolloff in downsampled data