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

nucleardatapy/fig/eos_setupCC_fig.py

@@ -0,0 +1,270 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def eos_setupCC_eos_fig( pname, band, crust_model, core_kind, core_model, core_param ):
+    """
+    Plot the EoS with crust and core parts.
+
+    The plot is 1x1.
+
+    [0]: pressure as a function of the energy density.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param band: object instantiated on the reference band.
+    :type band: object.
+    :param crust_model: the name of the crust model.
+    :type crust_model: str.
+    :param core_kind: the kind of the core model ('micro' or 'pheno').
+    :type core_kind: str.
+    :param core_model: the name of the core model.
+    :type core_model: str.
+    :param core_param: the parameters of the core model.
+    :type core_param: str.
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    p_den = 312.0
+    p_cen = 12.5
+    p_std = 11.0
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=0.95, top=0.95, wspace=0.05, hspace=0.05 )
+    #
+    axs.set_xlabel(r'$\epsilon_\text{tot}$ (MeV fm$^{-3}$)',fontsize='14')
+    axs.set_ylabel(r'$p_\text{tot}$ (MeV fm$^{-3}$)',fontsize='14')
+    axs.set_xlim([1.e-2, 2.e3])
+    axs.set_ylim([1.e-5, 2.e3])
+    #axs.set_tick_params('y', right=True)
+    #axs.set_tick_params('x', top=True)
+    axs.set_xscale('log')
+    axs.set_yscale('log')
+    #
+    # check the core EOS against the band in NM
+    #
+    if core_kind == 'micro':
+        eos = nuda.eos.setupAMBeq( model = core_model, kind = core_kind )
+    elif core_kind == 'pheno':
+        eos = nuda.eos.setupAMBeq( model = core_model, param = core_param, kind = core_kind )
+    else:
+        print('eos_setupCC_fig.py: core_kind ',core_kind,' is ill defined.')
+        print('eos_setupCC_fig.py: -- Exit the code --')
+        exit()
+    #        
+    check = nuda.matter.setupCheck( eos = eos, band = band )
+    #
+    if check.isInside:
+        lstyle = 'solid'
+    else:
+        lstyle = 'dashed'
+    #
+    # plot
+    #
+    emp = None
+    for k in range(0,3):
+        if k == 0:
+            connect = 'density'
+            #emp = 'simple'
+            #emp = 'Steiner'
+            boundaries = [ 0.016, 0.16 ] # in units of fm-3
+        elif k == 1:
+            connect = 'epsilon'
+            boundaries = [ 15.0, 150.0 ] # in units of MeV fm-3
+        elif k == 2:
+            connect = 'pressure'
+            boundaries = [ 0.1, 1.0 ] # in units of MeV fm-3
+        #
+        # create the crust+core EOS
+        #
+        eos_cc = nuda.eos.setupCC( crust_model = crust_model, core_kind=core_kind, core_model = core_model, core_param = core_param, connect = connect, boundaries = boundaries, emp = emp )
+        #
+        if eos_cc.pre is not None:
+            axs.plot( eos_cc.eps, eos_cc.pre, marker='.', linestyle=lstyle, markevery=eos_cc.every, label=eos_cc.label )
+    #
+    axs.errorbar( p_den, p_cen, yerr=p_std, color='k', linewidth = 3 )
+    #
+    fig.legend(loc='upper left',bbox_to_anchor=(0.2,0.9),fontsize='10')
+    #fig.legend(loc='upper left',bbox_to_anchor=(0.15,1.0),columnspacing=2,fontsize='8',ncol=3,frameon=False)
+    #
+    if pname is not None: 
+        plt.savefig(pname, dpi=200)
+        plt.close()
+    #
+
+def eos_setupCC_checkpre_fig( pname, band, crust_model, core_kind, core_model, core_param, connect, boundaries, emp ):
+    """
+    Check the EoS with crust and core parts.
+
+    The plot is 1x1.
+
+    [0]: pressure and band to check visualy.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param band: object instantiated on the reference band.
+    :type band: object.
+    :param crust_model: the name of the crust model.
+    :type crust_model: str.
+    :param core_kind: the kind of the core model ('micro' or 'pheno').
+    :type core_kind: str.
+    :param core_model: the name of the core model.
+    :type core_model: str.
+    :param core_param: the parameters of the core model.
+    :type core_param: str.
+    :param connect: The thermodynamic quantity employed to connect the crust and the core. Can be: 'density', 'epsilon' (energy density) or 'pressure'.
+    :type connect: str.
+    :param boundaries: array with lower and upper boundaries to consider. Connected to the variable `connect`. For instance boundaries = [ 0.016, 0.16 ] for 'density' (in units of fm-3), boundaries = [ 15.0, 150.0 ] for 'epsilon' (in units of MeV fm-3), or boundaries = [ 0.1, 1.0 ] for pressure (in units of MeV fm-3).
+    :type boundaries: array of real.
+    :param emp: way to connect the crust and the core. Can be: 'None', 'simple', 'Steiner'.
+    :type emp: str.
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    p_den = 312.0
+    p_cen = 12.5
+    p_std = 11.0
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=0.95, top=0.90, wspace=0.05, hspace=0.05 )
+    #
+    axs.set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs.set_ylabel(r'$p_\text{tot}$ (MeV fm$^{-3}$)')
+    axs.set_xlim([1.e-7, 2.e0])
+    axs.set_ylim([1.e-8, 2.e3])
+    #axs.set_tick_params('y', right=True)
+    #axs.set_tick_params('x', top=True)
+    axs.set_xscale('log')
+    axs.set_yscale('log')
+    #
+    # create the crust+core EOS
+    #
+    eos_cc = nuda.eos.setupCC( crust_model = crust_model, core_kind=core_kind, core_model = core_model, core_param = core_param, connect = connect, boundaries = boundaries, emp = emp )
+    #
+    # check the core EOS against the band in NM
+    #
+    if core_kind == 'micro':
+        eos = nuda.eos.setupAMBeq( model = core_model, kind = core_kind )
+    elif core_kind == 'pheno':
+        eos = nuda.eos.setupAMBeq( model = core_model, param = core_param, kind = core_kind )
+    else:
+        print('eos_setupCC_fig.py: core_kind ',core_kind,' is ill defined.')
+        print('eos_setupCC_fig.py: -- Exit the code --')
+        exit()
+    #        
+    check = nuda.matter.setupCheck( eos = eos, band = band )
+    #
+    if check.isInside:
+        lstyle = 'solid'
+    else:
+        lstyle = 'dashed'
+    #
+    # plot
+    #
+    if eos_cc.pre is not None:
+        axs.plot( eos_cc.den, eos_cc.pre, marker='o', linestyle=lstyle, markevery=eos_cc.every )
+        axs.plot( eos_cc.crust_den, eos_cc.crust_pre, marker='x', linestyle=lstyle, markevery=10 )
+        axs.plot( eos_cc.core_den, eos_cc.core_pre, marker='+', linestyle=lstyle, markevery=eos_cc.every )
+        axs.plot( eos_cc.crux_den, 0.4+eos_cc.crux_pre, marker='s', linestyle=lstyle, markevery=10, color='r' )
+        axs.plot( eos_cc.corx_den, 0.8+eos_cc.corx_pre, marker='s', linestyle=lstyle, markevery=eos_cc.every, color='k' )
+    #
+    #axs.errorbar( p_den, p_cen, yerr=p_std, color='k', linewidth = 3 )
+    #
+    #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()
+    #
+
+def eos_setupCC_checkeos_fig( pname, band, crust_model, core_kind, core_model, core_param, connect, boundaries, emp ):
+    """
+    Check the EoS with crust and core parts.
+
+    The plot is 1x1.
+
+    [0]: pressure and band to check visualy.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param band: object instantiated on the reference band.
+    :type band: object.
+    :param crust_model: the name of the crust model.
+    :type crust_model: str.
+    :param core_kind: the kind of the core model ('micro' or 'pheno').
+    :type core_kind: str.
+    :param core_model: the name of the core model.
+    :type core_model: str.
+    :param core_param: the parameters of the core model.
+    :type core_param: str.
+    :param connect: The thermodynamic quantity employed to connect the crust and the core. Can be: 'density', 'epsilon' (energy density) or 'pressure'.
+    :type connect: str.
+    :param boundaries: array with lower and upper boundaries to consider. Connected to the variable `connect`. For instance boundaries = [ 0.016, 0.16 ] for 'density' (in units of fm-3), boundaries = [ 15.0, 150.0 ] for 'epsilon' (in units of MeV fm-3), or boundaries = [ 0.1, 1.0 ] for pressure (in units of MeV fm-3).
+    :type boundaries: array of real.
+    :param emp: way to connect the crust and the core. Can be: 'None', 'simple', 'Steiner'.
+    :type emp: str.
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    p_den = 312.0
+    p_cen = 12.5
+    p_std = 11.0
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=0.95, top=0.90, wspace=0.05, hspace=0.05 )
+    #
+    axs.set_xlabel(r'$\epsilon_\text{tot}$ (MeV fm$^{-3}$)')
+    axs.set_ylabel(r'$p_\text{tot}$ (MeV fm$^{-3}$)')
+    axs.set_xlim([1.e-4, 3.e3])
+    axs.set_ylim([1.e-8, 2.e3])
+    #axs.set_tick_params('y', right=True)
+    #axs.set_tick_params('x', top=True)
+    axs.set_xscale('log')
+    axs.set_yscale('log')
+    #
+    # create the crust+core EOS
+    #
+    eos_cc = nuda.eos.setupCC( crust_model = crust_model, core_kind=core_kind, core_model = core_model, core_param = core_param, connect = connect, boundaries = boundaries, emp = emp )
+    #
+    # check the core EOS against the band in NM
+    #
+    if core_kind == 'micro':
+        eos = nuda.eos.setupAMBeq( model = core_model, kind = core_kind )
+    elif core_kind == 'pheno':
+        eos = nuda.eos.setupAMBeq( model = core_model, param = core_param, kind = core_kind )
+    else:
+        print('eos_setupCC_fig.py: core_kind ',core_kind,' is ill defined.')
+        print('eos_setupCC_fig.py: -- Exit the code --')
+        exit()
+    #        
+    check = nuda.matter.setupCheck( eos = eos, band = band )
+    #
+    if check.isInside:
+        lstyle = 'solid'
+    else:
+        lstyle = 'dashed'
+    #
+    # plot
+    #
+    if eos_cc.pre is not None:
+        axs.plot( eos_cc.eps, eos_cc.pre, marker='o', linestyle=lstyle, markevery=eos_cc.every )
+        axs.plot( eos_cc.crust_eps, eos_cc.crust_pre, marker='x', linestyle=lstyle, markevery=10 )
+        axs.plot( eos_cc.core_eps, eos_cc.core_pre, marker='+', linestyle=lstyle, markevery=eos_cc.every )
+        axs.plot( eos_cc.crux_eps, 0.4+eos_cc.crux_pre, marker='s', linestyle=lstyle, markevery=10, color='r' )
+        axs.plot( eos_cc.corx_eps, 0.8+eos_cc.corx_pre, marker='s', linestyle=lstyle, markevery=eos_cc.every, color='k' )
+    #
+    #axs.errorbar( p_den, p_cen, yerr=p_std, color='k', linewidth = 3 )
+    #
+    #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/hnuc_setupChart_fig.py

@@ -3,45 +3,48 @@ import matplotlib.pyplot as plt
 
 import nucleardatapy as nuda
 
-def hnuc_setupChart_fig( pname, table1L, table2L, table1Xi ):
+def hnuc_setupChart_fig( pname, tables1L, tables2L, tables1Xi ):
     """
-    Plot hyper-nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
+    Plot hyper-nuclear chart (N versus Z).
+
+    The plot is 1x1 with:
+
     [0]: nuclear chart.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
     :param table1L: table.
-    :type table1L: str.
+    :type table1L: array of str.
     :param table2L: table.
-    :type table2L: str.
+    :type table2L: array of str.
     :param table1Xi: table.
-    :type table1Xi: str.
-
+    :type table1Xi: array of str.
     """
     #
     print(f'Plot name: {pname}')
     #
-    hnuc1L = nuda.hnuc.setupRE1LExp( table = table1L )
-    hnuc2L = nuda.hnuc.setupRE2LExp( table = table2L )
-    hnuc1Xi = nuda.hnuc.setupRE1XiExp( table = table1Xi )
-    #
     # plot
     #
     fig, axs = plt.subplots(1,1)
     fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
     fig.subplots_adjust(left=0.12, bottom=0.15, right=None, top=0.95, wspace=0.3, hspace=0.3)
     #
-    axs.set_ylabel(r'Proton number $Z$')
-    axs.set_xlabel(r'Neutron number $N$')
+    axs.set_ylabel(r'Proton number $Z$ (log scale)',fontsize='14')
+    axs.set_xlabel(r'Neutron number $N$ (log scale)',fontsize='14')
     axs.set_yscale('log')
     axs.set_xscale('log')
     axs.set_xlim([0.8, 140])
     axs.set_ylim([0.8, 90.0])
     #
-    axs.scatter( hnuc1L.N,      hnuc1L.Z,      marker='s', s=12, color=nuda.param.col[0], label=r'1$\Lambda$ from table '+table1L )
-    axs.scatter( hnuc2L.N+0.05, hnuc2L.Z+0.05, marker='s', s=12, color=nuda.param.col[1], label=r'2$\Lambda$ from table '+table2L )
-    axs.scatter( hnuc1Xi.N+0.1, hnuc1Xi.Z+0.1, marker='s', s=12, color=nuda.param.col[2], label=r'1$\Xi^{-}$ from table '+table1Xi )
+    for table1L in tables1L:
+        hnuc1L = nuda.hnuc.setupRE1LExp( table = table1L )
+        axs.scatter( hnuc1L.N,      hnuc1L.Z,      marker='s', s=12, color=nuda.param.col[0], label=r'1$\Lambda$ from table '+table1L )
+    for table2L in tables2L:
+        hnuc2L = nuda.hnuc.setupRE2LExp( table = table2L )
+        axs.scatter( hnuc2L.N+0.05, hnuc2L.Z+0.05, marker='s', s=12, color=nuda.param.col[1], label=r'2$\Lambda$ from table '+table2L )
+    for table1Xi in tables1Xi:
+        hnuc1Xi = nuda.hnuc.setupRE1XiExp( table = table1Xi )
+        axs.scatter( hnuc1Xi.N+0.1, hnuc1Xi.Z+0.1, marker='s', s=12, color=nuda.param.col[2], label=r'1$\Xi^{-}$ from table '+table1Xi )
     #
     #axs.text(0.15,12,r'$K_{sym}$='+str(int(Ksym))+' MeV',fontsize='12')
     axs.legend(loc='lower right',fontsize='11',ncol=1)

nucleardatapy/fig/hnuc_setupRE1LExp_fig.py

@@ -4,15 +4,25 @@ import matplotlib.pyplot as plt
 import nucleardatapy as nuda
 
 def hnuc_setupRE1LExp_fig( pname, tables ):
-    #
-    # plot
+    """
+    Plot the removal energies as a function of A^{-2/3}.
+
+    The plot is 1x1 with:
+
+    [0]: nuclear chart.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param tables: table.
+    :type tables: str.
+    """
     #
     fig, axs = plt.subplots(1,1)
     fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.12, bottom=0.15, right=None, top=0.95, wspace=0.3, hspace=0.3)
+    fig.subplots_adjust(left=0.12, bottom=0.15, right=0.95, top=0.95, wspace=0.3, hspace=0.3)
     #
-    axs.set_ylabel(r'Binding Energy (MeV)')
-    axs.set_xlabel(r'$A^{-2/3}$')
+    axs.set_ylabel(r'Removal energy $B_\Lambda$ (MeV)',fontsize='14')
+    axs.set_xlabel(r'$A^{-2/3}$',fontsize='14')
     axs.set_xlim([0.0, 0.28])
     axs.set_ylim([-5.0, 32.0])
     #