nucleardatapy 0.2.0__py3-none-any.whl

nucleardatapy/fig/matter_setupMicroEsym_fig.py

@@ -0,0 +1,85 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_setupMicroEsym_fig( pname, models, band ):
+    """
+    Plot nuclear chart (N versus Z).\
+    The plot is 1x2 with:\
+    [0]: nuclear chart.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param table: table.
+    :type table: str.
+    :param version: version of table to run on.
+    :type version: str.
+    :param theo_tables: object instantiated on the reference band.
+    :type theo_tables: object.
+
+    """
+    #
+    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=None, top=0.8, wspace=0.2, hspace=0.2 )
+    #
+    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[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}$)')
+    axs[1,1].set_xlim([0.5, 2.0])
+    axs[1,1].set_ylim([1, 4])
+    #
+    for ind,model in enumerate(models):
+        #
+        print('in Sample: model',model)
+        esym = nuda.matter.setupMicroEsym( model = model )
+        #
+        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 )
+            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' )
+
+    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[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)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=300)
+    	plt.close()

nucleardatapy/fig/matter_setupMicro_LP_fig.py

@@ -0,0 +1,109 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_setupMicro_LP_fig( pname, models, matter = 'SM', ell = 0 ):
+    """
+    Plot nucleonic energy per particle E/A in matter.\
+    The plot is 2x2 with:\
+    [0,0]: E/A versus den.       [0,1]: E/A versus kfn.\
+    [1,0]: E/E_NRFFG versus den. [1,1]: E/E_NRFFG versus kfn.\
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param models: models to run on.
+    :type models: array of str.
+    :param matter: can be 'SM' (default) or 'NM'.
+    :type matter: str.
+    :param ell: Value of the angular momentum L of the Landau residual parameter. Can be 0 (default) or 1.
+    :type ell: int.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    if matter.lower() == 'sm':
+    	fig, axs = plt.subplots(2,2)
+    	fig.subplots_adjust(left=0.12, bottom=0.12, right=0.98, top=0.98, wspace=0.15, hspace=0.05 )
+    elif matter.lower() == 'nm':
+    	fig, axs = plt.subplots(1,2)
+    	fig.subplots_adjust(left=0.12, bottom=0.12, right=0.98, top=0.98, wspace=0.15, hspace=0.3 )
+    else:
+    	print('matter_setupMicro_LP_fig: issue with matter:',matter)
+    	print('matter_setupMicro_LP_fig: --- exit() ---')
+    	exit()
+    #
+    #
+    if matter.lower() == 'sm':
+    	axs[0,0].set_xlim([0, 2.0])
+    	axs[0,0].set_ylim([-2.1, 1.6])
+    	if ell == 0: axs[0,0].set_ylabel(r'$F_0$')
+    	if ell == 1: axs[0,0].set_ylabel(r'$F_1$')
+    	axs[0,0].tick_params('x', labelbottom=False)
+    	if ell == 0: axs[0,1].set_ylabel(r'$G_0$')
+    	if ell == 1: axs[0,1].set_ylabel(r'$G_1$')
+    	axs[0,1].set_xlim([0, 2.0])
+    	axs[0,1].set_ylim([-2.1, 1.6])
+    	axs[0,1].tick_params('x', labelbottom=False)
+    	axs[0,1].tick_params('y', labelleft=False)
+    	if ell == 0: axs[1,0].set_ylabel(r'$F_0^\prime$')
+    	if ell == 1: axs[1,0].set_ylabel(r'$F_1^\prime$')
+    	axs[1,0].set_xlabel(r'$k_F$ (fm$^{-1}$)')
+    	axs[1,0].set_xlim([0, 2.0])
+    	axs[1,0].set_ylim([-2.1, 1.6])
+    	if ell == 0: axs[1,1].set_ylabel(r'$G_0^\prime$')
+    	if ell == 1: axs[1,1].set_ylabel(r'$G_1^\prime$')
+    	axs[1,1].set_xlabel(r'$k_F$ (fm$^{-1}$)')
+    	axs[1,1].set_xlim([0, 2.0])
+    	axs[1,1].set_ylim([-2.1, 1.6])
+    	axs[1,1].tick_params('y', labelleft=False)
+    elif matter.lower() == 'nm':
+    	if ell == 0: axs[0].set_ylabel(r'$F_0$')
+    	if ell == 1: axs[0].set_ylabel(r'$F_1$')
+    	axs[0].set_xlabel(r'$k_{F_n}$ (fm$^{-1}$)')
+    	axs[0].set_xlim([0, 2.0])
+    	axs[0].set_ylim([-1.1, 1.1])
+    	if ell == 0: axs[1].set_ylabel(r'$G_0$')
+    	if ell == 1: axs[1].set_ylabel(r'$G_1$')
+    	axs[1].set_xlabel(r'$k_{F_n}$ (fm$^{-1}$)')
+    	axs[1].set_xlim([0, 2.0])
+    	axs[1].set_ylim([-1.1, 1.1])
+    	axs[1].tick_params('y', labelleft=False)
+    #
+    for model in models:
+    #
+    	mic = nuda.matter.setupMicroLP( model = model )
+    	if matter in model:
+    		print('\nmodel:',model,'\n')
+    		if mic.err:
+    			if matter.lower() == 'sm' and mic.sm_kfn is not None:
+    				axs[0,0].errorbar( mic.sm_kfn, mic.sm_LP['F'][ell], yerr=mic.sm_LP_F_err[ell], marker=mic.marker, linestyle='none', label=mic.label )
+    				axs[0,1].errorbar( mic.sm_kfn, mic.sm_LP['G'][ell], yerr=mic.sm_LP_G_err[ell], marker=mic.marker, linestyle='none', label=mic.label )
+    				axs[1,0].errorbar( mic.sm_kfn, mic.sm_LP['Fp'][ell], yerr=mic.sm_LP_Fp_err[ell], marker=mic.marker, linestyle='none', label=mic.label )
+    				axs[1,1].errorbar( mic.sm_kfn, mic.sm_LP['Gp'][ell], yerr=mic.sm_LP_Gp_err[ell], marker=mic.marker, linestyle='none', label=mic.label )
+    			elif matter.lower() == 'nm' and mic.nm_kfn is not None:
+    				axs[0].errorbar( mic.nm_kfn, mic.nm_LP['F'][ell], marker=mic.marker, linestyle='none', label=mic.label )
+    				axs[1].errorbar( mic.nm_kfn, mic.nm_LP['G'][ell], marker=mic.marker, linestyle='none', label=mic.label )
+    		else:
+    			if matter.lower() == 'sm' and mic.sm_kfn is not None:
+    				axs[0,0].plot( mic.sm_kfn, mic.sm_LP['F'][ell], marker=mic.marker, linestyle=mic.linestyle, markevery=mic.every, label=mic.label )
+    				axs[0,1].plot( mic.sm_kfn, mic.sm_LP['G'][ell], marker=mic.marker, linestyle=mic.linestyle, markevery=mic.every, label=mic.label )
+    				axs[1,0].plot( mic.sm_kfn, mic.sm_LP['Fp'][ell], marker=mic.marker, linestyle=mic.linestyle, markevery=mic.every, label=mic.label )
+    				axs[1,1].plot( mic.sm_kfn, mic.sm_LP['Gp'][ell], marker=mic.marker, linestyle=mic.linestyle, markevery=mic.every, label=mic.label )
+    			elif matter.lower() == 'nm' and mic.nm_kfn is not None:
+    				axs[0].plot( mic.nm_kfn, mic.nm_LP['F'][ell], marker=mic.marker, linestyle=mic.linestyle, markevery=mic.every, label=mic.label )
+    				axs[1].plot( mic.nm_kfn, mic.nm_LP['G'][ell], marker=mic.marker, linestyle=mic.linestyle, markevery=mic.every, label=mic.label )
+    	#
+    	if nuda.env.verb_output: mic.print_outputs( )
+    	#
+    #
+    if matter.lower() == 'sm':
+    	axs[1,0].legend(loc='lower left',fontsize='6', ncol=2)
+    elif matter.lower() == 'nm':
+    	axs[0].legend(loc='upper left',fontsize='8', ncol=1)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=300)
+    	plt.close()
+    #

nucleardatapy/fig/matter_setupMicro_band_fig.py

@@ -0,0 +1,85 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_setupMicro_band_fig( pname, models, den, matter ):
+    """
+    Plot the correlation between Esym and Lsym.\
+    The plot is 1x1 with:\
+    [0]: Esym - Lsym correlation plot
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param constraints: list of constraints to run on.
+    :type constraints: array of str.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    #den = 0.03*np.arange(13)
+    #
+    fig, axs = plt.subplots(1,2)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.15, bottom=0.12, right=None, top=0.98, wspace=0.38, hspace=0.3 )
+    #
+    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)',fontsize='14')
+    axs[0].set_xlim([0, 0.2])
+    #
+    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)',fontsize='14')
+    axs[1].set_xlim([0, 0.2])
+    #
+    if matter.lower() == 'nm':
+        axs[0].set_ylabel(r'$E_\text{NM}$',fontsize='14')
+        axs[1].set_ylabel(r'$E_\text{NM}/E_\text{NM,NRFFG}$',fontsize='14')
+        axs[0].set_ylim([0, 20])
+        axs[1].set_ylim([0.2, 0.9])
+    elif matter.lower() == 'sm':
+        axs[0].set_ylabel(r'$E_\text{SM}$',fontsize='14')
+        axs[1].set_ylabel(r'$E_\text{SM}/E_\text{SM,NRFFG}$',fontsize='14')
+        axs[0].set_ylim([-20, 0])
+        axs[1].set_ylim([-1.0, 0.0])
+    elif matter.lower() == 'esym':
+        axs[0].set_ylabel(r'$E_\text{sym}$',fontsize='14')
+        axs[1].set_ylabel(r'$E_\text{sym}/E_\text{sym,NRFFG}$',fontsize='14')
+        axs[0].set_ylim([0, 50])
+        axs[1].set_ylim([1.5, 2.8])
+    #
+    for model in models:
+        #
+        mic = nuda.matter.setupMicro( model = model )
+        if nuda.env.verb_output: mic.print_outputs( )
+        if mic.nm_e2a is not None: 
+            print('model:',model)
+            if matter.lower() == 'nm':
+                axs[0].errorbar( mic.nm_den, mic.nm_e2a, yerr= mic.nm_e2a_err, linestyle=mic.linestyle, label=mic.label, errorevery=mic.every )
+                axs[1].errorbar( mic.nm_den, mic.nm_e2a/nuda.effg_nr(mic.nm_kfn), yerr= mic.nm_e2a_err/nuda.effg_nr(mic.nm_kfn), linestyle=mic.linestyle, label=mic.label, errorevery=mic.every )
+            elif matter.lower() == 'sm':
+                axs[0].errorbar( mic.sm_den, mic.sm_e2a, yerr= mic.sm_e2a_err, linestyle=mic.linestyle, label=mic.label, errorevery=mic.every )
+                axs[1].errorbar( mic.sm_den, mic.sm_e2a/nuda.effg_nr(mic.sm_kfn), yerr= mic.sm_e2a_err/nuda.effg_nr(mic.sm_kfn), linestyle=mic.linestyle, label=mic.label, errorevery=mic.every )
+            elif matter.lower() == 'esym':
+                esym = nuda.matter.setupMicroEsym( model = model )
+                axs[0].errorbar( esym.den, esym.esym, yerr= esym.esym_err, linestyle=esym.linestyle, label=esym.label, errorevery=esym.every )
+                axs[1].errorbar( esym.den, esym.esym/nuda.esymffg_nr(esym.kf), yerr= esym.esym_err/nuda.esymffg_nr(esym.kf), linestyle=esym.linestyle, label=esym.label, errorevery=esym.every )
+    #
+    band = nuda.matter.setupMicroBand( models, den=den, matter=matter )
+    #
+    axs[0].fill_between( band.den, y1=band.e2a-band.e2a_std, y2=band.e2a+band.e2a_std, color=band.color, alpha=band.alpha )
+    if matter.lower() == 'nm':
+        axs[1].fill_between( band.den, y1=(band.e2a-band.e2a_std)/nuda.effg_nr(band.kfn), y2=(band.e2a+band.e2a_std)/nuda.effg_nr(band.kfn), color=band.color, alpha=band.alpha )
+    elif matter.lower() == 'sm':
+        axs[1].fill_between( band.den, y1=(band.e2a-band.e2a_std)/nuda.effg_nr(band.kf), y2=(band.e2a+band.e2a_std)/nuda.effg_nr(band.kf), color=band.color, alpha=band.alpha )
+    elif matter.lower() == 'esym':
+        axs[1].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 )
+    #
+    if matter.lower() == 'nm':
+        axs[1].legend(loc='upper left',fontsize='12', ncol=1)
+    elif matter.lower() == 'sm':
+        axs[1].legend(loc='upper left',fontsize='12', ncol=1)
+    elif matter.lower() == 'esym':
+        axs[0].legend(loc='upper left',fontsize='12', ncol=1)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()

nucleardatapy/fig/matter_setupMicro_effmass_fig.py

@@ -0,0 +1,70 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_setupMicro_effmass_fig( pname, models, matter = 'NM' ):
+    """
+    Plot the effective mass as function of the density and Fermi momentum.\
+    The plot is 1x2 with:\
+    [0]: effmass(den)
+    [1]: effmass(kF)
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param models: list of models to run on.
+    :type models: array of str.
+    :param matter: chose between 'SM' and 'NM' (default).
+    :type matter: str.
+    """
+    #
+    print(f'Plot name: {pname}')
+    print('models:',models)
+    print('matter:',matter)
+    #
+    fig, axs = plt.subplots(1,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.05, hspace=0.05 )
+    #
+    axs[0].set_ylabel(r'Landau effective mass $m^*/m$')
+    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs[0].set_xlim([0, 0.34])
+    axs[0].set_ylim([0.4, 1.2])
+    #axs[0].tick_params('x', labelbottom=False)
+    #
+    axs[1].set_xlim([0, 2.0])
+    axs[1].set_ylim([0.4, 1.2])
+    #axs[1].tick_params('x', labelbottom=False)
+    axs[1].tick_params('y', labelleft=False)
+    axs[1].set_xlabel(r'$k_{F_n}$ (fm$^{-1}$)')
+    #
+    #axs[1,0].set_ylabel(r'$\Delta_{1S0}/E_F$')
+    #
+    for model in models:
+        #
+        ms = nuda.matter.setupMicroEffmass( model = model, matter = matter )
+        #
+        if matter.lower() == 'nm':
+            if ms.nm_effmass is not None:
+                if ms.nm_effmass_err is not None:
+                    axs[0].errorbar( ms.nm_den, ms.nm_effmass, yerr=ms.nm_effmass_err, marker=ms.marker, markevery=ms.every, linestyle='none', label=ms.label )
+                    axs[1].errorbar( ms.nm_kfn, ms.nm_effmass, yerr=ms.nm_effmass_err, marker=ms.marker, markevery=ms.every, linestyle='none' )
+                else:
+                    axs[0].plot( ms.nm_den, ms.nm_effmass, marker=ms.marker, markevery=ms.every, linestyle='none', label=ms.label )
+                    axs[1].plot( ms.nm_kfn, ms.nm_effmass, marker=ms.marker, markevery=ms.every, linestyle='none' )
+        elif matter.lower() == 'sm':
+            if ms.sm_effmass is not None:
+                if ms.sm_effmass_err is not None:
+                    axs[0].errorbar( ms.sm_den, ms.sm_effmass, yerr=ms.sm_effmass_err, marker=ms.marker, markevery=ms.every, linestyle='none', label=ms.label )
+                    axs[1].errorbar( ms.sm_kfn, ms.sm_effmass, yerr=ms.sm_effmass_err, marker=ms.marker, markevery=ms.every, linestyle='none' )
+                else:
+                    axs[0].plot( ms.sm_den, ms.sm_effmass, marker=ms.marker, markevery=ms.every, linestyle='none', label=ms.label )
+                    axs[1].plot( ms.sm_kfn, ms.sm_effmass, marker=ms.marker, markevery=ms.every, linestyle='none' )
+        if nuda.env.verb_output: ms.print_outputs( )
+    #
+    #axs[1,0].legend(loc='upper right',fontsize='8')
+    fig.legend(loc='upper left',bbox_to_anchor=(0.08,1.0),columnspacing=2,fontsize='8',ncol=3,frameon=False)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=300)
+    	plt.close()

nucleardatapy/fig/matter_setupMicro_err_NM_fig.py

@@ -0,0 +1,47 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_setupMicro_err_NM_fig( pname, models ):
+    """
+    Plot uncertainties (err) estimated by different authors in NM.\
+    The plot is 1x1 with:\
+    [0]: uncertainty for E/A versus den.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param models: array of models.
+    :type models: array of str.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    den = 0.03*np.arange(13)
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.98, wspace=0.2, hspace=0.2 )
+    #
+    axs.set_xlabel(r'n (fm$^{-3}$)',fontsize='12')
+    axs.set_ylabel(r'$\Delta E_{NM}/E_{NM}$',fontsize='12')
+    axs.set_xlim([0, 0.4])
+    axs.set_ylim([0, 0.3])
+    #
+    for model in models:
+        #
+        mic = nuda.matter.setupMicro( model = model )
+        if nuda.env.verb_output: mic.print_outputs( )
+        if mic.nm_e2a is not None:
+            print('model:',model)
+            axs.plot( mic.nm_den, mic.nm_e2a_err/mic.nm_e2a, marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle, label=mic.label )
+    #axs.plot( den, nuda.matter.uncertainty_stat(den), linestyle='dashed', linewidth=3, label='fit (MBPT)' )
+    axs.plot( den, nuda.matter.uncertainty_stat(den,err='MBPT'), linestyle='dashed', linewidth=3, label='fit (MBPT)' )
+    axs.plot( den, nuda.matter.uncertainty_stat(den,err='QMC'),  linestyle='dashed', linewidth=3, label='fit (QMC)' )
+    #
+    axs.legend(loc='upper left',fontsize='12', ncol=3)
+    #
+    if pname is not None:
+        plt.savefig(pname, dpi=200)
+        plt.close()

nucleardatapy/fig/matter_setupMicro_fig.py

@@ -0,0 +1,165 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_setupMicro_E_fig( pname, mb, models, band, matter ):
+    """
+    Plot nucleonic energy per particle E/A in matter.\
+    The plot is 2x2 with:\
+    [0,0]: E/A versus den.       [0,1]: E/A versus kfn.\
+    [1,0]: E/E_NRFFG versus den. [1,1]: E/E_NRFFG versus kfn.\
+
+    :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}')
+    #
+    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.85, wspace=0.05, hspace=0.05 )
+    #
+    axs[0,0].set_ylabel(r'$E_\text{NM}/A$ (MeV)')
+    axs[0,0].set_xlim([0, 0.33])
+    axs[1,0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs[1,0].set_ylabel(r'$E_\text{NM}/E_\text{NRFFG}$')
+    axs[1,0].set_xlim([0, 0.33])
+    axs[1,1].set_xlabel(r'$k_{F_n}$ (fm$^{-1}$)')
+    axs[0,1].tick_params('y', labelleft=False)
+    axs[1,1].tick_params('y', labelleft=False)
+    axs[0,0].tick_params('x', labelbottom=False)
+    axs[0,1].tick_params('x', labelbottom=False)
+    if matter.lower() == 'nm':
+        axs[0,0].set_ylim([0, 35])
+        axs[0,1].set_xlim([0.5, 2.0])
+        axs[0,1].set_ylim([0, 35])
+        axs[1,0].set_ylim([0.3, 0.84])
+        axs[1,1].set_xlim([0.5, 2.0])
+        axs[1,1].set_ylim([0.3, 0.84])
+        delta = 1.0
+    elif matter.lower() == 'sm':
+        axs[0,0].set_ylim([-20, 10])
+        axs[0,1].set_xlim([0, 1.5])
+        axs[0,1].set_ylim([-20, 10])
+        axs[1,0].set_ylim([-2.0, 0.1])
+        axs[1,1].set_xlim([0, 1.5])
+        axs[1,1].set_ylim([-2.0, 0.1])
+        delta = 0.0
+    #
+    for model in models:
+        #
+        mic = nuda.matter.setupMicro( model = model, var2 = delta )
+        print('model:',model,' delta:',delta)
+        if mic.nm_e2a is not None and 'FIT' not in model:
+            if 'NLEFT' in model:
+                if matter.lower() == 'nm':
+                    axs[0,0].errorbar( mic.nm_den, mic.nm_e2adata, yerr=mic.nm_e2adata_err, linestyle = 'dotted', linewidth = 1, alpha=0.6 )
+                    axs[0,1].errorbar( mic.nm_kfn, mic.nm_e2adata, yerr=mic.nm_e2adata_err, linestyle = 'dotted', linewidth = 1, alpha=0.6 )
+                    axs[1,0].errorbar( mic.nm_den, mic.nm_e2adata/nuda.effg_nr(mic.nm_kfn), yerr=mic.nm_e2adata_err/nuda.effg_nr(mic.nm_kfn), linestyle = 'dotted', linewidth = 1, alpha=0.6 )
+                    axs[1,1].errorbar( mic.nm_kfn, mic.nm_e2adata/nuda.effg_nr(mic.nm_kfn), yerr=mic.nm_e2adata_err/nuda.effg_nr(mic.nm_kfn), linestyle = 'dotted', linewidth = 1, alpha=0.6 )
+                elif matter.lower() == 'sm':
+                    axs[0,0].errorbar( mic.sm_den, mic.sm_e2adata, yerr=mic.sm_e2adata_err, linestyle = 'dotted', linewidth = 1, alpha=0.6 )
+                    axs[0,1].errorbar( mic.sm_kfn, mic.sm_e2adata, yerr=mic.sm_e2adata_err, linestyle = 'dotted', linewidth = 1, alpha=0.6 )
+                    axs[1,0].errorbar( mic.sm_den, mic.sm_e2adata/nuda.effg_nr(mic.sm_kfn), yerr=mic.sm_e2adata_err/nuda.effg_nr(mic.sm_kfn), linestyle = 'dotted', linewidth = 1, alpha=0.6 )
+                    axs[1,1].errorbar( mic.sm_kfn, mic.sm_e2adata/nuda.effg_nr(mic.sm_kfn), yerr=mic.sm_e2adata_err/nuda.effg_nr(mic.sm_kfn), linestyle = 'dotted', linewidth = 1, alpha=0.6 )
+                #
+                if matter.lower() == 'nm':
+                    axs[0,0].fill_between( mic.nm_den, y1=(mic.nm_e2a-mic.nm_e2a_err), y2=(mic.nm_e2a+mic.nm_e2a_err), alpha=0.3 )
+                    axs[0,1].fill_between( mic.nm_kfn, y1=(mic.nm_e2a-mic.nm_e2a_err), y2=(mic.nm_e2a+mic.nm_e2a_err), alpha=0.3 )
+                    axs[1,0].fill_between( mic.nm_den, y1=(mic.nm_e2a-mic.nm_e2a_err)/nuda.effg_nr(mic.nm_kfn), y2=(mic.nm_e2a+mic.nm_e2a_err)/nuda.effg_nr(mic.nm_kfn), alpha=0.3 )
+                    axs[1,1].fill_between( mic.nm_kfn, y1=(mic.nm_e2a-mic.nm_e2a_err)/nuda.effg_nr(mic.nm_kfn), y2=(mic.nm_e2a+mic.nm_e2a_err)/nuda.effg_nr(mic.nm_kfn), alpha=0.3 )
+                elif matter.lower() == 'sm':
+                    axs[0,0].fill_between( mic.sm_den, y1=(mic.sm_e2a-mic.sm_e2a_err), y2=(mic.sm_e2a+mic.sm_e2a_err), alpha=0.3 )
+                    axs[0,1].fill_between( mic.sm_kfn, y1=(mic.sm_e2a-mic.sm_e2a_err), y2=(mic.sm_e2a+mic.sm_e2a_err), alpha=0.3 )
+                    axs[1,0].fill_between( mic.sm_den, y1=(mic.sm_e2a-mic.sm_e2a_err)/nuda.effg_nr(mic.sm_kfn), y2=(mic.sm_e2a+mic.sm_e2a_err)/nuda.effg_nr(mic.sm_kfn), alpha=0.3 )
+                    axs[1,1].fill_between( mic.sm_kfn, y1=(mic.sm_e2a-mic.sm_e2a_err)/nuda.effg_nr(mic.sm_kfn), y2=(mic.sm_e2a+mic.sm_e2a_err)/nuda.effg_nr(mic.sm_kfn), alpha=0.3 )
+            if mic.marker:
+                if mic.err:
+                    if matter.lower() == 'nm':
+                        axs[0,0].errorbar( mic.nm_den, mic.nm_e2a, yerr=mic.nm_e2a_err, marker=mic.marker, markevery=mic.every, linestyle='None', label=mic.label, errorevery=mic.every )
+                        axs[0,1].errorbar( mic.nm_kfn, mic.nm_e2a, yerr=mic.nm_e2a_err, marker=mic.marker, markevery=mic.every, linestyle='None', errorevery=mic.every )
+                        axs[1,0].errorbar( mic.nm_den, mic.nm_e2a/nuda.effg_nr(mic.nm_kfn), yerr=mic.nm_e2a_err/nuda.effg_nr(mic.nm_kfn), marker=mic.marker, markevery=mic.every, linestyle='None', errorevery=mic.every )
+                        axs[1,1].errorbar( mic.nm_kfn, mic.nm_e2a/nuda.effg_nr(mic.nm_kfn), yerr=mic.nm_e2a_err/nuda.effg_nr(mic.nm_kfn), marker=mic.marker, markevery=mic.every, linestyle='None', errorevery=mic.every )
+                    elif matter.lower() == 'sm':
+                        axs[0,0].errorbar( mic.sm_den, mic.sm_e2a, yerr=mic.sm_e2a_err, marker=mic.marker, markevery=mic.every, linestyle='None', label=mic.label, errorevery=mic.every )
+                        axs[0,1].errorbar( mic.sm_kfn, mic.sm_e2a, yerr=mic.sm_e2a_err, marker=mic.marker, markevery=mic.every, linestyle='None', errorevery=mic.every )
+                        axs[1,0].errorbar( mic.sm_den, mic.sm_e2a/nuda.effg_nr(mic.sm_kfn), yerr=mic.sm_e2a_err/nuda.effg_nr(mic.sm_kfn), marker=mic.marker, markevery=mic.every, linestyle='None', errorevery=mic.every )
+                        axs[1,1].errorbar( mic.sm_kfn, mic.sm_e2a/nuda.effg_nr(mic.sm_kfn), yerr=mic.sm_e2a_err/nuda.effg_nr(mic.sm_kfn), marker=mic.marker, markevery=mic.every, linestyle='None', errorevery=mic.every )
+                else:
+                    if matter.lower() == 'nm':
+                        axs[0,0].plot( mic.nm_den, mic.nm_e2a, marker=mic.marker, markevery=mic.every, linestyle='None', label=mic.label )
+                        axs[0,1].plot( mic.nm_kfn, mic.nm_e2a, marker=mic.marker, markevery=mic.every, linestyle='None' )
+                        axs[1,0].plot( mic.nm_den, mic.nm_e2a/nuda.effg_nr(mic.nm_kfn), marker=mic.marker, markevery=mic.every, linestyle='None' )
+                        axs[1,1].plot( mic.nm_kfn, mic.nm_e2a/nuda.effg_nr(mic.nm_kfn), marker=mic.marker, markevery=mic.every, linestyle='None' )
+                    elif matter.lower() == 'sm':
+                        axs[0,0].plot( mic.sm_den, mic.sm_e2a, marker=mic.marker, markevery=mic.every, linestyle='None', label=mic.label )
+                        axs[0,1].plot( mic.sm_kfn, mic.sm_e2a, marker=mic.marker, markevery=mic.every, linestyle='None' )
+                        axs[1,0].plot( mic.sm_den, mic.sm_e2a/nuda.effg_nr(mic.sm_kfn), marker=mic.marker, markevery=mic.every, linestyle='None' )
+                        axs[1,1].plot( mic.sm_kfn, mic.sm_e2a/nuda.effg_nr(mic.sm_kfn), marker=mic.marker, markevery=mic.every, linestyle='None' )
+            else:
+                if mic.err:
+                    if matter.lower() == 'nm':
+                        axs[0,0].errorbar( mic.nm_den, mic.nm_e2a, yerr=mic.nm_e2a_err, marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle, label=mic.label, errorevery=mic.every )
+                        axs[0,1].errorbar( mic.nm_kfn, mic.nm_e2a, yerr=mic.nm_e2a_err, marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle, errorevery=mic.every )
+                        axs[1,0].errorbar( mic.nm_den, mic.nm_e2a/nuda.effg_nr(mic.nm_kfn), yerr=mic.nm_e2a_err/nuda.effg_nr(mic.nm_kfn), marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle, errorevery=mic.every )
+                        axs[1,1].errorbar( mic.nm_kfn, mic.nm_e2a/nuda.effg_nr(mic.nm_kfn), yerr=mic.nm_e2a_err/nuda.effg_nr(mic.nm_kfn), marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle, errorevery=mic.every )
+                    elif matter.lower() == 'sm':
+                        axs[0,0].errorbar( mic.sm_den, mic.sm_e2a, yerr=mic.sm_e2a_err, marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle, label=mic.label, errorevery=mic.every )
+                        axs[0,1].errorbar( mic.sm_kfn, mic.sm_e2a, yerr=mic.sm_e2a_err, marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle, errorevery=mic.every )
+                        axs[1,0].errorbar( mic.sm_den, mic.sm_e2a/nuda.effg_nr(mic.sm_kfn), yerr=mic.sm_e2a_err/nuda.effg_nr(mic.sm_kfn), marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle, errorevery=mic.every )
+                        axs[1,1].errorbar( mic.sm_kfn, mic.sm_e2a/nuda.effg_nr(mic.sm_kfn), yerr=mic.sm_e2a_err/nuda.effg_nr(mic.sm_kfn), marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle, errorevery=mic.every )
+                else:
+                    if matter.lower() == 'nm':
+                        axs[0,0].plot( mic.nm_den, mic.nm_e2a, marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle, label=mic.label )
+                        axs[0,1].plot( mic.nm_kfn, mic.nm_e2a, marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle )
+                        axs[1,0].plot( mic.nm_den, mic.nm_e2a/nuda.effg_nr(mic.nm_kfn), marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle )
+                        axs[1,1].plot( mic.nm_kfn, mic.nm_e2a/nuda.effg_nr(mic.nm_kfn), marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle )
+                    elif matter.lower() == 'sm':
+                        axs[0,0].plot( mic.sm_den, mic.sm_e2a, marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle, label=mic.label )
+                        axs[0,1].plot( mic.sm_kfn, mic.sm_e2a, marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle )
+                        axs[1,0].plot( mic.sm_den, mic.sm_e2a/nuda.effg_nr(mic.sm_kfn), marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle )
+                        axs[1,1].plot( mic.sm_kfn, mic.sm_e2a/nuda.effg_nr(mic.sm_kfn), marker=mic.marker, markevery=mic.every, linestyle=mic.linestyle )
+        elif mic.e2a is not None: 
+            if 'fit' in model:
+                #print('den:',mic.den)
+                #print('e2a:',mic.e2a)                
+                axs[0,0].plot( mic.den, mic.e2a, marker=mic.marker, linestyle=mic.linestyle, markevery=mic.every, label=mic.label )
+                axs[0,1].plot( mic.kfn, mic.e2a, marker=mic.marker, linestyle=mic.linestyle )
+                axs[1,0].plot( mic.den, mic.e2a/nuda.effg_nr(mic.kfn), marker=mic.marker, linestyle=mic.linestyle )
+                axs[1,1].plot( mic.kfn, mic.e2a/nuda.effg_nr(mic.kfn), marker=mic.marker, linestyle=mic.linestyle )
+        #
+        if nuda.env.verb_output: mic.print_outputs( )
+        #
+    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' )
+    axs[0,0].plot( band.den, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
+    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' )
+    axs[0,1].plot( band.kfn, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
+    axs[1,0].fill_between( band.den, y1=(band.e2a-band.e2a_std)/nuda.effg_nr(band.kfn), y2=(band.e2a+band.e2a_std)/nuda.effg_nr(band.kfn), color=band.color, alpha=band.alpha, visible=True )
+    axs[1,0].plot( band.den, (band.e2a-band.e2a_std)/nuda.effg_nr(band.kfn), color='k', linestyle='dashed' )
+    axs[1,0].plot( band.den, (band.e2a+band.e2a_std)/nuda.effg_nr(band.kfn), color='k', linestyle='dashed' )
+    axs[1,1].fill_between( band.kfn, y1=(band.e2a-band.e2a_std)/nuda.effg_nr(band.kfn), y2=(band.e2a+band.e2a_std)/nuda.effg_nr(band.kfn), color=band.color, alpha=band.alpha, visible=True )
+    axs[1,1].plot( band.kfn, (band.e2a-band.e2a_std)/nuda.effg_nr(band.kfn), color='k', linestyle='dashed' )
+    axs[1,1].plot( band.kfn, (band.e2a+band.e2a_std)/nuda.effg_nr(band.kfn), color='k', linestyle='dashed' )
+    #
+    #axs[0,1].legend(loc='upper left',fontsize='8', ncol=2)
+    #if mb not in 'BHF':
+    #    axs[0,1].legend(loc='upper left',fontsize='8', ncol=2)
+    #
+    #plt.tight_layout(pad=3.0)
+    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()