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

simtools/applications/db_get_array_layouts_from_db.py

@@ -85,7 +85,7 @@ def _parse():
     )
 
 
-def _layout_from_db(args_dict, db_config):
+def _layout_from_db(args_dict):
     """
     Read array elements and their positions from data base using the layout name.
 
@@ -93,8 +93,6 @@ def _layout_from_db(args_dict, db_config):
     ----------
     args_dict : dict
         Dictionary with the command line arguments.
-    db_config : dict
-        Database configuration.
 
     Returns
     -------
@@ -102,7 +100,6 @@ def _layout_from_db(args_dict, db_config):
         Table with array element positions.
     """
     array_model = ArrayModel(
-        db_config=db_config,
         model_version=args_dict["model_version"],
         site=args_dict["site"],
         layout_name=args_dict.get("array_layout_name", None),
@@ -121,14 +118,13 @@ def main():
         if app_context.args.get("site", None) is None:
             raise ValueError("Site must be provided to list available layouts.")
         site_model = SiteModel(
-            db_config=app_context.db_config,
             model_version=app_context.args["model_version"],
             site=app_context.args["site"],
         )
         print(site_model.get_list_of_array_layouts())
     else:
         app_context.logger.info("Array layout: %s", app_context.args["array_layout_name"])
-        layout = _layout_from_db(app_context.args, app_context.db_config)
+        layout = _layout_from_db(app_context.args)
         layout.pprint()
 
         if not app_context.args.get("output_file_from_default", False):

simtools/applications/db_get_file_from_db.py

@@ -60,7 +60,7 @@ def main():
     """Get file from database."""
     app_context = startup_application(_parse)
 
-    db = db_handler.DatabaseHandler(db_config=app_context.db_config)
+    db = db_handler.DatabaseHandler()
     file_id = db.export_model_files(
         dest=app_context.io_handler.get_output_directory(),
         file_names=app_context.args["file_name"],

simtools/applications/db_get_parameter_from_db.py

@@ -108,7 +108,7 @@ def main():
     """Get a parameter entry from DB for a specific telescope or a site."""
     app_context = startup_application(_parse)
 
-    db = db_handler.DatabaseHandler(db_config=app_context.db_config)
+    db = db_handler.DatabaseHandler()
 
     pars = db.get_model_parameter(
         parameter=app_context.args["parameter"],

simtools/applications/db_inspect_databases.py

@@ -32,7 +32,7 @@ def main():
     """Inspect databases."""
     app_context = startup_application(_parse, setup_io_handler=False)
 
-    db = db_handler.DatabaseHandler(db_config=app_context.db_config)
+    db = db_handler.DatabaseHandler()
     # databases without internal databases we don't have rights to modify
     databases = [
         d

simtools/applications/db_upload_model_repository.py

@@ -87,7 +87,7 @@ def main():
     """Application main."""
     app_context = startup_application(_parse)
 
-    db = db_handler.DatabaseHandler(db_config=app_context.db_config)
+    db = db_handler.DatabaseHandler()
     db.print_connection_info()
 
     db_model_upload.add_complete_model(

simtools/applications/derive_ctao_array_layouts.py

@@ -92,7 +92,7 @@ def main():
         branch_name=app_context.args["repository_branch"],
     )
 
-    db = db_handler.DatabaseHandler(db_config=app_context.db_config)
+    db = db_handler.DatabaseHandler()
     db_array_layouts = db.get_model_parameter(
         parameter="array_layouts",
         site=app_context.args["site"],
@@ -107,7 +107,6 @@ def main():
     write_array_layouts(
         array_layouts=merge_array_layouts(db_array_layouts["array_layouts"], ctao_array_layouts),
         args_dict=app_context.args,
-        db_config=app_context.db_config,
     )
 
 

simtools/applications/derive_mirror_rnda.py

@@ -218,9 +218,7 @@ def main():
     """Derive mirror random reflection angle of a single mirror panel."""
     app_context = startup_application(_parse)
 
-    panel_psf = MirrorPanelPSF(
-        app_context.args.get("label"), app_context.args, app_context.db_config
-    )
+    panel_psf = MirrorPanelPSF(app_context.args.get("label"), app_context.args)
     panel_psf.derive_random_reflection_angle(save_figures=True)
     panel_psf.print_results()
     panel_psf.write_optimization_data()

simtools/applications/derive_psf_parameters.py

@@ -187,6 +187,11 @@ def _parse():
         type=float,
         default=0.8,
     )
+    config.parser.add_argument(
+        "--cleanup",
+        help="Remove intermediate *.log and *.lis* files after optimization.",
+        action="store_true",
+    )
     return config.initialize(
         db_config=True,
         simulation_model=["telescope", "model_version", "parameter_version"],
@@ -199,7 +204,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"],

simtools/applications/derive_pulse_shape_parameters.py

@@ -0,0 +1,194 @@
+#!/usr/bin/env python3
+r"""Derive Gaussian sigma and exponential tau from specified rise/fall widths.
+
+Solve (sigma, tau) for a Gaussian convolved with a causal exponential so the
+pulse matches user-provided rise and fall widths between fractional amplitude
+levels (e.g. 0.1-0.9 rise, 0.9-0.1 fall).
+
+Command line arguments
+----------------------
+site (str, required)
+        North or South.
+telescope (str, required)
+        Telescope model name.
+model_version (str, required)
+        Model version.
+parameter_version (str, required)
+    Parameter version.
+rise_width_ns (float, required)
+    Rising-edge width between rise_range fractions (ns).
+fall_width_ns (float, required)
+    Falling-edge width between fall_range fractions (ns).
+rise_range (float float, optional)
+    Fractional amplitudes (low high) for rise width (default: 0.1 0.9).
+fall_range (float float, optional)
+    Fractional amplitudes (high low) for fall width (default: 0.9 0.1).
+dt_ns (float, optional)
+    Time sampling step (ns). Default: 0.1.
+time_margin_ns (float, optional)
+    Margin added at both ends of readout window. Default: 5.
+
+
+Example
+-------
+Derive parameters for a pulse with 2.5 ns rise (10-90%) and
+ 5 ns fall (90-10%) for LSTN-01:
+
+.. code-block:: console
+
+    simtools-derive-pulse-shape-parameters \
+    --site North \
+    --telescope MSTx-NectarCam \
+    --model_version 7.0 \
+    --parameter_version 1.0.0 \
+    --rise_width_ns 2.5 \
+    --fall_width_ns 5.0 \
+    --rise_range 0.1 0.9 \
+    --fall_range 0.9 0.1 \
+    --dt_ns 0.1 \
+    --time_margin_ns 10
+"""
+
+import logging
+
+import simtools.data_model.model_data_writer as writer
+from simtools.application_control import get_application_label, startup_application
+from simtools.configuration import configurator
+from simtools.model.model_utils import initialize_simulation_models
+from simtools.simtel.pulse_shapes import solve_sigma_tau_from_rise_fall
+
+
+def _parse():
+    """Parse command line configuration for parameter derivation."""
+    config = configurator.Configurator(
+        label=get_application_label(__file__),
+        description=(
+            "Derive Gaussian sigma and exponential tau from rise/fall width specifications."
+        ),
+    )
+
+    config.parser.add_argument(
+        "--rise_width_ns",
+        help="Wdth on the rising edge in ns between rise_range fractions.",
+        type=float,
+        required=True,
+    )
+    config.parser.add_argument(
+        "--fall_width_ns",
+        help="Width on the falling edge in ns between fall_range fractions.",
+        type=float,
+        required=True,
+    )
+    config.parser.add_argument(
+        "--rise_range",
+        help="Fractional amplitudes (low high) for rise width, e.g. 0.1 0.9",
+        type=float,
+        nargs=2,
+        default=[0.1, 0.9],
+        required=False,
+    )
+    config.parser.add_argument(
+        "--fall_range",
+        help="Fractional amplitudes (high low) for fall width, e.g. 0.9 0.1",
+        type=float,
+        nargs=2,
+        default=[0.9, 0.1],
+        required=False,
+    )
+    config.parser.add_argument(
+        "--dt_ns",
+        help="Time sampling step in ns used by the solver.",
+        type=float,
+        default=0.1,
+        required=False,
+    )
+    config.parser.add_argument(
+        "--time_margin_ns",
+        help=(
+            "Margin (ns) added to both ends of the instrument readout window when deriving the "
+            "internal time window."
+        ),
+        type=float,
+        default=10.0,
+        required=False,
+    )
+
+    return config.initialize(
+        db_config=True,
+        simulation_model=["site", "telescope", "model_version", "parameter_version"],
+        output=True,
+    )
+
+
+def main():
+    """Run parameter derivation and write results."""
+    app_context = startup_application(_parse)
+    log = logging.getLogger(__name__)
+
+    rise_width_ns = app_context.args["rise_width_ns"]
+    fall_width_ns = app_context.args["fall_width_ns"]
+    rise_range = tuple(app_context.args["rise_range"])
+    fall_range = tuple(app_context.args["fall_range"])
+    dt_ns = app_context.args["dt_ns"]
+    time_margin_ns = app_context.args["time_margin_ns"]
+    site = app_context.args["site"]
+    label = app_context.args.get("label") or get_application_label(__file__)
+    telescope_model, _, _ = initialize_simulation_models(
+        label=label,
+        model_version=app_context.args["model_version"],
+        site=site,
+        telescope_name=app_context.args["telescope"],
+    )
+    fadc_sum_bins = telescope_model.get_parameter_value("fadc_sum_bins")
+
+    window_ns = fadc_sum_bins + time_margin_ns
+    t_start_ns = -window_ns
+    t_stop_ns = window_ns
+
+    sigma_ns, tau_ns = solve_sigma_tau_from_rise_fall(
+        rise_width_ns=rise_width_ns,
+        fall_width_ns=fall_width_ns,
+        dt_ns=dt_ns,
+        rise_range=rise_range,
+        t_start_ns=t_start_ns,
+        t_stop_ns=t_stop_ns,
+    )
+    # Apply reasonable rounding for output precision.
+    sigma_ns = round(sigma_ns, 4)
+    tau_ns = round(tau_ns, 4)
+
+    log.info(
+        f"Derived pulse parameters: sigma={sigma_ns:.6g} ns, tau={tau_ns:.6g} ns "
+        f"(rise={rise_width_ns} ns @ {rise_range}, fall={fall_width_ns} ns @ {fall_range})"
+    )
+
+    output_path = app_context.args.get("output_path")
+    instrument = app_context.args.get("telescope")
+    parameter_version = app_context.args.get("parameter_version")
+
+    writer.ModelDataWriter.dump_model_parameter(
+        parameter_name="flasher_pulse_width",
+        value=sigma_ns,
+        instrument=instrument,
+        parameter_version=parameter_version,
+        output_file="flasher_pulse_width.json",
+        output_path=output_path,
+        unit="ns",
+    )
+    writer.ModelDataWriter.dump_model_parameter(
+        parameter_name="flasher_pulse_exp_decay",
+        value=tau_ns,
+        instrument=instrument,
+        parameter_version=parameter_version,
+        output_file="flasher_pulse_exp_decay.json",
+        output_path=output_path,
+        unit="ns",
+    )
+    log.info(
+        f"Wrote model parameter files flasher_pulse_width.json and "
+        f"flasher_pulse_exp_decay.json (sigma={sigma_ns:.6g} ns, tau={tau_ns:.6g} ns)"
+    )
+
+
+if __name__ == "__main__":
+    main()

simtools/applications/derive_trigger_rates.py

@@ -81,7 +81,7 @@ def main():
     """Derive trigger rates for a single telescope or an array of telescopes."""
     app_context = startup_application(_parse)
 
-    telescope_trigger_rates(app_context.args, app_context.db_config)
+    telescope_trigger_rates(app_context.args)
 
 
 if __name__ == "__main__":

simtools/applications/docs_produce_array_element_report.py

@@ -61,18 +61,12 @@ def main():
             app_context.args.get("all_model_versions"),
         ]
     ):
-        ReportGenerator(
-            app_context.db_config,
-            app_context.args,
-            output_path,
-        ).auto_generate_array_element_reports()
+        ReportGenerator(app_context.args, output_path).auto_generate_array_element_reports()
 
     else:
         model_version = app_context.args["model_version"]
         ReadParameters(
-            app_context.db_config,
-            app_context.args,
-            Path(output_path / f"{model_version}"),
+            app_context.args, Path(output_path / f"{model_version}")
         ).produce_array_element_report()
 
         app_context.logger.info(

simtools/applications/docs_produce_calibration_reports.py

@@ -32,9 +32,7 @@ def main():
 
     output_path = app_context.io_handler.get_output_directory()
 
-    generator = ReportGenerator(
-        db_config=app_context.db_config, args=app_context.args, output_path=output_path
-    )
+    generator = ReportGenerator(args=app_context.args, output_path=output_path)
     generator.auto_generate_calibration_reports()
 
     if app_context.args.get("all_model_versions"):

simtools/applications/docs_produce_model_parameter_reports.py

@@ -41,14 +41,12 @@ def main():
 
     if any([app_context.args.get("all_telescopes"), app_context.args.get("all_sites")]):
         ReportGenerator(
-            app_context.db_config,
             app_context.args,
             output_path,
         ).auto_generate_parameter_reports()
 
     else:
         ReadParameters(
-            app_context.db_config,
             app_context.args,
             output_path,
         ).produce_model_parameter_reports()

simtools/applications/docs_produce_simulation_configuration_report.py

@@ -33,9 +33,7 @@ def main():
 
     output_path = app_context.io_handler.get_output_directory()
 
-    report_generator = ReportGenerator(
-        db_config=app_context.db_config, args=app_context.args, output_path=output_path
-    )
+    report_generator = ReportGenerator(args=app_context.args, output_path=output_path)
     report_generator.auto_generate_simulation_configuration_reports()
 
     app_context.logger.info(

simtools/applications/generate_array_config.py

@@ -46,7 +46,6 @@ def main():
     array_model = ArrayModel(
         label=app_context.args["label"],
         model_version=app_context.args["model_version"],
-        db_config=app_context.db_config,
         site=app_context.args.get("site"),
         layout_name=app_context.args.get("array_layout_name"),
         array_elements=app_context.args.get("array_elements"),