dummy changes

Library/config.py

@@ -23,8 +23,8 @@ def identity(param):
     return param
 param_functions = [np.cos, identity, identity, np.log, np.log, np.log, np.sin, identity, np.log]
 
-param_indices = [6]
-# param_indices = [0,1,2,3,4,5,6,7,8]
+# param_indices = [6]
+param_indices = [0,1,2,3,4,5,6,7,8]
 
 param_scales = [0,0,0,0,0,0,3.14e-03,0,0]
 

Library/set_injection_parameters.py

@@ -88,6 +88,7 @@ f_high = Compute_fmax(R_isco, m1, m2)
 # import IPython; IPython.embed();
 duration = tc_3p5PN + 2
 sampling_frequency = 4096
+# sampling_frequency = 4397
 # sampling_frequency = 2 * fmax
 start_time = injection_parameters['geocent_time'] - tc_3p5PN
 
@@ -96,7 +97,8 @@ meta_params = dict(
     tc_3p5PN = tc_3p5PN,
     f_high = f_high,
     duration = duration,
-    sampling_frequency = sampling_frequency
+    sampling_frequency = sampling_frequency,
+    start_time = start_time
 )
 
 injection_parameters['meta'] = meta_params

Tests/compare_waveform.py

@@ -117,7 +117,7 @@ print("overlap(h_wave_1_bilby,h_wave_1_bilby) = {}".format(ov_hbilby_hbilby))
 h_wave_bilby = [h_wave_1_bilby, h_wave_2_bilby, h_wave_3_bilby]
 h_wave_hmc = [h_wave_1_hmc, h_wave_2_hmc, h_wave_3_hmc]
 
-ifo_number = 2
+ifo_number = 1
 ifo_index = ifo_number - 1
 
 # freqs_for_plot = freqs[i_low:i_high+1]
@@ -130,8 +130,8 @@ f_min = 0
 f_max = freqs_for_plot[-1]
 f_min = max(f_min, sigpar[13])
 # f_max = min(f_max, sigpar[14])
-# f_min = 800
-# f_max = 804
+f_min = 30
+f_max = 50
 # f_min = 1500
 # f_max = 1504
 indices = np.where((f_min<freqs_for_plot) & (freqs_for_plot<f_max))[0]
@@ -153,6 +153,7 @@ print("")
 plot_fourier_domain = True
 plot_time_domain = False
 one_figure_per_plot = False
+plot_power_spectra = True
 
 if plot_fourier_domain or plot_time_domain: