gammasimtools 0.9.0__py3-none-any.whl → 0.10.0__py3-none-any.whl

simtools/applications/generate_regular_arrays.py

@@ -1,27 +1,27 @@
 #!/usr/bin/python3
 
 """
-    Make a regular array of telescopes and save it as astropy table.
+Make a regular array of telescopes and save it as astropy table.
 
-    The arrays consist of one telescope at the center of the array and or of 4 telescopes
-    in a square grid. These arrays are used for trigger rate simulations.
+The arrays consist of one telescope at the center of the array and or of 4 telescopes
+in a square grid. These arrays are used for trigger rate simulations.
 
-    The array layout files created will be available at the data/layout directory.
+The array layout files created will be available at the data/layout directory.
 
-    Command line arguments
-    ----------------------
-    site (str, required)
-        observatory site (e.g., North or South).
-    model_version (str, optional)
-        Model version to use (e.g., 6.0.0). If not provided, the latest version is used.
+Command line arguments
+----------------------
+site (str, required)
+    observatory site (e.g., North or South).
+model_version (str, optional)
+    Model version to use (e.g., 6.0.0). If not provided, the latest version is used.
 
-    Example
-    -------
-    Runtime < 10 s.
+Example
+-------
+Runtime < 10 s.
 
-    .. code-block:: console
+.. code-block:: console
 
-        simtools-generate-regular-arrays --site=North
+    simtools-generate-regular-arrays --site=North
 """
 
 import logging

simtools/applications/generate_simtel_array_histograms.py

@@ -13,8 +13,6 @@ r"""
         Name of the histogram files to be plotted.
         It can be given as the histogram file names (more than one option allowed) or as a text
         file with the names of the histogram files in it.
-    figure_name (str, required)
-        File name for the pdf output (without extension).
     pdf (bool, optional)
         If set, histograms are saved into pdf files.
         One pdf file contains all the histograms found in the file.
@@ -29,6 +27,8 @@ r"""
         If the output output_file_name.hdf5 file already exists and hdf5 is set, the tables
         associated to hdf5 will be overwritten. The remaining tables, if any, will stay
         untouched.
+    test: bool
+        Test option. Generate only two histograms for testing purposes.
 
     Raises
     ------
@@ -108,38 +108,6 @@ def _parse(label, description):
     return config_parser
 
 
-def build_histogram_files(config_parser, logger):
-    """
-    Build a list of histogram files from command line arguments.
-
-    Parameters
-    ----------
-    config_parser: dict
-        Parsed command line arguments.
-    logger: logging.Logger
-        Logger object for logging messages.
-
-    Returns
-    -------
-    list
-        List of histogram file paths.
-    """
-    histogram_files = []
-    for one_file in config_parser["hist_file_names"]:
-        try:
-            if Path(one_file).suffix in [".zst", ".simtel", ".hdata"]:
-                histogram_files.append(one_file)
-            else:
-                with open(one_file, encoding="utf-8") as file:
-                    for line in file:
-                        histogram_files.append(line.strip())
-        except FileNotFoundError as exc:
-            msg = f"{one_file} is not a file."
-            logger.error(msg)
-            raise FileNotFoundError from exc
-    return histogram_files
-
-
 def check_and_log_overwrite(config_parser, logger):
     """
     Check if the output hdf5 file already exists and log a warning if it does.
@@ -197,36 +165,31 @@ def create_pdf(simtel_histograms, output_file_name, config_parser, logger):
         logger.info(f"Wrote histograms to the pdf file {output_file_name}.pdf")
 
 
-def export_to_hdf5(simtel_histograms, output_file_name, overwrite, config_parser, logger):
-    """Export histograms to an HDF5 file."""
-    if config_parser["hdf5"]:
-        logger.info(f"Wrote histograms to the hdf5 file {output_file_name}.hdf5")
-        simtel_histograms.export_histograms(f"{output_file_name}.hdf5", overwrite=overwrite)
-
-
 def main():  # noqa: D103
     label = Path(__file__).stem
-    description = "Display the simtel_array histograms."
+    description = "Display simtel_array histograms and/or write them into hdf5 format."
     io_handler_instance = io_handler.IOHandler()
     config_parser = _parse(label, description)
     output_path = io_handler_instance.get_output_directory(label, sub_dir="application-plots")
     logger = logging.getLogger()
     logger.setLevel(gen.get_log_level_from_user(config_parser["log_level"]))
-    logger.info("Starting the application.")
 
-    histogram_files = build_histogram_files(config_parser, logger)
+    histogram_files = gen.get_list_of_files_from_command_line(
+        config_parser["hist_file_names"], [".zst", ".simtel", ".hdata"]
+    )
 
     # If no output name is passed, the tool gets the name of the first histogram of the list
     if config_parser["output_file_name"] is None:
         config_parser["output_file_name"] = Path(histogram_files[0]).absolute().name
     output_file_name = Path(output_path).joinpath(f"{config_parser['output_file_name']}")
 
-    # If the hdf5 output file already exists, it is overwritten
-    overwrite = check_and_log_overwrite(config_parser, logger)
-
     simtel_histograms = SimtelIOHistograms(histogram_files)
     create_pdf(simtel_histograms, output_file_name, config_parser, logger)
-    export_to_hdf5(simtel_histograms, output_file_name, overwrite, config_parser, logger)
+    if config_parser["hdf5"]:
+        simtel_histograms.export_histograms(
+            f"{output_file_name}.hdf5",
+            overwrite=check_and_log_overwrite(config_parser, logger),
+        )
 
 
 if __name__ == "__main__":

simtools/applications/production_generate_simulation_config.py

@@ -1,11 +1,24 @@
 #!/usr/bin/python3
 
 r"""
-Configure a simulation based on command-line arguments.
+Derive simulation configuration parameters for a simulation production.
 
-This application configures and
-generates simulation parameters for a specific grid point in a statistical uncertainty
-evaluation setup.
+Derived simulation configuration parameters include:
+
+* energy range
+* shower core scatter radius
+* view cone radius
+* total number of events to be simulated
+
+Configuration parameters depend on characteristics of the observations, especially elevation,
+azimuth, and night sky background.
+
+The configuration parameters are derived according to the required precision. The metrics are:
+
+* statistical uncertainty on the determination of the effective area as function of primary energy
+* fraction of lost events to the selected core scatter and view cone radius (to be implemented)
+* statistical uncertainty of the energy migration matrix as function of primary energy
+  (to be implemented)
 
 Command line arguments
 ----------------------
@@ -38,9 +51,9 @@ To run the simulation configuration, execute the script as follows:
         --nsb 0.3 --ctao_data_level "A" --science_case "high_precision" \
         --file_path tests/resources/production_dl2_fits/dl2_mc_events_file.fits \
         --file_type "point-like"    \
-        --metrics_file tests/resources/production_simulation_config_metrics.yaml --site North
+        --metrics_file tests/resources/production_simulation_config_metrics.yml --site North
 
-The output will show the configured simulation parameters.
+The output will show the derived simulation parameters.
 """
 
 import json
@@ -51,6 +64,7 @@ import astropy.units as u
 
 import simtools.utils.general as gen
 from simtools.configuration import configurator
+from simtools.data_model import schema
 from simtools.io_operations import io_handler
 from simtools.production_configuration.generate_simulation_config import (
     SimulationConfig,
@@ -60,7 +74,8 @@ from simtools.production_configuration.generate_simulation_config import (
 def _parse(label):
     """Parse command-line arguments."""
     config = configurator.Configurator(
-        label=label, description="Configure and run a simulation based on input parameters."
+        label=label,
+        description="Derive simulation configuration parameters for a simulation production.",
     )
     config.parser.add_argument(
         "--azimuth", type=float, required=True, help="Azimuth angle in degrees."
@@ -128,6 +143,9 @@ def main():
     metrics = (
         gen.collect_data_from_file(args_dict["metrics_file"]) if "metrics_file" in args_dict else {}
     )
+    schema.validate_dict_using_schema(
+        data=metrics, schema_file="production_configuration_metrics.schema.yml"
+    )
 
     simulation_config = SimulationConfig(
         grid_point=grid_point_config,

simtools/applications/production_scale_events.py

@@ -93,8 +93,8 @@ def _parse(label, description):
     config.parser.add_argument(
         "--metrics_file",
         type=str,
-        default="production_simulation_config_metrics.yaml",
-        help="Metrics definition file. (default: production_simulation_config_metrics.yaml)",
+        default="production_simulation_config_metrics.yml",
+        help="Metrics definition file. (default: production_simulation_config_metrics.yml)",
     )
     config.parser.add_argument(
         "--science_case", type=str, required=True, help="Science case for the simulation."

simtools/applications/simulate_prod.py

@@ -31,7 +31,7 @@ r"""
         Zenith angle in degrees.
     nshow (int, optional)
         Number of showers to simulate.
-        The Number of simulated events depends on the number of times a shower is re-used in the
+        The Number of simulated events depends on the number of times a shower is reused in the
         telescope simulation. The number provided here is before any reuse factors.
     start_run (int, required)
         Start run number such that the actual run number will be 'start_run' + 'run'.

simtools/applications/simulate_prod_htcondor_generator.py

@@ -1,40 +1,40 @@
 #!/usr/bin/python3
 
 r"""
-    Generate a run script and submit file for HT Condor job submission of a simulation production.
+Generate a run script and submit file for HT Condor job submission of a simulation production.
 
-    This tool facilitates the submission of multiple simulations to the HT Condor batch system,
-    enabling:
+This tool facilitates the submission of multiple simulations to the HT Condor batch system,
+enabling:
 
-    - Execution of simulations using the "simtools-simulate-prod" application.
-    - 'number_of_runs' jobs are submitted to the HT Condor batch system.
-    - Utilization of an Apptainer image containing the SimPipe simulation software and tools.
-    - Packaging of data and histogram files, and writing them to a specified directory.
+- Execution of simulations using the "simtools-simulate-prod" application.
+- 'number_of_runs' jobs are submitted to the HT Condor batch system.
+- Utilization of an Apptainer image containing the SimPipe simulation software and tools.
+- Packaging of data and histogram files, and writing them to a specified directory.
 
-    This tool is intended for use in an HT Condor environment. Jobs run in a container universe
-    using the Apptainer image specified in the command line ('--apptainer_image'). Output is written
-    to the 'output_path' directory, with 'simtools-output' and 'logs' subdirectories.
+This tool is intended for use in an HT Condor environment. Jobs run in a container universe
+using the Apptainer image specified in the command line ('--apptainer_image'). Output is written
+to the 'output_path' directory, with 'simtools-output' and 'logs' subdirectories.
 
-    Requirements for the 'simtools-simulate-prod-htcondor-generator' application:
+Requirements for the 'simtools-simulate-prod-htcondor-generator' application:
 
-    - Availability of an Apptainer image 'simtools-prod' (obtainable from the package registry on
-      GitHub, e.g., via 'apptainer pull --force docker://ghcr.io/gammasim/simtools-prod:latest').
-    - Environment parameters required to run CORSIKA and sim_telarray, as well as DB access
-      credentials.  These should be listed similarly to a '.env' file and copied to
-      'output_path/env.txt'.  Ensure that the path to the simulation software is correctly set to
-      'SIMTOOLS_SIMTEL_PATH=/workdir/sim_telarray'.
+- Availability of an Apptainer image 'simtools-prod' (obtainable from the package registry on
+  GitHub, e.g., via 'apptainer pull --force docker://ghcr.io/gammasim/simtools-prod:latest').
+- Environment parameters required to run CORSIKA and sim_telarray, as well as DB access
+  credentials.  These should be listed similarly to a '.env' file and copied to
+  'output_path/env.txt'.  Ensure that the path to the simulation software is correctly set to
+  'SIMTOOLS_SIMTEL_PATH=/workdir/sim_telarray'.
 
 
-    Command line arguments
-    ----------------------
-    output_path (str, required)
-        Directory where the output and the simulation data files will be written.
-    apptainer_image (str, optional)
-        Apptainer image to use for the simulation (full path).
-    priority (int, optional)
-        Job priority (default: 1).
+Command line arguments
+----------------------
+output_path (str, required)
+    Directory where the output and the simulation data files will be written.
+apptainer_image (str, optional)
+    Apptainer image to use for the simulation (full path).
+priority (int, optional)
+    Job priority (default: 1).
 
-    (all other command line arguments are identical to those of :ref:`simulate_prod`).
+(all other command line arguments are identical to those of :ref:`simulate_prod`).
 
 """
 

simtools/applications/submit_data_from_external.py

@@ -24,7 +24,7 @@ r"""
         simtools-submit-data-from-external \\
             --input_meta ./tests/resources/MLTdata-preproduction.meta.yml \\
             --input ./tests/resources/MLTdata-preproduction.ecsv \\
-            --schema ./tests/resources/schema_MST_mirror_2f_measurements.yml \\
+            --schema src/simtools/schemas/input/MST_mirror_2f_measurements.schema.yml \\
             --output_file TEST-submit_data_from_external.ecsv
 
     Expected final print-out message:
@@ -83,6 +83,12 @@ def _parse(label, description):
         type=str,
         required=False,
     )
+    config.parser.add_argument(
+        "--ignore_metadata",
+        help="Ignore metadata",
+        action="store_true",
+        required=False,
+    )
     return config.initialize(output=True)
 
 
@@ -95,16 +101,18 @@ def main():  # noqa: D103
     logger = logging.getLogger()
     logger.setLevel(gen.get_log_level_from_user(args_dict["log_level"]))
 
-    _metadata = MetadataCollector(args_dict=args_dict)
+    _metadata = None if args_dict.get("ignore_metadata") else MetadataCollector(args_dict)
 
     data_validator = validate_data.DataValidator(
-        schema_file=_metadata.get_data_model_schema_file_name(),
+        schema_file=(
+            _metadata.get_data_model_schema_file_name() if _metadata else args_dict.get("schema")
+        ),
         data_file=args_dict["input"],
     )
 
     writer.ModelDataWriter.dump(
         args_dict=args_dict,
-        metadata=_metadata.get_top_level_metadata(),
+        metadata=_metadata.get_top_level_metadata() if _metadata else None,
         product_data=data_validator.validate_and_transform(),
     )
 

simtools/applications/submit_model_parameter_from_external.py

@@ -16,8 +16,8 @@ r"""
         instrument name.
     site (str)
         site location.
-    model_version (str)
-        Model version.
+    parameter_version (str)
+        Parameter version.
     input_meta (str, optional)
         input meta data file (yml format)
 
@@ -33,7 +33,7 @@ r"""
             --value 2 \\
             --instrument LSTN-design \\
             --site North \\
-            --model_version 6.0.0 \\
+            --parameter_version 0.1.0 \\
             --input_meta num_gains.metadata.yml
 
 """
@@ -70,7 +70,9 @@ def _parse(label, description):
     )
     config.parser.add_argument("--instrument", type=str, required=True, help="Instrument name")
     config.parser.add_argument("--site", type=str, required=True, help="Site location")
-    config.parser.add_argument("--model_version", type=str, required=True, help="Model version")
+    config.parser.add_argument(
+        "--parameter_version", type=str, required=True, help="Parameter version"
+    )
 
     config.parser.add_argument(
         "--value",
@@ -102,7 +104,7 @@ def main():  # noqa: D103
     logger.setLevel(gen.get_log_level_from_user(args_dict["log_level"]))
 
     output_path = (
-        Path(args_dict["output_path"]) / args_dict["model_version"] / args_dict["instrument"]
+        Path(args_dict["output_path"]) / args_dict["parameter_version"] / args_dict["instrument"]
         if args_dict.get("output_path")
         else None
     )
@@ -110,7 +112,7 @@ def main():  # noqa: D103
         parameter_name=args_dict["parameter"],
         value=args_dict["value"],
         instrument=args_dict["instrument"],
-        model_version=args_dict["model_version"],
+        parameter_version=args_dict["parameter_version"],
         output_file=Path(args_dict["parameter"]).with_suffix(".json"),
         output_path=output_path,
         use_plain_output_path=args_dict.get("use_plain_output_path"),

simtools/applications/validate_camera_efficiency.py

@@ -47,7 +47,7 @@ import logging
 from pathlib import Path
 
 import simtools.utils.general as gen
-from simtools.camera_efficiency import CameraEfficiency
+from simtools.camera.camera_efficiency import CameraEfficiency
 from simtools.configuration import configurator
 
 
@@ -106,7 +106,7 @@ def main():  # noqa: D103
     ce = CameraEfficiency(
         db_config=_db_config,
         simtel_path=args_dict["simtel_path"],
-        label=label,
+        label=args_dict.get("label", label),
         config_data=args_dict,
     )
     ce.simulate()

simtools/applications/validate_file_using_schema.py

@@ -36,12 +36,10 @@ import logging
 import re
 from pathlib import Path
 
-import jsonschema
-
 import simtools.utils.general as gen
 from simtools.configuration import configurator
 from simtools.constants import MODEL_PARAMETER_SCHEMA_PATH
-from simtools.data_model import metadata_collector, metadata_model, validate_data
+from simtools.data_model import metadata_collector, schema, validate_data
 
 
 def _parse(label, description):
@@ -116,7 +114,20 @@ def _get_schema_file_name(args_dict, data_dict=None):
     return schema_file
 
 
-def validate_schema(args_dict, logger):
+def _get_json_file_list(file_directory=None, file_name=None):
+    """Return list of json files in a directory."""
+    file_list = []
+    if file_directory is not None:
+        file_list = list(Path(file_directory).rglob("*.json"))
+        if not file_list:
+            raise FileNotFoundError(f"No files found in {file_directory}")
+    elif file_name is not None:
+        file_list = [file_name]
+
+    return file_list
+
+
+def validate_dict_using_schema(args_dict, logger):
     """
     Validate a schema file (or several files) given in yaml or json format.
 
@@ -124,33 +135,26 @@ def validate_schema(args_dict, logger):
     the metadata section of the data dictionary.
 
     """
-    if args_dict.get("file_directory") is not None:
-        file_list = list(Path(args_dict["file_directory"]).rglob("*.json"))
-    else:
-        file_list = [args_dict["file_name"]]
-    for file_name in file_list:
+    for file_name in _get_json_file_list(
+        args_dict.get("file_directory"), args_dict.get("file_name")
+    ):
         try:
             data = gen.collect_data_from_file(file_name=file_name)
         except FileNotFoundError as exc:
             logger.error(f"Error reading schema file from {file_name}")
             raise exc
-        try:
-            metadata_model.validate_schema(data, _get_schema_file_name(args_dict, data))
-        except jsonschema.exceptions.ValidationError as exc:
-            logger.error(f"Failed validation of file {file_name}")
-            raise exc
+        schema.validate_dict_using_schema(data, _get_schema_file_name(args_dict, data))
         logger.info(f"Successful validation of file {file_name}")
 
 
 def validate_data_files(args_dict, logger):
     """Validate data files."""
-    file_directory = args_dict.get("file_directory")
-    if file_directory is not None:
+    if args_dict.get("file_directory") is not None:
         tmp_args_dict = {}
-        for file_name in Path(file_directory).rglob("*.json"):
+        for file_name in _get_json_file_list(args_dict.get("file_directory")):
             tmp_args_dict["file_name"] = file_name
             parameter_name = re.sub(r"-\d+\.\d+\.\d+", "", file_name.stem)
-            schema_file = MODEL_PARAMETER_SCHEMA_PATH / f"{parameter_name}.schema.yml"
+            schema_file = schema.get_model_parameter_schema_file(f"{parameter_name}")
             tmp_args_dict["schema"] = schema_file
             tmp_args_dict["data_type"] = "model_parameter"
             tmp_args_dict["require_exact_data_type"] = args_dict["require_exact_data_type"]
@@ -192,7 +196,7 @@ def main():  # noqa: D103
     if args_dict["data_type"].lower() == "metadata":
         validate_metadata(args_dict, logger)
     elif args_dict["data_type"].lower() == "schema":
-        validate_schema(args_dict, logger)
+        validate_dict_using_schema(args_dict, logger)
     else:
         validate_data_files(args_dict, logger)
 

simtools/camera/single_photon_electron_spectrum.py

@@ -0,0 +1,168 @@
+"""Single photon electron spectral analysis."""
+
+import logging
+import re
+import subprocess
+import tempfile
+from io import BytesIO
+from pathlib import Path
+
+from astropy.table import Table
+
+import simtools.data_model.model_data_writer as writer
+from simtools.constants import SCHEMA_PATH
+from simtools.data_model import validate_data
+from simtools.data_model.metadata_collector import MetadataCollector
+from simtools.io_operations import io_handler
+
+
+class SinglePhotonElectronSpectrum:
+    """
+    Single photon electron spectral analysis.
+
+    Parameters
+    ----------
+    args_dict: dict
+        Dictionary with input arguments.
+    """
+
+    prompt_column = "frequency (prompt)"
+    prompt_plus_afterpulse_column = "frequency (prompt+afterpulsing)"
+    afterpulse_column = "frequency (afterpulsing)"
+
+    input_schema = SCHEMA_PATH / "input" / "single_pe_spectrum.schema.yml"
+
+    def __init__(self, args_dict):
+        """Initialize SinglePhotonElectronSpectrum class."""
+        self._logger = logging.getLogger(__name__)
+        self._logger.debug("Initialize SinglePhotonElectronSpectrum class.")
+
+        self.args_dict = args_dict
+        # default output is of ecsv format
+        self.args_dict["output_file"] = str(
+            Path(self.args_dict["output_file"]).with_suffix(".ecsv")
+        )
+        self.io_handler = io_handler.IOHandler()
+        self.data = ""  # Single photon electron spectrum data (as string)
+        self.metadata = MetadataCollector(args_dict=self.args_dict)
+
+    def derive_single_pe_spectrum(self):
+        """Derive single photon electron spectrum."""
+        if self.args_dict.get("use_norm_spe"):
+            return self._derive_spectrum_norm_spe()
+
+        raise NotImplementedError(
+            "Derivation of single photon electron spectrum using a simtool is not yet implemented."
+        )
+
+    def write_single_pe_spectrum(self):
+        """
+        Write single photon electron spectrum plus metadata to disk.
+
+        Includes writing in simtel and simtools (ecsv) formats.
+
+        """
+        simtel_file = self.io_handler.get_output_directory() / Path(
+            self.args_dict["output_file"]
+        ).with_suffix(".dat")
+        self._logger.debug(f"norm_spe output file: {simtel_file}")
+        with open(simtel_file, "w", encoding="utf-8") as simtel:
+            simtel.write(self.data)
+
+        cleaned_data = re.sub(r"%%%.+", "", self.data)  # remove norm_spe row metadata
+        table = Table.read(
+            BytesIO(cleaned_data.encode("utf-8")),
+            format="ascii.no_header",
+            comment="#",
+            delimiter="\t",
+        )
+        table.rename_columns(
+            ["col1", "col2", "col3"],
+            ["amplitude", self.prompt_column, self.prompt_plus_afterpulse_column],
+        )
+
+        writer.ModelDataWriter.dump(
+            args_dict=self.args_dict,
+            metadata=self.metadata.top_level_meta,
+            product_data=table,
+            validate_schema_file=None,
+        )
+
+    def _derive_spectrum_norm_spe(self):
+        """
+        Derive single photon electron spectrum using sim_telarray tool 'norm_spe'.
+
+        Returns
+        -------
+        int
+            Return code of the executed command
+
+        Raises
+        ------
+        subprocess.CalledProcessError
+            If the command execution fails.
+        """
+        tmp_input_file = self._get_input_data(
+            input_file=self.args_dict["input_spectrum"],
+            frequency_column=self.prompt_column,
+        )
+        tmp_ap_file = self._get_input_data(
+            input_file=self.args_dict.get("afterpulse_spectrum"),
+            frequency_column=self.afterpulse_column,
+        )
+
+        command = [
+            f"{self.args_dict['simtel_path']}/sim_telarray/bin/norm_spe",
+            "-r",
+            f"{self.args_dict['step_size']},{self.args_dict['max_amplitude']}",
+            tmp_input_file.name,
+        ]
+        if tmp_ap_file:
+            command.insert(1, "-a")
+            command.insert(2, f"{tmp_ap_file.name}")
+
+        self._logger.debug(f"Running norm_spe command: {' '.join(command)}")
+        try:
+            result = subprocess.run(command, capture_output=True, text=True, check=True)
+        except subprocess.CalledProcessError as exc:
+            self._logger.error(f"Error running norm_spe: {exc}")
+            self._logger.error(f"stderr: {exc.stderr}")
+            raise exc
+        finally:
+            for tmp_file in [tmp_input_file, tmp_ap_file]:
+                try:
+                    Path(tmp_file.name).unlink()
+                except (AttributeError, FileNotFoundError):
+                    pass
+
+        self.data = result.stdout
+        return result.returncode
+
+    def _get_input_data(self, input_file, frequency_column):
+        """
+        Return input data for norm_spe command.
+
+        Input data need to be space separated values of the amplitude spectrum.
+        """
+        input_data = ""
+        if not input_file:
+            return None
+        input_file = Path(input_file)
+
+        if input_file.suffix == ".ecsv":
+            data_validator = validate_data.DataValidator(
+                schema_file=self.input_schema, data_file=input_file
+            )
+            table = data_validator.validate_and_transform()
+            input_data = "\n".join(f"{row['amplitude']} {row[frequency_column]}" for row in table)
+        else:
+            with open(input_file, encoding="utf-8") as f:
+                input_data = (
+                    f.read().replace(",", " ")
+                    if frequency_column == self.prompt_column
+                    else f.read()
+                )
+
+        with tempfile.NamedTemporaryFile(delete=False, mode="w", encoding="utf-8") as tmpfile:
+            tmpfile.write(input_data)
+        return tmpfile