Continuous Wave Analysis of the NANOGrav 15-Year Dataset

This repository contains the data used in the analysis detailed within "<em>The NANOGrav 15-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries</em>" (DOI arXiv:2306.16222). It also provides instruction on how one could reproduce those results on their own. <code>jar</code> folder: contains pickled <code>enterprise</code> pulsar objects that contain each pulsar's TOA data and timing model information <code>v1p1_all_dict.json</code>: dictionary with measured white noise and red noise parameters for each pulsar - used to fix white noise parameters when running QuickCW <code>emp_dist_15yr_v1p1_bence_my_run_v4.3.pkl</code>: <code>enterprise_extension</code> <code>EmpiricalDistribution2D</code> objects for the amplitude and spectral index of each pulsar's red noise - used to make informed red noise proposals when running QuickCW <code>pulsar_distances_15yr.pkl</code>: dictionary containing 3 element list for each pulsar specifying its distance [kpc], statistical error of distance [kpc], and method of distance measurement [PX for parallax or DM for dispersion measure] <code>15yr_quickCW_detection.h5</code>: HDF5 file containing results from QuickCW detection runs (log uniform amplitude prior). Includes the following data: MCMC samples (<code>samples_cold</code>) log likelihood values corresponding to those samples (<code>log_likelihood</code> ) temperature ladder used for the parallel tempering (<code>Ts</code>) names of parameters corresponding to each column in the samples array (<code>par_names</code>) rate of different MCMC jumps being accepted (<code>acc_fraction</code>) diagonal fisher matrix used for some proposals (<code>fisher_diag</code>) <code>15yr_quickCW_UL.h5</code>: Same as <code>15yr_quickCW_detection.h5</code> but for upper limits based on runs with uniform amplitude prior. <code>15yr_cw_3d_limits_v4.npz</code>: npz file containing limits as a function of frequency and sky location. Specifically it has the following arrays: <code>F_edges</code>: array specifying the edges of the frequency bins <code>strain_limit_skies</code>: 2D array containing the strain upper limits as a function of frequency bin and sky pixel <code>strain_limit_sky_sigmas</code>: 2D array containing 1-sigma statistical errors on the strain upper limits <code>dist_limit_skies</code>: 2D array containing the luminosity distance limits in Mpc as a function of frequency bin and sky pixel <code>dist_limit_sky_sigmas</code>: 2D array containing 1-sigma statistical errors on the distance limits in Mpc