gammasimtools 0.23.0__py3-none-any.whl → 0.25.0__py3-none-any.whl

simtools/simtel/simulator_light_emission.py

@@ -12,6 +12,7 @@ import numpy as np
 from simtools.io import io_handler
 from simtools.model.model_utils import initialize_simulation_models
 from simtools.runners.simtel_runner import SimtelRunner
+from simtools.simtel.simtel_config_writer import SimtelConfigWriter
 from simtools.utils.general import clear_default_sim_telarray_cfg_directories
 from simtools.utils.geometry import fiducial_radius_from_shape
 
@@ -335,18 +336,54 @@ class SimulatorLightEmission(SimtelRunner):
     def _add_flasher_command_options(self):
         """Add flasher options for all telescope types (ff-1m style)."""
         flasher_xyz = self.calibration_model.get_parameter_value_with_unit("flasher_position")
-        camera_radius = fiducial_radius_from_shape(
+        camera_diam_cm = (
             self.telescope_model.get_parameter_value_with_unit("camera_body_diameter")
             .to(u.cm)
-            .value,
-            self.telescope_model.get_parameter_value("camera_body_shape"),
+            .value
         )
+        camera_shape = self.telescope_model.get_parameter_value("camera_body_shape")
+        camera_radius = fiducial_radius_from_shape(camera_diam_cm, camera_shape)
         flasher_wavelength = self.calibration_model.get_parameter_value_with_unit(
             "flasher_wavelength"
         )
         dist_cm = self.calculate_distance_focal_plane_calibration_device().to(u.cm).value
         angular_distribution = self._get_angular_distribution_string_for_sim_telarray()
 
+        # Build pulse table for ff-1m using model width/exp parameters; else use token.
+        pulse_arg = self._get_pulse_shape_string_for_sim_telarray()
+        pulse_shape = self.calibration_model.get_parameter_value("flasher_pulse_shape")
+        width_q = self.calibration_model.get_parameter_value_with_unit("flasher_pulse_width")
+        exp_q = self.calibration_model.get_parameter_value_with_unit("flasher_pulse_exp_decay")
+        if (
+            isinstance(exp_q, u.Quantity)
+            and isinstance(width_q, u.Quantity)
+            and pulse_shape == "Gauss-Exponential"
+        ):
+            try:
+                base_dir = self.io_handler.get_output_directory("pulse_shapes")
+
+                def _sanitize_name(value):
+                    return "".join(
+                        ch if (ch.isalnum() or ch in ("-", "_")) else "_" for ch in str(value)
+                    )
+
+                tel = self.light_emission_config.get("telescope") or "telescope"
+                cal = self.light_emission_config.get("light_source") or "calibration"
+                fname = f"flasher_pulse_shape_{_sanitize_name(tel)}_{_sanitize_name(cal)}.dat"
+                table_path = base_dir / fname
+                fadc_bins = self.telescope_model.get_parameter_value("fadc_sum_bins")
+
+                SimtelConfigWriter.write_lightpulse_table_gauss_expconv(
+                    file_path=table_path,
+                    width_ns=width_q.to(u.ns).value,
+                    exp_decay_ns=exp_q.to(u.ns).value,
+                    fadc_sum_bins=fadc_bins,
+                    time_margin_ns=5.0,
+                )
+                pulse_arg = str(table_path)
+            except (ValueError, OSError) as err:
+                self._logger.warning(f"Failed to write pulse shape table, using token: {err}")
+
         return [
             f"--events {self.light_emission_config['number_of_events']}",
             f"--photons {self.light_emission_config['flasher_photons']}",
@@ -355,7 +392,7 @@ class SimulatorLightEmission(SimtelRunner):
             f"--distance {dist_cm}",
             f"--camera-radius {camera_radius}",
             f"--spectrum {int(flasher_wavelength.to(u.nm).value)}",
-            f"--lightpulse {self._get_pulse_shape_string_for_sim_telarray()}",
+            f"--lightpulse {pulse_arg}",
             f"--angular-distribution {angular_distribution}",
         ]
 
@@ -475,9 +512,8 @@ class SimulatorLightEmission(SimtelRunner):
         str
             The angular distribution string.
         """
-        option_string = self.calibration_model.get_parameter_value(
-            "flasher_angular_distribution"
-        ).lower()
+        opt = self.calibration_model.get_parameter_value("flasher_angular_distribution")
+        option_string = str(opt).lower() if opt is not None else ""
         width = self.calibration_model.get_parameter_value_with_unit(
             "flasher_angular_distribution_width"
         )
@@ -492,6 +528,7 @@ class SimulatorLightEmission(SimtelRunner):
         str
             The pulse shape string.
         """
-        option_string = self.calibration_model.get_parameter_value("flasher_pulse_shape").lower()
+        opt = self.calibration_model.get_parameter_value("flasher_pulse_shape")
+        option_string = str(opt).lower() if opt is not None else ""
         width = self.calibration_model.get_parameter_value_with_unit("flasher_pulse_width")
         return f"{option_string}:{width.to(u.ns).value}" if width is not None else option_string