Skip to content

Zstone19/pypetal

Repository files navigation

pyPETaL: A Pipeline for Estimating AGN Time Lags

Documentation Status pre-commit Workflow Status License codecov

pypi-ver python-ver pypi-downloads

ascl:2401.004

pyPETaL is a time-series data analysis pipeline for AGN reverberation mapping (RM) data. It combines multiple different popular softwares using for AGN RM analysis, including PyCCF, PyZDCF, JAVELIN, and PyROA. This pipeline also implements outlier rejection using Damped Random Walk Gaussian proces fitting, and detrending through the LinMix algorithm. pyPetal implements a weighting scheme (for all modules) in order to mitigate aliasing in peaks of time lag distributions between light curves.

pyPetal is very flexible, with almost every argument for each module allowing user input. pyPetal is designed to work with any combination of modules being run, allowing it to scale from the simplest to the most complex of projects.

Installation

pyPetal

pyPetal is available on PyPI and can be installed with pip:

    pip install pypetal

Or, if you want to install the latest development version:

    git clone https://github.com/Zstone19/pypetal.git
    cd pypetal
    pip install .

PLIKE is an optional algorithm that is used in pyPetal. There is a script available in the main directory to install and compile PLIKE (assuming that gfortran is installed). To install PLIKE, run the following command:

    sh build_plike.sh

Detrending

pyPetal offers the option to detrend the input light curves via the LinMix algorithm. This package is not offered in the base version of pyPetal, but can be installed with pip:

    pip install "linmix @ git+https://github.com/jmeyers314/linmix.git"

Or with poetry:

    poetry install --with extra

MICA2

pyPetal also offers MICA2 as an optional module. This package is more complex to install than the others - to find out more, read the MICA2 installation guide.

poetry can help to install some easily-installable python dependencies for MICA2:

    poetry install --with extra

The functionality and inputs of the MICA2 module are identical to the original software, so the best way to learn how this module functions is to read the MICA2 documentation.

NOTE: Seeing as this is an optional module, pyPetal will still compile and run if MICA2's dependencies are not installed. Each pyPetal run will assume run_mica2=False.

pyPetal and JAVELIN

The JAVELIN software used in pyPetal runs on Python 2, though the bulk of pyPetal software relies on Python >=3.8. To circumvent this issue, pyPetal has a JAVELIN "module" (pypetal-jav) that can be installed as a separate package and used in conjunction with pyPetal, in a Python 2 environment.

External requirements (not installed by pip or setup.py):

    A Fortran compiler (>F90)

pyPetal-jav is available on PyPI and can be installed with pip:

    pip install --no-deps pypetal-jav
    pip install pypetal-jav

Or, if you want to install the latest development version:

    git clone https://github.com/Zstone19/pypetal-jav.git
    cd pypetal-jav
    pip install .

NOTE: The user may need to install NumPy before installing pyPetal-jav through pip or setup.py. This is because JAVELIN requires NumPy in order to be installed. This can be done with:

    pip install numpy

Citing pyPetal

To cite the pyPetal code itself, use the ASCL reference:

    Stone Z., 2024, version xxxx, Astrophysics Source Code Library, record ascl:2401.004

Cit the paper pyPetal was used in: Shen et al. (2023)