nucleardatapy 0.2.0__py3-none-any.whl

nucleardatapy/fig/matter_ENM_fig.py

@@ -0,0 +1,119 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_ENM_fig( pname, micro_mbs, pheno_models, band ):
+    """
+    Plot nucleonic energy per particle E/A in NM.\
+    The plot is 1x2 with:\
+    [0,0]: E/A versus den (micro). [0,1]: E/A versus den (pheno).\
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param micro_mbs: many-body (mb) approach considered.
+    :type micro_mbs: str.
+    :param pheno_models: models to run on.
+    :type pheno_models: array of str.
+    :param band: object instantiated on the reference band.
+    :type band: object.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    fig, axs = plt.subplots(1,2)
+    #fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.10, bottom=0.12, right=None, top=0.9, wspace=0.05, hspace=0.3 )
+    #
+    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs[0].set_ylabel(r'$e_\text{NM}(n)$')
+    axs[0].set_xlim([0, 0.33])
+    axs[0].set_ylim([0, 30])
+    #
+    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    #axs[1].set_ylabel(r'$e_{sym}(n)$')
+    axs[1].set_xlim([0, 0.33])
+    axs[1].set_ylim([0, 30])
+    axs[1].tick_params('y', labelleft=False)
+    #
+    mb_check = []
+    k = 0
+    #
+    for mb in micro_mbs:
+        #
+        models, models_lower = nuda.matter.micro_models_mb( mb )
+        #
+        for model in models:
+            #
+            enm = nuda.matter.setupMicro( model = model )
+            #
+            if enm.nm_e2a is not None:
+                print('mb:',mb,'model:',model)
+                if mb in mb_check:
+                    if enm.marker:
+                        if enm.err:
+                            axs[0].errorbar( enm.nm_den, enm.nm_e2a, yerr=enm.nm_e2a_err, marker=enm.marker, linestyle=None, errorevery=enm.every, color=nuda.param.col[k] )
+                        else:
+                            axs[0].plot( enm.nm_den, enm.nm_e2a, marker=enm.marker, linestyle=None, markevery=enm.every, color=nuda.param.col[k] )
+                    else:
+                        if enm.err:
+                            axs[0].errorbar( enm.nm_den, enm.nm_e2a, yerr=enm.nm_e2a_err, marker=enm.marker, linestyle=enm.linestyle, errorevery=enm.every, color=nuda.param.col[k] )
+                        else:
+                            axs[0].plot( enm.nm_den, enm.nm_e2a, marker=enm.marker, linestyle=enm.linestyle, markevery=enm.every, color=nuda.param.col[k] )
+                else:
+                    mb_check.append(mb)
+                    k += 1
+                    if enm.marker:
+                        if enm.err:
+                            axs[0].errorbar( enm.nm_den, enm.nm_e2a, yerr=enm.nm_e2a_err, marker=enm.marker, linestyle=None, label=mb, errorevery=enm.every, color=nuda.param.col[k] )
+                        else:
+                            axs[0].plot( enm.nm_den, enm.nm_e2a, marker=enm.marker, linestyle=None, label=mb, markevery=enm.every, color=nuda.param.col[k] )
+                    else:
+                        if enm.err:
+                            axs[0].errorbar( enm.nm_den, enm.nm_e2a, yerr=enm.nm_e2a_err, marker=enm.marker, linestyle=enm.linestyle, label=mb, errorevery=enm.every, color=nuda.param.col[k] )
+                        else:
+                            axs[0].plot( enm.nm_den, enm.nm_e2a, marker=enm.marker, linestyle=enm.linestyle, label=mb, markevery=enm.every, color=nuda.param.col[k] )
+                    #axs[0].plot( esym.den, esym.esym, color=nuda.param.col[k], label=mb )
+            if nuda.env.verb: esm.print_outputs( )
+    axs[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].plot( band.den, (band.e2a-band.e2a_std), color='k', linestyle='dashed' )
+    axs[0].plot( band.den, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
+    axs[0].text(0.03,2,'microscopic models',fontsize='10')
+    #axs[0].legend(loc='upper left',fontsize='8', ncol=3)
+    #
+    model_check = []
+    k = 0
+    #
+    for model in pheno_models:
+        #
+        params, params_lower = nuda.matter.pheno_params( model = model )
+        #
+        for param in params:
+            #
+            enm = nuda.matter.setupPheno( model = model, param = param )
+            #
+            if enm.nm_e2a is not None: 
+                print('model:',model,' param:',param)
+                if model in model_check:
+                    axs[1].plot( enm.nm_den, enm.nm_e2a, color=nuda.param.col[k] )
+                else:
+                    model_check.append(model)
+                    k += 1
+                    axs[1].plot( enm.nm_den, enm.nm_e2a, color=nuda.param.col[k], label=model )
+                #pheno.label=None
+                #axs[1].plot( esym.den, esym.esym, label=esym.label )
+            if nuda.env.verb: esym.print_outputs( )
+    axs[1].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[1].plot( band.den, (band.e2a-band.e2a_std), color='k', linestyle='dashed' )
+    axs[1].plot( band.den, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
+    axs[1].text(0.03,2,'phenomenological models',fontsize='10')
+    #axs[1].legend(loc='upper left',fontsize='8', ncol=2)
+    #axs[0,1].legend(loc='upper left',fontsize='xx-small', ncol=2)
+    fig.legend(loc='upper left',bbox_to_anchor=(0.15,1.0),columnspacing=2,fontsize='8',ncol=5,frameon=False)
+    #
+    #plt.tight_layout()
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+    #

nucleardatapy/fig/matter_ESM_fig.py

@@ -0,0 +1,119 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_ESM_fig( pname, micro_mbs, pheno_models, band ):
+    """
+    Plot nucleonic energy per particle E/A in matter.\
+    The plot is 1x2 with:\
+    [0,0]: E/A versus den (micro). [0,1]: E/A versus den (pheno).\
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param micro_mbs: many-body (mb) approach considered.
+    :type micro_mbs: str.
+    :param pheno_models: models to run on.
+    :type pheno_models: array of str.
+    :param band: object instantiated on the reference band.
+    :type band: object.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    fig, axs = plt.subplots(1,2)
+    #fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.10, bottom=0.12, right=None, top=0.9, wspace=0.05, hspace=0.3 )
+    #
+    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs[0].set_ylabel(r'$e_\text{SM}(n)$')
+    axs[0].set_xlim([0, 0.33])
+    axs[0].set_ylim([-25, 5])
+    #
+    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    #axs[1].set_ylabel(r'$e_{sym}(n)$')
+    axs[1].set_xlim([0, 0.33])
+    axs[1].set_ylim([-25, 5])
+    axs[1].tick_params('y', labelleft=False)
+    #
+    mb_check = []
+    k = 0
+    #
+    for mb in micro_mbs:
+        #
+        models, models_lower = nuda.matter.micro_models_mb( mb )
+        #
+        for model in models:
+            #
+            esm = nuda.matter.setupMicro( model = model )
+            #
+            if esm.sm_e2a is not None:
+                print('mb:',mb,'model:',model)
+                if mb in mb_check:
+                    if esm.marker:
+                        if esm.err:
+                            axs[0].errorbar( esm.sm_den, esm.sm_e2a, yerr=esm.sm_e2a_err, marker=esm.marker, linestyle=None, errorevery=esm.every, color=nuda.param.col[k] )
+                        else:
+                            axs[0].plot( esm.sm_den, esm.sm_e2a, marker=esm.marker, linestyle=None, markevery=esm.every, color=nuda.param.col[k] )
+                    else:
+                        if esm.err:
+                            axs[0].errorbar( esm.sm_den, esm.sm_e2a, yerr=esm.sm_e2a_err, marker=esm.marker, linestyle=esm.linestyle, errorevery=esm.every, color=nuda.param.col[k] )
+                        else:
+                            axs[0].plot( esm.sm_den, esm.sm_e2a, marker=esm.marker, linestyle=esm.linestyle, markevery=esm.every, color=nuda.param.col[k] )
+                else:
+                    mb_check.append(mb)
+                    k += 1
+                    if esm.marker:
+                        if esm.err:
+                            axs[0].errorbar( esm.sm_den, esm.sm_e2a, yerr=esm.sm_e2a_err, marker=esm.marker, linestyle=None, label=mb, errorevery=esm.every, color=nuda.param.col[k] )
+                        else:
+                            axs[0].plot( esm.sm_den, esm.sm_e2a, marker=esm.marker, linestyle=None, label=mb, markevery=esm.every, color=nuda.param.col[k] )
+                    else:
+                        if esm.err:
+                            axs[0].errorbar( esm.sm_den, esm.sm_e2a, yerr=esm.sm_e2a_err, marker=esm.marker, linestyle=esm.linestyle, label=mb, errorevery=esm.every, color=nuda.param.col[k] )
+                        else:
+                            axs[0].plot( esm.sm_den, esm.sm_e2a, marker=esm.marker, linestyle=esm.linestyle, label=mb, markevery=esm.every, color=nuda.param.col[k] )
+                    #axs[0].plot( esym.den, esym.esym, color=nuda.param.col[k], label=mb )
+            if nuda.env.verb: esm.print_outputs( )
+    axs[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].plot( band.den, (band.e2a-band.e2a_std), color='k', linestyle='dashed' )
+    axs[0].plot( band.den, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
+    axs[0].text(0.03,2,'microscopic models',fontsize='10')
+    #axs[0].legend(loc='upper left',fontsize='8', ncol=3)
+    #
+    model_check = []
+    k = 0
+    #
+    for model in pheno_models:
+        #
+        params, params_lower = nuda.matter.pheno_params( model = model )
+        #
+        for param in params:
+            #
+            esm = nuda.matter.setupPheno( model = model, param = param )
+            #
+            if esm.sm_e2a is not None: 
+                print('model:',model,' param:',param)
+                if model in model_check:
+                    axs[1].plot( esm.sm_den, esm.sm_e2a, color=nuda.param.col[k] )
+                else:
+                    model_check.append(model)
+                    k += 1
+                    axs[1].plot( esm.sm_den, esm.sm_e2a, color=nuda.param.col[k], label=model )
+                #pheno.label=None
+                #axs[1].plot( esym.den, esym.esym, label=esym.label )
+            if nuda.env.verb: esym.print_outputs( )
+    axs[1].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[1].plot( band.den, (band.e2a-band.e2a_std), color='k', linestyle='dashed' )
+    axs[1].plot( band.den, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
+    axs[1].text(0.03,2,'phenomenological models',fontsize='10')
+    #axs[1].legend(loc='upper left',fontsize='8', ncol=2)
+    #axs[0,1].legend(loc='upper left',fontsize='xx-small', ncol=2)
+    fig.legend(loc='upper left',bbox_to_anchor=(0.15,1.0),columnspacing=2,fontsize='8',ncol=4,frameon=False)
+    #
+    #plt.tight_layout()
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+    #

nucleardatapy/fig/matter_Esym_fig.py

@@ -0,0 +1,122 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_Esym_fig( pname, micro_mbs, pheno_models, band ):
+    """
+    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(1,2)
+    #fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.10, bottom=0.12, right=None, top=0.9, wspace=0.05, hspace=0.3 )
+    #
+    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs[0].set_ylabel(r'$e_\text{sym}(n)$')
+    axs[0].set_xlim([0, 0.33])
+    axs[0].set_ylim([0, 60])
+    #
+    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    #axs[1].set_ylabel(r'$e_{sym}(n)$')
+    axs[1].set_xlim([0, 0.33])
+    axs[1].set_ylim([0, 60])
+    axs[1].tick_params('y', labelleft=False)
+    #
+    mb_check = []
+    k = 0
+    #
+    for mb in micro_mbs:
+        #
+        models, models_lower = nuda.matter.micro_esym_models_mb( mb )
+        #
+        for model in models:
+            #
+            esym = nuda.matter.setupMicroEsym( model = model )
+            #
+            if esym.esym is not None:
+                print('mb:',mb,'model:',model)
+                if mb in mb_check:
+                    if esym.marker:
+                        if esym.err:
+                            axs[0].errorbar( esym.den, esym.esym, yerr=esym.esym_err, marker=esym.marker, linestyle=None, errorevery=esym.every, color=nuda.param.col[k] )
+                        else:
+                            axs[0].plot( esym.den, esym.esym, marker=esym.marker, linestyle=None, markevery=esym.every, color=nuda.param.col[k] )
+                    else:
+                        if esym.err:
+                            axs[0].errorbar( esym.den, esym.esym, yerr=esym.esym_err, marker=esym.marker, linestyle=esym.linestyle, errorevery=esym.every, color=nuda.param.col[k] )
+                        else:
+                            axs[0].plot( esym.den, esym.esym, marker=esym.marker, linestyle=esym.linestyle, markevery=esym.every, color=nuda.param.col[k] )
+                else:
+                    mb_check.append(mb)
+                    k += 1
+                    if esym.marker:
+                        if esym.err:
+                            axs[0].errorbar( esym.den, esym.esym, yerr=esym.esym_err, marker=esym.marker, linestyle=None, label=mb, errorevery=esym.every, color=nuda.param.col[k] )
+                        else:
+                            axs[0].plot( esym.den, esym.esym, marker=esym.marker, linestyle=None, label=mb, markevery=esym.every, color=nuda.param.col[k] )
+                    else:
+                        if esym.err:
+                            axs[0].errorbar( esym.den, esym.esym, yerr=esym.esym_err, marker=esym.marker, linestyle=esym.linestyle, label=mb, errorevery=esym.every, color=nuda.param.col[k] )
+                        else:
+                            axs[0].plot( esym.den, esym.esym, marker=esym.marker, linestyle=esym.linestyle, label=mb, markevery=esym.every, color=nuda.param.col[k] )
+                    #axs[0].plot( esym.den, esym.esym, color=nuda.param.col[k], label=mb )
+            if nuda.env.verb: esym.print_outputs( )
+    axs[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].plot( band.den, (band.e2a-band.e2a_std), color='k', linestyle='dashed' )
+    axs[0].plot( band.den, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
+    axs[0].text(0.05,5,'microscopic models',fontsize='10')
+    #axs[0].legend(loc='upper left',fontsize='8', ncol=3)
+    #
+    model_check = []
+    k = 0
+    #
+    for model in pheno_models:
+        #
+        params, params_lower = nuda.matter.pheno_params( model = model )
+        #
+        for param in params:
+            #
+            esym = nuda.matter.setupPhenoEsym( model = model, param = param )
+            #
+            if esym.esym is not None: 
+                print('model:',model,' param:',param)
+                if model in model_check:
+                    axs[1].plot( esym.den, esym.esym, color=nuda.param.col[k] )
+                else:
+                    model_check.append(model)
+                    k += 1
+                    axs[1].plot( esym.den, esym.esym, color=nuda.param.col[k], label=model )
+                #pheno.label=None
+                #axs[1].plot( esym.den, esym.esym, label=esym.label )
+            if nuda.env.verb: esym.print_outputs( )
+    axs[1].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[1].plot( band.den, (band.e2a-band.e2a_std), color='k', linestyle='dashed' )
+    axs[1].plot( band.den, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
+    axs[1].text(0.05,5,'phenomenological models',fontsize='10')
+    #axs[1].legend(loc='upper left',fontsize='8', ncol=2)
+    #axs[0,1].legend(loc='upper left',fontsize='xx-small', ncol=2)
+    fig.legend(loc='upper left',bbox_to_anchor=(0.2,1.0),columnspacing=2,fontsize='8',ncol=4,frameon=False)
+    #
+    #plt.tight_layout()
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+    #

nucleardatapy/fig/matter_setupFFGNuc_fig.py

@@ -0,0 +1,172 @@
+
+import numpy as np
+import matplotlib.pyplot as plt
+
+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) ):
+    """
+    Plot nucleonic 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}')
+    #
+    denkfn = nuda.matter.den( kfn )
+    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 )
+        #
+        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 )
+            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)')
+            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 )
+            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].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)
+    #
+    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) ):
+    """
+    Plot nucleonic FFG EOS in NM and SM.\
+    The plot is 1x2 with:\
+    [0]: EOS (pre) versus energy density rho.\
+    [1]: Sound speed c_s^2 versus energy density rho.\
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param den: density.
+    :type den: float or numpy vector of real numbers.
+
+    """
+    #
+    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 )
+    #
+    axs[0].set_ylabel(r'$p_\text{FFG}$')
+    axs[0].set_xlim([0, 360])
+    axs[0].set_ylim([0, 20])
+    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_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 )
+        #
+        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 )
+            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)')
+            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 )
+            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[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)
+    #
+    if pname is not None: 
+        plt.savefig(pname, dpi=300)
+        plt.close()
+    #

nucleardatapy/fig/matter_setupHIC_fig.py

@@ -0,0 +1,101 @@
+import numpy as np
+import matplotlib.pyplot as plt
+#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']
+
+import nucleardatapy as nuda
+
+def matter_setupHIC_fig( pname, constraints ):
+    """
+    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.10, bottom=0.12, right=None, 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_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_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].tick_params('x', labelbottom=False)
+    axs[0,1].tick_params('x', labelbottom=False)
+    #
+    for constraint in constraints:
+        #
+        print('constraint:',constraint)
+        hic = nuda.matter.setupHIC( constraint = constraint )
+        #
+        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 )
+        #
+        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.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)
+        #
+        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.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)
+        #
+    #
+    #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()
+    #