def from_hdulist(cls, hdulist, name=""):
"""Create map dataset from list of HDUs.
Parameters
----------
hdulist : `~astropy.io.fits.HDUList`
List of HDUs.
Returns
-------
dataset : `MapDataset`
Map dataset.
"""
init_kwargs = {}
init_kwargs["name"] = name
if "COUNTS" in hdulist:
init_kwargs["counts"] = Map.from_hdulist(hdulist, hdu="counts")
if "EXPOSURE" in hdulist:
init_kwargs["exposure"] = Map.from_hdulist(hdulist, hdu="exposure")
if "BACKGROUND" in hdulist:
background_map = Map.from_hdulist(hdulist, hdu="background")
init_kwargs["background_model"] = BackgroundModel(background_map)
if "EDISP_MATRIX" in hdulist:
init_kwargs["edisp"] = EnergyDispersion.from_hdulist(
hdulist, hdu1="EDISP_MATRIX", hdu2="EDISP_MATRIX_EBOUNDS"
)
if "PSF_KERNEL" in hdulist:
psf_map = Map.from_hdulist(hdulist, hdu="psf_kernel")
init_kwargs["psf"] = PSFKernel(psf_map)
if "MASK_SAFE" in hdulist:
mask_safe_map = Map.from_hdulist(hdulist, hdu="mask_safe")
init_kwargs["mask_safe"] = mask_safe_map.data.astype(bool)
if "MASK_FIT" in hdulist:
mask_fit_map = Map.from_hdulist(hdulist, hdu="mask_fit")
init_kwargs["mask_fit"] = mask_fit_map.data.astype(bool)
if "GTI" in hdulist:
gti = GTI.from_hdulist(hdulist, hdu="GTI")
init_kwargs["gti"] = gti
return cls(**init_kwargs)
def _read_ogip_hdulist(hdulist, hdu1="SPECTRUM", hdu2="EBOUNDS", hdu3="GTI"):
"""Create from `~astropy.io.fits.HDUList`."""
counts_table = Table.read(hdulist[hdu1])
ebounds = Table.read(hdulist[hdu2])
emin = ebounds["E_MIN"].quantity
emax = ebounds["E_MAX"].quantity
if hdu3 in hdulist:
gti = GTI.from_hdulist(hdulist, hdu3)
else:
gti = None
# Check if column are present in the header
quality = None
areascal = None
backscal = None
if "QUALITY" in counts_table.colnames:
quality = counts_table["QUALITY"].data
if "AREASCAL" in counts_table.colnames:
areascal = counts_table["AREASCAL"].data
if "BACKSCAL" in counts_table.colnames:
backscal = counts_table["BACKSCAL"].data
return dict(
data=counts_table["COUNTS"],
backscal=backscal,
energy_lo=emin,
energy_hi=emax,
quality=quality,
areascal=areascal,
livetime=counts_table.meta["EXPOSURE"] * u.s,
obs_id=counts_table.meta["OBS_ID"],
is_bkg=False,
gti=gti,
)