gammasimtools 0.16.0__py3-none-any.whl → 0.18.0__py3-none-any.whl

simtools/schemas/model_parameters/discriminator_pulse_shape.schema.yml

@@ -31,3 +31,34 @@ source:
   - Calibration
 simulation_software:
   - name: sim_telarray
+plot_configuration:
+  - type: discriminator_pulse_shape
+    title: 'Discriminator pulse shape'
+    xtitle: 'Time [ns]'
+    ytitle: 'Amplitude [a.u.]'
+    no_markers: true
+    tables:
+      - parameter: discriminator_pulse_shape
+        column_x: 'time'
+        column_y: 'amplitude'
+  - type: discriminator_pulse_shape_low_gain
+    title: 'Discriminator pulse shape (low gain)'
+    xtitle: 'Time [ns]'
+    ytitle: 'Amplitude [a.u.]'
+    no_markers: true
+    tables:
+      - parameter: discriminator_pulse_shape
+        column_x: 'time'
+        column_y: 'amplitude (low gain)'
+  - type: discriminator_pulse_shape_high_and_low_gain
+    title: 'Discriminator pulse shape (high and low gain)'
+    xtitle: 'Time [ns]'
+    ytitle: 'Amplitude [a.u.]'
+    no_markers: true
+    tables:
+      - parameter: discriminator_pulse_shape
+        column_x: 'time'
+        column_y: 'amplitude'
+      - parameter: discriminator_pulse_shape
+        column_x: 'time'
+        column_y: 'amplitude (low gain)'

simtools/schemas/model_parameters/dsum_threshold.schema.yml

@@ -1,6 +1,47 @@
 %YAML 1.2
 ---
 title: Schema for dsum_threshold model parameter
+version: 0.2.0
+meta_schema: simpipe-schema
+meta_schema_url: https://raw.githubusercontent.com/gammasim/simtools/main/src/simtools/schemas/model_parameter_and_data_schema.metaschema.yml
+meta_schema_version: 0.1.0
+name: dsum_threshold
+description: |-
+  The amplitude level above pedestal sum above which a
+  digital sum leads to a telescope trigger.
+  Note that, like for discriminator/comparator and analog sum, the signal
+  must exceed (\'>\') the threshold here before we declare the telescope
+  triggered. The assigned threshold value would have to be one count lower
+  than in a camera-internal trigger implementation (like MSTx-FlashCam) where
+  reaching (\'>=\') the threshold is enough.
+short_description: Amplitude level above which a digital sum leads to a telescope
+  trigger.
+data:
+  - type: int64
+    unit: count
+    default: 0
+    allowed_range:
+      min: 0
+    condition: default_trigger==DigitalSum
+instrument:
+  class: Camera
+  type:
+    - MSTx-NectarCam
+    - MSTx-FlashCam
+activity:
+  setting:
+    - SetParameterFromExternal
+    - SetTriggerThresholdsFromRateScan
+  validation:
+    - ValidateParameterByExpert
+    - ValidateTriggerPerformance
+source:
+  - Observation execution
+simulation_software:
+  - name: sim_telarray
+...
+---
+title: Schema for dsum_threshold model parameter
 version: 0.1.0
 meta_schema: simpipe-schema
 meta_schema_url: https://raw.githubusercontent.com/gammasim/simtools/main/src/simtools/schemas/model_parameter_and_data_schema.metaschema.yml

simtools/schemas/model_parameters/fadc_pulse_shape.schema.yml

@@ -37,3 +37,15 @@ source:
   - Calibration
 simulation_software:
   - name: sim_telarray
+plot_configuration:
+  - type: fadc_pulse_shape
+    title: 'FADC pulse shape'
+    xtitle: 'Time [ns]'
+    ytitle: 'Amplitude [a.u.]'
+    xscale: 'linear'
+    yscale: 'linear'
+    no_markers: true
+    tables:
+      - parameter: fadc_pulse_shape
+        column_x: 'time'
+        column_y: 'amplitude'

simtools/schemas/model_parameters/lightguide_efficiency_vs_incidence_angle.schema.yml

@@ -39,3 +39,13 @@ source:
   - Calibration
 simulation_software:
   - name: simtools
+plot_configuration:
+  - type: lightguide_efficiency_vs_incidence_angle
+    title: 'Light guide efficiency'
+    xtitle: 'Incident angle [deg]'
+    ytitle: 'Efficiency'
+    no_markers: true
+    tables:
+      - parameter: lightguide_efficiency_vs_incidence_angle
+        column_x: 'angle'
+        column_y: 'efficiency'

simtools/schemas/model_parameters/mirror_reflectivity.schema.yml

@@ -38,3 +38,13 @@ source:
   - Calibration
 simulation_software:
   - name: sim_telarray
+plot_configuration:
+  - type: mirror_reflectivity
+    title: 'Mirror reflectivity'
+    xtitle: 'Wavelength [nm]'
+    ytitle: 'Reflectivity'
+    no_markers: true
+    tables:
+      - parameter: mirror_reflectivity
+        column_x: 'wavelength'
+        column_y: 'reflectivity'

simtools/schemas/model_parameters/nsb_reference_spectrum.schema.yml

@@ -32,3 +32,15 @@ source:
   - Calibration
 simulation_software:
   - name: simtools
+plot_configuration:
+  - type: nsb_reference_spectrum
+    title: 'NSB spectrum'
+    xtitle: 'Wavelength [nm]'
+    ytitle: 'Differential photon rate [10^9 / (nm s m^2 sr)]'
+    xscale: 'linear'
+    yscale: 'log'
+    no_markers: true
+    tables:
+      - parameter: nsb_reference_spectrum
+        column_x: 'wavelength'
+        column_y: 'differential photon rate'

simtools/schemas/model_parameters/pm_photoelectron_spectrum.schema.yml

@@ -39,3 +39,22 @@ source:
   - Calibration
 simulation_software:
   - name: sim_telarray
+plot_configuration:
+  - type: pm_photoelectron_spectrum
+    title: 'Single p.e. response'
+    xtitle: 'Amplitude'
+    ytitle: 'Response'
+    xscale: 'linear'
+    yscale: 'log'
+    xlim: [0., 40.]
+    ylim: [1.e-11, null]
+    no_markers: true
+    tables:
+      - parameter: pm_photoelectron_spectrum
+        label: 'prompt'
+        column_x: 'amplitude'
+        column_y: 'response'
+      - parameter: pm_photoelectron_spectrum
+        label: 'prompt + afterpulsing'
+        column_x: 'amplitude'
+        column_y: 'response_with_ap'

simtools/schemas/model_parameters/quantum_efficiency.schema.yml

@@ -40,3 +40,13 @@ source:
   - Calibration
 simulation_software:
   - name: sim_telarray
+plot_configuration:
+  - type: quantum_efficiency
+    title: 'Quantum/detection efficiency'
+    xtitle: 'Wavelength [nm]'
+    ytitle: 'Quantum/detection efficiency'
+    no_markers: true
+    tables:
+      - parameter: quantum_efficiency
+        column_x: 'wavelength'
+        column_y: 'efficiency'

simtools/schemas/plot_configuration.metaschema.yml

@@ -6,157 +6,146 @@ description: YAML representation of plot configuration metaschema
 version: 0.1.0
 name: plot_configuration.metaschema
 type: object
-additionalProperties: false
 
 definitions:
   SimtoolsPlotConfiguration:
     type: object
     additionalProperties: false
     properties:
-      CTA_SIMPIPE:
-        "$ref": "#/definitions/CTASIMPIPE"
-      SCHEMA_VERSION:
+      plot:
+        "$ref": "#/definitions/plot"
+      schema_version:
         type: string
         description: "Version of the schema."
-      SCHEMA_URL:
+      schema_url:
         type: string
         format: uri
         description: "URL of the schema."
-      SCHEMA_NAME:
+      schema_name:
         type: string
         description: "Name of the schema."
     required:
-      - CTA_SIMPIPE
+      - plot
     title: "SimtoolsPlotConfiguration"
-  CTASIMPIPE:
-    type: object
-    additionalProperties: false
-    properties:
-      PLOT:
-        "$ref": "#/definitions/Plot"
-    required:
-      - PLOT
-    title: "CTASimpipe"
-  Plot:
+  plot:
     description: "Plot configuration of simtools data visualization."
     type: object
     additionalProperties: false
     properties:
-      TYPE:
+      type:
         type: string
         description: "Type of the plot (e.g., mirror reflectivity, quantum_efficiency)."
-      TITLE:
+      title:
         type: string
         description: "Title of the plot."
-      XTITLE:
+      xtitle:
         type: string
         description: "Title of x-axis."
-      YTITLE:
+      ytitle:
         type: string
         description: "Title of y-axis."
-      XSCALE:
+      xscale:
         type: string
         enum: ["linear", "log"]
         description: "Scale of x-axis (linear or log)."
-      YSCALE:
+      yscale:
         type: string
         enum: ["linear", "log"]
         description: "Scale of y-axis (linear or log)."
-      XLIM:
+      xlim:
         type: array
         items:
           type: [number, "null"]
           description: "Limits for x-axis [min, max]. Use null for auto-limit."
         minItems: 2
         maxItems: 2
-      YLIM:
+      ylim:
         type: array
         items:
           type: [number, "null"]
           description: "Limits for y-axis [min, max]. Use null for auto-limit."
         minItems: 2
         maxItems: 2
-      ERROR_TYPE:
+      error_type:
         type: string
         enum: ["fill_between", "errorbar", "none"]
         description: "Type of errors."
-      NO_MARKERS:
+      no_markers:
         type: boolean
         description: "Whether to display markers on data points."
-      PLOT_RATIO:
+      plot_ratio:
         type: boolean
         description: "Whether to display ratio plot."
-      TABLES:
+      tables:
         type: array
         description: "List of tables to plot."
         items:
-          "$ref": "#/definitions/TableConfig"
+          "$ref": "#/definitions/table_config"
     required:
-      - TYPE
-      - TITLE
-      - XTITLE
-      - YTITLE
-      - TABLES
+      - type
+      - title
+      - xtitle
+      - ytitle
+      - tables
     title: "Plot"
-  TableConfig:
+  table_config:
     type: object
     description: "Configuration for a data table to plot."
     additionalProperties: false
     properties:
-      PARAMETER:
+      parameter:
         type: string
         description: "Parameter name to retrieve data for."
-      FILE_NAME:
+      file_name:
         type: string
         description: "Path to the data file."
-      TYPE:
+      type:
         type: string
         description: "Type of data file."
-      TELESCOPE:
+      telescope:
         type: string
         description: "Telescope descriptor to retrieve data for."
-      SITE:
+      site:
         type: string
         description: "Site name (North/South)."
-      MODEL_VERSION:
+      model_version:
         type: string
         description: "Model version to use."
-      PARAMETER_VERSION:
+      parameter_version:
         type: string
         description: "Parameter version to use."
-      LABEL:
+      label:
         type: string
         description: "Label for the plot legend."
-      COLUMN_X:
+      column_x:
         type: string
         description: "Column name to use for x-axis."
-      COLUMN_Y:
+      column_y:
         type: string
         description: "Column name to use for y-axis."
-      COLUMN_X_ERR:
+      column_x_err:
         type: string
         description: "Column name to use for x-axis error."
-      COLUMN_Y_ERR:
+      column_y_err:
         type: string
         description: "Column name to use for y-axis error."
-      NORMALIZE_Y:
+      normalize_y:
         type: boolean
         description: "Whether to normalize y values."
-      SELECT_VALUES:
+      select_values:
         type: object
         description: "Selection criteria for data filtering."
         properties:
-          COLUMN_NAME:
+          column_name:
             type: string
             description: "Column name to use for filtering."
-          VALUE:
+          value:
             type: [number, string]
             description: "Value to filter by."
         required:
-          - COLUMN_NAME
-          - VALUE
+          - column_name
+          - value
     required:
-      - COLUMN_X
-      - COLUMN_Y
-      - LABEL
+      - column_x
+      - column_y
     title: "TableConfig"

simtools/schemas/production_configuration_metrics.schema.yml

@@ -29,12 +29,12 @@ definitions:
       description:
         type: string
         description: Description of the error metric.
-      target_error:
+      target_uncertainty:
         $ref: '#/definitions/TargetError'
       energy_range:
         $ref: '#/definitions/EnergyRange'
     required:
-      - target_error
+      - target_uncertainty
       - energy_range
 
   TargetError:

simtools/simtel/simtel_config_writer.py

@@ -54,12 +54,20 @@ class SimtelConfigWriter:
         Layout name.
     label: str
         Instance label. Important for output file naming.
+    simtel_path: str or Path
+        Path to the sim_telarray installation directory.
     """
 
     TAB = " " * 3
 
     def __init__(
-        self, site, model_version, layout_name=None, telescope_model_name=None, label=None
+        self,
+        site,
+        model_version,
+        layout_name=None,
+        telescope_model_name=None,
+        label=None,
+        simtel_path=None,
     ):
         """Initialize SimtelConfigWriter."""
         self._logger = logging.getLogger(__name__)
@@ -70,10 +78,9 @@ class SimtelConfigWriter:
         self._label = label
         self._layout_name = layout_name
         self._telescope_model_name = telescope_model_name
+        self._simtel_path = simtel_path
 
-    def write_telescope_config_file(
-        self, config_file_path, parameters, telescope_name=None, write_dummy_config=False
-    ):
+    def write_telescope_config_file(self, config_file_path, parameters, telescope_name=None):
         """
         Write the sim_telarray config file for a single telescope.
 
@@ -85,8 +92,6 @@ class SimtelConfigWriter:
             Model parameters
         telescope_name: str
             Name of the telescope (use self._telescope_model_name if None)
-        write_dummy_config: bool
-            Flag to write a dummy telescope configuration file.
         """
         self._logger.debug(f"Writing telescope config file {config_file_path}")
 
@@ -97,8 +102,6 @@ class SimtelConfigWriter:
             file.write("#ifdef TELESCOPE\n")
             file.write(f"   echo Configuration for {telescope_name} - TELESCOPE $(TELESCOPE)\n")
             file.write("#endif\n\n")
-            if write_dummy_config:
-                file.write("#define DUMMY_CONFIG 1\n")
 
             for par, value in parameters.items():
                 simtel_name, value = self._convert_model_parameters_to_simtel_format(
@@ -252,6 +255,8 @@ class SimtelConfigWriter:
             file.write(self.TAB + f"echo ModelVersion: {self._model_version}\n")
             file.write(self.TAB + "echo *****************************\n\n")
 
+            self._write_simtools_parameters(file)
+
             self._write_site_parameters(
                 file,
                 site_model.parameters,
@@ -392,6 +397,23 @@ class SimtelConfigWriter:
         header += f"{comment_char}\n"
         file.write(header)
 
+    def _write_simtools_parameters(self, file):
+        """Write simtools-specific parameters."""
+        meta_items = {
+            "simtools_version": simtools.version.__version__,
+            "simtools_model_production_version": self._model_version,
+        }
+        try:
+            build_opts = gen.collect_data_from_file(Path(self._simtel_path) / "build_opts.yml")
+            for key, value in build_opts.items():
+                meta_items[f"simtools_{key}"] = value
+        except (FileNotFoundError, TypeError):
+            pass  # don't expect build_opts.yml to be present on all systems
+
+        file.write(f"{self.TAB}% Simtools parameters\n")
+        for key, value in meta_items.items():
+            file.write(f"{self.TAB}metaparam global set {key} = {value}\n")
+
     def _write_site_parameters(
         self, file, site_parameters, model_path, telescope_model, sim_telarray_seeds=None
     ):
@@ -556,9 +578,9 @@ class SimtelConfigWriter:
             "disc_bins": 10,
             "fadc_sum_bins": 10,
             "fadc_sum_offset": 0,
-            "asum_threshold": 0,
-            "dsum_threshold": 0,
-            "discriminator_threshold": 1,
+            "asum_threshold": 9999,
+            "dsum_threshold": 9999,
+            "discriminator_threshold": 9999,
             "fadc_amplitude": 1.0,
             "discriminator_amplitude": 1.0,
         }
@@ -567,9 +589,7 @@ class SimtelConfigWriter:
             if key in parameters:
                 parameters[key]["value"] = val
 
-        self.write_telescope_config_file(
-            config_file_path, parameters, telescope_name, write_dummy_config=True
-        )
+        self.write_telescope_config_file(config_file_path, parameters, telescope_name)
 
         config_file_directory = Path(config_file_path).parent
         self._write_dummy_mirror_list_files(config_file_directory, telescope_name)