paulthebaker/la_forge

La Forge

Pulsar Timing Array Bayesian Data Visualization

Graphic Credit: Stewart Vernon, via Deviant Art

Python package for conveniently plotting results from pulsar timing array bayesian analyses. Many of the functions are best used with enterprise outputs.

Not yet available on PyPI, please use

pip install git+https://github.com/Hazboun6/la_forge@master

to install or run the setup.py script once cloned to your pc.

• Free software: MIT license
• Documentation: https://la-forge.readthedocs.io.

Example code

from la_forge import rednoise
from la_forge.core import Core
from la_forge import utils

normal_ul_dir = '../BF_standard/DE436/'
free_spec_ul_dir = '../BF_free_spec/DE436/'

a = Core('plaw',chaindir=normal_ul_dir)
b = Core('free_spec',chaindir=free_spec_ul_dir)

tspan = 11.4*365.25*24*3600

a.set_rn_freqs(Tspan=tspan)
b.set_rn_freqs(Tspan=tspan)

compare = [a,b]
plot_filename = './noise_model_plots.png'
Colors = ['blue','red']
Labels = ['PTA PLaw', 'PTA Free Spec']

rednoise.plot_rednoise_spectrum(pulsar=psr, cores=compare, chaindir=chaindir,
                                show_figure=True, rn_type='', verbose=False,
                                Tspan=tspan, Colors=Colors, n_plaw_realizations=100,
                                labels=Labels, plotpath=plot_filename)

Features

• TODO

Credits

This package was created with Cookiecutter and the audreyr/cookiecutter-pypackage project template.