gammasimtools 0.20.0__tar.gz → 0.21.0__tar.gz

{gammasimtools-0.20.0 → gammasimtools-0.21.0}/.github/workflows/CI-integrationtests.yml

@@ -57,7 +57,7 @@ jobs:
     strategy:
       fail-fast: false
       matrix:
-        model_version: ['6.0.0', '5.0.0', '6.0.0,6.1.0']
+        model_version: ['6.0', '5.0', '6.0,6.1']
 
     defaults:
       run:

{gammasimtools-0.20.0 → gammasimtools-0.21.0}/CHANGELOG.md

@@ -7,6 +7,21 @@ This changelog is generated using [Towncrier](https://towncrier.readthedocs.io/)
 
 <!-- towncrier release notes start -->
 
+## [v0.21.0](https://github.com/gammasim/simtools/releases/tag/v0.21.0) - 2025-09-23
+
+### Bugfixes
+
+- Fix output file error when running application `simtools-print-version`. ([#1774](https://github.com/gammasim/simtools/pull/1774))
+- Bugfix in configuring hard-stereo array triggers (`hard_stereo` should be `hardstereo` in `array_trigger.dat`). ([#1780](https://github.com/gammasim/simtools/pull/1780))
+
+### New Features
+
+- Implemented gradient descent minimisation for psf parameter optimisation, KS statistic as a sanity check, and monte carlo error on the simulated psf. ([#1717](https://github.com/gammasim/simtools/pull/1717))
+- Add functionality to work with delta-type production tables of the simulation model for patch updates. ([#1758](https://github.com/gammasim/simtools/pull/1758))
+- Automatic retrieval of the latest patch version for the simulation model.
+  Allow `--model_version=6.0` being resolved to the newest version. ([#1768](https://github.com/gammasim/simtools/pull/1768))
+
+
 ## [v0.20.0](https://github.com/gammasim/simtools/releases/tag/v0.20.0) - 2025-09-12
 
 ### Bugfixes

{gammasimtools-0.20.0/src/gammasimtools.egg-info → gammasimtools-0.21.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: gammasimtools
-Version: 0.20.0
+Version: 0.21.0
 Summary: Tools for the Simulation System of the CTA Observatory
 Author-email: simtools developers <simtools-developer@desy.de>
 License-Expression: BSD-3-Clause

{gammasimtools-0.20.0 → gammasimtools-0.21.0}/docs/_static/switcher.json

@@ -2,5 +2,7 @@
   {"version": "latest", "url": "https://gammasim.github.io/simtools/", "preferred": true},
   {"version": "v0.16.0", "url": "https://gammasim.github.io/simtools/v0.16.0/"},
   {"version": "v0.17.0", "url": "https://gammasim.github.io/simtools/v0.17.0/"},
-  {"version": "v0.18.0", "url": "https://gammasim.github.io/simtools/v0.18.0/"}
+  {"version": "v0.18.0", "url": "https://gammasim.github.io/simtools/v0.18.0/"},
+  {"version": "v0.20.0", "url": "https://gammasim.github.io/simtools/v0.20.0/"},
+  {"version": "v0.21.0", "url": "https://gammasim.github.io/simtools/v0.21.0/"}
 ]

{gammasimtools-0.20.0 → gammasimtools-0.21.0}/docs/source/api-reference/corsika.md

@@ -29,12 +29,3 @@ Support modules for running CORSIKA.
    :members:
 
 ```
-
-## corsika_histograms_visualize
-
-(corsika-histograms-visualize-1)=
-
-```{eval-rst}
-.. automodule:: corsika.corsika_histograms_visualize
-   :members:
-```

{gammasimtools-0.20.0 → gammasimtools-0.21.0}/docs/source/api-reference/visualization.md

@@ -41,6 +41,15 @@ the visualization module.
    :members:
 ```
 
+## plot_corsika_histograms
+
+(corsika-histograms-visualize-1)=
+
+```{eval-rst}
+.. automodule:: visualization.plot_corsika_histograms
+   :members:
+```
+
 (plot-tables)=
 
 ## plot_tables
@@ -75,11 +84,19 @@ the visualization module.
    :members:
 ```
 
+## plot_psf
+
+```{eval-rst}
+.. automodule:: visualization.plot_psf
+   :members:
+```
+
 (plot-incident-angles)=
 
 ## plot_incident_angles
 
 ```{eval-rst}
 .. automodule:: visualization.plot_incident_angles
+
    :members:
 ```

gammasimtools-0.21.0/docs/source/user-guide/applications/simtools-maintain-simulation-model-add-production.rst

@@ -0,0 +1,6 @@
+
+simtools-maintain-simulation-model-add-production
+=================================================
+
+.. automodule:: simtools.applications.maintain_simulation_model_add_production
+   :members:

{gammasimtools-0.20.0 → gammasimtools-0.21.0}/docs/source/user-guide/applications.md

@@ -75,7 +75,7 @@ simtools-generate-corsika-histograms <applications/simtools-generate-corsika-his
 simtools-generate-default-metadata <applications/simtools-generate-default-metadata>
 simtools-generate-regular-arrays <applications/simtools-generate-regular-arrays>
 simtools-generate-simtel-event-data <applications/simtools-generate-simtel-event-data>
-simtools-maintain-simulation-model-add-production-table <applications/simtools-maintain-simulation-model-add-production-table>
+simtools-maintain-simulation-model-add-production <applications/simtools-maintain-simulation-model-add-production>
 simtools-maintain-simulation-model-compare-productions <applications/simtools-maintain-simulation-model-compare-productions>
 simtools-maintain-simulation-model-verify-production-tables <applications/simtools-maintain-simulation-model-verify-production-tables>
 simtools-merge-tables <applications/simtools-merge-tables>

{gammasimtools-0.20.0 → gammasimtools-0.21.0}/docs/source/user-guide/manage_simulation_models/simulation_models_database.md

@@ -34,5 +34,5 @@ Applications modifying database entries:
 Applications to manage the simulation model repository:
 
 * verify consistency and completeness of production tables and model parameters: [simtools-maintain-simulation-model-verify-production-tables](maintain_simulation_model_verify_production_tables)
-* generate a new simulation model production tables by copying existing table and apply modifications: [simtools-maintain-simulation-model-add-production-table](maintain_simulation_model_add_production_table)
+* generate a new simulation model production by copying existing table and apply modifications: [simtools-maintain-simulation-model-add-production](maintain_simulation_model_add_production)
 * compare two simulation model production directories and report differences: [simtools-maintain-simulation-model-compare-productions](maintain_simulation_model_compare_productions)

{gammasimtools-0.20.0 → gammasimtools-0.21.0}/pyproject.toml

@@ -98,7 +98,7 @@ scripts.simtools-generate-corsika-histograms = "simtools.applications.generate_c
 scripts.simtools-generate-default-metadata = "simtools.applications.generate_default_metadata:main"
 scripts.simtools-generate-regular-arrays = "simtools.applications.generate_regular_arrays:main"
 scripts.simtools-generate-simtel-event-data = "simtools.applications.generate_simtel_event_data:main"
-scripts.simtools-maintain-simulation-model-add-production-table = "simtools.applications.maintain_simulation_model_add_production_table:main"
+scripts.simtools-maintain-simulation-model-add-production = "simtools.applications.maintain_simulation_model_add_production:main"
 scripts.simtools-maintain-simulation-model-compare-productions = "simtools.applications.maintain_simulation_model_compare_productions:main"
 scripts.simtools-maintain-simulation-model-verify-production-tables = "simtools.applications.maintain_simulation_model_verify_production_tables:main"
 scripts.simtools-merge-tables = "simtools.applications.merge_tables:main"

{gammasimtools-0.20.0 → gammasimtools-0.21.0/src/gammasimtools.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: gammasimtools
-Version: 0.20.0
+Version: 0.21.0
 Summary: Tools for the Simulation System of the CTA Observatory
 Author-email: simtools developers <simtools-developer@desy.de>
 License-Expression: BSD-3-Clause

{gammasimtools-0.20.0 → gammasimtools-0.21.0}/src/gammasimtools.egg-info/SOURCES.txt

@@ -108,7 +108,7 @@ docs/source/user-guide/applications/simtools-generate-corsika-histograms.rst
 docs/source/user-guide/applications/simtools-generate-default-metadata.rst
 docs/source/user-guide/applications/simtools-generate-regular-arrays.rst
 docs/source/user-guide/applications/simtools-generate-simtel-event-data.rst
-docs/source/user-guide/applications/simtools-maintain-simulation-model-add-production-table.rst
+docs/source/user-guide/applications/simtools-maintain-simulation-model-add-production.rst
 docs/source/user-guide/applications/simtools-maintain-simulation-model-compare-productions.rst
 docs/source/user-guide/applications/simtools-maintain-simulation-model-verify-production-tables.rst
 docs/source/user-guide/applications/simtools-merge-tables.rst
@@ -138,6 +138,7 @@ docs/source/user-guide/applications/simtools-validate-optics.rst
 docs/source/user-guide/applications/images/derive_mirror_rnda_North-MST-FlashCam-D.png
 docs/source/user-guide/applications/images/derive_mirror_rnda_North-MST-FlashCam-D_D80-distributions.png
 docs/source/user-guide/applications/images/derive_psf_parameters.png
+docs/source/user-guide/applications/images/gradient_descent.png
 docs/source/user-guide/applications/images/incident_angles_multi_calculate_incident_angles_SSTS-04.png
 docs/source/user-guide/applications/images/incident_angles_primary_multi_calculate_incident_angles_SSTS-04.png
 docs/source/user-guide/applications/images/incident_angles_secondary_multi_calculate_incident_angles_SSTS-04.png
@@ -193,7 +194,7 @@ src/simtools/applications/generate_corsika_histograms.py
 src/simtools/applications/generate_default_metadata.py
 src/simtools/applications/generate_regular_arrays.py
 src/simtools/applications/generate_simtel_event_data.py
-src/simtools/applications/maintain_simulation_model_add_production_table.py
+src/simtools/applications/maintain_simulation_model_add_production.py
 src/simtools/applications/maintain_simulation_model_compare_productions.py
 src/simtools/applications/maintain_simulation_model_verify_production_tables.py
 src/simtools/applications/merge_tables.py
@@ -227,7 +228,6 @@ src/simtools/configuration/commandline_parser.py
 src/simtools/configuration/configurator.py
 src/simtools/corsika/corsika_config.py
 src/simtools/corsika/corsika_histograms.py
-src/simtools/corsika/corsika_histograms_visualize.py
 src/simtools/corsika/primary_particle.py
 src/simtools/data_model/data_reader.py
 src/simtools/data_model/format_checkers.py
@@ -543,8 +543,10 @@ src/simtools/utils/value_conversion.py
 src/simtools/visualization/legend_handlers.py
 src/simtools/visualization/plot_array_layout.py
 src/simtools/visualization/plot_camera.py
+src/simtools/visualization/plot_corsika_histograms.py
 src/simtools/visualization/plot_incident_angles.py
 src/simtools/visualization/plot_pixels.py
+src/simtools/visualization/plot_psf.py
 src/simtools/visualization/plot_simtel_event_histograms.py
 src/simtools/visualization/plot_simtel_events.py
 src/simtools/visualization/plot_tables.py

{gammasimtools-0.20.0 → gammasimtools-0.21.0}/src/gammasimtools.egg-info/entry_points.txt

@@ -25,7 +25,7 @@ simtools-generate-corsika-histograms = simtools.applications.generate_corsika_hi
 simtools-generate-default-metadata = simtools.applications.generate_default_metadata:main
 simtools-generate-regular-arrays = simtools.applications.generate_regular_arrays:main
 simtools-generate-simtel-event-data = simtools.applications.generate_simtel_event_data:main
-simtools-maintain-simulation-model-add-production-table = simtools.applications.maintain_simulation_model_add_production_table:main
+simtools-maintain-simulation-model-add-production = simtools.applications.maintain_simulation_model_add_production:main
 simtools-maintain-simulation-model-compare-productions = simtools.applications.maintain_simulation_model_compare_productions:main
 simtools-maintain-simulation-model-verify-production-tables = simtools.applications.maintain_simulation_model_verify_production_tables:main
 simtools-merge-tables = simtools.applications.merge_tables:main

{gammasimtools-0.20.0 → gammasimtools-0.21.0}/src/simtools/_version.py

@@ -28,7 +28,7 @@ version_tuple: VERSION_TUPLE
 commit_id: COMMIT_ID
 __commit_id__: COMMIT_ID
 
-__version__ = version = '0.20.0'
-__version_tuple__ = version_tuple = (0, 20, 0)
+__version__ = version = '0.21.0'
+__version_tuple__ = version_tuple = (0, 21, 0)
 
-__commit_id__ = commit_id = 'g660b0430e'
+__commit_id__ = commit_id = 'g05415aa4c'

{gammasimtools-0.20.0 → gammasimtools-0.21.0}/src/simtools/applications/db_generate_compound_indexes.py

@@ -46,7 +46,7 @@ def main():  # noqa: D103
     ]
     requested = db.get_db_name(
         db_name=args_dict["db_name"],
-        model_version=args_dict.get("db_simulation_model_version"),
+        db_simulation_model_version=args_dict.get("db_simulation_model_version"),
         model_name=args_dict.get("db_simulation_model"),
     )
     if requested != "all" and requested not in databases:

{gammasimtools-0.20.0 → gammasimtools-0.21.0}/src/simtools/applications/derive_psf_parameters.py

@@ -6,25 +6,22 @@ r"""
     This includes parameters mirror_reflection_random_angle, \
     mirror_align_random_horizontal and mirror_align_random_vertical.
 
-    The telescope zenith angle and the source distance can be set by command line arguments.
-
     The measured cumulative PSF should be provided by using the command line argument data. \
     A file name is expected, in which the file should contain 3 columns: radial distance in mm, \
     differential value of photon intensity and its integral value.
 
-    The derivation is performed through a random search. A number of random combination of the \
-    parameters are tested and the best ones are selected based on the minimum value of \
-    the Root Mean Squared Deviation between data and simulations. The range in which the \
-    parameter are drawn uniformly are defined based on the previous value on the telescope model.
+    The derivation is performed through gradient descent optimization that minimizes either the \
+    Root Mean Squared Deviation (RMSD) between measured and simulated PSF curves (default) or the \
+    Kolmogorov-Smirnov (KS) statistic when the --ks_statistic flag is used.
 
     The optimization workflow includes:
 
     * Loading and preprocessing PSF data from measurement files
-    * Generating random parameter combinations for optimization
-    * Running ray-tracing simulations for each parameter set
-    * Calculating RMSD between measured and simulated PSF curves
-    * Identifying the best-fit parameters with minimum RMSD
-    * Creating comprehensive plots and D80 vs off-axis angle analysis
+    * Running gradient descent optimization to minimize RMSD
+    * Generating cumulative PSF plots for each iteration showing optimization progression
+    * Logging parameter evolution through gradient descent steps
+    * Creating convergence plots showing RMSD and D80 evolution
+    * Automatically generating D80 vs off-axis angle analysis for best parameters
     * Optionally exporting optimized parameters as simulation model files
 
     The assumption are:
@@ -36,7 +33,7 @@ r"""
     One example of the plot generated by this applications are shown below.
 
     .. _derive_psf_parameters_plot:
-    .. image::  images/derive_psf_parameters.png
+    .. image::  images/gradient_descent.png
       :width: 49 %
 
     Command line arguments
@@ -63,47 +60,53 @@ r"""
         If activated, application will be faster by simulating fewer photons.
     write_psf_parameters (activation mode, optional)
         Write the optimized PSF parameters as simulation model parameter files.
-    random_seed (int, optional)
-        Random seed for parameter generation.
-    n_runs (int, optional)
-        Number of parameter combinations to test.
+    rmsd_threshold (float, optional)
+        RMSD threshold for gradient descent convergence (default: 0.007).
+    learning_rate (float, optional)
+        Learning rate for gradient descent optimization (default: 0.01).
+    monte_carlo_analysis (activation mode, optional)
+        Run Monte Carlo analysis to find statistical uncertainties.
 
     Example
     -------
-    LSTN-01 5.0.0
-
-    Runtime < 3 min.
-
-    Get PSF data from the DB:
+    --telescope LSTN-01 --model_version 6.0.0
 
-    .. code-block:: console
 
-        simtools-db-get-file-from-db --file_name PSFcurve_data_v2.txt
 
     Run the application:
 
     .. code-block:: console
 
         simtools-derive-psf-parameters --site North --telescope LSTN-01 \\
-            --model_version 6.0.0 --data tests/resources/PSFcurve_data_v2.txt --plot_all --test
+            --model_version 6.0.0 --data tests/resources/PSFcurve_data_v2.ecsv --plot_all --test
 
     Run with parameter export:
 
     .. code-block:: console
 
-        simtools-derive-psf-parameters --site North --telescope LSTN-01 \\
-            --model_version 6.0.0 \\
-            --data tests/resources/PSFcurve_data_v2.txt --write_psf_parameters
+        simtools-derive-psf-parameters --site North --telescope LSTN-01 --model_version 6.0.0 \\
+            --plot_all --test --rmsd_threshold 0.01 --learning_rate 0.001 \\
+            --data tests/resources/PSFcurve_data_v2.ecsv \\
+            --write_psf_parameters
 
-    The output is saved in simtools-output/tune_psf.
+    Run monte carlo analysis:
+
+    .. code-block:: console
+
+        simtools-derive-psf-parameters --site North --telescope LSTN-01 --model_version 6.0.0 \\
+            --plot_all --test --monte_carlo_analysis \\
+            --data tests/resources/PSFcurve_data_v2.ecsv \\
+            --write_psf_parameters
+
+    The output is saved in simtools-output/derive_psf_parameters.
 
     Output files include:
 
-    * Parameter optimization results in tested_psf_parameters.txt
-    * PSF comparison plots in tune_psf_[telescope].pdf
+    * Gradient descent progression log in psf_gradient_descent_[telescope].log
+    * Gradient descent convergence plots in gradient_descent_convergence_[telescope].png
+    * PSF progression plots showing evolution through iterations (if --plot_all is specified)
     * D80 vs off-axis angle plots (d80_vs_offaxis_cm.png, d80_vs_offaxis_deg.png)
     * Optimized simulation model parameter files (if --write_psf_parameters is specified)
-    * Cumulative PSF plots for all tested combinations (if --plot_all is specified)
 
 """
 
@@ -154,16 +157,32 @@ def _parse():
         required=False,
     )
     config.parser.add_argument(
-        "--random_seed",
-        help="Random seed for parameter generation.",
-        type=int,
-        default=None,
+        "--rmsd_threshold",
+        help=(
+            "RMSD threshold for gradient descent convergence "
+            "(not used with --monte_carlo_analysis)."
+        ),
+        type=float,
+        default=0.01,
+    )
+    config.parser.add_argument(
+        "--learning_rate",
+        help=(
+            "Learning rate for gradient descent optimization "
+            "(not used with --monte_carlo_analysis)."
+        ),
+        type=float,
+        default=0.01,
     )
     config.parser.add_argument(
-        "--n_runs",
-        help="Number of parameter combinations to test.",
-        type=int,
-        default=5,
+        "--monte_carlo_analysis",
+        help="Run analysis to find monte carlo uncertainties.",
+        action="store_true",
+    )
+    config.parser.add_argument(
+        "--ks_statistic",
+        help="Use KS statistic for monte carlo uncertainty analysis.",
+        action="store_true",
     )
     return config.initialize(
         db_config=True,

{gammasimtools-0.20.0 → gammasimtools-0.21.0}/src/simtools/applications/generate_corsika_histograms.py

@@ -165,15 +165,11 @@ r"""
 """
 
 import logging
-import re
 import time
 from pathlib import Path
 
-import numpy as np
-
 import simtools.utils.general as gen
 from simtools.configuration import configurator
-from simtools.corsika import corsika_histograms_visualize
 from simtools.corsika.corsika_histograms import CorsikaHistograms
 from simtools.io import io_handler
 
@@ -284,134 +280,6 @@ def _parse(label, description):
     return config_parser, _
 
 
-def _plot_figures(corsika_histograms_instance, test=False):
-    """
-    Auxiliary function to centralize the plotting functions.
-
-    Parameters
-    ----------
-    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)
-        if plotting_method.startswith("plot_")
-        and "event_header_distribution" not in plotting_method
-    ]
-    if test:
-        plot_function_names = plot_function_names[:2]
-
-    figure_list = []
-    for function_name in plot_function_names:
-        plot_function = getattr(corsika_histograms_visualize, function_name)
-        figures = plot_function(corsika_histograms_instance)
-        for fig in figures:
-            figure_list.append(fig)
-
-    figure_list = np.array(figure_list).flatten()
-    core_name = re.sub(r"\.hdf5$", "", corsika_histograms_instance.hdf5_file_name)
-    output_file_name = Path(corsika_histograms_instance.output_path).joinpath(f"{core_name}.pdf")
-    corsika_histograms_visualize.save_figs_to_pdf(figure_list, output_file_name)
-
-
-def _derive_event_1d_histograms(
-    corsika_histograms_instance, event_1d_header_keys, pdf, hdf5, overwrite=False
-):
-    """
-    Auxiliary function to derive the histograms for the arguments given by event_1d_histograms.
-
-    Parameters
-    ----------
-    corsika_histograms_instance: CorsikaHistograms instance.
-        The CorsikaHistograms instance created in main.
-    event_1d_header_keys: str
-        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.
-    hdf5: bool
-        If true, histograms are saved into hdf5 files.
-    overwrite: bool
-        If true, overwrites the current output hdf5 file.
-    """
-    figure_list = []
-    for event_header_element in event_1d_header_keys:
-        if pdf:
-            figure = corsika_histograms_visualize.plot_1d_event_header_distribution(
-                corsika_histograms_instance, event_header_element
-            )
-            figure_list.append(figure)
-        if hdf5:
-            corsika_histograms_instance.export_event_header_1d_histogram(
-                event_header_element, bins=50, hist_range=None, overwrite=overwrite
-            )
-    if pdf:
-        figures_list = np.array(figure_list).flatten()
-        output_file_name = Path(corsika_histograms_instance.output_path).joinpath(
-            f"{corsika_histograms_instance.hdf5_file_name}_event_1d_histograms.pdf"
-        )
-        corsika_histograms_visualize.save_figs_to_pdf(figures_list, output_file_name)
-
-
-def _derive_event_2d_histograms(
-    corsika_histograms_instance, event_2d_header_keys, pdf, hdf5, overwrite=False
-):
-    """
-    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.
-        The CorsikaHistograms instance created in main.
-    event_2d_header_keys: str
-        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.
-    hdf5: bool
-        If true, histograms are saved into hdf5 files.
-    overwrite: bool
-        If true, overwrites the current output hdf5 file.
-    """
-    figure_list = []
-    for i_event_header_element, _ in enumerate(event_2d_header_keys[::2]):
-        # [::2] to discard the last one in case an odd number of keys are passed
-
-        if len(event_2d_header_keys) % 2 == 1:  # if odd number of keys
-            msg = (
-                "An odd number of keys was passed to generate 2D histograms."
-                "The last key is being ignored."
-            )
-            logger.warning(msg)
-
-        if pdf:
-            figure = corsika_histograms_visualize.plot_2d_event_header_distribution(
-                corsika_histograms_instance,
-                event_2d_header_keys[i_event_header_element],
-                event_2d_header_keys[i_event_header_element + 1],
-            )
-            figure_list.append(figure)
-        if hdf5:
-            corsika_histograms_instance.export_event_header_2d_histogram(
-                event_2d_header_keys[i_event_header_element],
-                event_2d_header_keys[i_event_header_element + 1],
-                bins=50,
-                hist_range=None,
-                overwrite=overwrite,
-            )
-    if pdf:
-        figures_list = np.array(figure_list).flatten()
-        output_file_name = Path(corsika_histograms_instance.output_path).joinpath(
-            f"{corsika_histograms_instance.hdf5_file_name}_event_2d_histograms.pdf"
-        )
-        corsika_histograms_visualize.save_figs_to_pdf(figures_list, output_file_name)
-
-
 def main():  # noqa: D103
     label = Path(__file__).stem
     description = "Generate histograms for the Cherenkov photons saved in the CORSIKA IACT file."
@@ -427,62 +295,17 @@ def main():  # noqa: D103
     corsika_histograms_instance = CorsikaHistograms(
         args_dict["iact_file"], output_path=output_path, hdf5_file_name=args_dict["hdf5_file_name"]
     )
-    if args_dict["telescope_indices"] is not None:
-        try:
-            indices = np.array(args_dict["telescope_indices"]).astype(int)
-        except ValueError:
-            msg = (
-                f"{args_dict['telescope_indices']} not a valid input. "
-                f"Please use integer numbers for telescope_indices"
-            )
-            logger.error(msg)
-            raise
-    else:
-        indices = None
-    # If the hdf5 output file already exists, the results are appended to it.
-    if (Path(corsika_histograms_instance.hdf5_file_name).exists()) and (
-        args_dict["hdf5"] or args_dict["event_1d_histograms"] or args_dict["event_2d_histograms"]
-    ):
-        msg = (
-            f"Output hdf5 file {corsika_histograms_instance.hdf5_file_name} already exists. "
-            f"Overwriting it."
-        )
-        logger.warning(msg)
-        overwrite = True
-    else:
-        overwrite = False
-    corsika_histograms_instance.set_histograms(
-        telescope_indices=indices,
+    corsika_histograms_instance.run_export_pipeline(
         individual_telescopes=args_dict["individual_telescopes"],
         hist_config=args_dict["hist_config"],
+        indices_arg=args_dict["telescope_indices"],
+        write_pdf=args_dict["pdf"],
+        write_hdf5=args_dict["hdf5"],
+        event1d=args_dict["event_1d_histograms"],
+        event2d=args_dict["event_2d_histograms"],
+        test=args_dict["test"],
     )
 
-    # Cherenkov photons
-    if args_dict["pdf"]:
-        _plot_figures(
-            corsika_histograms_instance=corsika_histograms_instance, test=args_dict["test"]
-        )
-    if args_dict["hdf5"]:
-        corsika_histograms_instance.export_histograms(overwrite=overwrite)
-
-    # Event information
-    if args_dict["event_1d_histograms"] is not None:
-        _derive_event_1d_histograms(
-            corsika_histograms_instance,
-            args_dict["event_1d_histograms"],
-            args_dict["pdf"],
-            args_dict["hdf5"],
-            overwrite=not args_dict["hdf5"],
-        )
-    if args_dict["event_2d_histograms"] is not None:
-        _derive_event_2d_histograms(
-            corsika_histograms_instance,
-            args_dict["event_2d_histograms"],
-            args_dict["pdf"],
-            args_dict["hdf5"],
-            overwrite=not (args_dict["hdf5"] or args_dict["event_1d_histograms"]),
-        )
-
     final_time = time.time()
     logger.info(
         f"Finalizing the application. Total time needed: {round(final_time - initial_time)}s."