Single detector, O1 data
========================

Here we compare ``lalpulsar_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://gwosc.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:`\cos{\iota}`, all with uniform
priors.

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_corner.png

.. include:: scripts/outputs/single_detector_O1_data_compare.txt

.. _page:

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