GammaLib

#! /usr/bin/env python 

import gammalib as gl

from ctools 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_obs(pntdir,tstart=0.0,duration=obstime,deadc=0.95, \

                  emin=0.1, emax=100.0, rad=8.0, irf="DESY20140105_50h",\

                  caldb="aar")

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.GModelSpatialPointSource(crabdir)

crab_spectrum = gl.GModelSpectralLogParabola(3.2e-17,-2.47,gl.GEnergy(1,"TeV"),-0.24) 

crab = gl.GModelSky(crab_spatial, crab_spectrum)

crab.name("CrabNebula")

crab.ids("")

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

# Create GObservations and add models

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

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.fit(obssim,log=True, debug=False, edisp=False)