add zeroNoise PE example for LISA

gwastro · Sep 29, 2024 · 3ef908f · 3ef908f
1 parent d4d4fa0
commit 3ef908f
Showing 1 changed file with 16 additions and 3 deletions.
diff --git a/examples/inference/lisa_smbhb_inj/lisa_smbhb_relbin.ini b/examples/inference/lisa_smbhb_inj/lisa_smbhb_relbin.ini
@@ -5,6 +5,10 @@ analysis-start-time = -4800021
 analysis-end-time = 26735979
 pad-data = 0
 sample-rate = 0.2
+channel-name = LISA_A:LISA_A LISA_E:LISA_E LISA_T:LISA_T
+injection-file = injection_smbhb.hdf
+
+; with noise simulation
 fake-strain = LISA_A:analytical_psd_lisa_tdi_AE LISA_E:analytical_psd_lisa_tdi_AE LISA_T:analytical_psd_lisa_tdi_T
 ; fake-strain-extra-args = LISA_A:len_arm:2.5e9 LISA_A:acc_noise_level:2.4e-15 LISA_A:oms_noise_level:7.9e-12 LISA_A:tdi:1.5 LISA_E:len_arm:2.5e9 LISA_E:acc_noise_level:2.4e-15 LISA_E:oms_noise_level:7.9e-12 LISA_E:tdi:1.5 LISA_T:len_arm:2.5e9 LISA_T:acc_noise_level:2.4e-15 LISA_T:oms_noise_level:7.9e-12 LISA_T:tdi:1.5
 fake-strain-extra-args = len_arm:2.5e9 acc_noise_level:2.4e-15 oms_noise_level:7.9e-12 tdi:1.5
@@ -19,8 +23,16 @@ psd-segment-length = 267840
 psd-segment-stride = 133920
 psd-start-time = -4800021
 psd-end-time = 26735979
-channel-name = LISA_A:LISA_A LISA_E:LISA_E LISA_T:LISA_T
-injection-file = injection_smbhb.hdf
+
+; ; zero noise simulation, but still use PSD in likelihood
+; ; when using zero noise, please comment out above "with noise simulation" block and uncomment out this block
+; fake-strain = zeroNoise
+; fake-strain-flow = 0.0001
+; fake-strain-sample-rate = 0.2
+; fake-strain-filter-duration = 31536000
+; psd-model = LISA_A:analytical_psd_lisa_tdi_AE LISA_E:analytical_psd_lisa_tdi_AE LISA_T:analytical_psd_lisa_tdi_T
+; ; psd-extra-args = LISA_A:len_arm:2.5e9 LISA_A:acc_noise_level:2.4e-15 LISA_A:oms_noise_level:7.9e-12 LISA_A:tdi:1.5 LISA_E:len_arm:2.5e9 LISA_E:acc_noise_level:2.4e-15 LISA_E:oms_noise_level:7.9e-12 LISA_E:tdi:1.5 LISA_T:len_arm:2.5e9 LISA_T:acc_noise_level:2.4e-15 LISA_T:oms_noise_level:7.9e-12 LISA_T:tdi:1.5
+; psd-extra-args = len_arm:2.5e9 acc_noise_level:2.4e-15 oms_noise_level:7.9e-12 tdi:1.5
 
 [model]
 name = relative
@@ -40,7 +52,8 @@ tc =
 
 [static_params]
 ; Change it to "ref_frame = SSB", if you use SSB frame in injection file.
-ref_frame = LISA
+; ref_frame = LISA
+ref_frame = SSB
 approximant = BBHX_PhenomD
 ; You can use "1.5" or "2.0" for TDI.
 ; Please use the same TDI version for PSD and injection file.