def test_create_srcdep_dl3_energy_dependent_cuts(
temp_dir_observed_srcdep_files, observed_srcdep_dl2_file
):
"""
Generating a source-dependent DL3 file from a test DL2 files and test IRF file,
using energy-dependent cuts
"""
from lstchain.tools.lstchain_create_dl3_file import DataReductionFITSWriter
irf_file = temp_dir_observed_srcdep_files / "irf_edep.fits.gz"
assert (
run_tool(
DataReductionFITSWriter(),
argv=[
f"--input-dl2={observed_srcdep_dl2_file}",
f"--output-dl3-path={temp_dir_observed_srcdep_files}",
f"--input-irf={irf_file}",
"--source-name=Crab",
"--source-ra=83.633deg",
"--source-dec=22.01deg",
"--source-dep",
"--overwrite",
],
cwd=temp_dir_observed_srcdep_files,
)
== 0
)
def test_create_srcdep_dl3(
temp_dir_observed_srcdep_files, observed_srcdep_dl2_file, simulated_srcdep_irf_file
):
"""
Generating a source-dependent DL3 file from a test DL2 files and test IRF file
"""
from lstchain.tools.lstchain_create_dl3_file import DataReductionFITSWriter
from lstchain.paths import dl2_to_dl3_filename
assert (
run_tool(
DataReductionFITSWriter(),
argv=[
f"--input-dl2={observed_srcdep_dl2_file}",
f"--output-dl3-path={temp_dir_observed_srcdep_files}",
f"--input-irf={simulated_srcdep_irf_file}",
"--source-name=Crab",
"--source-ra=83.633deg",
"--source-dec=22.01deg",
"--source-dep",
"--overwrite",
],
cwd=temp_dir_observed_srcdep_files,
)
== 0
)
hdulist = fits.open(
temp_dir_observed_srcdep_files / dl2_to_dl3_filename(observed_srcdep_dl2_file)
)
ra = hdulist[1].data['RA']
dec = hdulist[1].data['DEC']
np.testing.assert_allclose(ra, 83.63, atol=1e-2)
np.testing.assert_allclose(dec, 22.01, atol=1e-2)
def test_create_dl3_with_config(temp_dir_observed_files, observed_dl2_file):
"""
Generating an DL3 file from a test DL2 files and test IRF file, using
a config file
"""
from lstchain.tools.lstchain_create_dl3_file import DataReductionFITSWriter
irf_file = temp_dir_observed_files / "fe_irf.fits.gz"
config_file = os.path.join(os.getcwd(), "docs/examples/dl3_tool_config.json")
assert (
run_tool(
DataReductionFITSWriter(),
argv=[
f"--input-dl2={observed_dl2_file}",
f"--output-dl3-path={temp_dir_observed_files}",
f"--input-irf={irf_file}",
"--source-name=Crab",
"--source-ra=83.633deg",
"--source-dec=22.01deg",
f"--config={config_file}",
"--overwrite",
],
cwd=temp_dir_observed_files,
)
== 0
)
def test_create_dl3(temp_dir_observed_files, observed_dl2_file, simulated_irf_file):
"""
Generating an DL3 file from a test DL2 files and test IRF file
"""
from lstchain.tools.lstchain_create_dl3_file import DataReductionFITSWriter
assert (
run_tool(
DataReductionFITSWriter(),
argv=[
f"--input-dl2={observed_dl2_file}",
f"--output-dl3-path={temp_dir_observed_files}",
f"--input-irf={simulated_irf_file}",
"--source-name=Crab",
"--source-ra=83.633deg",
"--source-dec=22.01deg",
"--overwrite",
],
cwd=temp_dir_observed_files,
)
== 0
)
def test_create_dl3_energy_dependent_cuts(
temp_dir_observed_files, observed_dl2_file
):
"""
Generating an DL3 file from a test DL2 files and test IRF file, using
energy dependent cuts. Here the previously created IRF is used.
"""
from lstchain.tools.lstchain_create_dl3_file import DataReductionFITSWriter
from gammapy.data import Observation
irf_file = temp_dir_observed_files / "pnt_irf.fits.gz"
dl2_name = observed_dl2_file.name
observed_dl3_file = temp_dir_observed_files / dl2_name.replace('dl2', 'dl3')
observed_dl3_file = observed_dl3_file.with_suffix(".fits")
assert (
run_tool(
DataReductionFITSWriter(),
argv=[
f"--input-dl2={observed_dl2_file}",
f"--output-dl3-path={temp_dir_observed_files}",
f"--input-irf={irf_file}",
"--source-name=Crab",
"--source-ra=83.633deg",
"--source-dec=22.01deg",
"--overwrite",
],
cwd=temp_dir_observed_files,
)
== 0
)
assert Observation.read(
event_file=observed_dl3_file, irf_file=irf_file
).obs_id == 2008