GammaLib

#! /usr/bin/env python

import os

import gammalib

import ctools

# Set globals

og_dol  = '/Users/jurgen/git/gammalib/gammalib/inst/spi/test/data/obs/og_spi.fits'

rsp_dol = '/project-data/integral/ic/spi/rsp/spi_irf_grp_0021.fits'

# ====== #

# ctlike #

# ====== #

def ctlike():

    """

    """

    # Create observations

    og  = gammalib.GSPIObservation(og_dol)

    rsp = gammalib.GSPIResponse(rsp_dol)

    rsp.set(og)

    og.response(rsp)

    obs = gammalib.GObservations()

    obs.append(og)

    # Set sky model

    spatial   = gammalib.GModelSpatialPointSource(83.6331, +22.0145)

    spatial['RA'].free()

    spatial['DEC'].free()

    spectral  = gammalib.GModelSpectralPlaw(1.0e-5, -2.0, gammalib.GEnergy(0.1, 'MeV'))

    sky_model = gammalib.GModelSky(spatial, spectral)

    sky_model.name('Crab')

    sky_model.tscalc(True)

    # Set background model

    bgm_model = gammalib.GSPIModelDataSpace(og, 'GEDSAT', 'orbit,dete', 0)

    # Create model and append it to observations

    models = gammalib.GModels()

    models.append(sky_model)

    models.append(bgm_model)

    obs.models(models)

    # Run ctlike

    like = ctools.ctlike(obs)

    like['outmodel'] = 'ctlike.xml'

    like['chatter']  = 3

    like['debug']    = True

    like.logFileOpen()

    like.execute()

    # Return

    return

# ======================== #

# Main routine entry point #

# ======================== #

if __name__ == '__main__':

    # Run ctlike

    ctlike()