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.