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

simtools/visualization/legend_handlers.py

@@ -3,9 +3,9 @@ import matplotlib.patches as mpatches
 import numpy as np
 
 from simtools.corsika.corsika_config import CorsikaConfig
-from simtools.io_handler import IOHandler
+from simtools.io_operations.io_handler import IOHandler
 from simtools.utils import names
-from simtools.utils.names import lst, mst, sct, sst
+from simtools.utils.names import hess, lst, magic, mst, sct, sst, veritas
 
 __all__ = [
     "EdgePixelObject",
@@ -50,9 +50,12 @@ class TelescopeHandler(object):
             mst: "dodgerblue",
             sct: "black",
             sst: "darkgreen",
+            hess: "grey",
+            magic: "grey",
+            veritas: "grey",
         }
-        for tel_type in names.all_telescope_class_names:
-            self.radius_dict[tel_type] = corsika_info["corsika_sphere_radius"][tel_type]["value"]
+        for key, value in corsika_info["corsika_sphere_radius"].items():
+            self.radius_dict[key] = value["value"]
 
 
 class PixelObject(object):
@@ -83,6 +86,18 @@ class SSTObject(object):
     """SST Object."""
 
 
+class HESSObject(object):
+    """HESS Object."""
+
+
+class MAGICObject(object):
+    """MAGIC Object."""
+
+
+class VERITASObject(object):
+    """VERITAS Object."""
+
+
 class MeanRadiusOuterEdgeObject(object):
     """Object for Mean radius outer edge."""
 
@@ -288,7 +303,7 @@ class SSTHandler(TelescopeHandler):
         return patch
 
 
-class SCTHandler(object):
+class SCTHandler(TelescopeHandler):
     """
     Legend handler class to plot a representation of an SCT in an array layout.
     """
@@ -308,6 +323,78 @@ class SCTHandler(object):
         return patch
 
 
+class HESSHandler(TelescopeHandler):
+    """
+    Legend handler class to plot a representation of an HESS in an array layout.
+    """
+
+    def legend_artist(self, legend, orig_handle, fontsize, handlebox):
+        x0, y0 = (
+            handlebox.xdescent + handlebox.width / 3,
+            handlebox.ydescent + handlebox.height / 3,
+        )
+        radius = handlebox.height
+        patch = mpatches.RegularPolygon(
+            (x0, y0),
+            numVertices=6,
+            radius=0.7 * radius * self.radius_dict[hess] / self.radius_dict[lst],
+            orientation=np.deg2rad(30),
+            facecolor=self.colors_dict[hess],
+            edgecolor=self.colors_dict[hess],
+            transform=handlebox.get_transform(),
+        )
+        handlebox.add_artist(patch)
+        return patch
+
+
+class MAGICHandler(TelescopeHandler):
+    """
+    Legend handler class to plot a representation of an MAGIC in an array layout.
+    """
+
+    def legend_artist(self, legend, orig_handle, fontsize, handlebox):
+        x0, y0 = (
+            handlebox.xdescent + handlebox.width / 3,
+            handlebox.ydescent + handlebox.height / 3,
+        )
+        radius = handlebox.height
+        patch = mpatches.RegularPolygon(
+            (x0, y0),
+            numVertices=6,
+            radius=0.7 * radius * self.radius_dict[magic] / self.radius_dict[lst],
+            orientation=np.deg2rad(30),
+            facecolor=self.colors_dict[magic],
+            edgecolor=self.colors_dict[magic],
+            transform=handlebox.get_transform(),
+        )
+        handlebox.add_artist(patch)
+        return patch
+
+
+class VERITASHandler(TelescopeHandler):
+    """
+    Legend handler class to plot a representation of an VERITAS in an array layout.
+    """
+
+    def legend_artist(self, legend, orig_handle, fontsize, handlebox):
+        x0, y0 = (
+            handlebox.xdescent + handlebox.width / 3,
+            handlebox.ydescent + handlebox.height / 3,
+        )
+        radius = handlebox.height
+        patch = mpatches.RegularPolygon(
+            (x0, y0),
+            numVertices=6,
+            radius=0.7 * radius * self.radius_dict[veritas] / self.radius_dict[lst],
+            orientation=np.deg2rad(30),
+            facecolor=self.colors_dict[veritas],
+            edgecolor=self.colors_dict[veritas],
+            transform=handlebox.get_transform(),
+        )
+        handlebox.add_artist(patch)
+        return patch
+
+
 class MeanRadiusOuterEdgeHandler(object):
     """
     Legend handler class to plot a the mean radius outer edge of the dish.
@@ -331,8 +418,27 @@ class MeanRadiusOuterEdgeHandler(object):
         return patch
 
 
-all_telescope_objects = {lst: LSTObject, mst: MSTObject, sct: SCTObject, sst: SSTObject}
-all_telescope_handlers = {lst: LSTHandler, mst: MSTHandler, sct: SCTHandler, sst: SSTHandler}
+all_telescope_objects = {
+    lst: LSTObject,
+    mst: MSTObject,
+    sct: SCTObject,
+    sst: SSTObject,
+    hess: HESSObject,
+    magic: MAGICObject,
+    veritas: VERITASObject,
+}
+all_telescope_handlers = {
+    lst: LSTHandler,
+    mst: MSTHandler,
+    sct: SCTHandler,
+    sst: SSTHandler,
+    hess: HESSHandler,
+    magic: MAGICHandler,
+    veritas: VERITASHandler,
+}
 legend_handler_map = {}
-for tel_type in names.all_telescope_class_names:
-    legend_handler_map[all_telescope_objects[tel_type]] = all_telescope_handlers[tel_type]
+try:
+    for tel_type in names.all_telescope_class_names:
+        legend_handler_map[all_telescope_objects[tel_type]] = all_telescope_handlers[tel_type]
+except KeyError:
+    pass

simtools/visualization/visualize.py

@@ -12,7 +12,7 @@ from cycler import cycler
 from matplotlib import gridspec
 from matplotlib.collections import PatchCollection
 
-from simtools.utils import general as gen
+from simtools.utils import geometry as transf
 from simtools.utils import names
 from simtools.utils.names import mst, sct
 from simtools.visualization import legend_handlers as leg_h
@@ -22,9 +22,9 @@ __all__ = [
     "get_lines",
     "get_markers",
     "get_telescope_patch",
-    "plot_1D",
+    "plot_1d",
     "plot_array",
-    "plot_hist_2D",
+    "plot_hist_2d",
     "plot_table",
     "set_style",
 ]
@@ -272,7 +272,7 @@ def get_lines():
     return LINES
 
 
-def plot_1D(data, **kwargs):
+def plot_1d(data, **kwargs):
     """
     Produce a high contrast one dimensional plot from multiple data sets. A ratio plot can be \
     added at the bottom to allow easy comparison. Additional options, such as plot title, plot
@@ -481,10 +481,10 @@ def plot_table(table, y_title, **kwargs):
     for column in table.keys()[1:]:
         data_dict[column] = QTable([table[x_axis], table[column]], names=[x_axis, y_title])
 
-    return plot_1D(data_dict, **kwargs)
+    return plot_1d(data_dict, **kwargs)
 
 
-def plot_hist_2D(data, **kwargs):
+def plot_hist_2d(data, **kwargs):
     """
     Produce a two dimensional histogram plot. Any option that can be changed after plotting (e.g.,\
     axes limits, log scale, etc.) should be done using the returned plt instance.
@@ -511,7 +511,7 @@ def plot_hist_2D(data, **kwargs):
     else:
         title = ""
 
-    # Set default style since the usual options do not affect 2D plots (for now).
+    # Set default style since the usual options do not affect 2d plots (for now).
     set_style()
 
     gs = gridspec.GridSpec(1, 1)
@@ -563,7 +563,7 @@ def get_telescope_patch(name, x, y, radius):
         Instance of mpatches.Circle.
     """
     tel_obj = leg_h.TelescopeHandler()
-    valid_name = names.get_telescope_type(name)
+    valid_name = names.get_telescope_class(name)
     fill_flag = False
 
     x = x.to(u.m)
@@ -627,13 +627,13 @@ def plot_array(telescopes, rotate_angle=0, show_tel_label=False):
     legend_labels = []
     tel_counters = {one_telescope: 0 for one_telescope in names.all_telescope_class_names}
     if rotate_angle != 0:
-        pos_x_rotated, pos_y_rotated = gen.rotate(
-            telescopes["pos_x"], telescopes["pos_y"], rotate_angle
+        pos_x_rotated, pos_y_rotated = transf.rotate(
+            telescopes["position_x"], telescopes["position_y"], rotate_angle
         )
     else:
-        pos_x_rotated, pos_y_rotated = telescopes["pos_x"], telescopes["pos_y"]
+        pos_x_rotated, pos_y_rotated = telescopes["position_x"], telescopes["position_y"]
 
-    pos_x_rotated, pos_y_rotated = gen.rotate(pos_x_rotated, pos_y_rotated, 90 * u.deg)
+    pos_x_rotated, pos_y_rotated = transf.rotate(pos_x_rotated, pos_y_rotated, 90 * u.deg)
 
     if len(pos_x_rotated) > 30:
         fontsize = 4
@@ -647,7 +647,7 @@ def plot_array(telescopes, rotate_angle=0, show_tel_label=False):
         for tel_type in tel_counters:
             if tel_type in tel_now["telescope_name"]:
                 tel_counters[tel_type] += 1
-        i_tel_name = names.get_telescope_type(telescopes[i_tel]["telescope_name"])
+        i_tel_name = names.get_telescope_class(telescopes[i_tel]["telescope_name"])
         patches.append(
             get_telescope_patch(
                 i_tel_name,

gammasimtools-0.5.1.dist-info/RECORD

@@ -1,83 +0,0 @@
-simtools/__init__.py,sha256=yvx0i5ygbllmMtF2Tzy7pohvhQNH_dBdj0Bt8bKFYpc,231
-simtools/_version.py,sha256=d4lWkB7Jtv-a_y1raRjfA2309LsV73tQCinVV-F9tv8,160
-simtools/camera_efficiency.py,sha256=Kpj0XOZ7bWUqR1TQY_7q0uHry_hrRicgsr8CgvIE3S4,17879
-simtools/db_handler.py,sha256=VBhYfoIx0VV-3z6dPm5o0_3vhH8pJwXtYjJJ5TEoSTo,46955
-simtools/io_handler.py,sha256=2fuOh119SdV74jHGiYHrsd_luhASYCq0LkFKpFzNWcw,4207
-simtools/psf_analysis.py,sha256=7w_fVKPGhut8ofp_RaJHwYiNh8FoRqOqKB93pSEWXTU,14475
-simtools/ray_tracing.py,sha256=ce52aoNxXAo5pcz8T4sP6Mx-_wRc8UkOtGxrlxbYv9Q,17199
-simtools/simulator.py,sha256=r2N9hWDCsrPbvvgO-83goWKemp_BlsZeb1nu3qjKxDM,26288
-simtools/version.py,sha256=QH3JDagS1ROImltyE_xRzRZlDxIr90poVpXToqNQUFI,595
-simtools/_dev_version/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-simtools/_dev_version/scm_version.py,sha256=o3HfXziUHq0iIHUQm1RBv_8HDEpckXHFl0kYgmmwdAA,375
-simtools/applications/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-simtools/applications/add_file_to_db.py,sha256=rbZ2MLynqXyWh9fqKquL9vkqnz-Ul3ImPIrr2ZiqOjE,4997
-simtools/applications/compare_cumulative_psf.py,sha256=yiFd2oldOiiz41Qh5fqs932cRJGB8GMebLj0sE3ZYPs,6772
-simtools/applications/derive_mirror_rnda.py,sha256=ocs1AoGQKBZ1CKaIkEl-N9jY0ErH_kIwtkIhtQU7ev8,15104
-simtools/applications/generate_corsika_histograms.py,sha256=6DaFUP2QIoKhnZgGj9ft-qkABk6EH2A4kHyZMpwKTIA,15189
-simtools/applications/get_file_from_db.py,sha256=dr66mSmbv0StsZph-lisPmW3Fdccv-rEWGPJTuC5JUY,2921
-simtools/applications/get_parameter.py,sha256=bmTbl9NAIEgXsWXYyYRNcozygTSWsiBphCNN0C_BSaE,2880
-simtools/applications/make_regular_arrays.py,sha256=VAEVoAfAl147pJe0yceQ4FQk8a8FMIF2VRl2cSSWWjc,5759
-simtools/applications/plot_layout_array.py,sha256=K-NJZv5nMl7lDSsF1RjD3Hap-hIQJb4siLZ-cqFzBwc,6326
-simtools/applications/plot_simtel_histograms.py,sha256=-7zcWYbW6RPJoDOGKO1ZMi8u3f1Ad_oanEORDkgn2kQ,3746
-simtools/applications/print_array_elements.py,sha256=1ecPKGtZVxWQMhx2ZHAwpGCMyhWi9b2B2sX7XF4fYSs,3989
-simtools/applications/produce_array_config.py,sha256=AnioMEsW7TZ6RxmPgsa1YOW1DAWmdWdnmPLHJv-Yl78,4142
-simtools/applications/production.py,sha256=gcMeOGQ3D3QaTAxmbwUW7tmb9pEpba9WzJxSDUrcLvQ,9172
-simtools/applications/sim_showers_for_trigger_rates.py,sha256=aQfs3fd5prGgFLeZ5Y8tLM2h5HphU2z5w6sC6C21274,6638
-simtools/applications/simulate_prod.py,sha256=CCkk_Pxk_sJorvGwNNdi603mNGcp4IBNJOP5OGDsOh8,13012
-simtools/applications/submit_data_from_external.py,sha256=FU6hPN5VyL8Ji-r2G0R8qfLepRZQR8ulqnO1U8d7mqA,3134
-simtools/applications/tune_psf.py,sha256=Ob6aH0OKrMYazcps8Gj46-gkBF8ciT7V3zDh7-rTt0s,11575
-simtools/applications/validate_camera_efficiency.py,sha256=jcD8-m2oMPxVEGzyhiKKE7P2nqiRS2tbkbUO68Vqy-I,4188
-simtools/applications/validate_camera_fov.py,sha256=H4gnj_nwFeYe5zJ1EJGdApmW6q2HuHHZ40V6_EDZ-bI,4899
-simtools/applications/validate_optics.py,sha256=MrcKAEPBfGlsh8ealE2zar0jRIqXedcNxG1vRUaGR5M,6617
-simtools/applications/validate_schema_files.py,sha256=HgpqqQVVMpoSVJoHA1sct_CNx3BgmKP7f-YNMHohuXQ,3013
-simtools/applications/db_development_tools/add_new_parameter_to_db.py,sha256=K4rdaXii-v7hOPveF6dh8zb5tZCPu6AzPeIkPS8NHJE,2565
-simtools/applications/db_development_tools/add_unit_to_parameter_in_db.py,sha256=eqcV_yFyaOl8GOEtQozY09J_IKFy9OtRNVekIcuX0m4,1809
-simtools/applications/db_development_tools/mark_non_optics_parameters_non_applicable.py,sha256=-0ImMFRyHmhB9XGA9pgXHe2laqQo1EoL9Zy13e92Ogs,3001
-simtools/configuration/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-simtools/configuration/commandline_parser.py,sha256=ofgOTFGm7tPf-xtg_LHxIGVcMxyMXPSwMDGzbbr6Lec,12519
-simtools/configuration/configurator.py,sha256=8VO6J8nw3NxIISjPgF4H_ajrFTAZqbFjrDeCra_3D3w,12160
-simtools/corsika/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-simtools/corsika/corsika_config.py,sha256=u-X5lu_GdkdrcaVxMEtAiFf2IuuIzBXtenIcp4cl4f0,19175
-simtools/corsika/corsika_histograms.py,sha256=dZVqVvx3uIkxcji6djL543DSM7eAo3aGU_iveMrUOKQ,65160
-simtools/corsika/corsika_histograms_visualize.py,sha256=9pQqy0iw5l9gEC46KsGGSNtUFGWNvqflmZgwNjLMa1I,18733
-simtools/corsika/corsika_output_visualize.py,sha256=Hi-K58pxMQ_DSII-35Ow6etx3cBo0SUcDiRTXJgi6LE,12080
-simtools/corsika/corsika_runner.py,sha256=_qypZytK4_WWw_d052hj-6twhPKwR_BG93FIG2cLS3I,16871
-simtools/corsika_simtel/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-simtools/corsika_simtel/corsika_simtel_runner.py,sha256=W_GXKulDpgXytaxjvmBtrJ8vTwA5Y9QLug75s-a7UtU,7804
-simtools/data_model/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-simtools/data_model/metadata_collector.py,sha256=oMjOVkVR7wv2NnvAc8Z2g4t9Gkz5qoAWW5Q_Ffi05aI,10245
-simtools/data_model/metadata_model.py,sha256=Po2k734pm9nMtqVRYaCLEnI73TKZHja-MFPZ9I0frc8,8151
-simtools/data_model/model_data_writer.py,sha256=p-yYgkQ8E9QxhSmnb7MWCMvE1GFwM4snTncNMX4mBhg,3882
-simtools/data_model/validate_data.py,sha256=4t2P6XaV3i3le0itsggfllac7TVLNPKN65_1SzVyt4A,16980
-simtools/data_model/validate_schema.py,sha256=Fqt3lpmwOulSjvDYlFNJIf-wXIG_BoExkBu7hxgJ8kE,8542
-simtools/job_execution/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-simtools/job_execution/job_manager.py,sha256=zj9cJLzY-I-_gNp42ClSBjoOzaczH8CYfdgJ5CCLoew,6048
-simtools/layout/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-simtools/layout/layout_array.py,sha256=pCSjENXIGNLTaz5qKGAkq7m8knuC4LDZYZYvLCgW0Z0,31305
-simtools/layout/telescope_position.py,sha256=_IMfwpKgWlp7ruQkqr_MZtETZMG6nTomAEwtXhS6rS4,17361
-simtools/model/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-simtools/model/array_model.py,sha256=qCl3SH07b1gtRQlrQU4_zKBfKwttdK7Wh4cDDhINLUQ,14179
-simtools/model/camera.py,sha256=ZOJhfMcRY5J8SQNZVAVFro4P7arHGzJkOGE1wi1Ic58,30775
-simtools/model/mirrors.py,sha256=jU13bG7kp5yQrxp6_hT5bVFVh0TbyCOAA4RtcjH32pw,5266
-simtools/model/model_utils.py,sha256=YrKl96I3pb3WzNkQLlejHMq5w-XwFMvfu5xlsED12P0,3748
-simtools/model/telescope_model.py,sha256=_yDh7XGfRsTd6EHBNOcWR9ZmR2lJlAXs-CMUxpth6rM,29204
-simtools/simtel/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-simtools/simtel/simtel_config_writer.py,sha256=DYVEt6-pZtFRgvbbh1YmEVz3wwuMMncptjQh07NbqV4,7903
-simtools/simtel/simtel_events.py,sha256=YzIvZ0p96LclPBCGfw1zzRGIMjgVnh4xIurzYokM8q4,9091
-simtools/simtel/simtel_histograms.py,sha256=0SDlmKnL_kQN-ETQC7-L3pCIxPRorxN0FOUZs39rKR4,7262
-simtools/simtel/simtel_runner.py,sha256=W16jA7epT4LXIDVoz76tBaHNln5haToY_Art7NscC_g,7489
-simtools/simtel/simtel_runner_array.py,sha256=f_AQvE4YwIg1j-RVjvPPqWg8-yofKgje8IUqHz6R_ww,10426
-simtools/simtel/simtel_runner_camera_efficiency.py,sha256=zUggzjYxcqjJBBjDJDH8DqtVvSWDhq9731jWEdSzm78,8244
-simtools/simtel/simtel_runner_ray_tracing.py,sha256=fkLyfxVeyEv49xGQWASsWmO110IAGGEA3lNbs3hs4Kg,10031
-simtools/utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-simtools/utils/general.py,sha256=bfNbGVpxpcNBXeWS00AkIB6NkcwESKe0VAC4CiuLnTs,28736
-simtools/utils/names.py,sha256=qiqgpQH3tgQa_LSnPU0tlqnoGO9DeDqds9odB_Imj2U,18192
-simtools/visualization/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-simtools/visualization/legend_handlers.py,sha256=ItWeTvwrXi7-ZMv3ZEyTdO832zFMPB7IGYDuK90LujE,9737
-simtools/visualization/visualize.py,sha256=9W-7l1uTtuPw_CQxrMvRnvyJWvZXXIs1PIR_G2JGczc,21783
-gammasimtools-0.5.1.dist-info/LICENSE,sha256=idBrtxX6eQcUlG0YBTx-O2RNKNKg1Gidb8OQf9q0X7s,1543
-gammasimtools-0.5.1.dist-info/METADATA,sha256=EqlVxeU00ZpLjrQ54ohGQP5KOXbt9BCcQJIh8ANO6Nk,6090
-gammasimtools-0.5.1.dist-info/WHEEL,sha256=5sUXSg9e4bi7lTLOHcm6QEYwO5TIF1TNbTSVFVjcJcc,92
-gammasimtools-0.5.1.dist-info/entry_points.txt,sha256=IXTMPrGSRbZg5Hy9Uqx09L_iPemB2L2uWyU5joAcAE4,1566
-gammasimtools-0.5.1.dist-info/top_level.txt,sha256=9Yi7Q7aQX6alvpP2YGQtIYo208wdGkXUApHboEdj-ys,9
-gammasimtools-0.5.1.dist-info/RECORD,,

simtools/applications/plot_simtel_histograms.py

@@ -1,120 +0,0 @@
-#!/usr/bin/python3
-
-"""
-    Summary
-    -------
-    This application plots and prints sim_telarray histograms into pdf file.
-    It accepts multiple lists of histograms files that are plot together,
-    side-by-side. Each histogram is plotted in a page of the pdf.
-
-
-    Command line arguments
-    ----------------------
-    file_lists (str, required)
-        Text file containing the list of sim_telarray histogram files to be plotted. \
-        Multiple text files can be given.
-    figure_name (str, required)
-        File name for the pdf output (without extension).
-    verbosity (str, optional)
-        Log level to print (default=INFO).
-
-    Example
-    -------
-    .. code-block:: console
-
-        simtools-plot-simtel-histograms \
-            --file_lists list_test1.txt list_test2.txt --figure_name histograms_comparison
-
-"""
-
-import logging
-from pathlib import Path
-
-import matplotlib.pyplot as plt
-from matplotlib.backends.backend_pdf import PdfPages
-
-import simtools.utils.general as gen
-from simtools.configuration import configurator
-from simtools.simtel.simtel_histograms import SimtelHistograms
-
-
-def main():
-    config = configurator.Configurator(
-        label=Path(__file__).stem,
-        description=("Plots sim_telarray histograms."),
-    )
-    config.parser.add_argument(
-        "--file_lists",
-        help="File containing the list of histogram files to be plotted.",
-        nargs="+",
-        type=str,
-        required=True,
-    )
-    config.parser.add_argument(
-        "--figure_name", help="File name for the pdf output.", type=str, required=True
-    )
-
-    args_dict, _ = config.initialize()
-
-    logger = logging.getLogger()
-    logger.setLevel(gen.get_log_level_from_user(args_dict["log_level"]))
-
-    n_lists = len(args_dict["file_lists"])
-    simtel_histograms = []
-    for this_list_of_files in args_dict["file_lists"]:
-        # Collecting hist files
-        histogram_files = []
-        with open(this_list_of_files, encoding="utf-8") as file:
-            for line in file:
-                # Removing '\n' from filename, in case it is left there.
-                histogram_files.append(line.replace("\n", ""))
-
-        # Building SimtelHistograms
-        sh = SimtelHistograms(histogram_files)
-        simtel_histograms.append(sh)
-
-    # Checking if number of histograms is consistent
-    number_of_hists = [sh.number_of_histograms for sh in simtel_histograms]
-    # Converting list to set will remove the duplicated entries.
-    # If all entries in the list are the same, len(set) will be 1
-    if len(set(number_of_hists)) > 1:
-        msg = (
-            "Number of histograms in different sets of simulations is inconsistent"
-            " - please make sure the simulations sets are consistent"
-        )
-        logger.error(msg)
-        raise IOError(msg)
-
-    # Plotting
-
-    # Checking if it is needed to add the pdf extension to the file name
-    if args_dict["figure_name"].split(".")[-1] == "pdf":
-        fig_name = args_dict["figure_name"]
-    else:
-        fig_name = args_dict["figure_name"] + ".pdf"
-
-    pdf_pages = PdfPages(fig_name)
-    for i_hist in range(number_of_hists[0]):
-        title = simtel_histograms[0].get_histogram_title(i_hist)
-
-        logger.debug(f"Processing: {title}")
-
-        fig, axs = plt.subplots(1, n_lists, figsize=(6 * n_lists, 6))
-
-        if n_lists == 1:
-            # If only one simulation set, axs is a single axes (not iterable)
-            sh.plot_one_histogram(i_hist, axs)
-        else:
-            for sh, ax in zip(simtel_histograms, axs):
-                sh.plot_one_histogram(i_hist, ax)
-
-        plt.tight_layout()
-        pdf_pages.savefig(fig)
-        plt.clf()
-
-    plt.close()
-    pdf_pages.close()
-
-
-if __name__ == "__main__":
-    main()

simtools/applications/validate_schema_files.py

@@ -1,135 +0,0 @@
-#!/usr/bin/python3
-"""
-    Summary
-    -------
-    Validate parameter description files using a json schema file.
-
-    Command line arguments
-    ----------------------
-    file_name (str)
-      input file to be validated
-    schema (str)
-      schema file (jsonschema format) used for validation
-
-    Example
-    -------
-
-    .. code-block:: console
-
-        simtools-validate-schema-file \
-         --file_name tests/resources/MST_mirror_2f_measurements.schema.yml \
-         --schema jsonschema.yml
-
-"""
-
-import logging
-from pathlib import Path
-
-import jsonschema
-import yaml
-
-import simtools.utils.general as gen
-from simtools.configuration import configurator
-
-
-def _parse(label, description):
-    """
-    Parse command line configuration
-
-    Parameters
-    ----------
-    label (str)
-        application label
-    description (str)
-        application description
-
-    Returns
-    -------
-    config (Configurator)
-        application configuration
-
-    """
-
-    config = configurator.Configurator(label=label, description=description)
-    config.parser.add_argument("--file_name", help="file to be validated", required=True)
-    config.parser.add_argument("--schema", help="json schema file", required=True)
-    return config.initialize()
-
-
-def load_schema(schema_file):
-    """
-    Load parameter schema from file.
-
-    Parameters
-    ----------
-    schema_file (str)
-        schema file
-
-    Returns
-    -------
-    parameter_schema (dict)
-        parameter schema
-
-    Raises
-    ------
-    FileNotFoundError
-        if schema file is not found
-
-    """
-
-    try:
-        with open(schema_file, "r", encoding="utf-8") as file:
-            parameter_schema = yaml.safe_load(file)
-    except FileNotFoundError:
-        logging.error(f"Schema file {schema_file} not found")
-        raise
-
-    return parameter_schema
-
-
-def validate_schema_file(input_file, schema, logger):
-    """
-    Validate parameter file against schema.
-
-    Parameters
-    ----------
-    input_file (str)
-        input file to be validated
-    schema (dict)
-        schema used for validation
-
-    Raises
-    ------
-    FileNotFoundError
-        if input file is not found
-
-    """
-
-    try:
-        with open(input_file, "r", encoding="utf-8") as file:
-            data = yaml.safe_load(file)
-    except FileNotFoundError:
-        logger.error(f"Input file {input_file} not found")
-        raise
-
-    try:
-        jsonschema.validate(data, schema=schema)
-    except jsonschema.exceptions.ValidationError:
-        logger.error(f"Schema validation failed for {input_file} using {schema}")
-        raise
-
-    logger.info(f"Schema validation successful for {input_file}")
-
-
-def main():
-    label = Path(__file__).stem
-    args_dict, _ = _parse(label, description="Parameter file schema checking")
-
-    logger = logging.getLogger()
-    logger.setLevel(gen.get_log_level_from_user(args_dict["log_level"]))
-
-    validate_schema_file(args_dict["file_name"], load_schema(args_dict["schema"]), logger)
-
-
-if __name__ == "__main__":
-    main()