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

nucleardatapy/fig/matter_Esym_fig.py

@@ -25,97 +25,109 @@ def matter_Esym_fig( pname, micro_mbs, pheno_models, band ):
     #
     print(f'Plot name: {pname}')
     #
+    matter = 'Esym'
+    #
     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_ylabel(r'$e_\text{sym}(n_\text{nuc})$')
+    axs[0].set_xlim([0, 0.34])
     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_xlim([0, 0.34])
     axs[1].set_ylim([0, 60])
     axs[1].tick_params('y', labelleft=False)
     #
     mb_check = []
-    k = 0
     #
-    for mb in micro_mbs:
+    for kmb,mb in enumerate(micro_mbs):
         #
         models, models_lower = nuda.matter.micro_esym_models_mb( mb )
         #
         for model in models:
             #
-            esym = nuda.matter.setupMicroEsym( model = model )
+            micro = nuda.matter.setupMicroEsym( model = model )
+            if nuda.env.verb: micro.print_outputs( )
+            #
+            check = nuda.matter.setupCheck( eos = micro, band = band )
+            #
+            if check.isInside:
+                lstyle = 'solid'
+            else:
+                lstyle = 'dashed'
             #
-            if esym.esym is not None:
+            if micro.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] )
+                    if micro.marker:
+                        if micro.err:
+                            axs[0].errorbar( micro.den, micro.esym, yerr=micro.esym_err, marker=micro.marker, linestyle=lstyle, errorevery=micro.every, color=nuda.param.col[kmb] )
                         else:
-                            axs[0].plot( esym.den, esym.esym, marker=esym.marker, linestyle=None, markevery=esym.every, color=nuda.param.col[k] )
+                            axs[0].plot( micro.den, micro.esym, marker=micro.marker, linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
                     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] )
+                        if micro.err:
+                            axs[0].errorbar( micro.den, micro.esym, yerr=micro.esym_err, marker=micro.marker, linestyle=lstyle, errorevery=micro.every, color=nuda.param.col[kmb] )
                         else:
-                            axs[0].plot( esym.den, esym.esym, marker=esym.marker, linestyle=esym.linestyle, markevery=esym.every, color=nuda.param.col[k] )
+                            axs[0].plot( micro.den, micro.esym, marker=micro.marker, linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
                 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] )
+                    if micro.marker:
+                        if micro.err:
+                            axs[0].errorbar( micro.den, micro.esym, yerr=micro.esym_err, marker=micro.marker, linestyle=lstyle, label=mb, errorevery=micro.every, color=nuda.param.col[kmb] )
                         else:
-                            axs[0].plot( esym.den, esym.esym, marker=esym.marker, linestyle=None, label=mb, markevery=esym.every, color=nuda.param.col[k] )
+                            axs[0].plot( micro.den, micro.esym, marker=micro.marker, linestyle=lstyle, label=mb, markevery=micro.every, color=nuda.param.col[kmb] )
                     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] )
+                        if micro.err:
+                            axs[0].errorbar( micro.den, micro.esym, yerr=micro.esym_err, marker=micro.marker, linestyle=lstyle, label=mb, errorevery=micro.every, color=nuda.param.col[kmb] )
                         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].plot( micro.den, micro.esym, marker=micro.marker, linestyle=lstyle, label=mb, markevery=micro.every, color=nuda.param.col[kmb] )
+            # end of model
+        # end of mb
     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:
+    for kmodel,model in enumerate(pheno_models):
         #
         params, params_lower = nuda.matter.pheno_params( model = model )
         #
         for param in params:
             #
-            esym = nuda.matter.setupPhenoEsym( model = model, param = param )
+            pheno = nuda.matter.setupPhenoEsym( model = model, param = param )
+            if nuda.env.verb: pheno.print_outputs( )
             #
-            if esym.esym is not None: 
+            check = nuda.matter.setupCheck( eos = pheno, band = band )
+            #
+            if check.isInside:
+                lstyle = 'solid'
+            else:
+                lstyle = 'dashed'
+            #
+            if pheno.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] )
+                    axs[1].plot( pheno.den, pheno.esym, linestyle=lstyle, color=nuda.param.col[kmodel] )
                 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].plot( pheno.den, pheno.esym, linestyle=lstyle, color=nuda.param.col[kmodel], label=model )
+            # end of param
+        # end of model
     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_cs2_fig.py

@@ -0,0 +1,83 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_cs2_fig( pname, micro_models, pheno_models, band ):
+    """
+    Plot nuclear chart (N versus Z).\
+    The plot is 1x1 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(1,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.3, hspace=0.3 )
+    #
+    axs[0].set_xlabel(r'n (fm$^{-3}$)')
+    axs[0].set_ylabel(r'$c_{s,NM}^2(n)$')
+    axs[0].set_xlim([0, 0.3])
+    axs[0].set_ylim([0, 0.5])
+    #
+    axs[1].set_xlabel(r'n (fm$^{-3}$)')
+    axs[1].set_ylabel(r'$c_{s,NM}^2(n)$')
+    axs[1].set_xlim([0, 0.3])
+    axs[1].set_ylim([0, 0.5])
+    #
+    for model in micro_models:
+        #
+        mic = nuda.matter.setupMicro( model = model )
+        if mic.nm_cs2 is not None: 
+            print('model:',model)
+            if mic.marker:
+                if mic.err:
+                    axs[0].errorbar( mic.nm_den, mic.nm_cs2, yerr=mic.nm_cs2_err, marker=mic.marker, linestyle=None, label=mic.label, errorevery=mic.every )
+                else:
+                    axs[0].plot( mic.nm_den, mic.nm_cs2, marker=mic.marker, linestyle=None, label=mic.label, markevery=mic.every )
+            else:
+                if mic.err:
+                    axs[0].errorbar( mic.nm_den, mic.nm_cs2, yerr=mic.nm_cs2_err, marker=mic.marker, linestyle=mic.linestyle, label=mic.label, errorevery=mic.every )
+                else:
+                    axs[0].plot( mic.nm_den, mic.nm_cs2, marker=mic.marker, linestyle=mic.linestyle, label=mic.label, markevery=mic.every )
+        if nuda.env.verb: mic.print_outputs( )
+    #axs[0].fill_between( band.den, y1=(band.pre-band.pre_std), y2=(band.pre+band.pre_std), color=band.color, alpha=band.alpha, visible=True )
+    #axs[0].plot( band.den, (band.pre-band.pre_std), color='k', linestyle='dashed' )
+    #axs[0].plot( band.den, (band.pre+band.pre_std), color='k', linestyle='dashed' )
+    axs[0].text(0.01,0.4,'microscopic models',fontsize='10')
+    axs[0].legend(loc='upper left',fontsize='8', ncol=3)
+    #
+    for model in pheno_models:
+        #
+        params, params_lower = nuda.matter.pheno_params( model = model )
+        #
+        for param in params:
+            #
+            pheno = nuda.matter.setupPheno( model = model, param = param )
+            if pheno.nm_pre is not None: 
+                print('model:',model,' param:',param)
+                #pheno.label=None
+                axs[1].plot( pheno.nm_den, pheno.nm_cs2, label=pheno.label )
+            if nuda.env.verb: pheno.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.01,0.4,'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)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+    #

nucleardatapy/fig/matter_preNM_fig.py

@@ -0,0 +1,146 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_preNM_fig( pname, micro_mbs, pheno_models, band ):
+    """
+    Plot nucleonic pressure 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}')
+    #
+    p_den = 0.32
+    p_cen = 23.0
+    p_std = 14.5 
+    p_micro_cen = 15.0
+    p_micro_std =  6.5
+    p_pheno_cen = 23.0
+    p_pheno_std = 14.5
+    #
+    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'$p_\text{NM}(n_\text{nuc})$')
+    axs[0].set_xlim([0, 0.35])
+    axs[0].set_ylim([-2, 45])
+    #
+    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    #axs[1].set_ylabel(r'$e_{sym}(n)$')
+    axs[1].set_xlim([0, 0.35])
+    axs[1].set_ylim([-2, 45])
+    axs[1].tick_params('y', labelleft=False)
+    #
+    mb_check = []
+    #
+    for kmb,mb in enumerate(micro_mbs): 
+        #
+        models, models_lower = nuda.matter.micro_models_mb( mb )
+        #
+        for model in models:
+            #
+            micro = nuda.matter.setupMicro( model = model )
+            if nuda.env.verb: micro.print_outputs( )
+            #
+            check = nuda.matter.setupCheck( eos = micro, band = band )
+            #
+            if check.isInside:
+                lstyle = 'solid'
+            else:
+                lstyle = 'dashed'
+                #continue
+            #
+            print('model:',model)
+            print('err:',micro.p_err)
+            print('den:',micro.nm_den)
+            print('pre:',micro.nm_pre)
+            print('pre_err:',micro.nm_pre_err)
+            if micro.nm_pre is not None:
+                print('mb:',mb,'model:',model)
+                if mb in mb_check:
+                    if micro.marker:
+                        if micro.p_err:
+                            axs[0].errorbar( micro.nm_den, micro.nm_pre, yerr=micro.nm_pre_err, marker=micro.marker, linestyle=lstyle, errorevery=micro.every, color=nuda.param.col[kmb] )
+                        else:
+                            axs[0].plot( micro.nm_den, micro.nm_pre, marker=micro.marker, linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
+                    else:
+                        if micro.p_err:
+                            axs[0].errorbar( micro.nm_den, micro.nm_pre, yerr=micro.nm_pre_err, marker=micro.marker, linestyle=lstyle, errorevery=micro.every, color=nuda.param.col[kmb] )
+                        else:
+                            axs[0].plot( micro.nm_den, micro.nm_pre, marker=micro.marker, linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
+                else:
+                    mb_check.append(mb)
+                    if micro.marker:
+                        if micro.p_err:
+                            axs[0].errorbar( micro.nm_den, micro.nm_pre, yerr=micro.nm_pre_err, marker=micro.marker, linestyle=lstyle, label=mb, errorevery=micro.every, color=nuda.param.col[kmb] )
+                        else:
+                            axs[0].plot( micro.nm_den, micro.nm_pre, marker=micro.marker, linestyle=lstyle, label=mb, markevery=micro.every, color=nuda.param.col[kmb] )
+                    else:
+                        if micro.p_err:
+                            axs[0].errorbar( micro.nm_den, micro.nm_pre, yerr=micro.nm_pre_err, marker=micro.marker, linestyle=lstyle, label=mb, errorevery=micro.every, color=nuda.param.col[kmb] )
+                        else:
+                            axs[0].plot( enm.nm_den, enm.nm_pre, marker=enm.marker, linestyle=lstyle, label=mb, markevery=micro.every, color=nuda.param.col[kmb] )
+            # end of model
+        # end of mb
+    axs[0].errorbar( p_den, p_cen, yerr=p_std, color='k' )
+    axs[0].errorbar( p_den+0.005, p_micro_cen, yerr=p_micro_std, color='r' )
+    axs[0].text(0.02,40,'microscopic models',fontsize='10')
+    #
+    model_check = []
+    #
+    for kmodel,model in enumerate(pheno_models):
+        #
+        params, params_lower = nuda.matter.pheno_params( model = model )
+        #
+        for param in params:
+            #
+            pheno = nuda.matter.setupPheno( model = model, param = param )
+            if nuda.env.verb: pheno.print_outputs( )
+            #
+            check = nuda.matter.setupCheck( eos = pheno, band = band )
+            #
+            if check.isInside:
+                lstyle = 'solid'
+            else:
+                lstyle = 'dashed'
+                #continue
+            #
+            if pheno.nm_pre is not None: 
+                print('model:',model,' param:',param)
+                if model in model_check:
+                    axs[1].plot( pheno.nm_den, pheno.nm_pre, linestyle=lstyle, color=nuda.param.col[kmodel] )
+                else:
+                    model_check.append(model)
+                    axs[1].plot( pheno.nm_den, pheno.nm_pre, linestyle=lstyle, color=nuda.param.col[kmodel], label=model )
+            # end of param
+        # end of model
+    axs[1].errorbar( p_den, p_cen, yerr=p_std, color='k' )
+    axs[1].errorbar( p_den+0.005, p_pheno_cen, yerr=p_pheno_std, color='r' )
+    #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.02,40,'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_preSM_fig.py

@@ -0,0 +1,144 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_preSM_fig( pname, micro_mbs, pheno_models, band ):
+    """
+    Plot nucleonic pressure in SM.\
+    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.
+    :param matter: variable `matter`employed to define the band.
+    :type matter: str.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    p_den = 0.32
+    p_cen = 11.5
+    p_std =  5.5
+    p_micro_cen =  9.0
+    p_micro_std =  3.0
+    p_pheno_cen = 14.5
+    p_pheno_std =  2.5
+    #
+    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'$p_\text{SM}(n_\text{nuc})$')
+    axs[0].set_xlim([0, 0.35])
+    axs[0].set_ylim([-2, 45])
+    #
+    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    #axs[1].set_ylabel(r'$e_{sym}(n)$')
+    axs[1].set_xlim([0, 0.35])
+    axs[1].set_ylim([-2, 45])
+    axs[1].tick_params('y', labelleft=False)
+    #
+    mb_check = []
+    #
+    for kmb,mb in enumerate(micro_mbs):
+        #
+        models, models_lower = nuda.matter.micro_models_mb( mb )
+        #
+        for model in models:
+            #
+            micro = nuda.matter.setupMicro( model = model )
+            if nuda.env.verb: micro.print_outputs( )
+            #
+            check = nuda.matter.setupCheck( eos = micro, band = band )
+            #
+            if check.isInside:
+                lstyle = 'solid'
+            else:
+                lstyle = 'dashed'
+                #continue
+            #
+            print('model:',model)
+            print('err:',micro.p_err)
+            print('den:',micro.sm_den)
+            print('pre:',micro.sm_pre)
+            print('pre_err:',micro.sm_pre_err)
+            if micro.sm_pre is not None:
+                print('mb:',mb,'model:',model)
+                if mb in mb_check:
+                    if micro.marker:
+                        if micro.p_err:
+                            axs[0].errorbar( micro.sm_den, micro.sm_pre, yerr=micro.sm_pre_err, marker=micro.marker, linestyle=lstyle, errorevery=micro.every, color=nuda.param.col[kmb] )
+                        else:
+                            axs[0].plot( micro.sm_den, micro.sm_pre, marker=micro.marker, linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
+                    else:
+                        if micro.p_err:
+                            axs[0].errorbar( micro.sm_den, micro.sm_pre, yerr=micro.sm_pre_err, marker=micro.marker, linestyle=lstyle, errorevery=micro.every, color=nuda.param.col[kmb] )
+                        else:
+                            axs[0].plot( micro.sm_den, micro.sm_pre, marker=micro.marker, linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
+                else:
+                    mb_check.append(mb)
+                    if micro.marker:
+                        if micro.p_err:
+                            axs[0].errorbar( esm.sm_den, micro.sm_pre, yerr=micro.sm_pre_err, marker=micro.marker, linestyle=lstyle, label=mb, errorevery=micro.every, color=nuda.param.col[kmb] )
+                        else:
+                            axs[0].plot( micro.sm_den, micro.sm_pre, marker=micro.marker, linestyle=lstyle, label=mb, markevery=micro.every, color=nuda.param.col[kmb] )
+                    else:
+                        if micro.p_err:
+                            axs[0].errorbar( micro.sm_den, micro.sm_pre, yerr=micro.sm_pre_err, marker=micro.marker, linestyle=lstyle, label=mb, errorevery=micro.every, color=nuda.param.col[kmb] )
+                        else:
+                            axs[0].plot( micro.sm_den, esm.sm_pre, marker=esm.marker, linestyle=lstyle, label=mb, markevery=micro.every, color=nuda.param.col[kmb] )
+            # end of model
+        # end of mb
+    axs[0].errorbar( p_den, p_cen, yerr=p_std, color='k' )
+    axs[0].errorbar( p_den+0.005, p_micro_cen, yerr=p_micro_std, color='r' )
+    axs[0].text(0.02,40,'microscopic models',fontsize='10')
+    #
+    model_check = []
+    #
+    for kmodel,model in enumerate(pheno_models):
+        #
+        params, params_lower = nuda.matter.pheno_params( model = model )
+        #
+        for param in params:
+            #
+            pheno = nuda.matter.setupPheno( model = model, param = param )
+            if nuda.env.verb: pheno.print_outputs( )
+            #
+            check = nuda.matter.setupCheck( eos = pheno, band = band )
+            #
+            if check.isInside:
+                lstyle = 'solid'
+            else:
+                lstyle = 'dashed'
+                #continue
+            #
+            if pheno.sm_pre is not None: 
+                print('model:',model,' param:',param)
+                if model in model_check:
+                    axs[1].plot( pheno.sm_den, pheno.sm_pre, linestyle=lstyle, color=nuda.param.col[kmodel] )
+                else:
+                    model_check.append(model)
+                    axs[1].plot( pheno.sm_den, pheno.sm_pre, linestyle=lstyle, color=nuda.param.col[kmodel], label=model )
+            # end of param
+        # end of model
+    axs[1].errorbar( p_den, p_cen, yerr=p_std, color='k' )
+    axs[1].errorbar( p_den+0.005, p_pheno_cen, yerr=p_pheno_std, color='r' )
+    axs[1].text(0.02,40,'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)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+    #

nucleardatapy/fig/matter_setupFFGNuc_fig.py

@@ -1,7 +1,8 @@
 
 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) ):

nucleardatapy/fig/matter_setupMicroEsym_fig.py

@@ -3,7 +3,7 @@ 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:\
@@ -34,7 +34,7 @@ def matter_setupMicroEsym_fig( pname, models, band ):
     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_xlabel(r'$n_\mathrm{nuc}$ (fm$^{-3}$)')
     axs[1,0].set_xlim([0, 0.3])
     axs[1,0].set_ylim([1, 4])
     #
@@ -42,27 +42,58 @@ def matter_setupMicroEsym_fig( pname, models, band ):
     axs[1,1].set_xlim([0.5, 2.0])
     axs[1,1].set_ylim([1, 4])
     #
-    for ind,model in enumerate(models):
+    mb_check = []
+    #
+    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( )
+                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] )
+
     #
-    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].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 )

nucleardatapy/fig/matter_setupMicro_band_fig.py

@@ -50,6 +50,7 @@ def matter_setupMicro_band_fig( pname, models, den, matter ):
         #
         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':