Configuration file

The default configuration file is mcmc_config.txt in the main directory. If you would like to use a different file, you can enter the absolute path of your configuration file as an argument when running blazar_run_mcmc.py, such as

$ python blazar_run_mcmc.py /home/myconfigfile.txt

There is a file named mcmc_config_template.txt as an example. Make a copy of this file named mcmc_config.txt. This file is automatically in .gitignore since it changes locally.

Notes:

• Labels must be exactly as listed in the configuration file.

• data_file is the relative path from the Bjet_MCMC directory to the file with data.

Example of a configuration file

(file named mcmc_config.txt, can be modified):

# Configuration file for running mcmc

# description:
description=J1010

# folder label:
folder_label=J1010

# eic (True/False):
eic=False

# Data file:
data_file=real_data/J1010_SED_reduced.dat

# Number of steps:
n_steps=1000

# Number of walkers:
n_walkers=100

# Discard number:
discard = 200

# Parallel processing (True/False):
parallel=True

# If parallel processing, # of cores (will use # of cores - 1 if not specified)
cores=15

# use tau variability (boolean):
use_variability=False

# tau variability (in hours):
tau_variability = 24

# redshift
redshift = 0.143

# Custom alpha2 limits (True/False, <val1>, <val2>) val1 and val2 optional
custom_alpha2_limits=False

#to freeze a parameter to a given value, replace null with the values you want it to be frozen
delta       = null    #       doppler factor                  linear
K           = null    #       particle density [cm^-3]        log10
n1          = null    #       first particle index            linear
n2          = null    #       second particle index           linear
gamma_min   = null    #       low-energy cutoff               log10
gamma_max   = null    #       high-energy cutoff              log10
gamma_break = null    #       energy break                    log10
B           = null    #       magnetic field strength [G]     log10
R           = null    #       blob radius [cm]                log10
#-----------------------------------------------------------------------#
#-----------------------Additional params for EIC-----------------------#
bb_temp     = null    #       Black body temp of disk [K]     log10
l_nuc       = null    #       Nucleus luminosity [ergs/s]     log10
tau         = null    #       Frac of luminosity scattered    log10
blob_dist   = null    #       Distance of blob from SMBH[cm]  log10