nucleardatapy 0.2.1__py3-none-any.whl → 1.0.0__py3-none-any.whl

nucleardatapy/fig/matter_setupCheck_fig.py

@@ -0,0 +1,103 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_setupCheck_fig( pname, mb, models, band ):
+    """
+    Plot nucleonic energy per particle E/A in matter.\
+    The plot is 1x1 with: E/A versus den.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param mb: many-body (mb) approach considered.
+    :type mb: str.
+    :param models: models to run on.
+    :type models: array of str.
+    :param band: object instantiated on the reference band.
+    :type band: object.
+    :param matter: can be 'SM' or 'NM'.
+    :type matter: str.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    matter = band.matter
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.10, bottom=0.12, right=0.95, top=0.85, wspace=0.05, hspace=0.05 )
+    #
+    axs.set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs.set_xlim([0, 0.33])
+    if matter.lower() == 'nm':
+        axs.set_ylabel(r'$E_\text{NM}^\text{int}/A$ (MeV)')
+        axs.set_ylim([0, 30])
+        delta = 1.0
+    elif matter.lower() == 'sm':
+        axs.set_ylabel(r'$E_\text{SM}^\text{int}/A$ (MeV)')
+        axs.set_ylim([-20, 10])
+        delta = 0.0
+    #
+    for model in models:
+        #
+        mic = nuda.matter.setupMicro( model = model, var2 = delta )
+        if nuda.env.verb_output: mic.print_outputs( )
+        #
+        print('model:',model,' delta:',delta)
+        #
+        check = nuda.matter.setupCheck( eos = mic, band = band )
+        #
+        if check.isInside:
+            lstyle = 'solid'
+        else:
+            lstyle = 'dashed'
+        #
+        if mic.e_err:
+            #
+            print('=> model (with err):',model,mic.e_err)
+            if 'NLEFT' in model:
+                #
+                print('   => model (NLEFT):',model)
+                if matter.lower() == 'nm':
+                    axs.errorbar( mic.nm_den, mic.nm_e2adata_int, yerr=mic.nm_e2adata_err, linestyle = 'dotted', markevery=mic.every, linewidth = 1, alpha=0.6 )
+                    axs.fill_between( mic.nm_den, y1=(mic.nm_e2a_int-mic.nm_e2a_err), y2=(mic.nm_e2a_int+mic.nm_e2a_err), alpha=0.3 )
+                elif matter.lower() == 'sm':
+                    axs.errorbar( mic.sm_den, mic.sm_e2adata_int, yerr=mic.sm_e2adata_err, linestyle = 'dotted', markevery=mic.every, linewidth = 1, alpha=0.6 )
+                    axs.fill_between( mic.sm_den, y1=(mic.sm_e2a_int-mic.sm_e2a_err), y2=(mic.sm_e2a_int+mic.sm_e2a_err), alpha=0.3 )
+                #
+            if mic.marker:
+                #
+                print('with marker:',mic.marker)
+                if matter.lower() == 'nm':
+                    axs.errorbar( mic.nm_den, mic.nm_e2a_int, yerr=mic.nm_e2a_err, marker=mic.marker, markevery=mic.every, linestyle=lstyle, label=mic.label, errorevery=mic.every )
+                elif matter.lower() == 'sm':
+                    axs.errorbar( mic.sm_den, mic.sm_e2a_int, yerr=mic.sm_e2a_err, marker=mic.marker, markevery=mic.every, linestyle=lstyle, label=mic.label, errorevery=mic.every )
+                #
+            else:
+                #
+                print('with no marker:',mic.marker)
+                if matter.lower() == 'nm':
+                    axs.errorbar( mic.nm_den, mic.nm_e2a_int, yerr=mic.nm_e2a_err, marker=mic.marker, markevery=mic.every, linestyle=lstyle, label=mic.label, errorevery=mic.every )
+                elif matter.lower() == 'sm':
+                    axs.errorbar( mic.sm_den, mic.sm_e2a_int, yerr=mic.sm_e2a_err, marker=mic.marker, markevery=mic.every, linestyle=lstyle, label=mic.label, errorevery=mic.every )
+        else:
+            print('=> model (no err):',model,mic.e_err)
+            if 'fit' in model:
+                axs.plot( mic.den, mic.e2a_int, marker=mic.marker, linestyle=lstyle, markevery=mic.every, label=mic.label )
+            else:
+                if matter.lower() == 'nm':
+                    axs.plot( mic.nm_den, mic.nm_e2a_int, marker=mic.marker, linestyle=lstyle, markevery=mic.every, label=mic.label )
+                elif matter.lower() == 'sm':
+                    axs.plot( mic.sm_den, mic.sm_e2a_int, marker=mic.marker, linestyle=lstyle, markevery=mic.every, label=mic.label )
+        #
+    axs.fill_between( band.den, y1=(band.e2a_int-band.e2a_std), y2=(band.e2a_int+band.e2a_std), color=band.color, alpha=band.alpha, visible=True )
+    axs.plot( band.den, (band.e2a_int-band.e2a_std), color='k', linestyle='dashed' )
+    axs.plot( band.den, (band.e2a_int+band.e2a_std), color='k', linestyle='dashed' )
+    #
+    fig.legend(loc='upper left',bbox_to_anchor=(0.1,1.0),fontsize='8',ncol=3,frameon=False)
+    #
+    if pname is not None: 
+        plt.savefig(pname, dpi=300)
+        plt.close()
+    #

nucleardatapy/fig/matter_setupFFGLep_fig.py

@@ -0,0 +1,70 @@
+import numpy as np
+import matplotlib as mpl
+import matplotlib.pyplot as plt
+import nucleardatapy as nuda
+
+mpl.use("Agg")  # Use a non-interactive backend
+
+def matter_setupFFGLep_fig( pname, den_el=None, den_mu1=None, den_mu2=None, den_mu3=None ):
+    """
+    Plot leptonic FFG energy per particle E/A and pressure in NM and SM.\
+    The plot is 2x2 with:\
+    [0,0]: E/A versus den.     [0,1]: E/A versus kfn.\
+    [1,0]: pre versus den.     [1,1]: pre versus kfn.\
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param den: density.
+    :type den: float or numpy vector of real numbers.
+    :param kfn: neutron Fermi momentum.
+    :type kfn: float or numpy vector of real numbers.
+
+    """
+    #
+    print(f"Plot name: {pname}")
+    #
+    if den_el is None: den_el=np.linspace(0.01, 0.1, num=20)
+    if den_mu1 is None: den_mu1 = 0.1*np.linspace(0.01, 0.1, num=20) 
+    if den_mu2 is None: den_mu2 = 0.2*np.linspace(0.01, 0.1, num=20) 
+    if den_mu3 is None: den_mu3 = 0.5*np.linspace(0.01, 0.1, num=20) 
+    #
+    fig, axs = plt.subplots(2, 1)
+    fig.tight_layout()  # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust( left=0.12, bottom=0.12, right=None, top=0.9, wspace=0.05, hspace=0.05 )
+    #
+    axs[0].set_ylabel(r"$e^\text{FFG}$ (MeV)", fontsize="14")
+    axs[0].set_xlim([0, 0.1])
+    axs[0].set_ylim([0, 200])
+    axs[0].tick_params("x", labelbottom=False)
+    #
+    axs[1].set_xlabel(r"$n_\text{el}$ (fm$^{-3}$)", fontsize="14")
+    axs[1].set_ylabel(r"$p^\text{FFG}$ (MeV fm$^{-3}$)", fontsize="14")
+    axs[1].set_xlim([0, 0.1])
+    axs[1].set_ylim([-0.01, 6])
+    #
+    lep1 = nuda.matter.setupFFGLep( den_el=den_el, den_mu=den_mu1 )
+    lep2 = nuda.matter.setupFFGLep( den_el=den_el, den_mu=den_mu2 )
+    lep3 = nuda.matter.setupFFGLep( den_el=den_el, den_mu=den_mu3 )
+    #
+    if any(lep1.e2n_el):
+        print(r"plot $\delta=0$ (SM)")
+        axs[0].plot( lep1.den_el, lep1.e2n_el, linestyle="solid",  color=nuda.param.col[0], label='electrons' )
+        axs[0].plot( lep1.den_el, lep1.e2n_mu, linestyle="dashed", color=nuda.param.col[1], label='muons (10%)' )
+        axs[0].plot( lep1.den_el, lep2.e2n_mu, linestyle="dashed", color=nuda.param.col[2], label='muons (20%)' )
+        axs[0].plot( lep1.den_el, lep3.e2n_mu, linestyle="dashed", color=nuda.param.col[3], label='muons (50%)' )
+        axs[1].plot( lep1.den_el, lep1.pre_el, linestyle="solid",  color=nuda.param.col[0] )
+        axs[1].plot( lep1.den_el, lep1.pre_mu, linestyle="dashed", color=nuda.param.col[1] )
+        axs[1].plot( lep1.den_el, lep2.pre_mu, linestyle="dashed", color=nuda.param.col[2] )
+        axs[1].plot( lep1.den_el, lep3.pre_mu, linestyle="dashed", color=nuda.param.col[3] )
+        if nuda.env.verb_output: lep.print_outputs()
+        #
+    #axs[0, 0].text(0.2, 16, r"$m=$" + str(mss[0]) + "$m_N$", rotation=8)
+    #axs[0, 0].text(0.2, 32, r"$m=$" + str(mss[1]) + "$m_N$", rotation=13)
+    #axs[0, 0].text(0.2, 50, r"$m=$" + str(mss[2]) + "$m_N$", rotation=20)
+    axs[0].legend(loc='lower right',fontsize='10')
+    #fig.legend( loc="upper left", bbox_to_anchor=(0.2, 0.97), fontsize="6", ncol=4, frameon=False )
+    #
+    if pname is not None:
+        plt.savefig(pname, dpi=300)
+        plt.close()
+    #

nucleardatapy/fig/matter_setupFFGNuc_fig.py

@@ -1,11 +1,13 @@
-
 import numpy as np
 import matplotlib as mpl
 import matplotlib.pyplot as plt
-mpl.use("Agg")  # Use a non-interactive backend
 import nucleardatapy as nuda
 
-def matter_setupFFGNuc_EP_fig( pname, mss = [ 1.0 ], den = np.linspace(0.01,0.35,10), kfn = np.linspace(0.5,2.0,10) ):
+mpl.use("Agg")  # Use a non-interactive backend
+
+def matter_setupFFGNuc_EP_fig(
+    pname, mss=[1.0], den=np.linspace(0.01, 0.35, 10), kf=np.linspace(0.5, 2.0, 10)
+):
     """
     Plot nucleonic FFG energy per particle E/A and pressure in NM and SM.\
     The plot is 2x2 with:\
@@ -21,87 +23,186 @@ def matter_setupFFGNuc_EP_fig( pname, mss = [ 1.0 ], den = np.linspace(0.01,0.35
 
     """
     #
-    print(f'Plot name: {pname}')
+    print(f"Plot name: {pname}")
     #
-    denkfn = nuda.matter.den( kfn )
+    denkf = nuda.matter.den(kf)
     delta0 = np.zeros(den.size)
     delta1 = np.ones(den.size)
     #
-    fig, axs = plt.subplots(2,2)
-    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.9, wspace=0.05, hspace=0.05 )
-    #
-    axs[0,0].set_ylabel(r'$e_\text{FFG,int}$')
-    axs[0,0].set_xlim([0, 0.33])
-    axs[0,0].set_ylim([0, 100])
-    axs[0,0].tick_params('x', labelbottom=False)
-    #
-    axs[1,0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    axs[1,0].set_ylabel(r'$p_\text{FFG}$ (MeV fm$^{-3}$)')
-    axs[1,0].set_xlim([0, 0.33])
-    axs[1,0].set_ylim([0, 23])
-    #
-    #axs[0,1].set_ylabel(r'$\Delta/E_F$')
-    axs[0,1].set_xlim([0.5, 2.0])
-    axs[0,1].set_ylim([0, 100])
-    axs[0,1].tick_params('y', labelleft=False)
-    axs[0,1].tick_params('x', labelbottom=False)
-    #
-    axs[1,1].set_xlabel(r'$k_{F_n}$ (fm$^{-1}$)')
-    #axs[1,1].set_ylabel(r'$\Delta$ (MeV)')
-    axs[1,1].set_xlim([0.5, 2.0])
-    axs[1,1].set_ylim([0, 23])
-    axs[1,1].tick_params('y', labelleft=False)
-    #
-    for ims,ms in enumerate(mss):
-        ffg0 = nuda.matter.setupFFGNuc( den, delta0, ms )
-        ffg0kf = nuda.matter.setupFFGNuc( denkfn, delta0, ms )
-        ffg1 = nuda.matter.setupFFGNuc( den, delta1, ms )
-        ffg1kf = nuda.matter.setupFFGNuc( denkfn, delta1, ms )
+    fig, axs = plt.subplots(2, 2)
+    fig.tight_layout()  # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust( left=0.12, bottom=0.12, right=None, top=0.9, wspace=0.05, hspace=0.05 )
+    #
+    axs[0, 0].set_ylabel(r"$e^\text{FFG}_\text{int}$ (MeV)")
+    axs[0, 0].set_xlim([0, 0.33])
+    axs[0, 0].set_ylim([0, 100])
+    axs[0, 0].tick_params("x", labelbottom=False)
+    #
+    axs[1, 0].set_xlabel(r"$n_\text{nuc}$ (fm$^{-3}$)")
+    axs[1, 0].set_ylabel(r"$p^\text{FFG}$ (MeV fm$^{-3}$)")
+    axs[1, 0].set_xlim([0, 0.33])
+    axs[1, 0].set_ylim([0, 23])
+    #
+    # axs[0,1].set_ylabel(r'$\Delta/E_F$')
+    #axs[0, 1].set_xlim([0.5, 2.0])
+    #axs[0, 1].set_ylim([0, 100])
+    axs[0, 1].tick_params("y", labelleft=False)
+    axs[0, 1].tick_params("x", labelbottom=False)
+    #
+    axs[1, 1].set_xlabel(r"$k_{F}$ (fm$^{-1}$)")
+    # axs[1,1].set_ylabel(r'$\Delta$ (MeV)')
+    #axs[1, 1].set_xlim([0.5, 2.0])
+    #axs[1, 1].set_ylim([0, 23])
+    axs[1, 1].tick_params("y", labelleft=False)
+    #
+    for ims, ms in enumerate(mss):
+        ffg0 = nuda.matter.setupFFGNuc(den, delta0, ms)
+        ffg0kf = nuda.matter.setupFFGNuc(denkf, delta0, ms)
+        ffg1 = nuda.matter.setupFFGNuc(den, delta1, ms)
+        ffg1kf = nuda.matter.setupFFGNuc(denkf, delta1, ms)
         #
-        if any(ffg0.e2a_int_nr): 
-            print(r'plot $\delta=0$ (SM)')
+        if any(ffg0.e2a_int_nr):
+            print(r"plot $\delta=0$ (SM)")
             if ims == 0:
-                axs[0,0].plot( ffg0.den, ffg0.e2a_int, linestyle='solid', color=nuda.param.col[0], label=ffg0.label )
-                axs[0,0].plot( ffg0.den, ffg0.e2a_int_nr, linestyle='None', marker='o', color=nuda.param.col[0], label='NR'+ffg0.label )
+                axs[0, 0].plot(
+                    ffg0.den,
+                    ffg0.e2a_int,
+                    linestyle="solid",
+                    color=nuda.param.col[0],
+                    label=ffg0.label,
+                )
+                axs[0, 0].plot(
+                    ffg0.den,
+                    ffg0.e2a_int_nr,
+                    linestyle="None",
+                    marker="o",
+                    color=nuda.param.col[0],
+                    label="NR" + ffg0.label,
+                )
             else:
-                axs[0,0].plot( ffg0.den, ffg0.e2a_int, linestyle='solid', color=nuda.param.col[0] )
-                axs[0,0].plot( ffg0.den, ffg0.e2a_int_nr, linestyle='None', marker='o', color=nuda.param.col[0] )
-            axs[1,0].plot( ffg0.den, ffg0.pre, linestyle='solid', color=nuda.param.col[0] )
-            axs[1,0].plot( ffg0.den, ffg0.pre_nr, linestyle='None', marker='o', color=nuda.param.col[0] )
-            axs[0,1].plot( ffg0kf.kf_n, ffg0kf.e2a_int, linestyle='solid', color=nuda.param.col[0] )
-            axs[0,1].plot( ffg0kf.kf_n, ffg0kf.e2a_int_nr, linestyle='None', marker='o', color=nuda.param.col[0] )
-            axs[1,1].plot( ffg0kf.kf_n, ffg0kf.pre, linestyle='solid', color=nuda.param.col[0] )
-            axs[1,1].plot( ffg0kf.kf_n, ffg0kf.pre_nr, linestyle='None', marker='o', color=nuda.param.col[0] )
-        if nuda.env.verb_output: ffg0.print_outputs( )
-        if any(ffg1.e2a_int_nr): 
-            print(r'plot $\delta=1$ (NM)')
+                axs[0, 0].plot(
+                    ffg0.den, ffg0.e2a_int, linestyle="solid", color=nuda.param.col[0]
+                )
+                axs[0, 0].plot(
+                    ffg0.den,
+                    ffg0.e2a_int_nr,
+                    linestyle="None",
+                    marker="o",
+                    color=nuda.param.col[0],
+                )
+            axs[1, 0].plot(
+                ffg0.den, ffg0.pre, linestyle="solid", color=nuda.param.col[0]
+            )
+            axs[1, 0].plot(
+                ffg0.den,
+                ffg0.pre_nr,
+                linestyle="None",
+                marker="o",
+                color=nuda.param.col[0],
+            )
+            axs[0, 1].plot(
+                ffg0kf.kf_nuc, ffg0kf.e2a_int, linestyle="solid", color=nuda.param.col[0]
+            )
+            axs[0, 1].plot(
+                ffg0kf.kf_nuc,
+                ffg0kf.e2a_int_nr,
+                linestyle="None",
+                marker="o",
+                color=nuda.param.col[0],
+            )
+            axs[1, 1].plot(
+                ffg0kf.kf_nuc, ffg0kf.pre, linestyle="solid", color=nuda.param.col[0]
+            )
+            axs[1, 1].plot(
+                ffg0kf.kf_nuc,
+                ffg0kf.pre_nr,
+                linestyle="None",
+                marker="o",
+                color=nuda.param.col[0],
+            )
+        if nuda.env.verb_output:
+            ffg0.print_outputs()
+        if any(ffg1.e2a_int_nr):
+            print(r"plot $\delta=1$ (NM)")
             if ims == 0:
-                axs[0,0].plot( ffg1.den, ffg1.e2a_int, linestyle='dashed', color=nuda.param.col[1], label=ffg1.label )
-                axs[0,0].plot( ffg1.den, ffg1.e2a_int_nr, linestyle='None', marker='o', color=nuda.param.col[1], label='NR'+ffg1.label )
+                axs[0, 0].plot(
+                    ffg1.den,
+                    ffg1.e2a_int,
+                    linestyle="dashed",
+                    color=nuda.param.col[1],
+                    label=ffg1.label,
+                )
+                axs[0, 0].plot(
+                    ffg1.den,
+                    ffg1.e2a_int_nr,
+                    linestyle="None",
+                    marker="o",
+                    color=nuda.param.col[1],
+                    label="NR" + ffg1.label,
+                )
             else:
-                axs[0,0].plot( ffg1.den, ffg1.e2a_int, linestyle='dashed', color=nuda.param.col[1] )
-                axs[0,0].plot( ffg1.den, ffg1.e2a_int_nr, linestyle='None', marker='o', color=nuda.param.col[1] )
-            axs[1,0].plot( ffg1.den, ffg1.pre, linestyle='dashed', color=nuda.param.col[1] )
-            axs[1,0].plot( ffg1.den, ffg1.pre_nr, linestyle='None', marker='o', color=nuda.param.col[1] )
-            axs[0,1].plot( ffg1kf.kf_n, ffg1kf.e2a_int, linestyle='dashed', color=nuda.param.col[1] )
-            axs[0,1].plot( ffg1kf.kf_n, ffg1kf.e2a_int_nr, linestyle='None', marker='o', color=nuda.param.col[1] )
-            axs[1,1].plot( ffg1kf.kf_n, ffg1kf.pre, linestyle='dashed', color=nuda.param.col[1] )
-            axs[1,1].plot( ffg1kf.kf_n, ffg1kf.pre_nr, linestyle='None', marker='o', color=nuda.param.col[1] )
-        if nuda.env.verb_output: ffg1.print_outputs( )
+                axs[0, 0].plot(
+                    ffg1.den, ffg1.e2a_int, linestyle="dashed", color=nuda.param.col[1]
+                )
+                axs[0, 0].plot(
+                    ffg1.den,
+                    ffg1.e2a_int_nr,
+                    linestyle="None",
+                    marker="o",
+                    color=nuda.param.col[1],
+                )
+            axs[1, 0].plot(
+                ffg1.den, ffg1.pre, linestyle="dashed", color=nuda.param.col[1]
+            )
+            axs[1, 0].plot(
+                ffg1.den,
+                ffg1.pre_nr,
+                linestyle="None",
+                marker="o",
+                color=nuda.param.col[1],
+            )
+            axs[0, 1].plot(
+                ffg1kf.kf_nuc, ffg1kf.e2a_int, linestyle="dashed", color=nuda.param.col[1]
+            )
+            axs[0, 1].plot(
+                ffg1kf.kf_nuc,
+                ffg1kf.e2a_int_nr,
+                linestyle="None",
+                marker="o",
+                color=nuda.param.col[1],
+            )
+            axs[1, 1].plot(
+                ffg1kf.kf_nuc, ffg1kf.pre, linestyle="dashed", color=nuda.param.col[1]
+            )
+            axs[1, 1].plot(
+                ffg1kf.kf_nuc,
+                ffg1kf.pre_nr,
+                linestyle="None",
+                marker="o",
+                color=nuda.param.col[1],
+            )
+        if nuda.env.verb_output:
+            ffg1.print_outputs()
         #
-    axs[0,0].text(0.2,16,r'$m^*/m=$'+str(mss[0]),rotation=8)
-    axs[0,0].text(0.2,32,r'$m^*/m=$'+str(mss[1]),rotation=13)
-    axs[0,0].text(0.2,50,r'$m^*/m=$'+str(mss[2]),rotation=20)
-    #axs[1,0].legend(loc='upper right',fontsize='xx-small')
-    fig.legend(loc='upper left',bbox_to_anchor=(0.2,0.97),fontsize='6',ncol=4,frameon=False)
+    axs[0, 0].text(0.2, 16, r"$m=$" + str(mss[0]) + "$m_N$", rotation=8)
+    axs[0, 0].text(0.2, 32, r"$m=$" + str(mss[1]) + "$m_N$", rotation=13)
+    axs[0, 0].text(0.2, 50, r"$m=$" + str(mss[2]) + "$m_N$", rotation=20)
+    # axs[1,0].legend(loc='upper right',fontsize='xx-small')
+    fig.legend(
+        loc="upper left",
+        bbox_to_anchor=(0.2, 0.97),
+        fontsize="6",
+        ncol=4,
+        frameon=False,
+    )
     #
-    if pname is not None: 
+    if pname is not None:
         plt.savefig(pname, dpi=300)
         plt.close()
     #
 
-def matter_setupFFGNuc_EOS_fig( pname, mss = [ 1.0 ], den = np.linspace(0.01,0.35,10) ):
+
+def matter_setupFFGNuc_EOS_fig(pname, mss=[1.0], den=np.linspace(0.01, 0.35, 10)):
     """
     Plot nucleonic FFG EOS in NM and SM.\
     The plot is 1x2 with:\
@@ -115,59 +216,121 @@ def matter_setupFFGNuc_EOS_fig( pname, mss = [ 1.0 ], den = np.linspace(0.01,0.3
 
     """
     #
-    print(f'Plot name: {pname}')
+    print(f"Plot name: {pname}")
     #
     delta0 = np.zeros(den.size)
     delta1 = np.ones(den.size)
     #
-    fig, axs = plt.subplots(2,1)
-    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.9, wspace=0.05, hspace=0.05 )
+    fig, axs = plt.subplots(2, 1)
+    fig.tight_layout()  # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust( left=0.12, bottom=0.12, right=None, top=0.9, wspace=0.05, hspace=0.05 )
     #
-    axs[0].set_ylabel(r'$p_\text{FFG}$')
+    axs[0].set_ylabel(r"$p^\text{FFG}$ (MeV fm$^{-3}$)", fontsize="14")
     axs[0].set_xlim([0, 360])
     axs[0].set_ylim([0, 20])
-    axs[0].tick_params('x', labelbottom=False)
+    axs[0].tick_params("x", labelbottom=False)
     #
-    axs[1].set_xlabel(r'$\rho$ (MeV fm$^{-3}$)')
-    axs[1].set_ylabel(r'$c_\text{s,FFG}^2$')
+    axs[1].set_xlabel(r"$\epsilon^\text{FFG}$ (MeV fm$^{-3}$)", fontsize="14")
+    axs[1].set_ylabel(r"$(c^\text{FFG}_\text{s}/c)^2$", fontsize="14")
     axs[1].set_xlim([0, 360])
     axs[1].set_ylim([0, 0.28])
     #
-    for ims,ms in enumerate(mss):
-        ffg0 = nuda.matter.setupFFGNuc( den, delta0, ms )
-        ffg1 = nuda.matter.setupFFGNuc( den, delta1, ms )
+    for ims, ms in enumerate(mss):
+        ffg0 = nuda.matter.setupFFGNuc(den, delta0, ms)
+        ffg1 = nuda.matter.setupFFGNuc(den, delta1, ms)
         #
-        if any(ffg0.e2a_int_nr): 
-            print(r'plot $\delta=0$ (SM)')
+        if any(ffg0.e2a_int_nr):
+            print(r"plot $\delta=0$ (SM)")
             if ims == 0:
-                axs[0].plot( ffg0.e2v, ffg0.pre, linestyle='solid', color=nuda.param.col[0], label=ffg0.label )
-                axs[0].plot( ffg0.e2v, ffg0.pre_nr, linestyle='None', marker='o', color=nuda.param.col[0], label='NR'+ffg0.label )
+                axs[0].plot(
+                    ffg0.eps,
+                    ffg0.pre,
+                    linestyle="solid",
+                    color=nuda.param.col[0],
+                    label=ffg0.label,
+                )
+                axs[0].plot(
+                    ffg0.eps,
+                    ffg0.pre_nr,
+                    linestyle="None",
+                    marker="o",
+                    color=nuda.param.col[0],
+                    label="NR" + ffg0.label,
+                )
             else:
-                axs[0].plot( ffg0.e2v, ffg0.pre, linestyle='solid', color=nuda.param.col[0] )
-                axs[0].plot( ffg0.e2v, ffg0.pre_nr, linestyle='None', marker='o', color=nuda.param.col[0] )
-            axs[1].plot( ffg0.e2v, ffg0.cs2, linestyle='solid', color=nuda.param.col[0] )
-            axs[1].plot( ffg0.e2v, ffg0.cs2_nr, linestyle='None', marker='o', color=nuda.param.col[0] )
-        if nuda.env.verb_output: ffg0.print_outputs( )
-        if any(ffg1.e2a_int_nr): 
-            print(r'plot $\delta=1$ (NM)')
+                axs[0].plot(
+                    ffg0.eps, ffg0.pre, linestyle="solid", color=nuda.param.col[0]
+                )
+                axs[0].plot(
+                    ffg0.eps,
+                    ffg0.pre_nr,
+                    linestyle="None",
+                    marker="o",
+                    color=nuda.param.col[0],
+                )
+            axs[1].plot(ffg0.eps, ffg0.cs2, linestyle="solid", color=nuda.param.col[0])
+            axs[1].plot(
+                ffg0.eps,
+                ffg0.cs2_nr,
+                linestyle="None",
+                marker="o",
+                color=nuda.param.col[0],
+            )
+        if nuda.env.verb_output:
+            ffg0.print_outputs()
+        if any(ffg1.e2a_int_nr):
+            print(r"plot $\delta=1$ (NM)")
             if ims == 0:
-                axs[0].plot( ffg1.e2v, ffg1.pre, linestyle='dashed', color=nuda.param.col[1], label=ffg1.label )
-                axs[0].plot( ffg1.e2v, ffg1.pre_nr, linestyle='None', marker='o', color=nuda.param.col[1], label='NR'+ffg1.label )
+                axs[0].plot(
+                    ffg1.eps,
+                    ffg1.pre,
+                    linestyle="dashed",
+                    color=nuda.param.col[1],
+                    label=ffg1.label,
+                )
+                axs[0].plot(
+                    ffg1.eps,
+                    ffg1.pre_nr,
+                    linestyle="None",
+                    marker="o",
+                    color=nuda.param.col[1],
+                    label="NR" + ffg1.label,
+                )
             else:
-                axs[0].plot( ffg1.e2v, ffg1.pre, linestyle='dashed', color=nuda.param.col[1] )
-                axs[0].plot( ffg1.e2v, ffg1.pre_nr, linestyle='None', marker='o', color=nuda.param.col[1] )
-            axs[1].plot( ffg1.e2v, ffg1.cs2, linestyle='dashed', color=nuda.param.col[1] )
-            axs[1].plot( ffg1.e2v, ffg1.cs2_nr, linestyle='None', marker='o', color=nuda.param.col[1] )
-        if nuda.env.verb_output: ffg1.print_outputs( )
+                axs[0].plot(
+                    ffg1.eps, ffg1.pre, linestyle="dashed", color=nuda.param.col[1]
+                )
+                axs[0].plot(
+                    ffg1.eps,
+                    ffg1.pre_nr,
+                    linestyle="None",
+                    marker="o",
+                    color=nuda.param.col[1],
+                )
+            axs[1].plot(ffg1.eps, ffg1.cs2, linestyle="dashed", color=nuda.param.col[1])
+            axs[1].plot(
+                ffg1.eps,
+                ffg1.cs2_nr,
+                linestyle="None",
+                marker="o",
+                color=nuda.param.col[1],
+            )
+        if nuda.env.verb_output:
+            ffg1.print_outputs()
         #
-    axs[1].text(300,0.07,r'$m^*/m=$'+str(mss[0]),rotation=3)
-    axs[1].text(240,0.12,r'$m^*/m=$'+str(mss[1]),rotation=4)
-    axs[1].text(180,0.2,r'$m^*/m=$'+str(mss[2]),rotation=5)
-    #axs[1,0].legend(loc='upper right',fontsize='xx-small')
-    fig.legend(loc='upper left',bbox_to_anchor=(0.15,0.97),fontsize='8',ncol=4,frameon=False)
+    axs[1].text(300, 0.07, r"$m=$" + str(mss[0]) + "$m_N$", rotation=3)
+    axs[1].text(240, 0.12, r"$m=$" + str(mss[1]) + "$m_N$", rotation=4)
+    axs[1].text(180, 0.20, r"$m=$" + str(mss[2]) + "$m_N$", rotation=5)
+    # axs[1,0].legend(loc='upper right',fontsize='xx-small')
+    fig.legend(
+        loc="upper left",
+        bbox_to_anchor=(0.15, 0.97),
+        fontsize="8",
+        ncol=4,
+        frameon=False,
+    )
     #
-    if pname is not None: 
+    if pname is not None:
         plt.savefig(pname, dpi=300)
         plt.close()
     #