ctools

import gammalib as gl

import cslike

import obsutils

import os

if __name__ == "__main__":

    dir = gl.GSkyDir()

    dir.radec_deg(83.6331,22.01)

    ob = obsutils.set(dir)

    obs = gl.GObservations()

    obs.append(ob)

    obs.models(os.environ["CTOOLS"]+"/test/data/crab.xml")

    filled_obs = obsutils.sim(obs)

    like = cslike.cslike(filled_obs)

    like.outfile = "crab_optmodel.xml"

    like["TSmethod"].string("UNBIASED")

    like.run()

    #print like.covar()