Esempio n. 1
0
def fixed_pointing_info():
    filename = "$GAMMAPY_DATA/cta-1dc/data/baseline/gps/gps_baseline_110380.fits"
    return FixedPointingInfo.read(filename)
Esempio n. 2
0
# Define map geometry
axis = MapAxis.from_edges(np.logspace(-1.0, 1.0, 10),
                          unit="TeV",
                          name="energy",
                          interp="log")
geom = WcsGeom.create(skydir=(0, 0),
                      binsz=0.02,
                      width=(5, 4),
                      coordsys="GAL",
                      axes=[axis])

# In[ ]:

# Define some observation parameters
# We read in the pointing info from one of the 1dc event list files as an example
pointing = FixedPointingInfo.read(
    "$GAMMAPY_DATA/cta-1dc/data/baseline/gps/gps_baseline_110380.fits")

livetime = 1 * u.hour
offset_max = 2 * u.deg
offset = Angle("2 deg")

# In[ ]:

exposure = make_map_exposure_true_energy(pointing=pointing.radec,
                                         livetime=livetime,
                                         aeff=irfs["aeff"],
                                         geom=geom)
exposure.slice_by_idx({
    "energy": 3
}).plot(add_cbar=True)
Esempio n. 3
0
 def setup_class(cls):
     filename = "$GAMMAPY_DATA/tests/pointing_table.fits.gz"
     cls.fpi = FixedPointingInfo.read(filename)
Esempio n. 4
0
 def setup_class(cls):
     filename = "$GAMMAPY_DATA/tests/hess_event_list.fits"
     cls.fpi = FixedPointingInfo.read(filename)