Commit 8c25d773 authored by Andrew Kim's avatar Andrew Kim Committed by Gregory Ashton

plot_skymap example update

parent 46af0329
......@@ -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 =
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 =
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,
# make some plots of the outputs
result.plot_corner(), nside=2**6)
# will require installation of ligo.skymap (pip install ligo.skymap)
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