GammaLib

#! /usr/bin/env python 

import gammalib as gl

import obsutils as obsu

obstime=0.5*3600

write=True

#-----------------------------------------------------

#Create one observation (run)

#-----------------------------------------------------

ra_degree=83.633212         #Crab position from NED

dec_degree=22.014460 

offset_degree=-1            #Observing offset in dec from Crab in deg

pntdir = gl.GSkyDir()

pntdir.radec_deg(ra_degree, dec_degree+offset_degree) 

ob = obsu.set(pntdir,tstart=0.0,duration=obstime,deadc=0.95, \

              emin=0.1, emax=100.0, rad=8.0, irf="irf_file.fits",\

              caldb="$PWD/caldb/data/cta/aar/bcf/DESY20140105_50h/")

ob.id("run_1")

#----------------------------------------------------

# Create background model

#----------------------------------------------------

spec = gl.GModelSpectralPlaw(1,0,gl.GEnergy(1.0,"TeV"))

filename="$PWD/caldb/data/cta/aar/bcf/DESY20140105_50h/irf_file.fits"

bck = gl.GCTAModelIrfBackground(spec)

bck.name("Background_run_1")

bck.instruments("CTA")

#----------------------------------------------------

# Create Crab nebula model

#----------------------------------------------------

crabdir = gl.GSkyDir()

crabdir.radec_deg(ra_degree, dec_degree)

crab_spatial = gl.GModelSpatialEllipticalDisk(crabdir,0.12,0.8,90.)

crab_spectrum = gl.GModelSpectralPlaw(2.83e-17, -2.62,gl.GEnergy(1,"TeV") )

crab = gl.GModelSky(crab_spatial, crab_spectrum)

crab.name("CrabNebula")

crab.ids("")

#----------------------------------------------------

# Create GObservations and add models

#----------------------------------------------------

degperpix=0.05

npix=200

obs = gl.GObservations()

obs.append(ob)

obs.models().append(bck)

obs.models().append(crab)

obssim = obsu.sim(obs, log=write, debug=False, edisp=False, seed=0)

obsu.modmap(obssim, eref=0.1, proj="TAN", coord="EQU", xval=ra_degree, yval=dec_degree+offset_degree, \

                    binsz=degperpix, nxpix=npix, nypix=npix, outname="crabextsim_modmap.fits")