Commit 74487037 authored by Gregory Ashton's avatar Gregory Ashton

Merge branch 'plot_skymap_update' into 'master'

plot_skymap example update

See merge request !527
parents 46af0329 8c25d773
Pipeline #67244 passed with stages
in 6 minutes and 4 seconds
......@@ -29,7 +29,7 @@ def vector_tensor_sine_gaussian(frequency_array, hrss, Q, frequency, epsilon):
Relative size of the vector modes compared to the tensor modes.
"""
waveform_polarizations = bilby.gw.source.sinegaussian(
frequency_array, hrss, Q, frequency, 0, 0, 0, 0)
frequency_array, hrss, Q, frequency)
waveform_polarizations['x'] = epsilon * waveform_polarizations['plus']
waveform_polarizations['y'] = epsilon * waveform_polarizations['cross']
......
......@@ -43,9 +43,10 @@ likelihood = bilby.gw.GravitationalWaveTransient(
result = bilby.run_sampler(
likelihood=likelihood, priors=priors, sampler='dynesty', npoints=1000,
injection_parameters=injection_parameters, outdir=outdir, label=label)
injection_parameters=injection_parameters, outdir=outdir, label=label, result_class=bilby.gw.result.CBCResult)
# make some plots of the outputs
result.plot_corner()
bilby.gw.utils.plot_skymap(result, nside=2**6)
# will require installation of ligo.skymap (pip install ligo.skymap)
result.plot_skymap(maxpts=5000)
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