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

nucleardatapy/fig/matter_ESM_fig.py

@@ -1,140 +0,0 @@
-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.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_\text{nuc})$')
-    axs[0].set_xlim([0, 0.34])
-    axs[0].set_ylim([-22, 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.34])
-    axs[1].set_ylim([-22, 5])
-    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'
-            #
-            if micro.sm_e2a is not None:
-                print('mb:',mb,'model:',model)
-                if mb in mb_check:
-                    if micro.marker:
-                        print('with marker 1:',micro.marker)
-                        if micro.e_err:
-                            print('with error',micro.e_err)
-                            axs[0].errorbar( micro.sm_den, micro.sm_e2a, yerr=micro.sm_e2a_err, marker=micro.marker, markevery=micro.every, linestyle=lstyle, errorevery=micro.every, color=nuda.param.col[kmb] )
-                        else:
-                            print('with no error',micro.e_err)
-                            axs[0].plot( micro.sm_den, micro.sm_e2a, marker=micro.marker, markevery=micro.every, linestyle=lstyle, color=nuda.param.col[kmb] )
-                    else:
-                        print('with no marker',micro.marker)
-                        if micro.e_err:
-                            print('with error',micro.e_err)
-                            axs[0].errorbar( micro.sm_den, micro.sm_e2a, yerr=micro.sm_e2a_err, marker=micro.marker, linestyle=lstyle, errorevery=micro.every, color=nuda.param.col[kmb] )
-                        else:
-                            print('with no error',micro.e_err)
-                            axs[0].plot( micro.sm_den, micro.sm_e2a, marker=micro.marker, linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
-                else:
-                    mb_check.append(mb)
-                    if micro.marker:
-                        print('with marker 2:',micro.marker)
-                        if micro.e_err:
-                            print('with error',micro.e_err)
-                            axs[0].errorbar( micro.sm_den, micro.sm_e2a, yerr=micro.sm_e2a_err, marker=micro.marker, markevery=micro.every, linestyle=lstyle, label=mb, errorevery=micro.every, color=nuda.param.col[kmb] )
-                        else:
-                            print('with no error',micro.e_err)
-                            axs[0].plot( micro.sm_den, micro.sm_e2a, marker=micro.marker, markevery=micro.every, linestyle=lstyle, label=mb, color=nuda.param.col[kmb] )
-                    else:
-                        print('with no marker',micro.marker)
-                        if micro.e_err:
-                            print('with error',micro.e_err)
-                            axs[0].errorbar( micro.sm_den, micro.sm_e2a, yerr=micro.sm_e2a_err, marker=micro.marker, linestyle=lstyle, label=mb, errorevery=micro.every, color=nuda.param.col[kmb] )
-                        else:
-                            print('with no error',micro.e_err)
-                            axs[0].plot( micro.sm_den, micro.sm_e2a, 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.03,2,'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'
-            #
-            if pheno.sm_e2a is not None: 
-                print('model:',model,' param:',param)
-                if model in model_check:
-                    axs[1].plot( pheno.sm_den, pheno.sm_e2a, linestyle=lstyle, color=nuda.param.col[kmodel] )
-                else:
-                    model_check.append(model)
-                    axs[1].plot( pheno.sm_den, pheno.sm_e2a, 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.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)
-    #
-    if pname is not None:
-    	plt.savefig(pname, dpi=200)
-    	plt.close()
-    #

nucleardatapy/fig/matter_Esym_fig.py

@@ -1,134 +0,0 @@
-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}')
-    #
-    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_\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.34])
-    axs[1].set_ylim([0, 60])
-    axs[1].tick_params('y', labelleft=False)
-    #
-    mb_check = []
-    #
-    for kmb,mb in enumerate(micro_mbs):
-        #
-        models, models_lower = nuda.matter.micro_esym_models_mb( mb )
-        #
-        for model in models:
-            #
-            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 micro.esym is not None:
-                print('mb:',mb,'model:',model)
-                if mb in mb_check:
-                    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( micro.den, micro.esym, marker=micro.marker, linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
-                    else:
-                        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( micro.den, micro.esym, marker=micro.marker, linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
-                else:
-                    mb_check.append(mb)
-                    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( micro.den, micro.esym, marker=micro.marker, linestyle=lstyle, label=mb, markevery=micro.every, color=nuda.param.col[kmb] )
-                    else:
-                        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( 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')
-    #
-    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.setupPhenoEsym( 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'
-            #
-            if pheno.esym is not None: 
-                print('model:',model,' param:',param)
-                if model in model_check:
-                    axs[1].plot( pheno.den, pheno.esym, linestyle=lstyle, color=nuda.param.col[kmodel] )
-                else:
-                    model_check.append(model)
-                    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)
-    #
-    if pname is not None:
-    	plt.savefig(pname, dpi=200)
-    	plt.close()
-    #

nucleardatapy/fig/matter_cs2_fig.py

@@ -1,83 +0,0 @@
-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

@@ -1,146 +0,0 @@
-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

@@ -1,144 +0,0 @@
-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()
-    #