xarpes 0.2.4__py3-none-any.whl → 0.6.0__py3-none-any.whl

xarpes/settings_parameters.py

@@ -0,0 +1,75 @@
+# Copyright (C) 2025 xARPES Developers
+# This program is free software under the terms of the GNU GPLv3 license.
+
+"""User-configurable numerical parameters for xARPES."""
+
+# Extend data range by this many Gaussian sigmas
+sigma_extend = 5.0
+
+# Gaussian confidence level expressed in "sigma"
+sigma_confidence = 2.0
+
+def parameter_settings(new_sigma_extend=None, new_sigma=None):
+    """
+    Configure global numerical parameters for xARPES.
+
+    Parameters
+    ----------
+    new_sigma_extend : float or None
+        Number of Gaussian sigmas used to extend arrays before convolution.
+    new_sigma : float or None
+        Gaussian confidence level expressed in units of sigma
+        (e.g. 1, 2, 3). Default is 2.
+    """
+    global sigma_extend, sigma_confidence
+
+    if new_sigma_extend is not None:
+        sigma_extend = float(new_sigma_extend)
+
+    if new_sigma is not None:
+        sigma_confidence = float(new_sigma)
+
+# ---------------- Defaults for MEM / chi2kink alpha-selection ----------------
+
+mem_defaults = {
+    "method": "chi2kink",
+    "parts": "both",
+    "iter_max": 1e4,
+    "alpha_min": 1.0,
+    "alpha_max": 9.0,
+    "alpha_num": 10,
+    "ecut_left": 0.0,
+    "ecut_right": None,
+    "f_chi_squared": None,
+    "W": None,
+    "power": 4,
+    "mu": 1.0,
+    "omega_S": 1.0,
+    "sigma_svd": 1e-4,
+    "t_criterion": 1e-8,
+    "a_guess": 1.0,
+    "b_guess": 2.5,
+    "c_guess": 3.0,
+    "d_guess": 1.5,
+    "h_n": None, 
+    "impurity_magnitude": 0.0,
+    "lambda_el": 0.0,
+}
+
+# ---------------- Defaults for bayesian_loop optimisation --------------------
+
+loop_defaults = {
+    "converge_iters": 50,
+    "tole": 1e-2,
+    "scale_vF": 1.0,
+    "scale_mb": 1.0,
+    "scale_imp": 1.0,
+    "scale_kF": 1.0,
+    "scale_lambda_el": 1.0,
+    "scale_hn": 1.0,
+    "opt_iter_max": 1e4,
+    "rollback_steps": 10,
+    "max_retries": 100,
+    "relative_best": 10.0,
+    "min_steps_for_regression": 25,
+}

xarpes/settings_plots.py

@@ -0,0 +1,54 @@
+# Copyright (C) 2025 xARPES Developers
+# This program is free software under the terms of the GNU GPLv3 license.
+
+"""Plotting and notebook behaviour settings for xARPES."""
+
+import matplotlib.pyplot as plt
+
+def plot_settings(name="default", register_pre_run=True):
+    """Configure default plotting style for xARPES.
+
+    Parameters
+    ----------
+    name : {"default", "large"}
+        Select a predefined style.
+    register_pre_run : bool
+        If True, register a Jupyter pre-run hook that closes figures.
+    """
+    import matplotlib as mpl
+
+    mpl.rc("xtick", labelsize=10, direction="in")
+    mpl.rc("ytick", labelsize=10, direction="in")
+    plt.rcParams["legend.frameon"] = False
+    lw = dict(default=2.0, large=4.0)[name]
+
+    mpl.rcParams.update({
+        "lines.linewidth": lw,
+        "lines.markersize": 3,
+        "xtick.major.size": 4,
+        "xtick.minor.size": 2,
+        "xtick.major.width": 0.8,
+        "font.size": 16,
+        "axes.ymargin": 0.15,
+    })
+
+    if register_pre_run:
+        _maybe_register_pre_run_close_all()
+
+
+def _maybe_register_pre_run_close_all():
+    """Register a pre_run_cell hook once, and only inside Jupyter."""
+    from IPython import get_ipython
+
+    if getattr(_maybe_register_pre_run_close_all, "_registered", False):
+        return
+
+    ip = get_ipython()
+    if ip is None or ip.__class__.__name__ != "ZMQInteractiveShell":
+        return
+
+    def _close_all(_info):
+        plt.close("all")
+
+    ip.events.register("pre_run_cell", _close_all)
+    _maybe_register_pre_run_close_all._registered = True

xarpes-0.6.0.dist-info/METADATA

@@ -0,0 +1,181 @@
+Metadata-Version: 2.1
+Name: xarpes
+Version: 0.6.0
+Summary: Extraction from angle resolved photoemission spectra
+Author: xARPES Developers
+Requires-Python: >=3.7.0
+Description-Content-Type: text/markdown
+Classifier: License :: OSI Approved :: GNU General Public License v3 (GPLv3)
+Classifier: Programming Language :: Python :: 3
+Requires-Dist: igor2>=0.5.8
+Requires-Dist: jupyterlab
+Requires-Dist: jupytext
+Requires-Dist: matplotlib<3.10.0
+Requires-Dist: numpy
+Requires-Dist: scipy
+Requires-Dist: lmfit
+Requires-Dist: pyqt5
+Requires-Dist: ipympl>=0.9.3
+Requires-Dist: ipywidgets>=8.1.5
+Requires-Dist: ipykernel<6.32.0
+Project-URL: Documentation, https://xarpes.github.io
+
+# xARPES 
+
+![xARPES](https://github.com/xARPES/xARPES/raw/main/logo/xarpes.svg)
+
+Repository for the code xARPES &ndash; extraction of the self-energy and Eliashberg function from angle-resolved photoemission spectroscopy. The documentation can be found on [Read the Docs](https://xarpes.readthedocs.io), while the code is maintained on [GitHub](https://github.com/xARPES/xARPES). Instructions for installing the code and downloading the code are found below. An extensive description of the functionalities and examples is found at the [arXiv preprint](https://arxiv.org/abs/2508.13845).
+
+# Warning
+
+This project is currently undergoing **beta testing**. Some of the functionalities are in the process of being implemented. If you encounter any bugs, you can open an issue.
+
+# Contributing
+
+Contributions to the code are most welcome. xARPES is intended to co-develop alongside the increasing complexity of experimental ARPES data sets. Contributions can be made by forking the code and creating a pull request. Importing of file formats from different beamlines is particularly encouraged. For development of the examples, the following scripts in `/dev_tools` could be useful:
+- The `Rmd2ipynb.py` file, to generate the `.ipynb` from which the examples can conveniently be developed; to be executed after cloning/pulling updated `.Rmd` and `.py` files, or if you prefer developing with `.Rmd` files.
+- The `ipynb2Rmd2py.py` file, to synchronise the `.Rmd` and `.py` files from the `.ipynb` file; to be executed right before pushing modifications to the repository. Note that this script resets some metadata in the `.Rmd` to prevent the tracking of their changes due to local virtual environments, etc.
+
+# Installation
+
+xARPES installation can be divided into graphical package manager instructions, and command-line instructions for either conda or pip. With command-line instructions, an editable installation of xARPES can be created; on Windows, we strongly recommend Windows Powershell to do so. Here is a summary for the three options:
+- via a graphical package manager (Anaconda Navigator, VS Code, PyCharm, Spyder, JupyterLab)  
+- via conda-forge, out-of-the-box or editable installation, sourcing the [conda-forge package](https://anaconda.org/conda-forge/xarpes).  
+- via Pip, out-of-the-box or editable installation, sourcing the [PyPI package](https://pypi.org/project/xarpes).
+
+We strongly recommend installing xARPES in a (conda/pip) virtual environment, and to activate the environment each time before activating xARPES.
+
+## Graphical package manager installation
+
+Most IDEs and scientific Python distributions include a GUI-based package manager.  
+These typically install from conda-forge (for conda environments) or PyPI (for venv/system Python).
+
+### Anaconda Navigator
+
+1. Open Anaconda Navigator  
+2. Select or create an environment  
+3. Set the package channel to conda-forge  
+4. Search for “xarpes”  
+5. Click Install
+
+This installs the latest stable release from conda-forge.
+
+### PyCharm, VS Code, Spyder, or JupyterLab
+
+These IDEs install from the active environment’s package source:
+- conda environment → installs from conda-forge  
+- venv/system Python → installs from PyPI
+
+### Installation steps (generic)
+
+1. Open your IDE  
+2. Select or create a Python environment  
+3. Open the environment/package manager panel  
+4. Search for “xarpes”  
+5. Click Install
+
+## Conda Forge installation
+
+Install xARPES inside a conda environment, either out of the box or as an editable.
+
+### Setting up a conda environment
+
+Download and install Miniconda (see the [Miniconda installation page](https://docs.anaconda.com/free/miniconda/)).
+
+Example for Linux:
+
+    wget https://repo.anaconda.com/miniconda/Miniconda3-latest-Linux-x86_64.sh  
+    bash Miniconda3-latest-Linux-x86_64.sh
+
+Answer `y` to questions. Create and activate a new environment:
+
+    conda create -n <my_env> -c conda-forge  
+    conda activate <my_env>
+
+Where `<my_env>` must be replaced by your desired name. Package compatibility issues may arise if conda installs from different channels. This can be prevented by appending `--strict-channel-priority` to the creation command.
+
+### Installing xARPES
+
+#### Option A — Out-of-the-box installation (from conda-forge)
+
+    conda install conda-forge::xarpes
+
+#### Option B — Editable installation (GitHub)
+
+First, clone the repository:
+
+    git clone git@github.com:xARPES/xARPES.git  
+    cd xARPES
+
+Then perform editable installation (this mixes conda and pip):
+
+    pip install -e .
+
+## Pip installation
+
+Install xARPES using pip, either out of the box or as an editable.
+
+### Setting up a virtual environment
+
+Install venv if necessary:
+
+    sudo apt install python3-venv
+
+Create and activate a virtual environment:
+
+    python3 -m venv <my_env>  
+    source <my_env>/bin/activate
+
+Upgrade pip:
+
+    python3 -m pip install --upgrade pip
+
+### Installing xARPES
+
+#### Option A — Out-of-the-box installation (PyPI)
+
+    python3 -m pip install xarpes
+
+#### Option B — Editable installation (GitHub)
+
+First, clone the repository:
+
+    git clone git@github.com:xARPES/xARPES.git  
+    cd xARPES
+
+Then perform editable installation:
+
+    pip install -e .
+
+# Examples
+
+After installation of xARPES, the `examples/` folder can be downloaded from the terminal into the current directory:
+
+    xarpes_download_examples
+
+This attempts to download the examples from the version corresponding encountered in `__init__.py`. If no corresponding tagged version can be downloaded, the code attempts to download the latest examples instead.
+
+The examples can also be installed by executing `.py`or `.ipynb` executable files containing the following:
+
+    import xarpes
+    xarpes.download_examples()
+
+# Execution
+
+It is recommended to use JupyterLab to analyse data. JupyterLab is launched using:
+
+    jupyter lab
+
+# Citation
+
+If you have used xARPES for your research, please cite the following preprint:  
+[arXiv preprint 2508.13845](https://arxiv.org/abs/2508.13845).
+
+# License
+
+Copyright (C) 2025 xARPES Developers
+
+This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License, version 3, as published by the Free Software Foundation.
+
+This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+

xarpes-0.6.0.dist-info/RECORD

@@ -0,0 +1,15 @@
+xarpes/__init__.py,sha256=dfT64reIAoDadP6Kz8T8s8lXi7jgcoau4Jy2Du_GzCY,756
+xarpes/bandmap.py,sha256=CGmXhY1ViZMlpDdqpgVUJsmhGRQdcPLa79l53QMP-Fc,34365
+xarpes/constants.py,sha256=XOdgSzyrmHr5xocHZfmcFHHoVAa1G05a305hm3XOTtY,504
+xarpes/distributions.py,sha256=pC8V5MlZDNFdooMonFREEASiN5QodHiyKc2ehnxMKvQ,23498
+xarpes/functions.py,sha256=_ecottx3vi1twu_7lUrO_G3UbhCd0YMIaKezYEAW0Sw,20837
+xarpes/mdcs.py,sha256=WRKSfGlRVKBssJp9FIHcAFsINVunPkmW9fBnFjqBHYI,42844
+xarpes/plotting.py,sha256=lGCReHcXhYLQXR5ns3EHFjCQjJ9Sc-HifV7n4BnWby4,5189
+xarpes/selfenergies.py,sha256=mal1SiezXq0Upg8xBL8PagUQN2ZAzCF84geRtY3gSDc,64611
+xarpes/settings_parameters.py,sha256=yOYvgEiDeDiLzzLkvysCTiVwqg6fKIkN48B-WSad728,1912
+xarpes/settings_plots.py,sha256=X-qteB2fIbBKOAcLMvMYDfQ8QdlUeA5xYQqF_Nyb4uA,1562
+xarpes-0.6.0.dist-info/entry_points.txt,sha256=917UR-cqFTMMI_vMqIbk7boYSuFX_zHwQlXKcj9vlCE,79
+xarpes-0.6.0.dist-info/LICENSE,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
+xarpes-0.6.0.dist-info/WHEEL,sha256=jPMR_Dzkc4X4icQtmz81lnNY_kAsfog7ry7qoRvYLXw,81
+xarpes-0.6.0.dist-info/METADATA,sha256=NlUZioy3lvyew1P0QlKuBbyxO0g_YkTC6A5MmzVZLdY,7154
+xarpes-0.6.0.dist-info/RECORD,,

{xarpes-0.2.4.dist-info → xarpes-0.6.0.dist-info}/WHEEL

@@ -1,4 +1,4 @@
 Wheel-Version: 1.0
-Generator: flit 3.9.0
+Generator: flit 3.6.0
 Root-Is-Purelib: true
 Tag: py3-none-any

xarpes-0.6.0.dist-info/entry_points.txt

@@ -0,0 +1,3 @@
+[console_scripts]
+xarpes_download_examples=xarpes.functions:download_examples
+

xarpes/.ipynb_checkpoints/__init__-checkpoint.py

@@ -1,8 +0,0 @@
-__version__ = '0.2.2'
-
-from . import plotting
-
-from .band_map import *
-from .distributions import *
-from .functions import *
-from .plotting import *

xarpes/band_map.py

@@ -1,302 +0,0 @@
-# Copyright (C) 2024 xARPES Developers
-# This program is free software under the terms of the GNU GPLv3 license.
-
-# get_ax_fig_plt and add_fig_kwargs originate from pymatgen/util/plotting.py.
-# Copyright (C) 2011-2024 Shyue Ping Ong and the pymatgen Development Team
-# Pymatgen is released under the MIT License.
-
-# See also abipy/tools/plotting.py.
-# Copyright (C) 2021 Matteo Giantomassi and the AbiPy Group
-# AbiPy is free software under the terms of the GNU GPLv2 license.
-
-"""The band map class and allowed operations on it."""
-
-import numpy as np
-
-from .plotting import get_ax_fig_plt, add_fig_kwargs
-from .functions import fit_leastsq
-from .distributions import fermi_dirac
-
-class MDCs():
-    r"""
-    """
-    def __init__(self, intensities, angles, ekin, angular_resolution):
-        self.intensities = intensities
-        self.angles = angles
-        self.ekin = ekin
-        self.angular_resolution = angular_resolution
-
-    # @add_fig_kwargs
-    # def fit(self):
-    #     r"""
-    #     """
-    #     return 0
-
-    # @add_fig_kwargs
-    # def fit_fermi_edge(self, hnuminphi_guess, background_guess=0.0,
-    #                    integrated_weight_guess=1.0, angle_min=-np.inf,
-    #                    angle_max=np.inf, ekin_min=-np.inf,
-    #                    ekin_max=np.inf, ax=None, **kwargs):
-
-    @add_fig_kwargs
-    def plot(self, ax=None, **kwargs):
-        r"""
-        """
-        ax, fig, plt = get_ax_fig_plt(ax=ax)
-
-        ax.scatter(self.angles, self.intensities)
-
-        ax.set_xlabel('Angle ($\degree$)')
-        ax.set_ylabel('Counts (-)')
-
-        return fig
-
-class band_map():
-    r"""Class for the band map from the ARPES experiment.
-
-    Parameters
-    ----------
-    intensities : ndarray
-        2D array of counts for given (E,k) or (E,angle) pairs [counts]
-    angles : ndarray
-        1D array of angular values for the abscissa [degrees]
-    ekin : ndarray
-        1D array of kinetic energy values for the ordinate [eV]
-    angular_resolution : float, None
-        Angular resolution of the detector [degrees]
-    energy_resolution : float, None
-        Energy resolution of the detector [eV]
-    temperature : float, None
-        Temperature of the sample [K]
-    hnuminphi : float, None
-        Kinetic energy minus the work function [eV]
-    hnuminphi_std : float, None
-        Standard deviation of kinetic energy minus work function [eV]
-    """
-    def __init__(self, intensities, angles, ekin, energy_resolution=None,
-                 angular_resolution=None, temperature=None, hnuminphi=None,
-                 hnuminphi_std=None):
-        self.intensities = intensities
-        self.angles = angles
-        self.ekin = ekin
-        self.energy_resolution = energy_resolution
-        self.angular_resolution = angular_resolution
-        self.temperature = temperature
-        self.hnuminphi = hnuminphi
-        self.hnuminphi_std = hnuminphi_std
-
-    @property
-    def hnuminphi(self):
-        r"""Returns the photon energy minus the work function in eV if it has
-        been set, either during instantiation, with the setter, or by fitting
-        the Fermi-Dirac distribution to the integrated weight.
-
-        Returns
-        -------
-        hnuminphi : float, None
-            Kinetic energy minus the work function [eV]
-        """
-        return self._hnuminphi
-
-    @hnuminphi.setter
-    def hnuminphi(self, hnuminphi):
-        r"""Manually sets the photon energy minus the work function in eV if it
-        has been set; otherwise returns None.
-
-        Parameters
-        ----------
-        hnuminphi : float, None
-            Kinetic energy minus the work function [eV]
-        """
-        self._hnuminphi = hnuminphi
-
-    @property
-    def hnuminphi_std(self):
-        r"""Returns standard deviation of the photon energy minus the work
-        function in eV.
-
-        Returns
-        -------
-        hnuminphi_std : float
-            Standard deviation of energy minus the work function [eV]
-        """
-        return self._hnuminphi_std
-
-    @hnuminphi_std.setter
-    def hnuminphi_std(self, hnuminphi_std):
-        r"""Manually sets the standard deviation of photon energy minus the
-        work function in eV.
-
-        Parameters
-        ----------
-        hnuminphi_std : float
-            Standard deviation of energy minus the work function [eV]
-        """
-        self._hnuminphi_std = hnuminphi_std
-
-    def shift_angles(self, shift):
-        r"""
-        Shifts the angles by the specified amount in degrees. Used to shift
-        from the detector angle to the material angle.
-
-        Parameters
-        ----------
-        shift : float
-            Angular shift [degrees]
-        """
-        self.angles = self.angles + shift
-
-    def slice(self, angle_min, angle_max, energy_value):
-        r"""
-        Parameters
-        ----------
-        angle_min : float
-            Minimum angle of integration interval [degrees]
-        angle_max : float
-            Maximum angle of integration interval [degrees]
-
-        Returns
-        -------
-        angle_range : ndarray
-            Array of size n containing the angular values
-        energy_range : ndarray
-            Array of size m containing the energy values
-        mdcs :
-            Array of size nxm containing the MDC intensities
-        """
-
-        energy_index = np.abs(self.ekin - energy_value).argmin()
-        angle_min_index = np.abs(self.angles - angle_min).argmin()
-        angle_max_index = np.abs(self.angles - angle_max).argmin()
-
-        angle_range = self.angles[angle_min_index:angle_max_index + 1]
-        energy_range = self.ekin[energy_index]
-        mdcs = self.intensities[energy_index,
-               angle_min_index:angle_max_index + 1]
-
-        return mdcs, angle_range, energy_range, self.angular_resolution
-
-    @add_fig_kwargs
-    def plot_band_map(self, ax=None, **kwargs):
-        r"""Plots the band map.
-
-        Parameters
-        ----------
-
-        Other parameters
-        ----------------
-        **kwargs : dict, optional
-            Additional arguments passed on to add_fig_kwargs. See the keyword
-            table below.
-
-        Returns
-        -------
-        fig : Matplotlib-Figure
-            Figure containing the Fermi edge fit
-        """
-        ax, fig, plt = get_ax_fig_plt(ax=ax)
-
-        Angl, Ekin = np.meshgrid(self.angles, self.ekin)
-        mesh = ax.pcolormesh(Angl, Ekin, self.intensities, shading='auto',
-               cmap=plt.get_cmap('bone').reversed())
-        cbar = plt.colorbar(mesh, ax=ax)
-
-        ax.set_xlabel('Angle ($\degree$)')
-        ax.set_ylabel('$E_{\mathrm{kin}}$ (eV)')
-
-        return fig
-
-    @add_fig_kwargs
-    def fit_fermi_edge(self, hnuminphi_guess, background_guess=0.0,
-                       integrated_weight_guess=1.0, angle_min=-np.inf,
-                       angle_max=np.inf, ekin_min=-np.inf,
-                       ekin_max=np.inf, ax=None, **kwargs):
-        r"""Fits the Fermi edge of the band map and plots the result.
-        Also sets hnuminphi, the kinetic energy minus the work function in eV.
-        The fitting includes an energy convolution with an abscissa range
-        expanded by 5 times the energy resolution standard deviation.
-
-        Parameters
-        ----------
-        hnuminphi_guess : float
-            Initial guess for kinetic energy minus the work function [eV]
-        background_guess : float
-            Initial guess for background intensity [counts]
-        integrated_weight_guess : float
-            Initial guess for integrated spectral intensity [counts]
-        angle_min : float
-            Minimum angle of integration interval [degrees]
-        angle_max : float
-            Maximum angle of integration interval [degrees]
-        ekin_min : float
-            Minimum kinetic energy of integration interval [eV]
-        ekin_max : float
-            Maximum kinetic energy of integration interval [eV]
-        ax : Matplotlib-Axes / NoneType
-            Axis for plotting the Fermi edge on. Created if not provided by
-            the user.
-
-        Other parameters
-        ----------------
-        **kwargs : dict, optional
-            Additional arguments passed on to add_fig_kwargs. See the keyword
-            table below.
-
-        Returns
-        -------
-        fig : Matplotlib-Figure
-            Figure containing the Fermi edge fit
-        """
-        ax, fig, plt = get_ax_fig_plt(ax=ax)
-
-        min_angle_index = np.argmin(np.abs(self.angles - angle_min))
-        max_angle_index = np.argmin(np.abs(self.angles - angle_max))
-
-        min_ekin_index = np.argmin(np.abs(self.ekin - ekin_min))
-        max_ekin_index = np.argmin(np.abs(self.ekin - ekin_max))
-
-        energy_range = self.ekin[min_ekin_index:max_ekin_index]
-
-        integrated_intensity = np.trapz(
-            self.intensities[min_ekin_index:max_ekin_index,
-                min_angle_index:max_angle_index], axis=1)
-
-        fdir_initial = fermi_dirac(temperature=self.temperature,
-                                   hnuminphi=hnuminphi_guess,
-                                   background=background_guess,
-                                   integrated_weight=integrated_weight_guess,
-                                   name='Initial guess')
-
-        parameters = np.array(
-            [hnuminphi_guess, background_guess, integrated_weight_guess])
-
-        extra_args = (self.energy_resolution)
-
-        popt, pcov = fit_leastsq(parameters, energy_range, integrated_intensity,
-                                 fdir_initial, extra_args)
-
-        fdir_final = fermi_dirac(temperature=self.temperature,
-                                 hnuminphi=popt[0], background=popt[1],
-                                 integrated_weight=popt[2],
-                                 name='Fitted result')
-
-        self.hnuminphi = popt[0]
-        self.hnuminphi_std = np.sqrt(np.diag(pcov))[0][0]
-
-        ax.set_xlabel(r'$E_{\mathrm{kin}}$ (-)')
-        ax.set_ylabel('Counts (-)')
-        ax.set_xlim([ekin_min, ekin_max])
-
-        ax.plot(energy_range, integrated_intensity, label='Data')
-
-        ax.plot(energy_range, fdir_initial.convolve(energy_range,
-                energy_resolution=self.energy_resolution),
-                label=fdir_initial.name)
-
-        ax.plot(energy_range, fdir_final.convolve(energy_range,
-                energy_resolution=self.energy_resolution),
-                label=fdir_final.name)
-
-        ax.legend()
-
-        return fig