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