def from_geoms(
cls,
geom,
geom_exposure,
geom_psf,
geom_edisp,
reference_time="2000-01-01",
name="",
**kwargs,
):
"""
Create a MapDatasetOnOff object with zero filled maps according to the specified geometries
Parameters
----------
geom : `Geom`
geometry for the counts, counts_off, acceptance and acceptance_off maps
geom_exposure : `Geom`
geometry for the exposure map
geom_psf : `Geom`
geometry for the psf map
geom_edisp : `Geom`
geometry for the energy dispersion map
reference_time : `~astropy.time.Time`
the reference time to use in GTI definition
name : str
Name of the dataset.
Returns
-------
empty_maps : `MapDatasetOnOff`
A MapDatasetOnOff containing zero filled maps
"""
maps = {}
for name in ["counts", "counts_off", "acceptance", "acceptance_off"]:
maps[name] = Map.from_geom(geom, unit="")
exposure = Map.from_geom(geom_exposure, unit="m2 s")
edisp = EDispMap.from_geom(geom_edisp)
psf = PSFMap.from_geom(geom_psf)
gti = GTI.create([] * u.s, [] * u.s, reference_time=reference_time)
mask_safe = Map.from_geom(geom, dtype=bool)
return cls(
counts=maps["counts"],
counts_off=maps["counts_off"],
acceptance=maps["acceptance"],
acceptance_off=maps["acceptance_off"],
exposure=exposure,
psf=psf,
edisp=edisp,
gti=gti,
mask_safe=mask_safe,
name=name,
**kwargs,
)
def from_geoms(
cls,
geom,
geom_exposure,
geom_psf,
geom_edisp,
reference_time="2000-01-01",
name="",
**kwargs,
):
"""
Create a MapDataset object with zero filled maps according to the specified geometries
Parameters
----------
geom : `Geom`
geometry for the counts and background maps
geom_exposure : `Geom`
geometry for the exposure map
geom_psf : `Geom`
geometry for the psf map
geom_edisp : `Geom`
geometry for the energy dispersion map
reference_time : `~astropy.time.Time`
the reference time to use in GTI definition
name : str
Name of the dataset.
Returns
-------
empty_maps : `MapDataset`
A MapDataset containing zero filled maps
"""
counts = Map.from_geom(geom, unit="")
background = Map.from_geom(geom, unit="")
background_model = BackgroundModel(background)
exposure = Map.from_geom(geom_exposure, unit="m2 s")
edisp = EDispMap.from_geom(geom_edisp)
psf = PSFMap.from_geom(geom_psf)
gti = GTI.create([] * u.s, [] * u.s, reference_time=reference_time)
mask_safe = Map.from_geom(geom, unit="", dtype=bool)
return cls(
counts=counts,
exposure=exposure,
psf=psf,
edisp=edisp,
background_model=background_model,
gti=gti,
mask_safe=mask_safe,
name=name,
**kwargs,
)