ctools

import gammalib as gl

import csbutterfly

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)

    bf = csbutterfly.csbutterfly(filled_obs)

    bf.source = "Crab"

    bf.m_emin = 0.1

    bf.m_emax = 100

    bf.outfile = "Crab.bf"

    bf.run()

    #ulimit = ulim.ulimit

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