从文件中设置的功率谱密度

Question

如何从文件执行一组的功率谱密度（PSD），并且有可能使用一个不同的PSD用于产生数据和似然的评价？

问题费雯丽雷蒙德提前通过电子邮件。

Answer 1

从设置文件中的PSD

从文件中设置的PSD，第一初始化干涉仪的列表，在这里我们只使用恒福：

>>> ifos = bilby.gw.detector.InterferometerList(['H1'])

列表中的每个元素都被初始化使用先进LIGO噪声曲线默认的PSD，以检查此

>>> ifos[0].power_spectral_density
PowerSpectralDensity(psd_file='/home/user1/miniconda3/lib/python3.6/site-packages/bilby-0.3.5-py3.6.egg/bilby/gw/noise_curves/aLIGO_ZERO_DET_high_P_psd.txt', asd_file='None')

请注意，尚未产生任何数据。要覆盖PSD，只需创建一个新PowerSpectralDensity对象，并为其分配（如果你有多个检测器，你就需要为列表中的每个元素做到这一点）

ifos[0].power_spectral_density = bilby.gw.detector.PowerSpectralDensity(psd_file=PATH_TO_FILE)

巢，生成从PSD的应变数据的一个实例：

ifos.set_strain_data_from_power_spectral_densities(
    sampling_frequency=4096, duration=4,
    start_time=-3)

您可以检查做什么数据的模样

ifos[0].plot_data()

请注意，您还可以使用注射的方法ifos.inject_signal信号。

使用可能性的评估不同的PSD

在ifos列表中的每个IFO既包含数据和PSD（或等效ASD）。为推论，我们通过该列表到bilby.gw.GravitationalWaveLikelihood对象作为第一个参数和用于列表的每个元素的PSD在计算似然性时使用。

所以，如果你想使用似然估计不同的PSD。第一生成数据（如上）。然后，您指定要用于采样到该对象进入的可能性，而不是ifos和传球的每个元素的PSD。这不会覆盖数据（只要你不叫，当然set_strain_data_from_power_spectral_densities）。