gammasimtools 0.24.0__py3-none-any.whl → 0.26.0__py3-none-any.whl

simtools/applications/plot_array_layout.py

@@ -62,6 +62,17 @@ highlighted_array_elements : list, optional
     List of array elements to plot with red circles around them.
 legend_location : str, optional
     Location of the legend (default "best").
+bounds : str, optional
+    Axis bounds mode. Use "symmetric" for +-R with padding (default) or "exact" for
+    per-axis min/max bounds.
+padding : float, optional
+    Fractional padding applied around computed extents in both modes (default 0.1).
+x_lim : tuple(float, float), optional
+    Explicit x-axis limits [xmin, xmax] in meters. When provided, overrides derived limits
+    and filters plotted elements by x.
+y_lim : tuple(float, float), optional
+    Explicit y-axis limits [ymin, ymax] in meters. When provided, overrides derived limits
+    and filters plotted elements by y.
 
 Examples
 --------
@@ -89,6 +100,20 @@ Plot layout from a file with a list of telescopes:
     simtools-plot-array-layout
         --array_layout_file tests/resources/telescope_positions-North-ground.ecsv
 
+Use exact bounds with default padding:
+
+.. code-block:: console
+
+    simtools-plot-array-layout --array_layout_name alpha \
+        --site North --model_version 6.0.0 --bounds exact
+
+Use symmetric bounds with custom padding:
+
+.. code-block:: console
+
+    simtools-plot-array-layout --array_layout_name alpha \
+        --site North --model_version 6.0.0 --bounds symmetric --padding 0.15
+
 Plot layout from a parameter file with a list of telescopes:
 
 .. code-block:: console
@@ -116,14 +141,10 @@ Plot layout with some telescopes grayed out and others highlighted:
                                --legend_location "upper right"
 """
 
-import matplotlib as mpl
-import matplotlib.pyplot as plt
-
 import simtools.layout.array_layout_utils as layout_utils
 from simtools.application_control import get_application_label, startup_application
 from simtools.configuration import configurator
-from simtools.visualization import visualize
-from simtools.visualization.plot_array_layout import plot_array_layout
+from simtools.visualization.plot_array_layout import plot_array_layouts
 
 
 def _parse():
@@ -173,6 +194,24 @@ def _parse():
         required=False,
         default=None,
     )
+    config.parser.add_argument(
+        "--x_lim",
+        help="Explicit x-axis limits [xmin xmax] in meters.",
+        type=float,
+        nargs=2,
+        required=False,
+        default=None,
+        metavar=("XMIN", "XMAX"),
+    )
+    config.parser.add_argument(
+        "--y_lim",
+        help="Explicit y-axis limits [ymin ymax] in meters.",
+        type=float,
+        nargs=2,
+        required=False,
+        default=None,
+        metavar=("YMIN", "YMAX"),
+    )
     config.parser.add_argument(
         "--array_layout_name_background",
         help="Name of the background layout array (e.g., test_layout, alpha, 4mst, etc.).",
@@ -201,12 +240,27 @@ def _parse():
         help=(
             "Location of the legend (e.g., 'best', 'upper right', 'upper left', "
             "'lower left', 'lower right', 'right', 'center left', 'center right', "
-            "'lower center', 'upper center', 'center')."
+            "'lower center', 'upper center', 'center', 'no_legend')."
         ),
         type=str,
         required=False,
         default="best",
     )
+    config.parser.add_argument(
+        "--bounds",
+        help=("Axis bounds mode: 'symmetric' uses +-R with padding, 'exact' uses per-axis min/max"),
+        type=str,
+        choices=["symmetric", "exact"],
+        required=False,
+        default="symmetric",
+    )
+    config.parser.add_argument(
+        "--padding",
+        help=("Fractional padding applied around computed extents (used for both modes)."),
+        type=float,
+        required=False,
+        default=0.1,
+    )
     return config.initialize(
         db_config=True,
         simulation_model=[
@@ -220,112 +274,14 @@ def _parse():
     )
 
 
-def read_layouts(args_dict, db_config, logger):
-    """
-    Read array layouts from the database or parameter file.
-
-    Parameters
-    ----------
-    args_dict : dict
-        Dictionary with command line arguments.
-    db_config : dict
-        Database configuration.
-    logger : logging.app_context.logger
-        app_context.logger instance.
-
-    Returns
-    -------
-    list
-        List of array layouts.
-    """
-    if args_dict["array_layout_name"] is not None or args_dict["plot_all_layouts"]:
-        logger.info("Plotting array from DB using layout array name(s).")
-        layouts = layout_utils.get_array_layouts_from_db(
-            args_dict["array_layout_name"],
-            args_dict["site"],
-            args_dict["model_version"],
-            db_config,
-            args_dict["coordinate_system"],
-        )
-        if isinstance(layouts, list):
-            return layouts
-        return [layouts]
-
-    if args_dict["array_layout_parameter_file"] is not None:
-        logger.info("Plotting array from parameter file(s).")
-        return layout_utils.get_array_layouts_from_parameter_file(
-            args_dict["array_layout_parameter_file"],
-            args_dict["model_version"],
-            db_config,
-            args_dict["coordinate_system"],
-        )
-
-    if args_dict["array_layout_file"] is not None:
-        logger.info("Plotting array from telescope table file(s).")
-        return layout_utils.get_array_layouts_from_file(args_dict["array_layout_file"])
-
-    if args_dict["array_element_list"] is not None:
-        logger.info("Plotting array from list of array elements.")
-        return layout_utils.get_array_layouts_using_telescope_lists_from_db(
-            [args_dict["array_element_list"]],
-            args_dict["site"],
-            args_dict["model_version"],
-            db_config,
-            args_dict["coordinate_system"],
-        )
-
-    return []
-
-
 def main():
     """Plot array layout application."""
     app_context = startup_application(_parse)
 
-    layouts = read_layouts(app_context.args, app_context.db_config, app_context.logger)
-
-    if app_context.args.get("array_layout_name_background"):
-        background_layout = layout_utils.get_array_layouts_from_db(
-            app_context.args["array_layout_name_background"],
-            app_context.args["site"],
-            app_context.args["model_version"],
-            app_context.db_config,
-            app_context.args["coordinate_system"],
-        )["array_elements"]
-    else:
-        background_layout = None
-
-    mpl.use("Agg")
-    for layout in layouts:
-        fig_out = plot_array_layout(
-            telescopes=layout["array_elements"],
-            show_tel_label=app_context.args["show_labels"],
-            axes_range=app_context.args["axes_range"],
-            marker_scaling=app_context.args["marker_scaling"],
-            background_telescopes=background_layout,
-            grayed_out_elements=app_context.args["grayed_out_array_elements"],
-            highlighted_elements=app_context.args["highlighted_array_elements"],
-            legend_location=app_context.args["legend_location"],
-        )
-        site_string = ""
-        if layout.get("site") is not None:
-            site_string = f"_{layout['site']}"
-        elif app_context.args["site"] is not None:
-            site_string = f"_{app_context.args['site']}"
-        coordinate_system_string = (
-            f"_{app_context.args['coordinate_system']}"
-            if app_context.args["coordinate_system"] not in layout["name"]
-            else ""
-        )
-        plot_file_name = app_context.args["figure_name"] or (
-            f"array_layout_{layout['name']}{site_string}{coordinate_system_string}"
-        )
-
-        visualize.save_figure(
-            fig_out,
-            app_context.io_handler.get_output_directory() / plot_file_name,
-            dpi=400,
-        )
-        plt.close()
+    layouts, background_layout = layout_utils.read_layouts(app_context.args)
+    plot_array_layouts(
+        app_context.args, app_context.io_handler.get_output_directory(), layouts, background_layout
+    )
 
 
 if __name__ == "__main__":

simtools/applications/plot_simulated_event_distributions.py

@@ -0,0 +1,57 @@
+#!/usr/bin/python3
+
+r"""
+Plot simulated event distributions for shower and/or triggered event data.
+
+Reads reduced event data files and generate histogram plots e.g. for energy or
+core distance distributions.
+
+Command line arguments
+----------------------
+input_file (str, required)
+    Input file path.
+output_path (str, required)
+    Output directory for the generated plots.
+
+Examples
+--------
+Generate plots from a given input file:
+
+.. code-block:: console
+
+    simtools-plot-simulated-event-distributions --input_file path/to/simtel_file.hdf5 \
+                                                --output_path simtools_output/
+
+
+"""
+
+from simtools.application_control import get_application_label, startup_application
+from simtools.configuration import configurator
+from simtools.sim_events.histograms import EventDataHistograms
+from simtools.visualization import plot_simtel_event_histograms
+
+
+def _parse():
+    """Parse command line configuration."""
+    config = configurator.Configurator(
+        label=get_application_label(__file__),
+        description="Plot simulated event distributions for shower and/or triggered event data.",
+    )
+    config.parser.add_argument("--input_file", type=str, required=True, help="Input file path")
+    return config.initialize(db_config=False, output=True)
+
+
+def main():
+    """Plot simulated event distributions."""
+    app_context = startup_application(_parse)
+    app_context.logger.info(f"Loading input file from: {app_context.args['input_file']}")
+
+    histograms = EventDataHistograms(app_context.args["input_file"])
+    histograms.fill()
+    plot_simtel_event_histograms.plot(
+        histograms.histograms, output_path=app_context.io_handler.get_output_directory()
+    )
+
+
+if __name__ == "__main__":
+    main()

simtools/applications/plot_tabular_data.py

@@ -78,7 +78,6 @@ def main():
     plot_tables.plot(
         config=plot_config["plot"],
         output_file=app_context.io_handler.get_output_file(app_context.args["output_file"]),
-        db_config=app_context.db_config,
         data_path=app_context.args.get("table_data_path"),
     )
 

simtools/applications/plot_tabular_data_for_model_parameter.py

@@ -81,15 +81,10 @@ def main():
         telescope=app_context.args.get("telescope"),
         output_path=app_context.io_handler.get_output_directory(),
         plot_type=app_context.args["plot_type"],
-        db_config=app_context.db_config,
     )
 
     for plot_config, output_file in zip(plot_configs, output_files):
-        plot_tables.plot(
-            config=plot_config,
-            output_file=output_file,
-            db_config=app_context.db_config,
-        )
+        plot_tables.plot(config=plot_config, output_file=output_file)
         MetadataCollector.dump(app_context.args, output_file=output_file, add_activity_name=True)
 
 

simtools/applications/production_derive_corsika_limits.py

@@ -134,7 +134,7 @@ def main():
     """Derive limits for energy, radial distance, and viewcone."""
     app_context = startup_application(_parse)
 
-    generate_corsika_limits_grid(app_context.args, app_context.db_config)
+    generate_corsika_limits_grid(app_context.args)
 
 
 if __name__ == "__main__":

simtools/applications/production_generate_grid.py

@@ -132,7 +132,6 @@ def main():
 
     axes = load_axes(app_context.args["axes"])
     site_model = SiteModel(
-        db_config=app_context.db_config,
         model_version=app_context.args["model_version"],
         site=app_context.args["site"],
     )

simtools/applications/run_application.py

@@ -101,7 +101,7 @@ def main():
     """Run several simtools applications using a configuration file."""
     app_context = startup_application(_parse, setup_io_handler=False)
 
-    simtools_runner.run_applications(app_context.args, app_context.db_config, app_context.logger)
+    simtools_runner.run_applications(app_context.args, app_context.logger)
 
 
 if __name__ == "__main__":

simtools/applications/simulate_flasher.py

@@ -29,8 +29,9 @@ Example Usage
     .. code-block:: console
 
         simtools-simulate-flasher --run_mode full_simulation \
-        --light_source MSFx-FlashCam --model_version 6.0.0 \
-        --telescope MSTS-04 --site South
+        --light_source MSFx-NectarCam --model_version 6.0 \
+        --telescope MSTS-04 --site South --run_number 1 \
+        --array_layout_name 1mst
 
 Command Line Arguments
 ----------------------
@@ -119,13 +120,11 @@ def main():
     if app_context.args["run_mode"] == "full_simulation":
         light_source = SimulatorLightEmission(
             light_emission_config=app_context.args,
-            db_config=app_context.db_config,
             label=app_context.args.get("label"),
         )
     elif app_context.args["run_mode"] == "direct_injection":
         light_source = Simulator(
             args_dict=app_context.args,
-            db_config=app_context.db_config,
             label=app_context.args.get("label"),
         )
     else:

simtools/applications/simulate_illuminator.py

@@ -120,7 +120,6 @@ def main():
 
     light_source = SimulatorLightEmission(
         light_emission_config=app_context.args,
-        db_config=app_context.db_config,
         label=app_context.args.get("label"),
     )
     light_source.simulate()

simtools/applications/simulate_pedestals.py

@@ -69,7 +69,7 @@ def _parse():
         help="Calibration run mode",
         type=str,
         required=True,
-        choices=["pedestals", "dark_pedestals", "nsb_only_pedestals"],
+        choices=["pedestals", "pedestals_dark", "pedestals_nsb_only"],
     )
     config.parser.add_argument(
         "--number_of_events",
@@ -108,11 +108,7 @@ def main():
     """Simulate pedestal events."""
     app_context = startup_application(_parse)
 
-    simulator = Simulator(
-        label=app_context.args.get("label"),
-        args_dict=app_context.args,
-        db_config=app_context.db_config,
-    )
+    simulator = Simulator(label=app_context.args.get("label"), args_dict=app_context.args)
     simulator.simulate()
 
 

simtools/applications/simulate_prod.py

@@ -39,10 +39,6 @@ r"""
         It allows the transformation system to keep using sequential run numbers without repetition.
     run (int, required)
         Run number (actual run number will be 'start_run' + 'run').
-    data_directory (str, optional)
-        The location of the output directories corsika-data and simtel-data.
-        the label is added to the data_directory, such that the output
-        will be written to data_directory/label/simtel-data.
     pack_for_grid_register (str, optional)
         Set whether to prepare a tarball for registering the output files on the grid.
         The files are written to the specified directory.
@@ -58,13 +54,6 @@ r"""
         simtools-simulate-prod \\
         --model_version 5.0.0 --site north --primary gamma --azimuth_angle north \\
         --zenith_angle 20 --start_run 0 --run 1
-
-    By default the configuration is saved in simtools-output/test-production
-    together with the actual simulation output in corsika-data and simtel-data within.
-    The location of the latter directories can be set
-    to a different location via the option --data_directory,
-    but the label is always added to the data_directory, such that the output
-    will be written to data_directory/label/simtel-data.
 """
 
 from simtools.application_control import startup_application
@@ -76,15 +65,12 @@ def _parse():
     """Parse command line configuration."""
     config = configurator.Configurator(description="Run simulations for productions")
     config.parser.add_argument(
-        "--data_directory",
+        "--corsika_file",
         help=(
-            "The directory where to save the corsika-data and simtel-data output directories."
-            "the label is added to the data_directory, such that the output"
-            "will be written to data_directory/label/simtel-data."
+            "Path to the CORSIKA input file (only relevant for simulation software 'sim_telarray')."
         ),
-        type=str.lower,
+        type=str,
         required=False,
-        default="./simtools-output/",
     )
     config.parser.add_argument(
         "--pack_for_grid_register",
@@ -141,28 +127,23 @@ def main():
     """Run simulations for productions."""
     app_context = startup_application(_parse, setup_io_handler=False)
 
-    simulator = Simulator(
-        label=app_context.args.get("label"),
-        args_dict=app_context.args,
-        db_config=app_context.db_config,
-    )
+    simulator = Simulator(label=app_context.args.get("label"), args_dict=app_context.args)
 
     simulator.simulate()
     simulator.validate_metadata()
 
-    app_context.logger.info(
-        f"Production run complete for primary {app_context.args['primary']} showers "
-        f"from {app_context.args['azimuth_angle']} azimuth and "
-        f"{app_context.args['zenith_angle']} zenith "
-        f"at {app_context.args['site']} site, using {app_context.args['model_version']} model."
-    )
     if app_context.args["save_reduced_event_lists"]:
         simulator.save_reduced_event_lists()
+
+    simulator.verify_simulations()
+
     if app_context.args.get("pack_for_grid_register"):
         simulator.pack_for_register(app_context.args["pack_for_grid_register"])
     if app_context.args["save_file_lists"]:
         simulator.save_file_lists()
 
+    simulator.report()
+
 
 if __name__ == "__main__":
     main()

simtools/applications/simulate_prod_htcondor_generator.py

@@ -22,7 +22,7 @@ Requirements for the 'simtools-simulate-prod-htcondor-generator' application:
 - Environment parameters required to run CORSIKA and sim_telarray, as well as DB access
   credentials.  These should be listed similarly to a '.env' file and copied to
   'output_path/env.txt'.  Ensure that the path to the simulation software is correctly set to
-  'SIMTOOLS_SIMTEL_PATH=/workdir/sim_telarray'.
+  'SIMTOOLS_SIM_TELARRAY_PATH=/workdir/sim_telarray'.
 
 To submit jobs, change to the output directory and run:
 
@@ -56,6 +56,13 @@ def _parse():
         label=get_application_label(__file__),
         description="Prepare simulations production for HT Condor job submission",
     )
+    config.parser.add_argument(
+        "--number_of_runs",
+        help="Number of runs to be simulated.",
+        type=int,
+        required=True,
+        default=1,
+    )
     config.parser.add_argument(
         "--apptainer_image",
         help="Apptainer image to use for the simulation (full path).",

simtools/applications/submit_array_layouts.py

@@ -1,9 +1,11 @@
 #!/usr/bin/python3
 
 r"""
-    Submit array-layouts definition and corresponding metadata and validate list of telescopes.
+    Read array layouts from file, validate with telescopes in DB, and prepare for submission.
 
-    Includes validation that all defined telescope exists.
+    Validates that all telescope defined in the array layouts exist in the database for the
+    specified model version. Prepares both JSON-style model parameters and corresponding
+    metadata for submission.
 
     Command line arguments
     ----------------------
@@ -19,7 +21,7 @@ r"""
     Example
     -------
 
-    Submit a new array layout dictionary:
+    Submit and validate a new array layout dictionary:
 
     .. code-block:: console
 
@@ -73,7 +75,7 @@ def main():
     """Submit and validate array layouts."""
     app_context = startup_application(_parse)
 
-    db = db_handler.DatabaseHandler(db_config=app_context.db_config)
+    db = db_handler.DatabaseHandler()
 
     array_layouts = validate_array_layouts_with_db(
         production_table=db.read_production_table_from_db(
@@ -82,9 +84,7 @@ def main():
         array_layouts=ascii_handler.collect_data_from_file(app_context.args["array_layouts"]),
     )
 
-    write_array_layouts(
-        array_layouts=array_layouts, args_dict=app_context.args, db_config=app_context.db_config
-    )
+    write_array_layouts(array_layouts=array_layouts, args_dict=app_context.args)
 
 
 if __name__ == "__main__":

simtools/applications/submit_model_parameter_from_external.py

@@ -106,9 +106,7 @@ def main():
         ),
         output_path=output_path,
         metadata_input_dict=app_context.args,
-        db_config=app_context.db_config
-        if app_context.args.get("check_parameter_version")
-        else None,
+        check_db_for_existing_parameter=app_context.args.get("check_parameter_version", False),
     )
 
 

simtools/applications/validate_camera_efficiency.py

@@ -107,7 +107,6 @@ def main():
     app_context = startup_application(_parse)
 
     ce = CameraEfficiency(
-        db_config=app_context.db_config,
         label=app_context.args.get("label"),
         config_data=app_context.args,
     )

simtools/applications/validate_camera_fov.py

@@ -92,7 +92,6 @@ def main():
     tel_model = TelescopeModel(
         site=app_context.args["site"],
         telescope_name=app_context.args["telescope"],
-        db_config=app_context.db_config,
         model_version=app_context.args["model_version"],
         label=label,
     )

simtools/applications/validate_cumulative_psf.py

@@ -130,7 +130,6 @@ def main():
 
     tel_model, site_model, _ = initialize_simulation_models(
         label=app_context.args.get("label"),
-        db_config=app_context.db_config,
         site=app_context.args["site"],
         telescope_name=app_context.args["telescope"],
         model_version=app_context.args["model_version"],
@@ -142,7 +141,6 @@ def main():
     ray = RayTracing(
         telescope_model=tel_model,
         site_model=site_model,
-        simtel_path=app_context.args["simtel_path"],
         zenith_angle=app_context.args["zenith"] * u.deg,
         source_distance=app_context.args["src_distance"] * u.km,
         off_axis_angle=[0.0] * u.deg,