gammasimtools 0.17.0__py3-none-any.whl → 0.19.0__py3-none-any.whl

simtools/simtel/simulator_array.py

@@ -2,7 +2,7 @@
 
 import logging
 
-from simtools.io_operations import io_handler
+from simtools.io import io_handler
 from simtools.runners.simtel_runner import InvalidOutputFileError, SimtelRunner
 from simtools.utils.general import clear_default_sim_telarray_cfg_directories
 
@@ -68,6 +68,73 @@ class SimulatorArray(SimtelRunner):
         str
             Command to run sim_telarray.
         """
+        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
+            ),
+        )
+
+        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
+    ):
+        """
+        Build and return the command to run sim_telarray for calibration 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
+
+        Returns
+        -------
+        str
+            Command to run sim_telarray for pedestal simulations.
+        """
+        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 += super().get_config_option(
+            "nsb_scaling_factor", calibration_runner_args["nsb_scaling_factor"]
+        )
+        if calibration_runner_args.get("stars"):
+            command += super().get_config_option("stars", calibration_runner_args["stars"])
+
+        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)
+
+        command += f" {input_file}"
+        command += f" | gzip > {self._log_file} 2>&1 || exit"
+
+        return clear_default_sim_telarray_cfg_directories(command)
+
+    def _common_run_command(self, run_number, weak_pointing=None):
+        """Build generic run command for sim_telarray."""
         config_dir = self.corsika_config.array_model.get_config_directory()
         self._log_file = self.get_file_name(file_type="log", run_number=run_number)
         histogram_file = self.get_file_name(file_type="histogram", run_number=run_number)
@@ -83,12 +150,6 @@ class SimulatorArray(SimtelRunner):
         command += super().get_config_option(
             "telescope_phi", self.corsika_config.azimuth_angle, weak_pointing
         )
-        command += super().get_config_option(
-            "power_law",
-            SimulatorArray.get_power_law_for_sim_telarray_histograms(
-                self.corsika_config.primary_particle
-            ),
-        )
         command += super().get_config_option("histogram_file", histogram_file)
         command += super().get_config_option("random_state", "none")
         if self.sim_telarray_seeds and self.sim_telarray_seeds.get("random_instrument_instances"):
@@ -100,10 +161,77 @@ class SimulatorArray(SimtelRunner):
             command += super().get_config_option("random_seed", self.sim_telarray_seeds["seed"])
         command += super().get_config_option("show", "all")
         command += super().get_config_option("output_file", output_file)
-        command += f" {input_file}"
-        command += f" | gzip > {self._log_file} 2>&1 || exit"
 
-        return clear_default_sim_telarray_cfg_directories(command)
+        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.
+
+        Returns
+        -------
+        str
+            Command to run sim_telarray.
+        """
+        null_values = [
+            "fadc_noise",
+            "fadc_lg_noise",
+            "qe_variation",
+            "gain_variation",
+            "fadc_var_pedestal",
+            "fadc_err_pedestal",
+            "fadc_sysvar_pedestal",
+            "fadc_dev_pedestal",
+        ]
+        null_command_parts = []
+        for param in null_values:
+            null_command_parts.append(super().get_config_option(param, 0.0))
+        command = " ".join(null_command_parts)
+
+        one_values = [
+            "fadc_lg_var_pedestal",
+            "fadc_lg_err_pedestal",
+            "fadc_lg_dev_pedestal",
+            "fadc_lg_sysvar_pedestal",
+        ]
+        one_command_parts = []
+        for param in one_values:
+            one_command_parts.append(super().get_config_option(param, -1.0))
+        command += " " + " ".join(one_command_parts)
+        return command
 
     def _check_run_result(self, run_number=None):
         """

simtools/simtel/simulator_camera_efficiency.py

@@ -3,7 +3,7 @@
 import logging
 from pathlib import Path
 
-from simtools.model.model_parameter import InvalidModelParameterError
+from simtools.io import ascii_handler
 from simtools.runners.simtel_runner import SimtelRunner
 from simtools.utils import general
 
@@ -82,41 +82,27 @@ class SimulatorCameraEfficiency(SimtelRunner):
         """Prepare the command used to run testeff."""
         self._logger.debug("Preparing the command to run testeff")
 
-        # Processing camera pixel features
         pixel_shape = self._telescope_model.camera.get_pixel_shape()
         pixel_shape_cmd = "-hpix" if pixel_shape in [1, 3] else "-spix"
         pixel_diameter = self._telescope_model.camera.get_pixel_diameter()
 
-        # Processing focal length
         focal_length = self._telescope_model.get_telescope_effective_focal_length("m", True)
 
-        # Processing mirror class
-        mirror_class = 1
-        try:
-            mirror_class = self._telescope_model.get_parameter_value("mirror_class")
-        except InvalidModelParameterError:
-            pass
-
-        # Processing camera transmission
-        camera_transmission = 1
-        try:
-            camera_transmission = self._telescope_model.get_parameter_value("camera_transmission")
-        except KeyError:
-            pass
-
-        # Processing camera filter
-        # A special case is testeff does not support 2D distributions
+        mirror_class = self._telescope_model.get_parameter_value("mirror_class")
+        curvature_radius = self._get_curvature_radius(mirror_class)
+        camera_transmission = self._telescope_model.get_parameter_value("camera_transmission")
+
         camera_filter_file = self._telescope_model.get_parameter_value("camera_filter")
+        # testeff does not support 2D distributions
         if self._telescope_model.is_file_2d("camera_filter"):
             camera_filter_file = self._get_one_dim_distribution(
                 "camera_filter", "camera_filter_incidence_angle"
             )
 
-        # Processing mirror reflectivity
-        # A special case is testeff does not support 2D distributions
         mirror_reflectivity = self._telescope_model.get_parameter_value("mirror_reflectivity")
         if mirror_class == 2:
             mirror_reflectivity_secondary = mirror_reflectivity
+        # testeff does not support 2D distributions
         if self._telescope_model.is_file_2d("mirror_reflectivity"):
             mirror_reflectivity = self._get_one_dim_distribution(
                 "mirror_reflectivity", "primary_mirror_incidence_angle"
@@ -136,6 +122,7 @@ class SimulatorCameraEfficiency(SimtelRunner):
         command += f" -alt {self._site_model.get_parameter_value('corsika_observation_level')}"
         command += f" -fatm {self._site_model.get_parameter_value('atmospheric_transmission')}"
         command += f" -flen {focal_length}"
+        command += f" -fcur {curvature_radius:.3f}"
         command += f" {pixel_shape_cmd} {pixel_diameter}"
         if mirror_class == 0:
             command += f" -fmir {self._telescope_model.get_parameter_value('mirror_list')}"
@@ -227,7 +214,8 @@ class SimulatorCameraEfficiency(SimtelRunner):
             two_dim_distribution, incidence_angle_distribution
         )
         new_file_name = (
-            f"weighted_average_1D_{self._telescope_model.get_parameter_value(two_dim_parameter)}"
+            f"weighted_average_1D_{weighting_distribution_parameter}"
+            f"_{self._telescope_model.get_parameter_value(two_dim_parameter)}"
         )
         return self._telescope_model.export_table_to_model_directory(
             new_file_name, distribution_to_export
@@ -257,7 +245,7 @@ class SimulatorCameraEfficiency(SimtelRunner):
             self._telescope_model.config_file_directory / Path(nsb_spectrum_file).name
         )
 
-        lines = general.read_file_encoded_in_utf_or_latin(nsb_spectrum_file)
+        lines = ascii_handler.read_file_encoded_in_utf_or_latin(nsb_spectrum_file)
 
         with open(validated_nsb_spectrum_file, "w", encoding="utf-8") as file:
             for line in lines:
@@ -271,3 +259,24 @@ class SimulatorCameraEfficiency(SimtelRunner):
                 else:
                     file.write(line)
         return validated_nsb_spectrum_file
+
+    def _get_curvature_radius(self, mirror_class=1):
+        """Get radius of curvature of dish."""
+        if mirror_class == 2:
+            return (
+                self._telescope_model.get_parameter_value_with_unit("primary_mirror_diameter")
+                .to("m")
+                .value
+            )
+
+        if self._telescope_model.get_parameter_value("parabolic_dish"):
+            return (
+                2.0
+                * self._telescope_model.get_parameter_value_with_unit("dish_shape_length")
+                .to("m")
+                .value
+            )
+
+        return (
+            self._telescope_model.get_parameter_value_with_unit("dish_shape_length").to("m").value
+        )