ctools

import gammalib

import ctools

def overall():

    model_name           = "models/model_ics.xml" 

    events_name          = "events-ics.fits"

    caldb                = "prod2"    

    irf                  = "South_50h" 

    ra                   =   258.37678

    dec                  =   -39.781175

    rad_sim              =   5.0

    tstart               =    0.0

    tstop                =  32400.0 # 9hr

    emin                 =    0.1 

    emax                 =  100.0

    rad_select           =    5.0

    # Simulate events

    print "Generating event file... "

    sim = ctools.ctobssim()

    sim.logFileOpen()

    sim["inmodel"].filename(model_name)

    sim["outevents"].filename(events_name)

    sim["caldb"].string(caldb)

    sim["irf"].string(irf)

    sim["ra"].real(ra)

    sim["dec"].real(dec)

    sim["rad"].real(rad_sim)

    sim["tmin"].real(tstart)

    sim["tmax"].real(tstop)

    sim["emin"].real(emin)

    sim["emax"].real(emax)

    sim.execute()

    print(sim.obs()[0].events())

    return

if __name__ == '__main__':

    overall()