ctools

import gammalib as gl

import csupperlimit

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)

    ulim = csupperlimit.csupperlimit(filled_obs)

    ulim.source = "Crab"

    ulim.parname="Prefactor"

    ulim.run()

    ulimit = ulim.ulimit

    print "Retrieved 95% U.L. to "+str("%.2e" % ulimit)+" 1/cm2/s/MeV"