pyADfit 0.0.1

A nested sampling approach to quasi-stellar object (QSO) accretion disc fitting.

pyADfit

A nested sampling approach to quasi-stellar object (QSO) accretion disc fitting.

This repository contains a Python module for modelling accretion discs around astrophysical objects. The module provides functions to calculate physical quantities related to accretion disks and perform parameter estimation using observational data. The accretion disc model is the alpha-disc model (see Shakura & Sunyaev 1976), while the parameter estimation can be performed either with Nessai, Raynest or CPnest.

Dependencies

• numpy
• scipy
• matplotlib
• raynest
• CPNest
• nessai
• h5py
• pandas

Installation

Clone this repository to your local machine:

git clone https://github.com/FabioRigamonti/pyADfit.git

Install the required dependencies using pip:

pip install -r requirements.txt

Or install it directly with pip:

pip install pyADfit

Usage

To fit quasar accretion disc data, follow these steps:

1. Define your input data in a text file with three columns: x-data [nu, i.e. frequency], y-data [log10 nu*Lnu], and y-errors.
2. Create a YAML configuration file specifying the hyperparameters, see "config.yaml" in the example directory, fitting parameters, and other settings.
3. Define your own "read_data" function to read and the path to the configuration file
4. Import the "read_config_and_launch" function from "disc_launch"
5. Run the parameter estimation by calling the "read_config_and_launch"

Example

from disc_launch import read_config_and_launch
import matplotlib.pyplot as plt 

def read_data(file_path):
  your function here

  return xdata,ydata,yerr

if __name__=='__main__':
  config_path = './config.yaml'

  post_df,best = read_config_and_launch(config_path,read_data)
  plt.show()