ridge/scripts/README.md

Don't forget to activate environment: source <MY PATH>/venv/bin/activate

The script number indicates the sequence of execution. The results are placed in ridge/products/.

The corresponding shell scripts run Python code with different parameters.

01_bgdmodel.py/01_bgdmodel.sh

Calibrates IBIS/ISGRI background model

01_crabmodel.py/01_crabmodel.sh

Calibrates Crab detector count rate model.

01_crabmodel_plot_poly.py

Plot long-term Crab detector count rate approximated by cubic polynomial function (Fig. B8 in paper).

01_crabmodel_plot_sys.py

Plots the normalized distribution of the residuals between the IBIS/ISGRI Crab count rate and the corresponding polynomial fit (Fig. B9 in the paper). The distribution is approximated with a Gaussian function.

02_grxe_resid.py/02_grxe_resid.sh

Calculates difference between detector count rate and that predicted by background model in mCrab units for GAL (Galaxy), BKG (Extragalactic), or ALL (all-sky) regions.

02_grxe_resid_plot.py

Plots normalized distribution of the relative residuals of the background model obtained in three energy bands for BKG region (Figs. A6 and A7 in the paper).

03_grxe_map.py

Makes the map of the residuals in mCrab units (not shown in the paper).

03_grxe_flux.py/03_grxe_flux.sh

Estimates flux in regions defined in config.py, and places result in products/Flux, e.g.:

cat ../products/Flux/GB.dat
GB A01  25.00   60.00 154.886019 3.087652
GB A02  60.00   80.00 96.472798 9.218161
GB A03  80.00  200.00 102.207302 13.291950

where GB means Galactic Bulge region (see Table 1 in the paper).

03_grxe_galprof.py/03_grxe_galprof.sh and 03_grxe_galprof_plot.py

Makes Galactic longitude profile and plots it.

03_grxe_spec.py

Extract spectrum of selected sky region, e.g.:

./03_grxe_spec.py GB or ./03_grxe_spec.py L+20

If you do not specify the sky region, the script makes spectra for all regions.

Sky regions are defined in config.py.

Spectra are stored in ../products/Spectra.

To load the spectrum, e.g. Galactic Bulge, do the following steps:

cd ../products/Spectra
cp ../../models/xspec/* .

Then data and spectral models can be loaded by the following scripts in xspec:

#xspec

@load_gb_cutoffpl.xcm
@load_gb_polar.xcm

@load_l+20_cutoffpl.xcm
@load_l+20_polar.xcm

@load_l-20_cutoffpl.xcm
@load_l-20_polar.xcm

03_grxe_syserr.py/03_plot_bkg.py/03_plot_skyreg.py

Prints statistics and plots distribution of residuals for extragalactic data, and for selected sky regions.