gammasimtools 0.21.0__py3-none-any.whl → 0.22.0__py3-none-any.whl

simtools/simtel/simtel_config_writer.py

@@ -96,6 +96,8 @@ class SimtelConfigWriter:
         """
         self._logger.debug(f"Writing telescope config file {config_file_path}")
 
+        simtel_par = self._get_parameters_for_sim_telarray(parameters, config_file_path)
+
         with open(config_file_path, "w", encoding="utf-8") as file:
             self._write_header(file, "TELESCOPE CONFIGURATION FILE")
 
@@ -104,22 +106,94 @@ class SimtelConfigWriter:
             file.write(f"   echo Configuration for {telescope_name} - TELESCOPE $(TELESCOPE)\n")
             file.write("#endif\n\n")
 
-            for par, value in parameters.items():
-                simtel_name, value = self._convert_model_parameters_to_simtel_format(
-                    names.get_simulation_software_name_from_parameter_name(
-                        par, software_name="sim_telarray"
-                    ),
-                    value["value"],
-                    config_file_path,
-                    None,
-                )
-                if simtel_name:
-                    file.write(f"{simtel_name} = {self._get_value_string_for_simtel(value)}\n")
-            if "stars" not in parameters:  # sim_telarray requires 'stars' to be set
-                file.write("stars = none\n")
+            for simtel_name, simtel_value in simtel_par.items():
+                file.write(f"{simtel_name} = {self._get_value_string_for_simtel(simtel_value)}\n")
             for meta in self._get_sim_telarray_metadata("telescope", parameters, telescope_name):
                 file.write(f"{meta}\n")
 
+    def _get_parameters_for_sim_telarray(self, parameters, config_file_path):
+        """
+        Convert parameter dictionary to sim_telarray configuration file format.
+
+        Accounts for differences between the data models for sim_telarray configuration
+        and the simulation models.
+
+        Parameters
+        ----------
+        parameters: dict
+            Model parameters.
+        config_file_path: str or Path
+            Path of the file to write on.
+
+        Returns
+        -------
+        dict
+            Model parameters in sim_telarray format.
+        """
+        simtel_par = {}
+        for par, value in parameters.items():
+            simtel_name, simtel_value = self._convert_model_parameters_to_simtel_format(
+                names.get_simulation_software_name_from_parameter_name(
+                    par, software_name="sim_telarray"
+                ),
+                value["value"],
+                config_file_path,
+                None,
+            )
+            if simtel_name:
+                simtel_par[simtel_name] = simtel_value
+        if "stars" not in parameters:  # sim_telarray requires 'stars' to be set
+            simtel_par["stars"] = None
+
+        return self._get_flasher_parameters_for_sim_telarray(parameters, simtel_par)
+
+    def _get_flasher_parameters_for_sim_telarray(self, parameters, simtel_par):
+        """
+        Combine flasher pulse time parameters into a single parameter.
+
+        Takes into account that sim_telarray expects a single parameter with a specific name.
+
+        Parameters
+        ----------
+        parameters: dict
+            Model parameters.
+        simtel_par: dict
+            Model parameters in sim_telarray format.
+
+        Returns
+        -------
+        dict
+            Model parameters in sim_telarray format including flasher parameters.
+
+        """
+        if "flasher_pulse_shape" not in parameters and "flasher_pulse_width" not in parameters:
+            return simtel_par
+
+        mapping = {
+            "gauss": "laser_pulse_sigtime",
+            "tophat": "laser_pulse_twidth",
+        }
+
+        shape = parameters.get("flasher_pulse_shape", {}).get("value", "").lower()
+        if "exponential" in shape:
+            simtel_par["laser_pulse_exptime"] = parameters.get("flasher_pulse_exp_decay", {}).get(
+                "value", 0.0
+            )
+        else:
+            simtel_par["laser_pulse_exptime"] = 0.0
+
+        width = parameters.get("flasher_pulse_width", {}).get("value", 0.0)
+
+        simtel_par.update(dict.fromkeys(mapping.values(), 0.0))
+        if shape == "gauss-exponential":
+            simtel_par["laser_pulse_sigtime"] = width
+        elif shape in mapping:
+            simtel_par[mapping[shape]] = width
+        else:
+            self._logger.warning(f"Flasher pulse shape '{shape}' without width definition")
+
+        return simtel_par
+
     def _get_value_string_for_simtel(self, value):
         """
         Return a value string for simtel.

simtools/simtel/simulator_array.py

@@ -34,6 +34,7 @@ class SimulatorArray(SimtelRunner):
         label=None,
         use_multipipe=False,
         sim_telarray_seeds=None,
+        calibration_config=None,
     ):
         """Initialize SimulatorArray."""
         self._logger = logging.getLogger(__name__)
@@ -47,6 +48,7 @@ class SimulatorArray(SimtelRunner):
 
         self.sim_telarray_seeds = sim_telarray_seeds
         self.corsika_config = corsika_config
+        self.calibration_config = calibration_config
         self.io_handler = io_handler.IOHandler()
         self._log_file = None
 
@@ -70,68 +72,69 @@ class SimulatorArray(SimtelRunner):
         """
         command = self._common_run_command(run_number, weak_pointing)
 
-        command += f" {input_file}"
-        command += f" | gzip > {self._log_file} 2>&1 || exit"
-        command += super().get_config_option(
-            "power_law",
-            SimulatorArray.get_power_law_for_sim_telarray_histograms(
-                self.corsika_config.primary_particle
-            ),
-        )
+        if self.calibration_config:
+            command += self._make_run_command_for_calibration_simulations(input_file)
+        else:
+            command += self._make_run_command_for_shower_simulations(input_file)
 
         return clear_default_sim_telarray_cfg_directories(command)
 
-    def make_run_command_for_calibration_simulations(
-        self, run_number=None, input_file=None, calibration_runner_args=None
-    ):
+    def _make_run_command_for_shower_simulations(self, input_file=None):
         """
-        Build and return the command to run sim_telarray for calibration simulations.
+        Build and return the command to run sim_telarray shower simulations.
 
         Parameters
         ----------
-        calibration_runner_args: dict
-            Dictionary with calibration runner arguments.
         input_file: str
             Full path of the input CORSIKA file
         run_number: int (optional)
             run number
+        weak_pointing: bool (optional)
+            Specify weak pointing option for sim_telarray.
 
         Returns
         -------
         str
-            Command to run sim_telarray for pedestal simulations.
+            Command to run sim_telarray.
         """
-        command = self._common_run_command(run_number)
-
-        command += super().get_config_option(
-            "Altitude",
-            self.corsika_config.array_model.site_model.get_parameter_value_with_unit(
-                "reference_point_altitude"
-            ).to_value("m"),
-        )
+        command = f" {input_file}"
+        command += f" | gzip > {self._log_file} 2>&1 || exit"
         command += super().get_config_option(
-            "nsb_scaling_factor", calibration_runner_args["nsb_scaling_factor"]
+            "power_law",
+            SimulatorArray.get_power_law_for_sim_telarray_histograms(
+                self.corsika_config.primary_particle
+            ),
         )
-        if calibration_runner_args.get("stars"):
-            command += super().get_config_option("stars", calibration_runner_args["stars"])
+        return command
 
-        if calibration_runner_args.get("run_mode") in ("pedestals", "nsb_only_pedestals"):
-            command += super().get_config_option(
-                "pedestal_events", calibration_runner_args["number_of_events"]
-            )
-        if calibration_runner_args.get("run_mode") == "nsb_only_pedestals":
-            command += self._nsb_only_pedestals_command()
-        if calibration_runner_args.get("run_mode") == "dark_pedestals":
-            command += super().get_config_option(
-                "dark_events", calibration_runner_args["number_of_events"]
-            )
-        if calibration_runner_args.get("run_mode") == "flasher":
-            command += self._flasher_command(calibration_runner_args)
+    def _make_run_command_for_calibration_simulations(self, input_file=None):
+        """Build sim_telarray command for calibration simulations."""
+        cfg = self.calibration_config
+        altitude = self.corsika_config.array_model.site_model.get_parameter_value_with_unit(
+            "reference_point_altitude"
+        ).to_value("m")
 
-        command += f" {input_file}"
-        command += f" | gzip > {self._log_file} 2>&1 || exit"
+        command = super().get_config_option("Altitude", altitude)
 
-        return clear_default_sim_telarray_cfg_directories(command)
+        for key in ("nsb_scaling_factor", "stars"):
+            if cfg.get(key):
+                command += super().get_config_option(key, cfg[key])
+
+        run_mode = cfg.get("run_mode")
+        if run_mode in ("pedestals", "nsb_only_pedestals"):
+            n_events = cfg.get("number_of_pedestal_events", cfg["number_of_events"])
+            command += super().get_config_option("pedestal_events", n_events)
+        if run_mode == "nsb_only_pedestals":
+            command += self._nsb_only_pedestals_command()
+        if run_mode == "dark_pedestals":
+            n_events = cfg.get("number_of_dark_events", cfg["number_of_events"])
+            command += super().get_config_option("dark_events", n_events)
+        if run_mode == "direct_injection":
+            n_events = cfg.get("number_of_flasher_events", cfg["number_of_events"])
+            command += super().get_config_option("laser_events", n_events)
+
+        command += f" {input_file} | gzip > {self._log_file} 2>&1 || exit"
+        return command
 
     def _common_run_command(self, run_number, weak_pointing=None):
         """Build generic run command for sim_telarray."""
@@ -164,39 +167,6 @@ class SimulatorArray(SimtelRunner):
 
         return command
 
-    def _flasher_command(self, calibration_runner_args):
-        """
-        Generate the command to run sim_telarray for flasher simulations.
-
-        Parameters
-        ----------
-        calibration_runner_args: dict
-            Dictionary with calibration runner arguments.
-
-        Returns
-        -------
-        str
-            Command to run sim_telarray for flasher simulations.
-        """
-        command = super().get_config_option(
-            "laser_events", calibration_runner_args["number_of_events"]
-        )
-        command += super().get_config_option(
-            "laser_photons", calibration_runner_args["flasher_photons"]
-        )
-        command += super().get_config_option(
-            "laser_var_photons", calibration_runner_args["flasher_var_photons"]
-        )
-        command += super().get_config_option(
-            "laser_pulse_exptime", calibration_runner_args["flasher_exp_time"]
-        )
-        command += super().get_config_option(
-            "laser_pulse_sigtime", calibration_runner_args["flasher_sig_time"]
-        )
-        command += super().get_config_option("laser_external_trigger", 1)
-
-        return command
-
     def _nsb_only_pedestals_command(self):
         """
         Generate the command to run sim_telarray for nsb-only pedestal simulations.