nanograv / 15yr_stochastic_analysis

Companion tutorials for the NANOGrav 15-year GWB analysis
MIT License
12 stars 6 forks source link

NANOGrav_15yr_tutorial

Tutorial to go with the 15 year isotropic GWB analysis

Generic badge MIT license made-with-python Generic badge

Authors: Aaron Johnson, Pat Meyers, Stephen Taylor, Sarah Vigeland for the NANOGrav Collaboration

Please send questions about this tutorial to aaron.johnson (at) nanograv.org

Important Note About Data!

If you want to use our data for publications, the full data products are available for download under the data_release section.

Data used in the tutorials section of this repository have been reduced to make them available on GitHub and may not reproduce the results of the 15-year analysis exactly.

Installing PTA software

  1. Install Miniconda

  2. To install a new environment: conda create -n enterprise -c conda-forge enterprise_extensions la_forge h5pulsar jupyterlab pyarrow pytables seaborn

    • Newer Macs (M1/M2/M3) can make a conda environment and install enterprise by first following the instructions here to create an environment and then using conda install -c conda-forge enterprise_extensions la_forge h5pulsar jupyterlab pyarrow pytables seaborn
  3. This will create a conda environment that can be activated by conda activate enterprise

  4. Next run jupyter notebook

  5. Set the Kernel

    • when opening a new notebook: click New and select Python [conda env:enterprise]
    • when opening an existing notebook (like these tutorials): click Kernel --> Change Kernel --> Python [conda env:enterprise]

Tutorials

These tutorials are split into several different files. The topic of each tutorial is shown below. These are roughly in the order that they should be viewed in to get a complete picture of how the isotropic GWB analysis is performed in NANOGrav.

Exploring NANOGrav Data

NANOGrav uses data files that may be unfamiliar to users that are new to pulsar timing or data analysis. Here, we investigate what information exists inside each par and tim file, how to load them into enterprise, and what information enterprise can use.

PTA GWB Analysis

This tutorial is split into two files, one for parameter estimation and one for model selection. We work through the Bayesian analysis of the full NANOGrav 15 year data to show what values each of the searched-over parameters prefers. Next, we show how to compare models and compute Bayes factors.

Optimal Statistic Analysis

This tutorial gives an introduction to frequentist methods we can use to look for an isotropic GWB.