Use abs(C) in fit_psd() (to match breit_wigner()).

diff --git a/calibcant/vib_analyze.py b/calibcant/vib_analyze.py
index d226d358504a819847a84ddabd6661649f10167c..59c83170daddabf5a30ec0e342b0a16fa944047c 100644 (file)
--- a/calibcant/vib_analyze.py
+++ b/calibcant/vib_analyze.py
@@ -436,7 +436,8 @@ def fit_psd(freq_axis, psd_data, min_frequency=500, max_frequency=25000,
         A,B,C = p
         D = 0
     A=abs(A) # A and B only show up as squares in f(x)
-    B=abs(B) # so ensure we get positive values
+    B=abs(B) # so ensure we get positive values.
+    C=abs(C) # Only abs(C) is used in breit_wigner().
     return (A, B, C, D)
 
 def breit_wigner_area(A, B, C):