fix f_ref for SEOBNRv5

bhooshan-gadre · Jul 25, 2023 · 236763f · 236763f
1 parent 6b15e0f
commit 236763f
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/pycbc/waveform/gwsignal_utils.py b/pycbc/waveform/gwsignal_utils.py
@@ -20,10 +20,17 @@
         'f_final': 'f_max',
         }
 
+def fix_SEOBNRv5_fref(par):
+    '''SEOBNRv5 expects f_ref != 0. Here we set it to f_lower if not given
+    '''
+    if 'SEOBNRv5' in par['approximant'] and not par['f_ref']:
+        par['f_ref'] = par['f_lower']
+
 
 def to_gwsignal_dict(par):
     '''convert param dict to gws dict
     '''
+    fix_SEOBNRv5_fref(par) # this is in-place
     params = {pycbc_to_gws.get(k, k): v for k, v in par.items() if (v is not
         None and pycbc_to_gws.get(k, k) in gws_units)}