Commit 6e7eb42a authored by moritz's avatar moritz
Browse files

Moritz Huebner: Made it work with new changes

parent ff03ed17
......@@ -58,7 +58,8 @@ IFOs = [H1, L1]
source = peyote.source.BinaryBlackHole('BBH', sampling_frequency, time_duration, spin_1=[0, 0, 0], spin_2=[0, 0, 0],
luminosity_distance=410., iota=2.97305, phase=1.145,
waveform_approximant='IMRPhenomPv2', reference_frequency=50., ra=1.375,
dec=-1.2108, geocent_time=1126259642.413, psi=2.659, mass_1=36, mass_2=28)
dec=-1.2108, geocent_time=1126259642.413, psi=2.659, mass_1=32, mass_2=32)
# ignore the fact that I hardcoded in some masses
likelihood = peyote.likelihood.Likelihood(IFOs, source)
prior = source.copy()
......
import peyote
import os
import numpy as np
np.random.seed(10)
time_duration = 1.
sampling_frequency = 4096.
simulation_parameters = dict(
mass_1=36.,
mass_2=29.,
spin_1=[0, 0, 0],
spin_2=[0, 0, 0],
luminosity_distance=410.,
iota=0.,
phase=0.,
waveform_approximant='IMRPhenomPv2',
reference_frequency=50.,
ra=0,
dec=1,
geocent_time=0,
psi=1
)
source = peyote.source.BinaryBlackHole(
'BBH', sampling_frequency, time_duration)
hf_signal = source.frequency_domain_strain(simulation_parameters)
source = peyote.source.BinaryBlackHole('BBH', sampling_frequency, time_duration, mass_1=36., mass_2=29.,
spin_1=[0, 0, 0], spin_2=[0, 0, 0], luminosity_distance=410., iota=0., phase=0.,
waveform_approximant='IMRPhenomPv2', reference_frequency=50., ra=0, dec=1,
geocent_time=0, psi=1)
hf_signal = source.frequency_domain_strain()
IFO = peyote.detector.H1
IFO.set_data(from_power_spectral_density=True, sampling_frequency=sampling_frequency, duration=time_duration)
IFO.inject_signal(source, simulation_parameters)
IFO.inject_signal(source)
hf_signal_and_noise = IFO.data
frequencies = peyote.utils.create_fequency_series(sampling_frequency=sampling_frequency, duration=time_duration)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment