gammasimtools 0.11.0__py3-none-any.whl → 0.13.0__py3-none-any.whl

simtools/applications/convert_all_model_parameters_from_simtel.py

@@ -1,8 +1,8 @@
 #!/usr/bin/python3
 r"""
-    Convert all simulation model parameters exported from sim_telarray format using schema files.
+    Convert all simulation model parameters exported from sim_telarray format.
 
-    Check value, type, and range, convert units, and write json files
+    Check value, type, and range, convert units using schema files. Write json files
     ready to be submitted to the model database. Prints out parameters which are not found
     in simtel configuration file and parameters which are not found in simtools schema files.
 
@@ -19,22 +19,15 @@ r"""
     telescope (str, optional)
         Telescope model name (e.g. LST-1, SST-D, ...)
 
+    skip_parameter (str, optional)
+        List of parameters to be skipped (use sim_telarray names).
+
     Example
     -------
 
-    Extract model parameters with schema files from simtel configuration file
-    (requires access to the model parameter repository)
-
-    .. code-block:: console
-
-       simtools-convert-all-model-parameters-from-simtel \\
-          --schema_directory ../model_parameters/schema\\
-          --simtel_cfg_file all_telescope_config_la_palma.cfg\\
-          --simtel_telescope_name CT1\\
-          --telescope LSTN-01\\
-          --parameter_version "1.0.0"
-
-    The export of the model parameters from sim_telarray for 6.0.0 can be done e.g., as follows:
+    To export the model parameters from sim_telarray, first copy and unpack the configuration
+    tar ball from sim_telarray (usually called 'sim_telarray_config.tar.gz') to the sim_telarray
+    working directory. Extract the configuration using the following command:
 
     .. code-block:: console
 
@@ -42,18 +35,32 @@ r"""
             -C limits=no-internal -C initlist=no-internal -C list=no-internal \\
             -C typelist=no-internal -C maximum_telescopes=30 -DNSB_AUTOSCALE \\
             -DNECTARCAM -DHYPER_LAYOUT -DNUM_TELESCOPES=30 /dev/null \\
-            2>|/dev/null | grep '(@cfg)'  >| all_telescope_config_la_palma.cfg
+            2>|/dev/null | grep '(@cfg)' | sed 's/^(@cfg) //' >| all_telescope_config_la_palma.cfg
 
         ./sim_telarray/bin/sim_telarray -c sim_telarray/cfg/CTA/CTA-PROD6-Paranal.cfg \\
             -C limits=no-internal -C initlist=no-internal -C list=no-internal \\
             -C typelist=no-internal -C maximum_telescopes=87 -DNSB_AUTOSCALE \\
             -DFLASHCAM -DHYPER_LAYOUT -DNUM_TELESCOPES=87 /dev/null \\
-            2>|/dev/null | grep '(@cfg)'  >| all_telescope_config_paranal.cfg
+            2>|/dev/null | grep '(@cfg)' | sed 's/^(@cfg) //' >| all_telescope_config_paranal.cfg
+
+
+    Extract then model parameters from the simtel_array configuration file for LSTN-01
+    (telescopes are named CT1, CT2, ..., in the sim_telarray configuration file and must be
+    provided in the "simtel_telescope_name" command line argument)
+    and write json files in the same format as the model parameter database:
+
+    .. code-block:: console
+
+       simtools-convert-all-model-parameters-from-simtel \\
+          --simtel_cfg_file all_telescope_config_la_palma.cfg\\
+          --simtel_telescope_name CT1\\
+          --telescope LSTN-01\\
+          --parameter_version "1.0.0"\\
+          --output_path /path/to/output
 
 """
 
 import logging
-import re
 from pathlib import Path
 
 import numpy as np
@@ -63,7 +70,7 @@ import simtools.utils.general as gen
 from simtools.configuration import configurator
 from simtools.data_model import schema
 from simtools.io_operations.io_handler import IOHandler
-from simtools.simtel.simtel_config_reader import SimtelConfigReader
+from simtools.simtel import simtel_config_reader
 
 
 def _parse(label=None, description=None):
@@ -85,14 +92,6 @@ def _parse(label=None, description=None):
     """
     config = configurator.Configurator(label=label, description=description)
 
-    config.parser.add_argument(
-        "--schema_directory",
-        help=(
-            "Directory with schema files for model parameter validation "
-            "(default: simtools schema directory)"
-        ),
-        required=False,
-    )
     config.parser.add_argument(
         "--simtel_cfg_file",
         help="File name for simtel_array configuration",
@@ -105,36 +104,17 @@ def _parse(label=None, description=None):
         type=str,
         required=True,
     )
+    config.parser.add_argument(
+        "--skip_parameter",
+        help="List of parameters to be skipped.",
+        type=str,
+        nargs="*",
+        default=[],
+    )
     return config.initialize(simulation_model=["telescope", "parameter_version"])
 
 
-def get_list_of_simtel_parameters(simtel_config_file, logger):
-    """
-    Return list of simtel parameters found in simtel configuration file.
-
-    Parameters
-    ----------
-    simtel_config_file: str
-        File name for sim_telarray configuration
-    logger: logging.Logger
-        Logger object
-
-    Returns
-    -------
-    list
-        List of parameters found in simtel configuration file.
-
-    """
-    simtel_parameter_set = set()
-    with open(simtel_config_file, encoding="utf-8") as file:
-        for line in file:
-            parts_of_lines = re.split(r",\s*|\s+", line.strip())
-            simtel_parameter_set.add(parts_of_lines[1].lower())
-    logger.info(f"Found {len(simtel_parameter_set)} parameters in simtel configuration file.")
-    return list(simtel_parameter_set)
-
-
-def read_simtel_config_file(args_dict, logger, schema_file, camera_pixels=None):
+def read_simtel_config_file(args_dict, schema_file, camera_pixels=None):
     """
     Read the simtel configuration file.
 
@@ -142,25 +122,26 @@ def read_simtel_config_file(args_dict, logger, schema_file, camera_pixels=None):
     ----------
     args_dict: dict
         Dictionary with command line arguments.
-    logger: logging.Logger
-        Logger object
     schema_file: str
         Schema path name.
     camera_pixels: int
         Number of camera pixels.
 
+    Returns
+    -------
+    SimtelConfigReader, None
+        SimtelConfigReader object (None if parameter not found)
+
     """
-    simtel_config_reader = SimtelConfigReader(
+    config_reader = simtel_config_reader.SimtelConfigReader(
         schema_file=schema_file,
         simtel_config_file=args_dict["simtel_cfg_file"],
         simtel_telescope_name=args_dict["simtel_telescope_name"],
         camera_pixels=camera_pixels,
     )
-    logger.info(f"Simtel parameter: {simtel_config_reader.parameter_dict}")
-
-    if simtel_config_reader.parameter_dict is None or len(simtel_config_reader.parameter_dict) == 0:
+    if config_reader.parameter_dict is None or len(config_reader.parameter_dict) == 0:
         return None
-    return simtel_config_reader
+    return config_reader
 
 
 def get_number_of_camera_pixel(args_dict, logger):
@@ -184,12 +165,10 @@ def get_number_of_camera_pixel(args_dict, logger):
 
     """
     try:
-        simtel_config_reader = SimtelConfigReader(
-            schema_file=schema.get_model_parameter_schema_file("camera_pixels"),
-            simtel_config_file=args_dict["simtel_cfg_file"],
-            simtel_telescope_name=args_dict["simtel_telescope_name"],
+        config_reader = read_simtel_config_file(
+            args_dict, schema.get_model_parameter_schema_file("camera_pixels")
         )
-        _camera_pixel = simtel_config_reader.parameter_dict.get(args_dict["simtel_telescope_name"])
+        _camera_pixel = config_reader.parameter_dict.get(args_dict["simtel_telescope_name"])
     except (FileNotFoundError, AttributeError):
         logger.warning("Failed to read camera pixel parameter.")
         _camera_pixel = None
@@ -219,13 +198,14 @@ def read_and_export_parameters(args_dict, logger):
         List of simtools parameter not found in simtel configuration file.
 
     """
-    _parameters, _schema_files = schema.get_get_model_parameter_schema_files(
-        args_dict.get("schema_directory")
+    _parameters, _schema_files = schema.get_get_model_parameter_schema_files()
+    _simtel_parameters = simtel_config_reader.get_list_of_simtel_parameters(
+        args_dict["simtel_cfg_file"]
     )
-    _simtel_parameters = get_list_of_simtel_parameters(args_dict["simtel_cfg_file"], logger)
+    logger.info(f"Found {len(_simtel_parameters)} parameters in simtel configuration file.")
 
     io_handler = IOHandler()
-    io_handler.set_paths(output_path=args_dict["output_path"], use_plain_output_path=True)
+    io_handler.set_paths(output_path=args_dict["output_path"])
 
     _camera_pixel = get_number_of_camera_pixel(args_dict, logger)
 
@@ -233,27 +213,34 @@ def read_and_export_parameters(args_dict, logger):
 
     for _parameter, _schema_file in zip(_parameters, _schema_files):
         logger.info(f"Parameter: {_parameter} Schema file: {_schema_file}")
-        simtel_config_reader = read_simtel_config_file(
-            args_dict, logger, _schema_file, _camera_pixel
-        )
+        if _parameter in args_dict["skip_parameter"]:
+            logger.info(f"Skipping {_parameter}")
+            continue
+        config_reader = read_simtel_config_file(args_dict, _schema_file, _camera_pixel)
 
-        if simtel_config_reader is None:
+        if config_reader is None:
             logger.info("Parameter not found in sim_telarray configuration file.")
             _parameters_not_in_simtel.append(_parameter)
             continue
 
+        logger.info(f"Simtel parameter: {config_reader.parameter_dict}")
+
         _json_dict = writer.ModelDataWriter.dump_model_parameter(
             parameter_name=_parameter,
-            value=simtel_config_reader.parameter_dict.get(args_dict["simtel_telescope_name"]),
+            value=config_reader.parameter_dict.get(args_dict["simtel_telescope_name"]),
             instrument=args_dict["telescope"],
             parameter_version=args_dict["parameter_version"],
-            output_file=io_handler.get_output_file(f"{_parameter}.json"),
+            output_file=io_handler.get_output_file(
+                f"{_parameter}-{args_dict['parameter_version']}.json",
+                label=f"{args_dict['telescope']}",
+                sub_dir=f"{_parameter}",
+            ),
         )
 
-        simtel_config_reader.compare_simtel_config_with_schema()
+        config_reader.compare_simtel_config_with_schema()
 
-        if simtel_config_reader.simtel_parameter_name.lower() in _simtel_parameters:
-            _simtel_parameters.remove(simtel_config_reader.simtel_parameter_name.lower())
+        if config_reader.simtel_parameter_name.lower() in _simtel_parameters:
+            _simtel_parameters.remove(config_reader.simtel_parameter_name.lower())
 
         if _json_dict["file"]:
             logger.info(f"File name for {_parameter} is {_json_dict['value']}")
@@ -263,7 +250,7 @@ def read_and_export_parameters(args_dict, logger):
 
 def print_parameters_not_found(_parameters_not_in_simtel, _simtel_parameters, args_dict, logger):
     """
-    Print simtel/simtools parameter not found in schema and ocnfiguration files.
+    Print simtel/simtools parameter not found in schema and configuration files.
 
     For simtel parameters not found, check if the setting for the chose
     telescope is different from the default values.
@@ -289,29 +276,31 @@ def print_parameters_not_found(_parameters_not_in_simtel, _simtel_parameters, ar
     logger.info(f"Simtel parameters not found in schema files ({len(_simtel_parameters)}):")
     for para in sorted(_simtel_parameters):
         logger.info(f"Simtel parameter: {para}")
-        simtel_config_reader = SimtelConfigReader(
+        config_reader = simtel_config_reader.SimtelConfigReader(
             schema_file=None,
             simtel_config_file=args_dict["simtel_cfg_file"],
             simtel_telescope_name=args_dict["simtel_telescope_name"],
             parameter_name=para,
         )
-        _default = simtel_config_reader.parameter_dict.get("default")
-        _tel_value = simtel_config_reader.parameter_dict.get(args_dict["simtel_telescope_name"])
+        _default = config_reader.parameter_dict.get("default")
+        _tel_value = config_reader.parameter_dict.get(args_dict["simtel_telescope_name"])
         # simple comparison of default value and telescope values, does not work for lists
         try:
             if _default == _tel_value or np.isclose(_default, _tel_value):
-                logger.info(f"    Default and telescope values are equal: {_default}")
+                logger.info(f"    Default and telescope values for {para} are equal: {_default}")
                 continue
         except (ValueError, TypeError):
             pass
         if isinstance(_default, np.ndarray):
-            logger.warning(f"    Default value: {_default} (length: {_default.size})")
+            logger.warning(f"    Default value ({para}): {_default} (length: {_default.size})")
         else:
-            logger.warning(f"    Default value: {_default}")
+            logger.warning(f"    Default value ({para}): {_default}")
         if isinstance(_tel_value, np.ndarray):
-            logger.warning(f"    Telescope value: {_tel_value} (length: {_tel_value.size})")
+            logger.warning(
+                f"    Telescope value ({para}): {_tel_value} (length: {_tel_value.size})"
+            )
         else:
-            logger.warning(f"    Telescope value: {_tel_value}")
+            logger.warning(f"    Telescope value ({para}): {_tel_value}")
 
 
 def print_list_of_files(args_dict, logger):

simtools/applications/convert_geo_coordinates_of_array_elements.py

@@ -178,7 +178,7 @@ def main():
         )
         writer.ModelDataWriter.dump(
             args_dict=args_dict,
-            metadata=metadata.get_top_level_metadata() if metadata else None,
+            metadata=metadata,
             product_data=product_data,
             validate_schema_file=validate_schema_file,
         )

simtools/applications/db_get_parameter_from_db.py

@@ -3,22 +3,36 @@
 r"""
     Get a parameter entry from DB for a specific telescope or a site.
 
-    The application receives a parameter name, a site, a telescope (if applicable) and \
-    a version. It then prints out the parameter entry.
+    The application receives a parameter name, a site, a telescope (if applicable) and
+    a version. Allow to print the parameter entry to screen or save it to a file.
+    Parameter describing a table file can be written to disk or exported as an astropy table
+    (if available).
 
     Command line arguments
     ----------------------
     parameter (str, required)
         Parameter name
 
+    parameter_version (str, optional)
+        Parameter version
+
+    model_version (str, required)
+        Model version
+
     site (str, required)
         South or North.
 
     telescope (str, optional)
         Telescope model name (e.g. LST-1, SST-D, ...)
 
-    log_level (str, optional)
-        Log level to print.
+    output_file (str, optional)
+        Output file name. If not given, print to stdout.
+
+    export_model_file (bool, optional)
+        Export model file (if parameter describes a file).
+
+    export_model_file_as_table (bool, optional)
+        Export model file as astropy table (if parameter describes a file).
 
     Raises
     ------
@@ -35,12 +49,14 @@ r"""
                 --model_version 5.0.0
 
     Get the mirror_list parameter using the parameter_version from the DB.
+    Write the mirror list to disk.
 
     .. code-block:: console
 
         simtools-db-get-parameter-from-db --parameter mirror_list \\
                 --site North --telescope LSTN-01 \\
-                --parameter_version 1.0.0
+                --parameter_version 1.0.0 \\
+                --export_model_file
 
 """
 
@@ -53,7 +69,6 @@ import simtools.utils.general as gen
 from simtools.configuration import configurator
 from simtools.db import db_handler
 from simtools.io_operations import io_handler
-from simtools.utils import names
 
 
 def _parse():
@@ -72,6 +87,18 @@ def _parse():
         type=str,
         required=False,
     )
+    config.parser.add_argument(
+        "--export_model_file",
+        help="Export model file (if parameter describes a file)",
+        action="store_true",
+        required=False,
+    )
+    config.parser.add_argument(
+        "--export_model_file_as_table",
+        help="Export model file as astropy table (if parameter describes a file)",
+        action="store_true",
+        required=False,
+    )
 
     return config.initialize(
         db_config=True, simulation_model=["telescope", "parameter_version", "model_version"]
@@ -86,26 +113,29 @@ def main():  # noqa: D103
 
     db = db_handler.DatabaseHandler(mongo_db_config=db_config)
 
-    # get parameter using 'parameter_version'
-    if args_dict["parameter_version"] is not None:
-        pars = db.get_model_parameter(
+    pars = db.get_model_parameter(
+        parameter=args_dict["parameter"],
+        site=args_dict["site"],
+        array_element_name=args_dict.get("telescope"),
+        parameter_version=args_dict.get("parameter_version"),
+        model_version=args_dict.get("model_version"),
+    )
+    if args_dict["export_model_file"] or args_dict["export_model_file_as_table"]:
+        table = db.export_model_file(
             parameter=args_dict["parameter"],
-            parameter_version=args_dict["parameter_version"],
             site=args_dict["site"],
             array_element_name=args_dict["telescope"],
+            parameter_version=args_dict.get("parameter_version"),
+            model_version=args_dict.get("model_version"),
+            export_file_as_table=args_dict["export_model_file_as_table"],
         )
-    # get parameter using 'model_version'
-    elif args_dict["model_version"] is not None:
-        pars = db.get_model_parameters(
-            site=args_dict["site"],
-            array_element_name=args_dict.get("telescope"),
-            model_version=args_dict["model_version"],
-            collection=names.get_collection_name_from_parameter_name(args_dict["parameter"]),
-        )
-    else:
-        raise ValueError("Either 'parameter_version' or 'model_version' must be provided.")
-    if args_dict["parameter"] not in pars:
-        raise KeyError(f"The requested parameter, {args_dict['parameter']}, does not exist.")
+        param_value = pars[args_dict["parameter"]]["value"]
+        table_file = Path(io_handler.IOHandler().get_output_directory()) / f"{param_value}"
+        logger.info(f"Exported model file {param_value} to {table_file}")
+        if table and table_file.suffix != ".ecsv":
+            table.write(table_file.with_suffix(".ecsv"), format="ascii.ecsv", overwrite=True)
+            logger.info(f"Exported model file {param_value} to {table_file.with_suffix('.ecsv')}")
+
     if args_dict["output_file"] is not None:
         _output_file = (
             Path(io_handler.IOHandler().get_output_directory()) / args_dict["output_file"]

simtools/applications/derive_photon_electron_spectrum.py

@@ -81,7 +81,7 @@ def _parse(label):
         required=False,
     )
 
-    return config.initialize(db_config=False, output=True)
+    return config.initialize(db_config=False, output=True, simulation_model=["telescope"])
 
 
 def main():  # noqa: D103

simtools/applications/docs_produce_array_element_report.py

@@ -12,7 +12,6 @@ import logging
 
 from simtools.configuration import configurator
 from simtools.io_operations import io_handler
-from simtools.model.telescope_model import TelescopeModel
 from simtools.reporting.docs_read_parameters import ReadParameters
 from simtools.utils import general as gen
 
@@ -46,17 +45,9 @@ def main():  # noqa: D103
     logger = logging.getLogger()
     logger.setLevel(gen.get_log_level_from_user(args["log_level"]))
 
-    telescope_model = TelescopeModel(
-        site=args["site"],
-        telescope_name=args["telescope"],
-        model_version=args["model_version"],
-        label=label_name,
-        mongo_db_config=db_config,
-    )
-
     ReadParameters(
         db_config,
-        telescope_model,
+        args,
         output_path,
     ).produce_array_element_report()
 

simtools/applications/docs_produce_model_parameter_reports.py

@@ -1,7 +1,7 @@
 #!/usr/bin/python3
 
 r"""
-Produce one file per model parameter for a given array element.
+Produce a model parameter report per array element.
 
 The markdown reports include detailed information on each parameter,
 comparing their values over various model versions.
@@ -12,7 +12,6 @@ import logging
 
 from simtools.configuration import configurator
 from simtools.io_operations import io_handler
-from simtools.model.telescope_model import TelescopeModel
 from simtools.reporting.docs_read_parameters import ReadParameters
 from simtools.utils import general as gen
 
@@ -24,9 +23,7 @@ def _parse(label):
         description=("Produce a markdown report for model parameters."),
     )
 
-    return config.initialize(
-        db_config=True, simulation_model=["site", "telescope", "model_version"]
-    )
+    return config.initialize(db_config=True, simulation_model=["site", "telescope"])
 
 
 def main():  # noqa: D103
@@ -38,24 +35,14 @@ def main():  # noqa: D103
     logger = logging.getLogger()
     logger.setLevel(gen.get_log_level_from_user(args["log_level"]))
 
-    telescope_model = TelescopeModel(
-        site=args["site"],
-        telescope_name=args["telescope"],
-        model_version=args["model_version"],
-        label=label_name,
-        mongo_db_config=db_config,
-    )
-
     ReadParameters(
         db_config,
-        telescope_model,
+        args,
         output_path,
     ).produce_model_parameter_reports()
 
     logger.info(
-        f"Markdown report generated for {args['site']}"
-        f" Telescope {args['telescope']} (v{args['model_version']}):"
-        f" {output_path}"
+        f"Markdown report generated for {args['site']} Telescope {args['telescope']}: {output_path}"
     )
 
 

simtools/applications/plot_tabular_data.py

@@ -17,6 +17,9 @@ from pathlib import Path
 
 import simtools.utils.general as gen
 from simtools.configuration import configurator
+from simtools.constants import PLOT_CONFIG_SCHEMA
+from simtools.data_model import schema
+from simtools.data_model.metadata_collector import MetadataCollector
 from simtools.io_operations import io_handler
 from simtools.visualization import plot_tables
 
@@ -54,7 +57,7 @@ def _parse(label, description, usage):
         type=str,
         required=True,
     )
-    return config.initialize(db_config=True)
+    return config.initialize(db_config=True, simulation_model=["telescope"])
 
 
 def main():
@@ -70,7 +73,10 @@ def main():
     io_handler_instance = io_handler.IOHandler()
 
     plot_config = gen.convert_keys_in_dict_to_lowercase(
-        gen.collect_data_from_file(args_dict["plot_config"])
+        schema.validate_dict_using_schema(
+            gen.collect_data_from_file(args_dict["plot_config"]),
+            PLOT_CONFIG_SCHEMA,
+        )
     )
 
     plot_tables.plot(
@@ -79,6 +85,12 @@ def main():
         db_config=db_config_,
     )
 
+    MetadataCollector.dump(
+        args_dict,
+        io_handler_instance.get_output_file(args_dict["output_file"]),
+        add_activity_name=True,
+    )
+
 
 if __name__ == "__main__":
     main()

simtools/applications/{production_derive_limits.py → production_derive_corsika_limits.py}

@@ -11,6 +11,10 @@ event_data_file (str, required)
     Path to the file containing the event data.
 loss_fraction (float, required)
     Fraction of events to be lost.
+telescope_ids (list of int, optional)
+    List of telescope IDs to filter the events (default is None).
+    Definition of the telescope IDs can be found
+    in the telescope definition file used for simulations.
 
 
 Example
@@ -19,17 +23,17 @@ Derive limits for a given file with a specified loss fraction.
 
 .. code-block:: console
 
-    simtools-production-derive-limits\\
+    simtools-production-derive-corsika-limits\\
         --event_data_file path/to/event_data_file.hdf5 \\
-        --loss_fraction 1e-6
+        --loss_fraction 1e-6 \\
+        --telescope_ids 1 2 3
 """
 
 import logging
 
 import simtools.utils.general as gen
 from simtools.configuration import configurator
-from simtools.io_operations.hdf5_handler import read_hdf5
-from simtools.production_configuration.limits_calculation import LimitCalculator
+from simtools.production_configuration.derive_corsika_limits import LimitCalculator
 
 _logger = logging.getLogger(__name__)
 
@@ -44,6 +48,10 @@ def _parse():
         The event data file.
     loss_fraction: float
         Loss fraction of events for limit derivation.
+    telescope_ids: list of int, optional
+        List of telescope IDs to filter the events (default is None).
+        Definition of the telescope IDs can be found in the telescope
+        definition file used for simulations.
 
     Returns
     -------
@@ -63,6 +71,12 @@ def _parse():
     config.parser.add_argument(
         "--loss_fraction", type=float, required=True, help="Fraction of events to be lost."
     )
+    config.parser.add_argument(
+        "--telescope_ids",
+        type=int,
+        nargs="*",
+        help="List of telescope IDs to filter the events as used in sim_telarray.",
+    )
     return config.initialize(db_config=False)
 
 
@@ -75,11 +89,9 @@ def main():
 
     event_data_file_path = args_dict["event_data_file"]
     loss_fraction = args_dict["loss_fraction"]
+    telescope_ids = args_dict.get("telescope_ids")
 
-    _logger.info(f"Loading event data file: {event_data_file_path}")
-    tables = read_hdf5(event_data_file_path)
-
-    calculator = LimitCalculator(tables)
+    calculator = LimitCalculator(event_data_file_path, telescope_list=telescope_ids)
 
     lower_energy_limit = calculator.compute_lower_energy_limit(loss_fraction)
     _logger.info(f"Lower energy threshold: {lower_energy_limit}")