|
1
|
|
|
from asgardpy.analysis import AsgardpyAnalysis |
|
2
|
|
|
|
|
3
|
|
|
|
|
4
|
|
|
def test_dataset1d(base_config_1d): |
|
5
|
|
|
"""Test for creating 1D stacked DL4 dataset.""" |
|
6
|
|
|
|
|
7
|
|
|
from gammapy.datasets import SpectrumDatasetOnOff |
|
8
|
|
|
|
|
9
|
|
|
analysis = AsgardpyAnalysis(base_config_1d) |
|
10
|
|
|
|
|
11
|
|
|
analysis.config.dataset1d.instruments[0].dataset_info.safe_mask.methods = ["custom-mask"] |
|
12
|
|
|
analysis.config.dataset1d.instruments[0].dataset_info.safe_mask.parameters = {} |
|
13
|
|
|
analysis.config.dataset1d.instruments[0].dataset_info.safe_mask.parameters["min"] = "200 GeV" |
|
14
|
|
|
analysis.config.dataset1d.instruments[0].dataset_info.safe_mask.parameters["max"] = "10 TeV" |
|
15
|
|
|
|
|
16
|
|
|
analysis.get_1d_datasets() |
|
17
|
|
|
|
|
18
|
|
|
assert isinstance(analysis.datasets[0], SpectrumDatasetOnOff) |
|
19
|
|
|
assert round(analysis.datasets[0].energy_range[-1].data[0][0]) == 10 |
|
20
|
|
|
|
|
21
|
|
|
analysis_1 = AsgardpyAnalysis(base_config_1d) |
|
22
|
|
|
analysis_1.config.dataset1d.instruments[0].input_dl3[0].glob_pattern.pop("dl3_files", None) |
|
23
|
|
|
analysis_1.config.dataset1d.instruments[0].input_dl3[0].glob_pattern["dl3_files"] = "data/hess_*fits.gz" |
|
24
|
|
|
analysis_1.run(["datasets-1d"]) |
|
25
|
|
|
|
|
26
|
|
|
assert len(analysis_1.datasets) == 2 |
|
27
|
|
|
|
|
28
|
|
|
|
|
29
|
|
|
def test_dataset1d_no_stack(base_config_1d): |
|
30
|
|
|
"""Test for creating 1D unstacked DL4 dataset.""" |
|
31
|
|
|
|
|
32
|
|
|
analysis = AsgardpyAnalysis(base_config_1d) |
|
33
|
|
|
|
|
34
|
|
|
analysis.config.general.stacked_dataset = False |
|
35
|
|
|
analysis.config.dataset1d.instruments[0].dataset_info.background.region_finder_method = "reflected" |
|
36
|
|
|
analysis.config.dataset1d.instruments[0].dataset_info.background.parameters = {} |
|
37
|
|
|
analysis.config.dataset1d.instruments[0].dataset_info.safe_mask.methods = [] |
|
38
|
|
|
|
|
39
|
|
|
analysis.get_1d_datasets() |
|
40
|
|
|
|
|
41
|
|
|
assert len(analysis.datasets) == 4 |
|
42
|
|
|
assert int(round(analysis.datasets[0].energy_range[-1].data[0][0])) == 100 |
|
43
|
|
|
|
chaimain /
asgardpy
