GammaLib

#! /usr/bin/env python

import gammalib

import ctools

import cscripts

from cscripts import obsutils

def test(irf, emin, emax):

    """

    """

    # Set pointing direction

    pntdir = gammalib.GSkyDir()

    pntdir.radec_deg(83.63, 22.01)

    # Set CTA Observation

    run = obsutils.set_obs(pntdir, emin=emin, emax=emax, rad=5.0, duration=21232, irf=irf, caldb='.')

    # Setup observation container

    obs = gammalib.GObservations()

    obs.append(run)

    # Simulate events

    sim = ctools.ctobssim(obs)

    sim['inmodel']   = 'model_src+bkg_CTA.xml'

    sim['outevents'] = 'events.fits'

    sim['chatter']   = 4

    sim['debug']     = True

    sim.execute()

    # Return

    return

if __name__ == '__main__':

    # Run tests

    #test('IRF_ASTRI_cut_mod.fits',   0.631, 398.0)

    test('IRF_ASTRI_cut_mod_2.fits', 0.631, 398.0)