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

nucleardatapy/fig/eos_setupCC_fig.py

@@ -0,0 +1,240 @@
+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.
+
+    :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}')
+    #
+    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 ):
+    """
+    Plot the EoS with crust and core parts.\
+    The plot is 1x1.
+
+    :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}')
+    #
+    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 ):
+    """
+    Plot the EoS with crust and core parts.\
+    The plot is 1x1.
+
+    :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}')
+    #
+    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

@@ -32,8 +32,8 @@ def hnuc_setupChart_fig( pname, table1L, table2L, table1Xi ):
     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])

nucleardatapy/fig/hnuc_setupRE1LExp_fig.py

@@ -3,16 +3,16 @@ import matplotlib.pyplot as plt
 
 import nucleardatapy as nuda
 
-def hnuc_setupRE1LExp_plot( pname, tables ):
+def hnuc_setupRE1LExp_fig( pname, tables ):
     #
     # 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)
+    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])
     #