def test_create_irf_point_like_srcdep(
temp_dir_observed_srcdep_files, simulated_srcdep_dl2_file
):
"""
Generating point-like source-dependent IRF file from a test DL2 files
"""
from lstchain.tools.lstchain_create_irf_files import IRFFITSWriter
irf_file = temp_dir_observed_srcdep_files / "irf.fits.gz"
assert (
run_tool(
IRFFITSWriter(),
argv=[
f"--input-gamma-dl2={simulated_srcdep_dl2_file}",
f"--output-irf-file={irf_file}",
"--point-like",
"--source-dep",
"--overwrite",
],
cwd=temp_dir_observed_srcdep_files,
)
== 0
)
with fits.open(irf_file) as hdul:
for hdu in hdul[1:]:
assert 'AL_CUT' in hdu.header
assert isinstance(hdu.header['AL_CUT'], float)
def test_create_irf_full_enclosure_with_config(temp_dir_observed_files, simulated_dl2_file):
"""
Generating full enclosure IRF file from a test DL2 files, using
a config file
"""
from lstchain.tools.lstchain_create_irf_files import IRFFITSWriter
irf_file = temp_dir_observed_files / "fe_irf.fits.gz"
config_file = os.path.join(os.getcwd(), "./docs/examples/irf_tool_config.json")
assert (
run_tool(
IRFFITSWriter(),
argv=[
f"--input-gamma-dl2={simulated_dl2_file}",
f"--input-proton-dl2={simulated_dl2_file}",
f"--input-electron-dl2={simulated_dl2_file}",
f"--output-irf-file={irf_file}",
f"--config={config_file}",
"--overwrite",
],
cwd=temp_dir_observed_files,
)
== 0
)
def test_create_irf_point_like(temp_dir_observed_files, simulated_dl2_file):
"""
Generating point-like IRF file from a test DL2 files
"""
from lstchain.tools.lstchain_create_irf_files import IRFFITSWriter
irf_file = temp_dir_observed_files / "pnt_irf.fits.gz"
assert (
run_tool(
IRFFITSWriter(),
argv=[
f"--input-gamma-dl2={simulated_dl2_file}",
f"--input-proton-dl2={simulated_dl2_file}",
f"--input-electron-dl2={simulated_dl2_file}",
f"--output-irf-file={irf_file}",
"--point-like",
"--overwrite",
],
cwd=temp_dir_observed_files,
)
== 0
)
with fits.open(irf_file) as hdul:
for hdu in hdul[1:]:
assert 'RAD_MAX' in hdu.header
assert isinstance(hdu.header['RAD_MAX'], float)
def test_create_irf_point_like_srcdep_energy_dependent_cuts(
temp_dir_observed_srcdep_files, simulated_srcdep_dl2_file
):
"""
Generating point-like source-dependent IRF file from a test DL2 files,
using energy-dependent cuts
"""
from lstchain.tools.lstchain_create_irf_files import IRFFITSWriter
from astropy.table import QTable
irf_file = temp_dir_observed_srcdep_files / "irf_edep.fits.gz"
assert (
run_tool(
IRFFITSWriter(),
argv=[
f"--input-gamma-dl2={simulated_srcdep_dl2_file}",
f"--output-irf-file={irf_file}",
"--point-like",
"--source-dep",
"--energy-dependent-gh",
"--energy-dependent-alpha",
"--overwrite",
],
cwd=temp_dir_observed_srcdep_files,
)
== 0
)
gh_cuts = QTable.read(irf_file, hdu="GH_CUTS")
assert isinstance(gh_cuts.meta["GH_EFF"], float)
al_cuts = QTable.read(irf_file, hdu="AL_CUTS")
assert isinstance(al_cuts.meta["AL_CONT"], float)
def test_create_irf_point_like_energy_dependent_cuts(
temp_dir_observed_files, simulated_dl2_file
):
"""
Generating point-like IRF file from a test DL2 files, using
energy-dependent cuts
"""
from lstchain.tools.lstchain_create_irf_files import IRFFITSWriter
from gammapy.irf import RadMax2D
irf_file = temp_dir_observed_files / "pnt_irf.fits.gz"
assert (
run_tool(
IRFFITSWriter(),
argv=[
f"--input-gamma-dl2={simulated_dl2_file}",
f"--input-proton-dl2={simulated_dl2_file}",
f"--input-electron-dl2={simulated_dl2_file}",
f"--output-irf-file={irf_file}",
"--overwrite",
"--energy-dependent-gh",
"--point-like",
"--energy-dependent-theta",
"--DL3Cuts.max_theta_cut=1",
"--DL3Cuts.fill_theta_cut=1"
],
cwd=temp_dir_observed_files,
)
== 0
)
assert RadMax2D.read(irf_file, hdu="RAD_MAX")
def test_create_irf_full_enclosure(temp_dir_observed_files, simulated_dl2_file):
"""
Generating full enclosure IRF file from a test DL2 files
"""
from lstchain.tools.lstchain_create_irf_files import IRFFITSWriter
irf_file = temp_dir_observed_files / "fe_irf.fits.gz"
assert (
run_tool(
IRFFITSWriter(),
argv=[
f"--input-gamma-dl2={simulated_dl2_file}",
f"--input-proton-dl2={simulated_dl2_file}",
f"--input-electron-dl2={simulated_dl2_file}",
f"--output-irf-file={irf_file}",
"--overwrite",
],
cwd=temp_dir_observed_files,
)
== 0
)