gammasimtools 0.6.1__py3-none-any.whl → 0.8.1__py3-none-any.whl

simtools/model/telescope_model.py

@@ -1,44 +1,36 @@
+"""MC model of a telescope."""
+
 import logging
-import shutil
-from copy import copy
-from pydoc import locate
 
 import astropy.io.ascii
 import numpy as np
-from astropy import units as u
 from astropy.table import Table
 
 import simtools.utils.general as gen
-from simtools import db_handler
-from simtools.io_operations import io_handler
 from simtools.model.camera import Camera
 from simtools.model.mirrors import Mirrors
-from simtools.simtel.simtel_config_writer import SimtelConfigWriter
+from simtools.model.model_parameter import InvalidModelParameterError, ModelParameter
 from simtools.utils import names
 
-__all__ = ["InvalidParameter", "TelescopeModel"]
-
-
-class InvalidParameter(Exception):
-    """Exception for invalid parameter."""
+__all__ = ["TelescopeModel"]
 
 
-class TelescopeModel:
+class TelescopeModel(ModelParameter):
     """
-    TelescopeModel represents the MC model of an individual telescope. It contains the list of \
-    parameters that can be read from the DB. A set of methods are available to manipulate \
-    parameters (changing, adding, removing etc).
+    TelescopeModel represents the MC model of an individual telescope.
+
+    TelescopeModel contains parameter names and values for a specific telescope model.
 
     Parameters
     ----------
     site: str
-        South or North.
-    telescope_model_name: str
-        Telescope name (ex. LST-1, ...).
+        Site name (e.g., South or North).
+    telescope_name: str
+        Telescope name (ex. LSTN-01, LSTN-design, ...).
     mongo_db_config: dict
         MongoDB configuration.
     model_version: str
-        Version of the model (ex. prod5).
+        Model version.
     label: str
         Instance label. Important for output file naming.
     """
@@ -46,561 +38,42 @@ class TelescopeModel:
     def __init__(
         self,
         site,
-        telescope_model_name,
-        mongo_db_config=None,
-        model_version="Released",
-        db=None,
+        telescope_name,
+        mongo_db_config,
+        model_version,
         label=None,
     ):
-        """
-        Initialize TelescopeModel.
-        """
+        """Initialize TelescopeModel."""
         self._logger = logging.getLogger(__name__)
-        self._logger.debug("Init TelescopeModel")
-
-        self.site = names.validate_site_name(site)
-        self.name = names.validate_telescope_model_name(telescope_model_name)
-        self.model_version = names.validate_model_version_name(model_version)
-        self.label = label
-        self._extra_label = None
-        self._added_parameter_files = None
+        self._logger.debug("Init TelescopeModel %s %s", site, telescope_name)
+        ModelParameter.__init__(
+            self,
+            site=site,
+            array_element_name=telescope_name,
+            mongo_db_config=mongo_db_config,
+            model_version=model_version,
+            db=None,
+            label=label,
+        )
+
         self._single_mirror_list_file_paths = None
-        self.simtel_config_writer = None
         self._mirrors = None
-        self._reference_data = None
-        self._derived = None
         self._camera = None
 
-        self.io_handler = io_handler.IOHandler()
-        self.db = None
-        if db is not None:
-            self.db = db
-        elif mongo_db_config is not None:
-            self._logger.debug("Connecting to DB")
-            self.db = db_handler.DatabaseHandler(mongo_db_config=mongo_db_config)
-
-        self._parameters = {}
-
-        self._load_parameters_from_db()
-
-        self._set_config_file_directory_and_name()
-        self._is_config_file_up_to_date = False
-        self._is_exported_model_files_up_to_date = False
-
     @property
     def mirrors(self):
-        """
-        Load the mirror information if the class instance hasn't done it yet.
-        """
+        """Load the mirror information if the class instance hasn't done it yet."""
         if self._mirrors is None:
             self._load_mirrors()
         return self._mirrors
 
     @property
     def camera(self):
-        """
-        Load the camera information if the class instance hasn't done it yet.
-        """
+        """Load the camera information if the class instance hasn't done it yet."""
         if self._camera is None:
             self._load_camera()
         return self._camera
 
-    @property
-    def reference_data(self):
-        """
-        Load the reference data information if the class instance hasn't done it yet.
-        """
-        if self._reference_data is None:
-            self._load_reference_data()
-        return self._reference_data
-
-    @property
-    def derived(self):
-        """
-        Load the derived values and export them if the class instance hasn't done it yet.
-        """
-        if self._derived is None:
-            self._load_derived_values()
-            self.export_derived_files()
-        return self._derived
-
-    @property
-    def extra_label(self):
-        """
-        Return the extra label if defined, if not return ''.
-        """
-        return self._extra_label if self._extra_label is not None else ""
-
-    @classmethod
-    def from_config_file(cls, config_file_name, site, telescope_model_name, label=None):
-        """
-        Create a TelescopeModel from a sim_telarray config file.
-
-        Notes
-        -----
-        This function does not deal with ifdef/indef etc., it just keeps the last version
-        of the parameters in the file. This is fine for now since we do not
-        expect to read from sim_telarray config files in the future.
-
-        Parameters
-        ----------
-        config_file_name: str or Path
-            Path to the input config file.
-        site: str
-            South or North.
-        telescope_model_name: str
-            Telescope model name for the base set of parameters (ex. LST-1, ...).
-        label: str
-            Instance label. Important for output file naming.
-
-        Returns
-        -------
-        TelescopeModel
-            Instance of TelescopeModel.
-        """
-        parameters = {}
-        tel = cls(
-            site=site,
-            telescope_model_name=telescope_model_name,
-            mongo_db_config=None,
-            label=label,
-        )
-
-        def _process_line(words):
-            """
-            Process a line of the input config file that contains a parameter.
-
-            Parameters
-            ----------
-            words: list of str
-                List of str from the split of a line from the file.
-
-            Returns
-            -------
-            (par_name, par_value)
-            """
-            i_comment = len(words)  # Index of any possible comment
-            for w in words:
-                if "%" in w:
-                    i_comment = words.index(w)
-                    break
-            words = words[0:i_comment]  # Removing comment
-            par_name = words[0].replace("=", "")
-            par_value = ""
-            for w in words[1:]:
-                w = w.replace("=", "")
-                w = w.replace(",", " ")
-                par_value += w + " "
-            par_value = par_value.rstrip().lstrip()  # Removing trailing spaces (left and right)
-            return par_name, par_value
-
-        with open(config_file_name, "r", encoding="utf-8") as file:
-            for line in file:
-                words = line.split()
-                if len(words) == 0:
-                    continue
-                if "%" in words[0] or "echo" in words:
-                    continue
-                if "#" not in line and len(words) > 0:
-                    par, value = _process_line(words)
-                    parameters[par] = value
-
-        for par, value in parameters.items():
-            tel.add_parameter(par, value)
-
-        tel._is_exported_model_files_up_to_date = True
-        return tel
-
-    def set_extra_label(self, extra_label):
-        """
-        Set an extra label for the name of the config file.
-
-        Notes
-        -----
-        The config file directory name is not affected by the extra label. Only the file name is \
-        changed. This is important for the ArrayModel class to export multiple config files in the\
-        same directory.
-
-        Parameters
-        ----------
-        extra_label: str
-            Extra label to be appended to the original label.
-        """
-
-        self._extra_label = extra_label
-        self._set_config_file_directory_and_name()
-
-    def _set_config_file_directory_and_name(self):
-        """Define the variable _config_file_directory and create directories, if needed."""
-
-        self._config_file_directory = self.io_handler.get_output_directory(
-            label=self.label, sub_dir="model"
-        )
-
-        # Setting file name and the location
-        config_file_name = names.simtel_telescope_config_file_name(
-            self.site, self.name, self.model_version, self.label, self._extra_label
-        )
-        self._config_file_path = self._config_file_directory.joinpath(config_file_name)
-
-    def _load_parameters_from_db(self):
-        """Read parameters from DB and store them in _parameters."""
-
-        if self.db is None:
-            return
-
-        self._logger.debug("Reading telescope parameters from DB")
-
-        self._set_config_file_directory_and_name()
-        self._parameters = self.db.get_model_parameters(
-            self.site, self.name, self.model_version, only_applicable=True
-        )
-
-        self._logger.debug("Reading site parameters from DB")
-        _site_pars = self.db.get_site_parameters(
-            self.site, self.model_version, only_applicable=True
-        )
-        self._parameters.update(_site_pars)
-
-    def has_parameter(self, par_name):
-        """
-        Verify if the parameter is in the model.
-
-        Parameters
-        ----------
-        par_name: str
-            Name of the parameter.
-
-        Returns
-        -------
-        bool
-            True if parameter is in the model.
-        """
-        return par_name in self._parameters
-
-    def get_parameter(self, par_name):
-        """
-        Get an existing parameter of the model, including derived parameters.
-
-        Parameters
-        ----------
-        par_name: str
-            Name of the parameter.
-
-        Returns
-        -------
-        Value of the parameter
-
-        Raises
-        ------
-        InvalidParameter
-            If par_name does not match any parameter in this model.
-        """
-        try:
-            return self._parameters[par_name]
-        except KeyError:
-            pass  # search in the derived parameters
-        try:
-            return self.derived[par_name]
-        except KeyError as e:
-            msg = f"Parameter {par_name} was not found in the model"
-            self._logger.error(msg)
-            raise InvalidParameter(msg) from e
-
-    def get_parameter_value(self, par_name):
-        """
-        Get the value of an existing parameter of the model.
-
-        Parameters
-        ----------
-        par_name: str
-            Name of the parameter.
-
-        Returns
-        -------
-        Value of the parameter.
-
-        Raises
-        ------
-        InvalidParameter
-            If par_name does not match any parameter in this model.
-        """
-        par_info = self.get_parameter(par_name)
-        return par_info["Value"]
-
-    def get_parameter_value_with_unit(self, par_name):
-        """
-        Get the value of an existing parameter of the model as an Astropy Quantity with its unit.\
-        If no unit is provided in the model, the value is returned without a unit.
-
-        Parameters
-        ----------
-        par_name: str
-            Name of the parameter.
-
-        Returns
-        -------
-        Astropy quantity with the value of the parameter multiplied by its unit. If no unit is \
-        provided in the model, the value is returned without a unit.
-
-        Raises
-        ------
-        InvalidParameter
-            If par_name does not match any parameter in this model.
-        """
-        par_info = self.get_parameter(par_name)
-        if "units" in par_info:
-            return par_info["Value"] * u.Unit(par_info["units"])
-        return par_info["Value"]
-
-    def add_parameter(self, par_name, value, is_file=False, is_aplicable=True):
-        """
-        Add a new parameters to the model. This function does not modify the DB, it affects only \
-        the current instance.
-
-        Parameters
-        ----------
-        par_name: str
-            Name of the parameter.
-        value:
-            Value of the parameter.
-        is_file: bool
-            Indicates whether the new parameter is a file or not.
-        is_aplicable: bool
-            Indicates whether the new parameter is applicable or not.
-
-        Raises
-        ------
-        InvalidParameter
-            If an existing parameter is tried to be added.
-        """
-        if par_name in self._parameters:
-            msg = f"Parameter {par_name} already in the model, use change_parameter instead"
-            self._logger.error(msg)
-            raise InvalidParameter(msg)
-
-        self._logger.info(f"Adding {par_name}={value} to the model")
-        self._parameters[par_name] = {}
-        self._parameters[par_name]["Value"] = value
-        self._parameters[par_name]["Type"] = type(value)
-        self._parameters[par_name]["Applicable"] = is_aplicable
-        self._parameters[par_name]["File"] = is_file
-
-        self._is_config_file_up_to_date = False
-        if is_file:
-            self._is_exported_model_files_up_to_date = False
-
-    def change_parameter(self, par_name, value):
-        """
-        Change the value of an existing parameter to the model. This function does not modify the \
-        DB, it affects only the current instance.
-
-        Parameters
-        ----------
-        par_name: str
-            Name of the parameter.
-        value:
-            Value of the parameter.
-
-        Raises
-        ------
-        InvalidParameter
-            If the parameter to be changed does not exist in this model.
-        """
-        if par_name not in self._parameters:
-            msg = f"Parameter {par_name} not in the model, use add_parameters instead"
-            self._logger.error(msg)
-            raise InvalidParameter(msg)
-
-        type_of_par_name = locate(self._parameters[par_name]["Type"])
-        if not isinstance(value, type_of_par_name):
-            self._logger.warning(
-                f"The type of the provided value ({value}, {type(value)}) "
-                f"is different from the type of {par_name} "
-                f"({self._parameters[par_name]['Type']}). "
-                f"Attempting to cast to the correct type."
-            )
-            try:
-                value = type_of_par_name(value)
-            except ValueError:
-                self._logger.error(
-                    f"Could not cast {value} to {self._parameters[par_name]['Type']}."
-                )
-                raise
-
-        self._logger.debug(
-            f"Changing parameter {par_name} "
-            f"from {self._parameters[par_name]['Value']} to {value}"
-        )
-        self._parameters[par_name]["Value"] = value
-
-        # In case parameter is a file, the model files will be outdated
-        if self._parameters[par_name]["File"]:
-            self._is_exported_model_files_up_to_date = False
-
-        self._is_config_file_up_to_date = False
-
-    def change_multiple_parameters(self, **kwargs):
-        """
-        Change the value of multiple existing parameters in the model. This function does not \
-        modify the DB, it affects only the current instance.
-
-        Parameters
-        ----------
-        **kwargs
-            Parameters should be passed as parameter_name=value.
-
-        Raises
-        ------
-        InvalidParameter
-            If at least one of the parameters to be changed does not exist in this model.
-        """
-        for par, value in kwargs.items():
-            if par in self._parameters:
-                self.change_parameter(par, value)
-            else:
-                self.add_parameter(par, value)
-
-        self._is_config_file_up_to_date = False
-
-    def remove_parameters(self, *args):
-        """
-        Remove a set of parameters from the model.
-
-        Parameters
-        ----------
-        *args
-            Each parameter to be removed has to be passed as args.
-
-        Raises
-        ------
-        InvalidParameter
-            If at least one of the parameter to be removed is not in this model.
-        """
-        for par in args:
-            if par in self._parameters:
-                self._logger.info(f"Removing parameter {par}")
-                del self._parameters[par]
-            else:
-                msg = f"Could not remove parameter {par} because it does not exist"
-                self._logger.error(msg)
-                raise InvalidParameter(msg)
-        self._is_config_file_up_to_date = False
-
-    def add_parameter_file(self, par_name, file_path):
-        """
-        Add a file to the config file directory.
-
-        Parameters
-        ----------
-        par_name: str
-            Name of the parameter.
-        file_path: str
-            Path of the file to be added to the config file directory.
-        """
-        if self._added_parameter_files is None:
-            self._added_parameter_files = []
-        self._added_parameter_files.append(par_name)
-        shutil.copy(file_path, self._config_file_directory)
-
-    def export_model_files(self):
-        """Exports the model files into the config file directory."""
-
-        # Removing parameter files added manually (which are not in DB)
-        pars_from_db = copy(self._parameters)
-        if self._added_parameter_files is not None:
-            for par in self._added_parameter_files:
-                pars_from_db.pop(par)
-
-        self.db.export_model_files(pars_from_db, self._config_file_directory)
-        self._is_exported_model_files_up_to_date = True
-
-    def print_parameters(self):
-        """Print parameters and their values for debugging purposes."""
-        for par, info in self._parameters.items():
-            print(f"{par} = {info['Value']}")
-
-    def export_config_file(self):
-        """Export the config file used by sim_telarray."""
-
-        # Exporting model file
-        if not self._is_exported_model_files_up_to_date:
-            self.export_model_files()
-
-        # Using SimtelConfigWriter to write the config file.
-        self._load_simtel_config_writer()
-        self.simtel_config_writer.write_telescope_config_file(
-            config_file_path=self._config_file_path, parameters=self._parameters
-        )
-
-    def export_derived_files(self):
-        """Write to disk a file from the derived values DB."""
-
-        for par_now in self.derived.values():
-            if par_now["File"]:
-                self.db.export_file_db(
-                    db_name=self.db.DB_DERIVED_VALUES,
-                    dest=self.io_handler.get_output_directory(label=self.label, sub_dir="derived"),
-                    file_name=par_now["Value"],
-                )
-
-    def get_config_file(self, no_export=False):
-        """
-        Get the path of the config file for sim_telarray. The config file is produced if the file\
-        is not updated.
-
-        Parameters
-        ----------
-        no_export: bool
-            Turn it on if you do not want the file to be exported.
-
-        Returns
-        -------
-        Path
-            Path of the exported config file for sim_telarray.
-        """
-        if not self._is_config_file_up_to_date and not no_export:
-            self.export_config_file()
-        return self._config_file_path
-
-    def get_config_directory(self):
-        """
-        Get the path where all the configuration files for sim_telarray are written to.
-
-        Returns
-        -------
-        Path
-            Path where all the configuration files for sim_telarray are written to.
-        """
-        return self._config_file_directory
-
-    def get_derived_directory(self):
-        """
-        Get the path where all the files with derived values for are written to.
-
-        Returns
-        -------
-        Path
-            Path where all the files with derived values are written to.
-        """
-        return self._config_file_directory.parents[0].joinpath("derived")
-
-    def get_telescope_transmission_parameters(self):
-        """
-        Get tel. transmission pars as a list of floats.
-
-        Returns
-        -------
-        list of floats
-            List of 4 parameters that describe the tel. transmission vs off-axis.
-        """
-
-        telescope_transmission = self.get_parameter_value("telescope_transmission")
-        if isinstance(telescope_transmission, str):
-            return [float(v) for v in self.get_parameter_value("telescope_transmission").split()]
-
-        return [float(telescope_transmission), 0, 0, 0]
-
     def export_single_mirror_list_file(self, mirror_number, set_focal_length_to_zero):
         """
         Export a mirror list file with a single mirror in it.
@@ -621,7 +94,7 @@ class TelescopeModel:
         )
         if self._single_mirror_list_file_paths is None:
             self._single_mirror_list_file_paths = {}
-        self._single_mirror_list_file_paths[mirror_number] = self._config_file_directory.joinpath(
+        self._single_mirror_list_file_paths[mirror_number] = self.config_file_directory.joinpath(
             file_name
         )
 
@@ -655,46 +128,69 @@ class TelescopeModel:
 
     def _load_mirrors(self):
         """Load the attribute mirrors by creating a Mirrors object with the mirror list file."""
-        mirror_list_file_name = self._parameters["mirror_list"]["Value"]
+        mirror_list_file_name = self.get_parameter_value("mirror_list")
+        self._logger.debug(f"Reading mirror list from {mirror_list_file_name}")
         try:
-            mirror_list_file = gen.find_file(mirror_list_file_name, self._config_file_directory)
+            mirror_list_file = gen.find_file(mirror_list_file_name, self.config_file_directory)
         except FileNotFoundError:
             mirror_list_file = gen.find_file(mirror_list_file_name, self.io_handler.model_path)
             self._logger.warning(
                 "Mirror_list_file was not found in the config directory - "
                 "Using the one found in the model_path"
             )
+        except TypeError as exc:
+            raise TypeError("Undefined mirror list") from exc
+
         self._mirrors = Mirrors(mirror_list_file, parameters=self._parameters)
 
-    def _load_reference_data(self):
-        """Load the reference data for this telescope from the DB."""
-        self._logger.debug("Reading reference data from DB")
-        self._reference_data = self.db.get_reference_data(
-            self.site, self.model_version, only_applicable=True
-        )
+    def get_telescope_effective_focal_length(self, unit="m", return_focal_length_if_zero=False):
+        """
+        Return effective focal length.
 
-    def _load_derived_values(self):
-        """Load the derived values for this telescope from the DB."""
-        self._logger.debug("Reading derived data from DB")
-        self._derived = self.db.get_derived_values(
-            self.site,
-            self.name,
-            self.model_version,
-        )
+        The function ensures backwards compatibility with older sim-telarray versions.
+
+        Parameters
+        ----------
+        unit: str
+            Unit of the effective focal length. Default is 'm'.
+        return_focal_length_if_zero: bool
+            If True, return the focal length if the effective focal length is 0.
+
+        Returns
+        -------
+        float:
+            Effective focal length.
+
+        """
+        try:
+            eff_focal_length = self.get_parameter_value_with_unit("effective_focal_length")[0]
+        except TypeError:
+            eff_focal_length = self.get_parameter_value_with_unit("effective_focal_length")
+        try:
+            eff_focal_length = eff_focal_length.to(unit).value
+        except AttributeError:
+            eff_focal_length = 0.0
+        if return_focal_length_if_zero and (eff_focal_length is None or eff_focal_length == 0.0):
+            self._logger.warning("Using focal_length because effective_focal_length is 0")
+            return self.get_parameter_value_with_unit("focal_length").to(unit).value
+        return eff_focal_length
 
     def _load_camera(self):
-        """Loading camera attribute by creating a Camera object with the camera config file."""
+        """Load camera attribute by creating a Camera object with the camera config file."""
         camera_config_file = self.get_parameter_value("camera_config_file")
-        focal_length = self.get_parameter_value("effective_focal_length")
-        if focal_length == 0.0:
-            self._logger.warning("Using focal_length because effective_focal_length is 0.")
-            focal_length = self.get_parameter_value("focal_length")
+        focal_length = self.get_telescope_effective_focal_length("cm", True)
         try:
-            camera_config_file_path = gen.find_file(camera_config_file, self._config_file_directory)
+            camera_config_file_path = gen.find_file(camera_config_file, self.config_file_directory)
+        except TypeError as exc:
+            self._logger.error(
+                f"Camera config file {camera_config_file} or "
+                f"config file directory ({self.config_file_directory}) is None"
+            )
+            raise TypeError from exc
         except FileNotFoundError:
             self._logger.warning(
-                "The camera_config_file was not found in the config directory - "
-                "Using the one found in the model_path"
+                f"Camera config file {camera_config_file} not found in the config directory "
+                f"{self.config_file_directory}. Using the one found in the model_path"
             )
             camera_config_file_path = gen.find_file(camera_config_file, self.io_handler.model_path)
 
@@ -704,19 +200,6 @@ class TelescopeModel:
             focal_length=focal_length,
         )
 
-    def _load_simtel_config_writer(self):
-        """
-        Load the SimtelConfigWriter object.
-
-        """
-        if self.simtel_config_writer is None:
-            self.simtel_config_writer = SimtelConfigWriter(
-                site=self.site,
-                telescope_model_name=self.name,
-                model_version=self.model_version,
-                label=self.label,
-            )
-
     def is_file_2d(self, par):
         """
         Check if the file referenced by par is a 2D table.
@@ -731,21 +214,22 @@ class TelescopeModel:
         bool:
             True if the file is a 2D map type, False otherwise.
         """
-        if not self.has_parameter(par):
+        try:
+            file_name = self.get_parameter_value(par)
+        except KeyError:
             logging.error(f"Parameter {par} does not exist")
             return False
 
-        file_name = self.get_parameter_value(par)
-        file = self.get_config_directory().joinpath(file_name)
-        with open(file, "r", encoding="utf-8") as f:
-            is_2d = "@RPOL@" in f.read()
-        return is_2d
+        file = self.config_file_directory.joinpath(file_name)
+        with open(file, encoding="utf-8") as f:
+            return "@RPOL@" in f.read()
 
     def read_two_dim_wavelength_angle(self, file_name):
         """
-        Read a two dimensional distribution of wavelngth and angle (z-axis can be anything). Return\
-        a dictionary with three arrays, wavelength, angles, z (can be transmission, reflectivity,\
-        etc.)
+        Read a two dimensional distribution of wavelength and angle (z-axis can be anything).
+
+        Return a dictionary with three arrays, wavelength, angles, z (can be transmission,
+        reflectivity, etc.)
 
         Parameters
         ----------
@@ -757,10 +241,10 @@ class TelescopeModel:
         dict:
             dict of three arrays, wavelength, degrees, z.
         """
-
-        _file = self.get_config_directory().joinpath(file_name)
+        _file = self.config_file_directory.joinpath(file_name)
+        self._logger.debug("Reading two dimensional distribution from %s", _file)
         line_to_start_from = 0
-        with open(_file, "r", encoding="utf-8") as f:
+        with open(_file, encoding="utf-8") as f:
             for i_line, line in enumerate(f):
                 if line.startswith("ANGLE"):
                     degrees = np.array(line.strip().split("=")[1].split(), dtype=np.float16)
@@ -777,14 +261,13 @@ class TelescopeModel:
 
     def get_on_axis_eff_optical_area(self):
         """Return the on-axis effective optical area (derived previously for this telescope)."""
-
         ray_tracing_data = astropy.io.ascii.read(
-            self.get_derived_directory().joinpath(self.derived["ray_tracing"]["Value"])
+            self.config_file_directory.joinpath(self.get_parameter_value("optics_properties"))
         )
         if not np.isclose(ray_tracing_data["Off-axis angle"][0], 0):
             self._logger.error(
                 f"No value for the on-axis effective optical area exists."
-                f" The minumum off-axis angle is {ray_tracing_data['Off-axis angle'][0]}"
+                f" The minimum off-axis angle is {ray_tracing_data['Off-axis angle'][0]}"
             )
             raise ValueError
         return ray_tracing_data["eff_area"][0]
@@ -803,17 +286,19 @@ class TelescopeModel:
         incidence_angle_dist: astropy.table.Table
             Instance of astropy.table.Table with the incidence angle distribution.
         """
-
-        incidence_angle_dist = astropy.io.ascii.read(
-            self.get_derived_directory().joinpath(incidence_angle_dist_file)
+        self._logger.debug(
+            "Reading incidence angle distribution from %s",
+            self.config_file_directory.joinpath(incidence_angle_dist_file),
         )
-        return incidence_angle_dist
+        return astropy.io.ascii.read(self.config_file_directory.joinpath(incidence_angle_dist_file))
 
     @staticmethod
     def calc_average_curve(curves, incidence_angle_dist):
         """
-        Calculate an average curve from a set of curves, using as weights the distribution of \
-        incidence angles provided in incidence_angle_dist.
+        Calculate an average curve from a set of curves.
+
+        The calculation uses weights the distribution of incidence angles provided in
+        incidence_angle_dist.
 
         Parameters
         ----------
@@ -830,7 +315,6 @@ class TelescopeModel:
         average_curve: astropy.table.Table
             Instance of astropy.table.Table with the averaged curve.
         """
-
         weights = []
         for angle_now in curves["Angle"]:
             weights.append(
@@ -839,13 +323,11 @@ class TelescopeModel:
                 ]
             )
 
-        average_curve = Table(
+        return Table(
             [curves["Wavelength"], np.average(curves["z"], weights=weights, axis=0)],
             names=("Wavelength", "z"),
         )
 
-        return average_curve
-
     def export_table_to_model_directory(self, file_name, table):
         """
         Write out a file with the provided table to the model directory.
@@ -862,7 +344,31 @@ class TelescopeModel:
         Path:
             Path to the file exported.
         """
-
-        file_to_write_to = self._config_file_directory.joinpath(file_name)
+        file_to_write_to = self.config_file_directory.joinpath(file_name)
         table.write(file_to_write_to, format="ascii.commented_header", overwrite=True)
         return file_to_write_to.absolute()
+
+    def position(self, coordinate_system="ground"):
+        """
+        Get coordinates in the given system.
+
+        Parameters
+        ----------
+        coordinate_system: str
+            Coordinates system. Default is 'ground'.
+
+        Returns
+        -------
+        list :
+            List of telescope position in the requested coordinate system.
+
+        Raises
+        ------
+        KeyError
+            If the coordinate system is not found.
+        """
+        try:
+            return self.get_parameter_value_with_unit(f"array_element_position_{coordinate_system}")
+        except InvalidModelParameterError as exc:
+            self._logger.error(f"Coordinate system {coordinate_system} not found.")
+            raise exc