gammasimtools 0.6.0__py3-none-any.whl → 0.8.1__py3-none-any.whl

simtools/applications/generate_corsika_histograms.py

@@ -1,21 +1,21 @@
 #!/usr/bin/python3
 
-"""
-    Summary
-    -------
-    This application produces a set of histograms of the distribution of Cherenkov photons on the
-    ground (at observation level) read from the CORSIKA IACT output file provided as input.
+r"""
+    Generates a set of histograms Cherenkov photon distributions from CORSIKA output.
+
+    The Cherenkov photons (from observation level) are read from a CORSIKA IACT
+    output file provided as input.
 
     The histograms can be saved both into pdfs and in a hdf5 file.
 
-    The following 2D histograms are produced:
+    The following 2D histograms are generated:
         - Number of Cherenkov photons on the ground;
         - Density of Cherenkov photons on the ground;
         - Incoming direction (directive cosines) of the Cherenkov photons;
         - Time of arrival (ns) vs altitude of production (km);
         - Number of Cherenkov photons per event per telescope.
 
-    The following 1D histograms are produced:
+    The following 1D histograms are generated:
         - Wavelength;
         - Counts;
         - Density;
@@ -25,7 +25,7 @@
         - Number of photons per event.
 
     Histograms for the distribution of CORSIKA event header elements can also be generated by using
-    the `--event_1d_histograms` and `--event_2d_histograms` arguments. The accepted arguments (keys)
+    the --event_1d_histograms and --event_2d_histograms arguments. The accepted arguments (keys)
     are to be found in the CORSIKA manual, e.g., "total_energy", "zenith", "azimuth".
 
     Command line arguments
@@ -65,7 +65,7 @@
                 bins: 100
                 scale: linear
                 start: !astropy.units.Quantity
-                    unit: &id001 !astropy.units.Unit {unit: m}
+                    unit: &id001 m
                     value: -1000.0
                 stop: &id002 !astropy.units.Quantity
                     unit: *id001
@@ -81,7 +81,7 @@
                 bins: 80
                 scale: linear
                 start: !astropy.units.Quantity
-                    unit: &id003 !astropy.units.Unit {unit: nm}
+                    unit: nm
                     value: 200.0
                 stop: !astropy.units.Quantity
                     unit: *id003
@@ -91,7 +91,7 @@
                 bins: 100
                 scale: linear
                 start: !astropy.units.Quantity
-                    unit: &id004 !astropy.units.Unit {unit: ns}
+                    unit: ns
                     value: -2000.0
                 stop: !astropy.units.Quantity
                     unit: *id004
@@ -100,7 +100,7 @@
                 bins: 100
                 scale: linear
                 start: !astropy.units.Quantity
-                    unit: &id005 !astropy.units.Unit {unit: km}
+                    unit: km
                     value: 120.0
                 stop: !astropy.units.Quantity
                     unit: *id005
@@ -110,10 +110,10 @@
     pdf (bool, optional)
         If set, histograms are saved into pdf files.
         One pdf file contains all the histograms for the Cherenkov photons.
-        The name of the file is controlled via `hdf5_file_name`.
+        The name of the file is controlled via hdf5_file_name.
         If event_1d_histograms and event_2d_histograms are used, two separate pdf files might be
         created to accommodate the histograms for the CORSIKA event header elements. The core names
-        of these output pdf files are also given by `hdf5_file_name` argument with the addition of
+        of these output pdf files are also given by hdf5_file_name argument with the addition of
         'event_1d_histograms' and 'event_2d_histograms'.
 
 
@@ -122,28 +122,28 @@
 
     hdf5_file_name (str, optional)
         The name of the output hdf5 data (without the path).
-        It requires the `--hdf5` flag.
-        If not given, `hdf5_file_name` takes the name from the input IACT file (`input_file`).
-        If the output `hdf5_file_name` file already exists, the tables associated to the chosen
-        flags (e.g. `hdf5`, `event_1d_histograms`, `event_2d_histograms`) will be overwritten. The
+        It requires the --hdf5 flag.
+        If not given, hdf5_file_name takes the name from the input IACT file (input_file).
+        If the output hdf5_file_name file already exists, the tables associated to the chosen
+        flags (e.g. hdf5, event_1d_histograms, event_2d_histograms) will be overwritten. The
         remaining tables, if any, will stay untouched.
 
 
     event_1d_histograms (str, optional)
-        Produce 1D histograms for elements given in `--event_1d_histograms` from the CORSIKA event
+        Generate 1D histograms for elements given in --event_1d_histograms from the CORSIKA event
         header and save into hdf5/pdf files.
         It allows more than one argument, separated by simple spaces.
-        Usage: `--event_1d_histograms first_interaction_height total_energy`.
+        Usage: --event_1d_histograms first_interaction_height total_energy.
 
     event_2d_histograms (str, optional)
-        Produce 2D histograms for elements given in `--event_2d_histograms` from the CORSIKA event
+        Generate 2D histograms for elements given in --event_2d_histograms from the CORSIKA event
         header and save into hdf5/pdf files.
         It allows more than one argument, separated by simple spaces.
-        The elements are grouped into pairs and the 2D histograms are produced always for two
+        The elements are grouped into pairs and the 2D histograms are generated always for two
         subsequent elements.
-        For example, `--event_2d_histograms first_interaction_height total_energy zenith azimuth`
-        will produce one 2D histogram for `first_interaction_height` `total_energy` and another 2D
-        histogram for `zenith` and `azimuth`.
+        For example, --event_2d_histograms first_interaction_height total_energy zenith azimuth
+        will generate one 2D histogram for first_interaction_height total_energy and another 2D
+        histogram for zenith and azimuth.
 
     Example
     -------
@@ -151,11 +151,11 @@
 
      .. code-block:: console
 
-        simtools-generate-corsika-histograms --iact_file /workdir/external/simtools/tests/\
-            resources/tel_output_10GeV-2-gamma-20deg-CTAO-South.corsikaio --pdf --hdf5
+        simtools-generate-corsika-histograms --iact_file /workdir/external/simtools/\\
+        testsresources/tel_output_10GeV-2-gamma-20deg-CTAO-South.corsikaio \\
+            --pdf --hdf5 \\
             --event_2d_histograms zenith azimuth --event_1d_histograms total_energy
 
-
     Expected final print-out message:
 
     .. code-block:: console
@@ -182,7 +182,7 @@ logger = logging.getLogger()
 
 def _parse(label, description):
     """
-    Parse command line configuration
+    Parse command line configuration.
 
     Parameters
     ----------
@@ -251,7 +251,7 @@ def _parse(label, description):
     config.parser.add_argument(
         "--event_1d_histograms",
         help="The keys from the CORSIKA event header to be used for the generation of 1D "
-        "histograms. The available choices can been found in the `all_event_keys` attribute of"
+        "histograms. The available choices can been found in the all_event_keys attribute of"
         "the CorsikaHistograms.",
         required=False,
         default=None,
@@ -261,7 +261,7 @@ def _parse(label, description):
     config.parser.add_argument(
         "--event_2d_histograms",
         help="The keys from the CORSIKA event header to be used for the generation of 2D "
-        "histograms. The available choices can been found in the `all_event_keys` attribute of"
+        "histograms. The available choices can been found in the all_event_keys attribute of"
         "the CorsikaHistograms.",
         required=False,
         default=None,
@@ -270,16 +270,16 @@ def _parse(label, description):
 
     config_parser, _ = config.initialize(db_config=False, paths=True)
 
-    if not config_parser["pdf"]:
-        if (
-            not config_parser["hdf5"]
-            and not config_parser["event_1d_histograms"]
-            and not config_parser["event_2d_histograms"]
-        ):
-            config.parser.error(
-                "At least one argument is required: `--pdf`, `--hdf5`, `--event_1d_histograms`, or "
-                "`--event_2d_histograms`."
-            )
+    if (
+        not config_parser["pdf"]
+        and not config_parser["hdf5"]
+        and not config_parser["event_1d_histograms"]
+        and not config_parser["event_2d_histograms"]
+    ):
+        config.parser.error(
+            "At least one argument is required: --pdf, --hdf5, --event_1d_histograms, or "
+            "--event_2d_histograms."
+        )
 
     return config_parser, _
 
@@ -290,12 +290,11 @@ def _plot_figures(corsika_histograms_instance, test=False):
 
     Parameters
     ----------
-    corsika_histograms_instance: `CorsikaHistograms` instance.
+    corsika_histograms_instance: CorsikaHistograms instance.
         The CorsikaHistograms instance created in main.
     test: bool
         If true plots the figures for the first two functions only.
     """
-
     plot_function_names = [
         plotting_method
         for plotting_method in dir(corsika_histograms_visualize)
@@ -307,8 +306,8 @@ def _plot_figures(corsika_histograms_instance, test=False):
 
     figure_list = []
     for function_name in plot_function_names:
-        function = getattr(corsika_histograms_visualize, function_name)
-        figures = function(corsika_histograms_instance)
+        plot_function = getattr(corsika_histograms_visualize, function_name)
+        figures = plot_function(corsika_histograms_instance)
         for fig in figures:
             figure_list.append(fig)
 
@@ -326,10 +325,10 @@ def _derive_event_1d_histograms(
 
     Parameters
     ----------
-    corsika_histograms_instance: `CorsikaHistograms` instance.
+    corsika_histograms_instance: CorsikaHistograms instance.
         The CorsikaHistograms instance created in main.
     event_1d_header_keys: str
-        Produce 1D histograms for elements given in `event_1d_header_keys` from the CORSIKA event
+        Generate 1D histograms for elements given in event_1d_header_keys from the CORSIKA event
         header and save into hdf5/pdf files.
     pdf: bool
         If true, histograms are saved into a pdf file.
@@ -362,14 +361,15 @@ def _derive_event_2d_histograms(
 ):
     """
     Auxiliary function to derive the histograms for the arguments given by event_1d_histograms.
+
     If an odd number of event header keys are given, the last one is discarded.
 
     Parameters
     ----------
-    corsika_histograms_instance: `CorsikaHistograms` instance.
+    corsika_histograms_instance: CorsikaHistograms instance.
         The CorsikaHistograms instance created in main.
     event_2d_header_keys: str
-        Produce 1D histograms for elements given in `event_1d_header_keys` from the CORSIKA event
+        Generate 1D histograms for elements given in event_1d_header_keys from the CORSIKA event
         header and save into hdf5/pdf files.
     pdf: bool
         If true, histograms are saved into a pdf file.
@@ -384,7 +384,7 @@ def _derive_event_2d_histograms(
 
         if len(event_2d_header_keys) % 2 == 1:  # if odd number of keys
             msg = (
-                "An odd number of keys was passed to produce 2D histograms."
+                "An odd number of keys was passed to generate 2D histograms."
                 "The last key is being ignored."
             )
             logger.warning(msg)
@@ -412,7 +412,7 @@ def _derive_event_2d_histograms(
         corsika_histograms_visualize.save_figs_to_pdf(figures_list, output_file_name)
 
 
-def main():
+def main():  # noqa: D103
     label = Path(__file__).stem
     description = "Generate histograms for the Cherenkov photons saved in the CORSIKA IACT file."
     io_handler_instance = io_handler.IOHandler()
@@ -433,7 +433,7 @@ def main():
         except ValueError:
             msg = (
                 f"{args_dict['telescope_indices']} not a valid input. "
-                f"Please use integer numbers for `telescope_indices`"
+                f"Please use integer numbers for telescope_indices"
             )
             logger.error(msg)
             raise

simtools/applications/generate_default_metadata.py

@@ -1,7 +1,5 @@
 #!/usr/bin/python3
-"""
-    Summary
-    -------
+r"""
     Generate a default simtools metadata file from a json schema.
 
     Command line arguments
@@ -16,8 +14,8 @@
     -------
     .. code-block:: console
 
-        simtools-generate-default-metadata
-            --schema simtools/schemas/metadata.metaschema.yml
+        simtools-generate-default-metadata \\
+            --schema simtools/schemas/metadata.metaschema.yml \\
             --output_file default_metadata.yml
 
 
@@ -37,7 +35,7 @@ from simtools.io_operations import io_handler
 
 def _parse(label, description):
     """
-    Parse command line configuration
+    Parse command line configuration.
 
     Parameters
     ----------
@@ -52,7 +50,6 @@ def _parse(label, description):
         Command line parser object
 
     """
-
     config = configurator.Configurator(label=label, description=description)
 
     config.parser.add_argument(
@@ -71,7 +68,7 @@ def _parse(label, description):
     return config.initialize(output=False, require_command_line=True)
 
 
-def main():
+def main():  # noqa: D103
     label = Path(__file__).stem
     args_dict, _ = _parse(
         label, description="Generate a default simtools metadata file from a json schema."

simtools/applications/generate_regular_arrays.py

@@ -0,0 +1,117 @@
+#!/usr/bin/python3
+
+"""
+    Make a regular array of telescopes and save it as astropy table.
+
+    The arrays consist of one telescope at the center of the array and or of 4 telescopes
+    in a square grid. These arrays are used for trigger rate simulations.
+
+    The array layout files created will be available at the data/layout directory.
+
+    Command line arguments
+    ----------------------
+    site (str, required)
+        observatory site (e.g., North or South).
+    model_version (str, optional)
+        Model version to use (e.g., 6.0.0). If not provided, the latest version is used.
+
+    Example
+    -------
+    Runtime < 10 s.
+
+    .. code-block:: console
+
+        simtools-generate-regular-arrays --site=North
+"""
+
+import logging
+import os
+from pathlib import Path
+
+import astropy.units as u
+from astropy.table import QTable
+
+import simtools.data_model.model_data_writer as writer
+import simtools.utils.general as gen
+from simtools.configuration import configurator
+from simtools.utils import names
+
+# Telescope distances for 4 tel square arrays
+# !HARDCODED
+telescope_distance = {"LST": 57.5 * u.m, "MST": 70 * u.m, "SST": 80 * u.m}
+
+
+def _parse():
+    config = configurator.Configurator(
+        label=Path(__file__).stem,
+        description=(
+            "Generate a regular array of telescope and save as astropy table.\n"
+            "Default telescope distances for 4 telescope square arrays are: \n"
+            f"  LST: {telescope_distance['LST']}\n"
+            f"  MST: {telescope_distance['MST']}\n"
+            f"  SST: {telescope_distance['SST']}\n"
+        ),
+    )
+    return config.initialize(db_config=False, simulation_model="site", output=True)
+
+
+def main():
+    """Create layout array files (ecsv) of regular arrays."""
+    args_dict, _ = _parse()
+
+    logger = logging.getLogger()
+    logger.setLevel(gen.get_log_level_from_user(args_dict["log_level"]))
+
+    if args_dict["site"] == "South":
+        array_list = ["1SST", "4SST", "1MST", "4MST", "1LST", "4LST"]
+    else:
+        array_list = ["1MST", "4MST", "1LST", "4LST"]
+
+    for array_name in array_list:
+        logger.info(f"Processing array {array_name}")
+
+        tel_name, pos_x, pos_y, pos_z = [], [], [], []
+        tel_size = array_name[1:4]
+
+        # Single telescope at the center
+        if array_name[0] == "1":
+            tel_name.append(
+                names.get_array_element_name_from_type_site_id(tel_size, args_dict["site"], "01")
+            )
+            pos_x.append(0 * u.m)
+            pos_y.append(0 * u.m)
+            pos_z.append(0 * u.m)
+        # 4 telescopes in a regular square grid
+        else:
+            for i in range(1, 5):
+                tel_name.append(
+                    names.get_array_element_name_from_type_site_id(
+                        tel_size, args_dict["site"], f"0{i}"
+                    )
+                )
+                pos_x.append(telescope_distance[tel_size] * (-1) ** (i // 2))
+                pos_y.append(telescope_distance[tel_size] * (-1) ** (i % 2))
+                pos_z.append(0 * u.m)
+
+        table = QTable(meta={"array_name": array_name, "site": args_dict["site"]})
+        table["telescope_name"] = tel_name
+        table["position_x"] = pos_x
+        table["position_y"] = pos_y
+        table["position_z"] = pos_z
+        table.sort("telescope_name")
+        table.pprint()
+
+        output_file = args_dict.get("output_file", None)
+        if output_file is not None:
+            base_name, file_extension = os.path.splitext(output_file)
+            output_file = f"{base_name}-{args_dict['site']}-{array_name}{file_extension}"
+        writer.ModelDataWriter.dump(
+            args_dict=args_dict,
+            output_file=output_file,
+            metadata=None,
+            product_data=table,
+        )
+
+
+if __name__ == "__main__":
+    main()

simtools/applications/generate_simtel_array_histograms.py

@@ -1,9 +1,8 @@
 #!/usr/bin/python3
 
-"""
-    Summary
-    -------
-    This application allows to write sim_telarray histograms into pdf and hdf5 files.
+r"""
+    Write sim_telarray histograms into pdf and hdf5 files.
+
     It accepts multiple lists of histograms files, a single list or a histogram file.
     Each histogram is plotted in a page of the pdf file if the --pdf option is activated.
 
@@ -19,31 +18,29 @@
     pdf (bool, optional)
         If set, histograms are saved into pdf files.
         One pdf file contains all the histograms found in the file.
-        The name of the file is controlled via `output_file_name`.
+        The name of the file is controlled via output_file_name.
     hdf5: bool
         If true, histograms are saved into hdf5 files.
-        At least one of `pdf` and `hdf5` has to be activated.
+        At least one of pdf and hdf5 has to be activated.
     output_file_name (str, optional)
         The name of the output hdf5 (and/or pdf) files (without the path).
-        If not given, `output_file_name` takes the name from the (first) input file
-        (`hist_file_names`).
-        If the output `output_file_name.hdf5` file already exists and `hdf5` is set, the tables
-        associated to `hdf5` will be overwritten. The remaining tables, if any, will stay
+        If not given, output_file_name takes the name from the (first) input file
+        (hist_file_names).
+        If the output output_file_name.hdf5 file already exists and hdf5 is set, the tables
+        associated to hdf5 will be overwritten. The remaining tables, if any, will stay
         untouched.
-    verbosity (str, optional)
-        Log level to print.
 
     Raises
     ------
     TypeError:
-        if argument passed through `hist_file_names` is not a file.
+        if argument passed through hist_file_names is not a file.
 
     Example
     -------
     .. code-block:: console
 
-        simtools-generate-simtel-array-histograms --hist_file_names tests/resources/
-            run2_gamma_za20deg_azm0deg-North-Prod5_test-production-5.hdata.zst
+        simtools-generate-simtel-array-histograms --hist_file_names tests/resources/ \\
+            run2_gamma_za20deg_azm0deg-North-Prod5_test-production-5.hdata.zst \\
             --output_file_name test_hist_hdata --hdf5 --pdf
 
 """
@@ -57,12 +54,12 @@ from matplotlib.backends.backend_pdf import PdfPages
 import simtools.utils.general as gen
 from simtools.configuration import configurator
 from simtools.io_operations import io_handler
-from simtools.simtel.simtel_histograms import SimtelHistograms
+from simtools.simtel.simtel_io_histograms import SimtelIOHistograms
 
 
 def _parse(label, description):
     """
-    Parse command line configuration
+    Parse command line configuration.
 
     Parameters
     ----------
@@ -105,88 +102,132 @@ def _parse(label, description):
     )
 
     config_parser, _ = config.initialize(db_config=False, paths=True)
-    if not config_parser["pdf"]:
-        if not config_parser["hdf5"]:
-            config.parser.error("At least one argument is required: `--pdf` or `--hdf5`.")
+    if not config_parser["pdf"] and not config_parser["hdf5"]:
+        config.parser.error("At least one argument is required: --pdf or --hdf5.")
 
     return config_parser
 
 
-def main():
-    label = Path(__file__).stem
-    description = "Display the simtel_array histograms."
-    io_handler_instance = io_handler.IOHandler()
-    config_parser = _parse(label, description)
-    output_path = io_handler_instance.get_output_directory(label, sub_dir="application-plots")
-    logger = logging.getLogger()
-    logger.setLevel(gen.get_log_level_from_user(config_parser["log_level"]))
-    logger.info("Starting the application.")
+def build_histogram_files(config_parser, logger):
+    """
+    Build a list of histogram files from command line arguments.
 
-    # Building list of histograms from the input files
+    Parameters
+    ----------
+    config_parser: dict
+        Parsed command line arguments.
+    logger: logging.Logger
+        Logger object for logging messages.
+
+    Returns
+    -------
+    list
+        List of histogram file paths.
+    """
     histogram_files = []
     for one_file in config_parser["hist_file_names"]:
         try:
             if Path(one_file).suffix in [".zst", ".simtel", ".hdata"]:
                 histogram_files.append(one_file)
             else:
-                # Collecting hist files
                 with open(one_file, encoding="utf-8") as file:
                     for line in file:
-                        # Removing '\n' from filename, in case it is left there.
-                        histogram_files.append(line.replace("\n", ""))
+                        histogram_files.append(line.strip())
         except FileNotFoundError as exc:
             msg = f"{one_file} is not a file."
             logger.error(msg)
             raise FileNotFoundError from exc
+    return histogram_files
 
-    # If no output name is passed, the tool gets the name of the first histogram of the list
-    if config_parser["output_file_name"] is None:
-        config_parser["output_file_name"] = Path(histogram_files[0]).absolute().name
-    output_file_name = Path(output_path).joinpath(f"{config_parser['output_file_name']}")
 
-    # If the hdf5 output file already exists, it is overwritten
-    if (Path(f"{config_parser['output_file_name']}.hdf5").exists()) and (config_parser["hdf5"]):
+def check_and_log_overwrite(config_parser, logger):
+    """
+    Check if the output hdf5 file already exists and log a warning if it does.
+
+    Parameters
+    ----------
+    config_parser: dict
+        Parsed command line arguments.
+    logger: logging.Logger
+        Logger object for logging messages.
+
+    Returns
+    -------
+    bool
+        True if the hdf5 file exists and should be overwritten, False otherwise.
+    """
+    if Path(f"{config_parser['output_file_name']}.hdf5").exists() and config_parser["hdf5"]:
         msg = (
             f"Output hdf5 file {config_parser['output_file_name']}.hdf5 already exists. "
             f"Overwriting it."
         )
         logger.warning(msg)
-        overwrite = True
-    else:
-        overwrite = False
+        return True
+    return False
 
-    # Building SimtelHistograms
-    simtel_histograms = SimtelHistograms(histogram_files)
 
+def create_pdf(simtel_histograms, output_file_name, config_parser, logger):
+    """
+    Create a PDF file containing histograms.
+
+    Parameters
+    ----------
+    simtel_histograms: SimtelIOHistograms
+        SimtelIOHistograms object containing histograms to plot.
+    output_file_name: str
+        Base name for the output PDF file.
+    config_parser: dict
+        Parsed command line arguments.
+    logger: logging.Logger
+        Logger object for logging messages.
+    """
     if config_parser["pdf"]:
         logger.debug(f"Creating the pdf file {output_file_name}.pdf")
         pdf_pages = PdfPages(f"{output_file_name}.pdf")
-
-        if config_parser["test"]:
-            number_of_histograms = 2
-        else:
-            number_of_histograms = simtel_histograms.number_of_histograms
-
+        number_of_histograms = 2 if config_parser["test"] else len(simtel_histograms.combined_hists)
         for i_hist in range(number_of_histograms):
-            title = simtel_histograms.get_histogram_title(i_hist)
-
-            logger.debug(f"Processing: {title}")
-
+            logger.debug(f"Processing: {i_hist + 1} histogram.")
             fig, ax = plt.subplots(1, 1, figsize=(6, 6))
             simtel_histograms.plot_one_histogram(i_hist, ax)
-
             plt.tight_layout()
             pdf_pages.savefig(fig)
             plt.clf()
-
         plt.close()
         pdf_pages.close()
         logger.info(f"Wrote histograms to the pdf file {output_file_name}.pdf")
 
+
+def export_to_hdf5(simtel_histograms, output_file_name, overwrite, config_parser, logger):
+    """Export histograms to an HDF5 file."""
     if config_parser["hdf5"]:
         logger.info(f"Wrote histograms to the hdf5 file {output_file_name}.hdf5")
         simtel_histograms.export_histograms(f"{output_file_name}.hdf5", overwrite=overwrite)
 
 
+def main():  # noqa: D103
+    label = Path(__file__).stem
+    description = "Display the simtel_array histograms."
+    io_handler_instance = io_handler.IOHandler()
+    config_parser = _parse(label, description)
+    output_path = io_handler_instance.get_output_directory(label, sub_dir="application-plots")
+    logger = logging.getLogger()
+    logger.setLevel(gen.get_log_level_from_user(config_parser["log_level"]))
+    logger.info("Starting the application.")
+
+    histogram_files = build_histogram_files(config_parser, logger)
+
+    # If no output name is passed, the tool gets the name of the first histogram of the list
+    if config_parser["output_file_name"] is None:
+        config_parser["output_file_name"] = Path(histogram_files[0]).absolute().name
+    output_file_name = Path(output_path).joinpath(f"{config_parser['output_file_name']}")
+
+    # If the hdf5 output file already exists, it is overwritten
+    overwrite = check_and_log_overwrite(config_parser, logger)
+
+    simtel_histograms = SimtelIOHistograms(histogram_files)
+    create_pdf(simtel_histograms, output_file_name, config_parser, logger)
+    export_to_hdf5(simtel_histograms, output_file_name, overwrite, config_parser, logger)
+
+
 if __name__ == "__main__":
     main()