Single detector, O1 data (restricted prior)
===========================================

Here we compare ``lalapps_pulsar_parameter_estimation_nested`` with ``cwinpy`` in the case of real
gravitational-wave data from a single detector (H1 in this case). This data (which can be downloaded
:download:`here <scripts/data/O1_Crab_H1.txt.gz>`) comes from the `O1 run
<https://www.gw-openscience.org/O1/>`_ and has been heterodyned using the parameters for the Crab
pulsar given in this file: :download:`J0534+2200.par <scripts/data/J0534+2200.par>`. The parameters
being estimated are :math:`h_0`, :math:`\phi_0`, :math:`\psi` and :math:`\iota`. In this case the
:math:`h_0` and :math:`\phi_0` parameter have uniform priors, while the :math:`\psi` and
:math:`\iota` parameters have "restricted" priors. The :math:`\psi` prior is a unimodal Gaussian
distribution, while the :math:`\iota` prior is a bimodal Gaussian distribution, with values dictated
by fitting the orientation of the Crab's pulsar wind nebula. 

The script for this comparison, using the `dynesty <https://dynesty.readthedocs.io/en/latest/>`_
nested sampling algorithm, is shown at the bottom of the page_. It produces the following comparison
data:

.. image:: scripts/outputs/single_detector_O1_data_restricted_prior_corner.png

.. include:: scripts/outputs/single_detector_O1_data_restricted_prior_compare.txt

.. _page:

.. literalinclude:: scripts/single_detector_O1_data_restricted_prior.py
   :language: python
