def dataset():
energy_axis = MapAxis.from_bounds(1,
10,
nbin=3,
unit="TeV",
name="energy",
interp="log")
geom = WcsGeom.create(skydir=(0, 0),
binsz=0.05,
width="5 deg",
frame="galactic",
axes=[energy_axis])
geom_true = geom.copy()
geom_true.axes[0].name = "energy_true"
dataset = get_map_dataset(geom=geom,
geom_etrue=geom_true,
edisp="edispmap",
name="test")
dataset.background /= 400
dataset.gti = GTI.create(start=0 * u.s,
stop=1000 * u.s,
reference_time="2000-01-01")
return dataset
def dataset():
position = SkyCoord(0.0, 0.0, frame="galactic", unit="deg")
energy_axis = MapAxis.from_bounds(1,
10,
nbin=3,
unit="TeV",
name="energy",
interp="log")
spatial_model = GaussianSpatialModel(lon_0="0 deg",
lat_0="0 deg",
sigma="0.2 deg",
frame="galactic")
spectral_model = PowerLawSpectralModel(amplitude="1e-11 cm-2 s-1 TeV-1")
t_min = 0 * u.s
t_max = 1000 * u.s
time = np.arange(t_max.value) * u.s
tau = u.Quantity("2e2 s")
norm = np.exp(-time / tau)
table = Table()
table["TIME"] = time
table["NORM"] = norm / norm.max()
t_ref = Time("2000-01-01")
table.meta = dict(MJDREFI=t_ref.mjd, MJDREFF=0, TIMEUNIT="s")
temporal_model = LightCurveTemplateTemporalModel(table)
skymodel = SkyModel(
spatial_model=spatial_model,
spectral_model=spectral_model,
temporal_model=temporal_model,
)
geom = WcsGeom.create(skydir=position,
binsz=1,
width="5 deg",
frame="galactic",
axes=[energy_axis])
gti = GTI.create(start=t_min, stop=t_max, reference_time=t_ref)
geom_true = geom.copy()
geom_true.axes[0].name = "energy_true"
dataset = get_map_dataset(sky_model=skymodel,
geom=geom,
geom_etrue=geom_true,
edisp="edispmap")
dataset.gti = gti
return dataset
def dataset():
position = SkyCoord(0.0, 0.0, frame="galactic", unit="deg")
energy_axis = MapAxis.from_bounds(1,
10,
nbin=3,
unit="TeV",
name="energy",
interp="log")
spatial_model = GaussianSpatialModel(lon_0="0 deg",
lat_0="0 deg",
sigma="0.2 deg",
frame="galactic")
spectral_model = PowerLawSpectralModel(amplitude="1e-11 cm-2 s-1 TeV-1")
skymodel = SkyModel(spatial_model=spatial_model,
spectral_model=spectral_model)
geom = WcsGeom.create(skydir=position,
binsz=1,
width="5 deg",
frame="galactic",
axes=[energy_axis])
t_min = 0 * u.s
t_max = 30000 * u.s
gti = GTI.create(start=t_min, stop=t_max)
geom_true = geom.copy()
geom_true.axes[0].name = "energy_true"
dataset = get_map_dataset(sky_model=skymodel,
geom=geom,
geom_etrue=geom_true,
edisp=True)
dataset.gti = gti
return dataset