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

nucleardatapy/fig/matter_setupHIC_fig.py

@@ -1,16 +1,18 @@
 import numpy as np
 import matplotlib.pyplot as plt
-#plt.rcParams.update({'font.size': 16})
+
+# plt.rcParams.update({'font.size': 16})
 # Set clean font settings
 # Set 'DejaVu Sans' as the font
 plt.rcParams["font.family"] = "Times New Roman"
-plt.rcParams['font.family'] = 'serif'
-plt.rcParams['mathtext.fontset'] = 'stixsans'
-plt.rcParams['font.serif'] = ['Times New Roman'] + plt.rcParams['font.serif']
+plt.rcParams["font.family"] = "serif"
+plt.rcParams["mathtext.fontset"] = "stixsans"
+plt.rcParams["font.serif"] = ["Times New Roman"] + plt.rcParams["font.serif"]
 
 import nucleardatapy as nuda
 
-def matter_setupHIC_fig( pname, constraints ):
+
+def matter_setupHIC_fig(pname, inferences):
     """
     Plot nuclear chart (N versus Z).\
     The plot is 1x2 with:\
@@ -27,75 +29,113 @@ def matter_setupHIC_fig( pname, constraints ):
 
     """
     #
-    print(f'Plot name: {pname}')
+    print(f"Plot name: {pname}")
     #
-    fig, axs = plt.subplots(2,2)
-    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.10, bottom=0.12, right=None, top=0.98, wspace=0.25, hspace=0.05)
+    fig, axs = plt.subplots(2, 2)
+    fig.tight_layout()  # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust( left=0.10, bottom=0.12, right=0.95, top=0.98, wspace=0.25, hspace=0.05 )
     #
-    #axs[0,0].set_xlabel(r'$n$ (fm$^{-3}$)',fontsize='12')
-    #axs[0,1].set_xlabel(r'$n$ (fm$^{-3}$)',fontsize='12')
-    axs[1,0].set_xlabel(r'$n$ (fm$^{-3}$)',fontsize='12')
-    axs[1,1].set_xlabel(r'$n$ (fm$^{-3}$)',fontsize='12')
+    # axs[0,0].set_xlabel(r'$n$ (fm$^{-3}$)',fontsize='12')
+    # axs[0,1].set_xlabel(r'$n$ (fm$^{-3}$)',fontsize='12')
+    axs[1, 0].set_xlabel(r"$n_\text{nuc}$ (fm$^{-3}$)", fontsize="12")
+    axs[1, 1].set_xlabel(r"$n_\text{nuc}$ (fm$^{-3}$)", fontsize="12")
     #
-    axs[0,0].set_ylabel(r'$p_\mathrm{SM}$ (MeV fm$^{-3}$)',fontsize='12')
-    axs[0,1].set_ylabel(r'$e_\mathrm{SM}$ (MeV)',fontsize='12')
-    axs[1,0].set_ylabel(r'$p_\mathrm{NM}$ (MeV fm$^{-3}$)',fontsize='12')
-    axs[1,1].set_ylabel(r'Symmetry energy (MeV)',fontsize='12')
+    axs[0, 0].set_ylabel(r"$p_\mathrm{SM}(n_\text{nuc})$ (MeV fm$^{-3}$)", fontsize="12")
+    axs[0, 1].set_ylabel(r"$e_\mathrm{SM}(n_\text{nuc})$ (MeV)", fontsize="12")
+    axs[1, 0].set_ylabel(r"$p_\mathrm{NM}(n_\text{nuc})$ (MeV fm$^{-3}$)", fontsize="12")
+    axs[1, 1].set_ylabel(r"$e_\text{sym}(n_\text{nuc})$ (MeV)", fontsize="12")
     #
-    axs[0,0].set_xlim([0.16, 0.8])
-    axs[0,0].set_ylim([0.5, 400])
-    axs[0,1].set_xlim([0.0, 0.44])
-    axs[0,1].set_ylim([-18, 20])
-    axs[1,0].set_xlim([0.16, 0.8])
-    axs[1,0].set_ylim([0.5, 400])
-    axs[1,1].set_xlim([0.0, 0.44])
-    axs[1,1].set_ylim([0.0, 80])
+    axs[0, 0].set_xlim([0.16, 0.8])
+    axs[0, 0].set_ylim([0.5, 400])
+    axs[0, 1].set_xlim([0.0, 0.44])
+    axs[0, 1].set_ylim([-18, 20])
+    axs[1, 0].set_xlim([0.16, 0.8])
+    axs[1, 0].set_ylim([0.5, 400])
+    axs[1, 1].set_xlim([0.0, 0.44])
+    axs[1, 1].set_ylim([0.0, 80])
     #
-    axs[0,0].set_yscale('log')
-    axs[1,0].set_yscale('log')
+    axs[0, 0].set_yscale("log")
+    axs[1, 0].set_yscale("log")
     #
-    axs[0,0].tick_params('x', labelbottom=False)
-    axs[0,1].tick_params('x', labelbottom=False)
+    axs[0, 0].tick_params("x", labelbottom=False)
+    axs[0, 1].tick_params("x", labelbottom=False)
     #
-    for constraint in constraints:
+    for inference in inferences:
         #
-        print('constraint:',constraint)
-        hic = nuda.matter.setupHIC( constraint = constraint )
+        print("inference:", inference)
+        hic = nuda.matter.setupHIC( inference = inference )
         #
         if hic.sm_pre is not None:
-            axs[0,0].fill_between( hic.den, y1=hic.sm_pre_lo, y2=hic.sm_pre_up, label=hic.label, alpha=hic.alpha*0.8, color = hic.color )
+            axs[0, 0].fill_between(
+                hic.den_pre,
+                y1=hic.sm_pre_lo,
+                y2=hic.sm_pre_up,
+                label=hic.label,
+                alpha=hic.alpha * 0.8,
+                color=hic.color,
+            )
         #
-        if hic.sm_e2a is not None:
-            axs[0,1].fill_between( hic.den_e2a, y1=hic.sm_e2a_lo, y2=hic.sm_e2a_up, label=hic.label, alpha=hic.alpha, color ='magenta' )
+        if hic.sm_e2a_int_lo is not None:
+            axs[0, 1].fill_between(
+                hic.den_e2a,
+                y1=hic.sm_e2a_int_lo,
+                y2=hic.sm_e2a_int_up,
+                label=hic.label,
+                alpha=hic.alpha,
+                color="magenta",
+            )
         #
         if hic.nm_pre is not None:
-            axs[1,0].fill_between( hic.den, y1=hic.nm_pre_lo, y2=hic.nm_pre_up, label=hic.label_so, alpha=0.2, color ='b' )
-            axs[1,0].fill_between( hic.den, y1=hic.nm_pre_st_lo, y2=hic.nm_pre_st_up, label=hic.label_st, alpha=0.2, color ='g' )
-        #
-        if hic.sym_enr_isodiff is not None:
-            axs[1,1].errorbar( hic.den_isodiff, hic.sym_enr_isodiff, xerr = hic.den_isodiff_err, yerr=hic.sym_enr_isodiff_err, fmt='o', 
-                               label=hic.label, color = hic.color, capsize=2, capthick=1, elinewidth=1 , markersize =3)
+            axs[1, 0].fill_between(
+                hic.den_pre,
+                y1=hic.nm_pre_lo,
+                y2=hic.nm_pre_up,
+                label=hic.label_so,
+                alpha=0.2,
+                color="b",
+            )
+            axs[1, 0].fill_between(
+                hic.den_pre,
+                y1=hic.nm_pre_st_lo,
+                y2=hic.nm_pre_st_up,
+                label=hic.label_st,
+                alpha=0.2,
+                color="g",
+            )
         #
-        if hic.sym_enr is not None:
-            axs[1,1].fill_between( hic.den, y1=hic.sym_enr_lo, y2=hic.sym_enr_up, label=hic.label, alpha=hic.alpha*0.7, color = hic.color )
+        if hic.esym is not None and hic.den_err:
+            axs[1, 1].errorbar(
+                hic.den_esym,
+                hic.esym,
+                xerr=hic.den_esym_err,
+                yerr=hic.esym_err,
+                fmt="o",
+                label=hic.label,
+                color=hic.color,
+                capsize=2,
+                capthick=1,
+                elinewidth=1,
+                markersize=3,
+            )
         #
-        if hic.sym_enr_np is not None:
-            axs[1,1].errorbar( hic.den_np, hic.sym_enr_np, xerr = hic.den_np_err, yerr=hic.sym_enr_np_err, fmt='o', 
-                               label=hic.label, color = hic.color, capsize=2, capthick=1, elinewidth=1, markersize =3 )
-        #
-        if hic.sym_enr_spirit is not None:
-            axs[1,1].errorbar( hic.den_spirit, hic.sym_enr_spirit, xerr = hic.den_spirit_err, yerr=hic.sym_enr_spirit_err, fmt='o', 
-                               label=hic.label, color = hic.color, capsize=2, capthick=1, elinewidth=1 , markersize =3)
+        if hic.esym is not None and not hic.den_err:
+            axs[1, 1].fill_between(
+                hic.den_esym,
+                y1=hic.esym_lo,
+                y2=hic.esym_up,
+                label=hic.label,
+                alpha=hic.alpha * 0.7,
+                color=hic.color,
+            )
         #
     #
-    #axs.text(0.15,12,r'$K_{sym}$='+str(int(Ksym))+' MeV',fontsize='12')
-    axs[0,0].legend(loc='lower right',fontsize='8')
-    axs[0,1].legend(loc='upper left',fontsize='8')
-    axs[1,0].legend(loc='lower right',fontsize='8')
-    axs[1,1].legend(loc='lower right',fontsize='8')
+    # axs.text(0.15,12,r'$K_{sym}$='+str(int(Ksym))+' MeV',fontsize='12')
+    axs[0, 0].legend(loc="lower right", fontsize="8")
+    axs[0, 1].legend(loc="upper left", fontsize="8")
+    axs[1, 0].legend(loc="lower right", fontsize="8")
+    axs[1, 1].legend(loc="lower right", fontsize="8")
     #
     if pname is not None:
-    	plt.savefig(pname, dpi=300)
-    	plt.close()
-    #
+        plt.savefig(pname, dpi=300)
+        plt.close()
+    #

nucleardatapy/fig/matter_setupMicroEsym_fig.py

@@ -3,7 +3,8 @@ import matplotlib.pyplot as plt
 
 import nucleardatapy as nuda
 
-def matter_setupMicroEsym_fig( pname, models, band ):
+
+def matter_setupMicroEsym_fig(pname, mbs, band):
     """
     Plot nuclear chart (N versus Z).\
     The plot is 1x2 with:\
@@ -20,66 +21,281 @@ def matter_setupMicroEsym_fig( pname, models, band ):
 
     """
     #
-    print(f'Plot name: {pname}')
+    print(f"Plot name: {pname}")
+    #
+    fig, axs = plt.subplots(2, 2)
+    fig.tight_layout()  # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust( left=0.12, bottom=0.12, right=0.95, top=0.9, wspace=0.05, hspace=0.05 )
     #
-    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.8, wspace=0.2, hspace=0.2 )
+    axs[0, 0].set_ylabel(r"$E_\mathrm{sym}$ (MeV)", fontsize="14")
+    axs[0, 0].set_xlim([0, 0.33])
+    axs[0, 0].set_ylim([0, 50])
+    axs[0, 0].tick_params("x", labelbottom=False)
     #
-    axs[0,0].set_ylabel(r'$E_\mathrm{sym}$ (MeV)')
-    axs[0,0].set_xlim([0, 0.3])
-    axs[0,0].set_ylim([0, 50])
+    axs[0, 1].set_xlim([0.5, 2.0])
+    axs[0, 1].set_ylim([0, 50])
+    axs[0, 1].tick_params("x", labelbottom=False)
+    axs[0, 1].tick_params("y", labelleft=False)
     #
-    axs[0,1].set_xlim([0.5, 2.0])
-    axs[0,1].set_ylim([0, 50])
+    axs[1, 0].set_ylabel(r"$E_\mathrm{sym}/E_\mathrm{sym, FFG, NR}$", fontsize="14")
+    axs[1, 0].set_xlabel(r"$n_\mathrm{nuc}$ (fm$^{-3}$)", fontsize="14")
+    axs[1, 0].set_xlim([0,   0.33])
+    axs[1, 0].set_ylim([1.5, 3.6])
     #
-    axs[1,0].set_ylabel(r'$E_\mathrm{sym}/E_\mathrm{sym, FFG, NR}$')
-    axs[1,0].set_xlabel(r'n (fm$^{-3}$)')
-    axs[1,0].set_xlim([0, 0.3])
-    axs[1,0].set_ylim([1, 4])
+    axs[1, 1].set_xlabel(r"$k_F$ (fm$^{-1}$)", fontsize="14")
+    axs[1, 1].set_xlim([0.5, 2.0])
+    axs[1, 1].set_ylim([1.5, 3.6])
+    axs[1, 1].tick_params("y", labelleft=False)
     #
-    axs[1,1].set_xlabel(r'$k_F$ (fm$^{-1}$)')
-    axs[1,1].set_xlim([0.5, 2.0])
-    axs[1,1].set_ylim([1, 4])
+    mb_check = []
     #
-    for ind,model in enumerate(models):
+    for kmb, mb in enumerate(mbs):
         #
-        print('in Sample: model',model)
-        esym = nuda.matter.setupMicroEsym( model = model )
+        models, models_lower = nuda.matter.micro_models_mb(mb)
         #
-        if esym.esym is not None:
-            if '2024-BHF' in model and (ind % 4 != 0.0): continue
-            if esym.esym_err is not None:
-                axs[0,0].errorbar( esym.den, esym.esym, yerr=esym.esym_err, marker=esym.marker, markevery=esym.every, linestyle='none', errorevery=esym.every, label=esym.label )
-                axs[0,1].errorbar( esym.kf,  esym.esym, yerr=esym.esym_err, marker=esym.marker, markevery=esym.every, linestyle='none', errorevery=esym.every )
-                axs[1,0].errorbar( esym.den, esym.esym/nuda.esymffg_nr(esym.kf), yerr=esym.esym_err/nuda.esymffg_nr(esym.kf), marker=esym.marker, markevery=esym.every, linestyle='none', errorevery=esym.every )
-                axs[1,1].errorbar( esym.kf,  esym.esym/nuda.esymffg_nr(esym.kf), yerr=esym.esym_err/nuda.esymffg_nr(esym.kf), marker=esym.marker, markevery=esym.every, linestyle='none', errorevery=esym.every )
+        for model in models:
+            #
+            print("in Sample: model", model)
+            #
+            micro = nuda.matter.setupMicroEsym(model=model)
+            if nuda.env.verb:
+                micro.print_outputs()
+            #
+            micro = nuda.matter.setupMicroEsym(model=model)
+            if nuda.env.verb_output:
+                micro.print_outputs()
+            #
+            check = nuda.matter.setupCheck(eos=micro, band=band)
+            #
+            if check.isInside:
+                lstyle = "solid"
             else:
-                axs[0,0].plot( esym.den, esym.esym, marker=esym.marker, markevery=esym.every, linestyle='none', label=esym.label )
-                axs[0,1].plot( esym.kf,  esym.esym, marker=esym.marker, markevery=esym.every, linestyle='none' )
-                axs[1,0].plot( esym.den, esym.esym/nuda.esymffg_nr(esym.kf), marker=esym.marker, markevery=esym.every, linestyle='none' )
-                axs[1,1].plot( esym.kf,  esym.esym/nuda.esymffg_nr(esym.kf), marker=esym.marker, markevery=esym.every, linestyle='none' )
-        if nuda.env.verb_output: esym.print_outputs( )
-    #
-    axs[0,0].plot( esym.den, nuda.esymffg_nr(esym.kf), linestyle='dashed' )
-    axs[0,1].plot( esym.kf,  nuda.esymffg_nr(esym.kf), linestyle='dashed' )
+                lstyle = "dashed"
+            #
+            if micro.esym is not None:
+                # if '2024-BHF' in model and (kmb % 4 != 0.0): continue
+                if mb in mb_check:
+                    if micro.esym_err is None:
+                        axs[0, 0].plot(
+                            micro.den,
+                            micro.esym,
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            color=nuda.param.col[kmb],
+                        )
+                        axs[0, 1].plot(
+                            micro.kf,
+                            micro.esym,
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            color=nuda.param.col[kmb],
+                        )
+                        axs[1, 0].plot(
+                            micro.den,
+                            micro.esym / nuda.esymffg_nr(micro.kf),
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            color=nuda.param.col[kmb],
+                        )
+                        axs[1, 1].plot(
+                            micro.kf,
+                            micro.esym / nuda.esymffg_nr(micro.kf),
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            color=nuda.param.col[kmb],
+                        )
+                    else:
+                        axs[0, 0].errorbar(
+                            micro.den,
+                            micro.esym,
+                            yerr=micro.esym_err,
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            errorevery=micro.every,
+                            color=nuda.param.col[kmb],
+                        )
+                        axs[0, 1].errorbar(
+                            micro.kf,
+                            micro.esym,
+                            yerr=micro.esym_err,
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            errorevery=micro.every,
+                            color=nuda.param.col[kmb],
+                        )
+                        axs[1, 0].errorbar(
+                            micro.den,
+                            micro.esym / nuda.esymffg_nr(micro.kf),
+                            yerr=micro.esym_err / nuda.esymffg_nr(micro.kf),
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            errorevery=micro.every,
+                            color=nuda.param.col[kmb],
+                        )
+                        axs[1, 1].errorbar(
+                            micro.kf,
+                            micro.esym / nuda.esymffg_nr(micro.kf),
+                            yerr=micro.esym_err / nuda.esymffg_nr(micro.kf),
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            errorevery=micro.every,
+                            color=nuda.param.col[kmb],
+                        )
+                else:
+                    mb_check.append(mb)
+                    if micro.esym_err is None:
+                        axs[0, 0].plot(
+                            micro.den,
+                            micro.esym,
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            color=nuda.param.col[kmb],
+                            label=mb,
+                        )
+                        axs[0, 1].plot(
+                            micro.kf,
+                            micro.esym,
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            color=nuda.param.col[kmb],
+                        )
+                        axs[1, 0].plot(
+                            micro.den,
+                            micro.esym / nuda.esymffg_nr(micro.kf),
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            color=nuda.param.col[kmb],
+                        )
+                        axs[1, 1].plot(
+                            micro.kf,
+                            micro.esym / nuda.esymffg_nr(micro.kf),
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            color=nuda.param.col[kmb],
+                        )
+                    else:
+                        axs[0, 0].errorbar(
+                            micro.den,
+                            micro.esym,
+                            yerr=micro.esym_err,
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            errorevery=micro.every,
+                            color=nuda.param.col[kmb],
+                            label=mb,
+                        )
+                        axs[0, 1].errorbar(
+                            micro.kf,
+                            micro.esym,
+                            yerr=micro.esym_err,
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            errorevery=micro.every,
+                            color=nuda.param.col[kmb],
+                        )
+                        axs[1, 0].errorbar(
+                            micro.den,
+                            micro.esym / nuda.esymffg_nr(micro.kf),
+                            yerr=micro.esym_err / nuda.esymffg_nr(micro.kf),
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            errorevery=micro.every,
+                            color=nuda.param.col[kmb],
+                        )
+                        axs[1, 1].errorbar(
+                            micro.kf,
+                            micro.esym / nuda.esymffg_nr(micro.kf),
+                            yerr=micro.esym_err / nuda.esymffg_nr(micro.kf),
+                            marker=micro.marker,
+                            markevery=micro.every,
+                            linestyle=lstyle,
+                            errorevery=micro.every,
+                            color=nuda.param.col[kmb],
+                        )
+
+    # FFG symmetry energy
+    #axs[0, 0].plot(micro.den, nuda.esymffg_nr(micro.kf), linestyle="dotted")
+    #axs[0, 1].plot(micro.kf, nuda.esymffg_nr(micro.kf), linestyle="dotted")
 
-    axs[0,0].fill_between( band.den, y1=(band.e2a-band.e2a_std), y2=(band.e2a+band.e2a_std), color=band.color, alpha=band.alpha, visible=True )
-    axs[0,0].plot( band.den, (band.e2a-band.e2a_std), color='k', linestyle='dashed', visible=True )
-    axs[0,0].plot( band.den, (band.e2a+band.e2a_std), color='k', linestyle='dashed', visible=True )
-    axs[0,1].fill_between( band.kfn, y1=(band.e2a-band.e2a_std), y2=(band.e2a+band.e2a_std), color=band.color, alpha=band.alpha, visible=True )
-    axs[0,1].plot( band.kfn, (band.e2a-band.e2a_std), color='k', linestyle='dashed', visible=True )
-    axs[0,1].plot( band.kfn, (band.e2a+band.e2a_std), color='k', linestyle='dashed', visible=True )
-    axs[1,0].fill_between( band.den, y1=(band.e2a-band.e2a_std)/nuda.esymffg_nr(band.kf), y2=(band.e2a+band.e2a_std)/nuda.esymffg_nr(band.kf), color=band.color, alpha=band.alpha, visible=True )
-    axs[1,0].plot( band.den, (band.e2a-band.e2a_std)/nuda.esymffg_nr(band.kf), color='k', linestyle='dashed', visible=True )
-    axs[1,0].plot( band.den, (band.e2a+band.e2a_std)/nuda.esymffg_nr(band.kf), color='k', linestyle='dashed', visible=True )
-    axs[1,1].fill_between( band.kfn, y1=(band.e2a-band.e2a_std)/nuda.esymffg_nr(band.kf), y2=(band.e2a+band.e2a_std)/nuda.esymffg_nr(band.kf), color=band.color, alpha=band.alpha, visible=True )
-    axs[1,1].plot( band.kfn, (band.e2a-band.e2a_std)/nuda.esymffg_nr(band.kf), color='k', linestyle='dashed', visible=True )
-    axs[1,1].plot( band.kfn, (band.e2a+band.e2a_std)/nuda.esymffg_nr(band.kf), color='k', linestyle='dashed', visible=True )
+    axs[0, 0].fill_between(
+        band.den,
+        y1=(band.e2a - band.e2a_std),
+        y2=(band.e2a + band.e2a_std),
+        color=band.color,
+        alpha=band.alpha,
+        visible=True,
+    )
+    axs[0, 0].plot(
+        band.den, (band.e2a - band.e2a_std), color="k", linestyle="dashed", visible=True, zorder = 100 )
+    axs[0, 0].plot(
+        band.den, (band.e2a + band.e2a_std), color="k", linestyle="dashed", visible=True, zorder = 100 )
+    axs[0, 1].fill_between(
+        band.kfn,
+        y1=(band.e2a - band.e2a_std),
+        y2=(band.e2a + band.e2a_std),
+        color=band.color,
+        alpha=band.alpha,
+        visible=True,
+    )
+    axs[0, 1].plot(
+        band.kfn, (band.e2a - band.e2a_std), color="k", linestyle="dashed", visible=True, zorder = 100 )
+    axs[0, 1].plot(
+        band.kfn, (band.e2a + band.e2a_std), color="k", linestyle="dashed", visible=True, zorder = 100 )
+    axs[1, 0].fill_between(
+        band.den,
+        y1=(band.e2a - band.e2a_std) / nuda.esymffg_nr(band.kf),
+        y2=(band.e2a + band.e2a_std) / nuda.esymffg_nr(band.kf),
+        color=band.color,
+        alpha=band.alpha,
+        visible=True,
+    )
+    axs[1, 0].plot(
+        band.den, (band.e2a - band.e2a_std) / nuda.esymffg_nr(band.kf),
+        color="k", linestyle="dashed", visible=True, zorder = 100 )
+    axs[1, 0].plot(
+        band.den, (band.e2a + band.e2a_std) / nuda.esymffg_nr(band.kf),
+        color="k", linestyle="dashed", visible=True, zorder = 100 )
+    axs[1, 1].fill_between(
+        band.kfn,
+        y1=(band.e2a - band.e2a_std) / nuda.esymffg_nr(band.kf),
+        y2=(band.e2a + band.e2a_std) / nuda.esymffg_nr(band.kf),
+        color=band.color,
+        alpha=band.alpha,
+        visible=True,
+    )
+    axs[1, 1].plot(
+        band.kfn, (band.e2a - band.e2a_std) / nuda.esymffg_nr(band.kf),
+        color="k", linestyle="dashed", visible=True, zorder = 100 )
+    axs[1, 1].plot(
+        band.kfn, (band.e2a + band.e2a_std) / nuda.esymffg_nr(band.kf),
+        color="k", linestyle="dashed", visible=True, zorder = 100 )
 
-    #axs[1,0].legend(loc='upper right',fontsize='8')
-    fig.legend(loc='upper left',bbox_to_anchor=(0.1,1.0),columnspacing=2,fontsize='8',ncol=3,frameon=False)
+    # axs[1,0].legend(loc='upper right',fontsize='8')
+    fig.legend(
+        loc="upper left",
+        bbox_to_anchor=(0.1, 1.0),
+        columnspacing=2,
+        fontsize="8",
+        ncol=4,
+        frameon=False,
+    )
     #
     if pname is not None:
-    	plt.savefig(pname, dpi=300)
-    	plt.close()
+        plt.savefig(pname, dpi=300)
+        plt.close()
+    #