Skip to content
Snippets Groups Projects

Resolve "issues at first running gaussian_example.py"

Merged Resolve "issues at first running gaussian_example.py"
205-add-shebangs into master
Merged Gregory Ashton requested to merge 205-add-shebangs into master
Compare
and
Show latest version
4 files
+ 0
− 358
Compare changes
  • Side-by-side
  • Inline
Files
4
examples/injection_examples/make_a_timeseries.py deleted 100644 → 0
+ 0
− 47
#!/usr/bin/env python
"""
"""
from __future__ import division, print_function
import numpy as np
import tupak
import gwpy
duration = 4.
sampling_frequency = 2048.
start_time = 0
# Specify the output directory and the name of the simulation.
outdir = 'outdir'
label = 'timeseries'
injection_parameters = dict(mass_1=36., mass_2=29., a_1=0.4, a_2=0.3,
tilt_1=0.5, tilt_2=1.0, phi_12=1.7, phi_jl=0.3,
luminosity_distance=2000., iota=0.4, psi=2.659,
phase=1.3, geocent_time=1126259642.413,
ra=1.375, dec=-1.2108)
waveform_arguments = dict(waveform_approximant='IMRPhenomPv2',
reference_frequency=50.)
waveform_generator = tupak.gw.WaveformGenerator(
duration=duration, sampling_frequency=sampling_frequency,
frequency_domain_source_model=tupak.gw.source.lal_binary_black_hole,
parameters=injection_parameters, waveform_arguments=waveform_arguments)
hf_signal = waveform_generator.frequency_domain_strain()
H1 = tupak.gw.detector.get_empty_interferometer('H1')
H1.set_strain_data_from_power_spectral_density(
sampling_frequency=sampling_frequency, duration=duration,
start_time=start_time)
H1.inject_signal(waveform_generator=waveform_generator,
parameters=injection_parameters)
time_domain_data = H1.strain_data.time_domain_strain
timeseries = gwpy.timeseries.TimeSeries(time_domain_data, t0=start_time,
sample_rate=sampling_frequency)
# This line will drop you to the ipython terminal. Useful for seeing what
# different objects are. I recommand setting up a hotkey to paste this in.
# Comment it out to stop it.
import IPython; IPython.embed()