最新消息:雨落星辰是一个专注网站SEO优化、网站SEO诊断、搜索引擎研究、网络营销推广、网站策划运营及站长类的自媒体原创博客

How can I get the full detector response in a parameter space point (Bilby)? - Stack Overflow

programmeradmin4浏览0评论

I defined my waveform in the following way

waveform_arguments = dict(
    waveform_approximant='IMRPhenomD', reference_frequency=10, catch_waveform_errors=True, minimum_frequency=10.,
maximum_frequency = 10000  )

waveform_generator = bilby.gw.WaveformGenerator(
    frequency_domain_source_model=bilby.gw.source.lal_binary_black_hole,
    waveform_arguments=waveform_arguments,
    parameter_conversion=convert_to_lal_binary_black_hole_parameters)

and I can't calculate hplus and hcross in this way, with some random point I defined:

dictionary = waveform_generator.frequency_domain_strain(random_points['points'][0])

and I want the full detector waveform (what would be subtracted from the observed strain). Given that I defined my detector as

det = bilby.gw.detector.Interferometer(
    power_spectral_density= bilby.gw.detector.PowerSpectralDensity(
        frequency_array=ASD[:,0], asd_array=ASD[:, 1]
    ),
    name="det",
    length=10,
    minimum_frequency=10,
    maximum_frequency=10000,
    latitude=40.5,
    longitude=9.4,
    elevation=0.0,
    xarm_azimuth=0,
    yarm_azimuth=60,
)

I am getting warnings about duration and sampling_frequency, which I don have, since there is no real signal. I wanted to do the obvious, just use the frequency array I provided in the detector. Anyone knows how to do that?

I tried to provide the array in the waveform generator as

waveform_arguments = dict(
    waveform_approximant='IMRPhenomD', reference_frequency=10, catch_waveform_errors=True, minimum_frequency=10.,
maximum_frequency = 10000, frequency_array=ASD[:,0]  )

but I still get the same problem.

发布评论

评论列表(0)

  1. 暂无评论