nucleardatapy 0.2.0__py3-none-any.whl

nucleardatapy/fig/__init__.py

@@ -0,0 +1,43 @@
+"""
+This module provides microscopic, phenomenological and experimental data constraints.
+"""
+from nucleardatapy.fig.matter_setupFFGNuc_fig        import *
+from nucleardatapy.fig.matter_setupMicro_fig         import *
+from nucleardatapy.fig.matter_setupMicro_err_NM_fig  import *
+from nucleardatapy.fig.matter_setupMicro_band_fig    import *
+from nucleardatapy.fig.matter_setupMicro_effmass_fig import *
+from nucleardatapy.fig.matter_setupMicro_gap_fig     import *
+from nucleardatapy.fig.matter_setupMicro_LP_fig      import *
+from nucleardatapy.fig.matter_setupMicroEsym_fig     import *
+#
+from nucleardatapy.fig.matter_setupPheno_fig        import *
+from nucleardatapy.fig.matter_setupPhenoEsym_fig    import *
+#
+from nucleardatapy.fig.matter_setupHIC_fig          import *
+#
+from nucleardatapy.fig.matter_setupNEPModelDist_fig import *
+#
+from nucleardatapy.fig.matter_Esym_fig              import *
+from nucleardatapy.fig.matter_ESM_fig               import *
+from nucleardatapy.fig.matter_ENM_fig               import *
+#
+from nucleardatapy.fig.corr_setupKsatQsat_fig       import *
+from nucleardatapy.fig.corr_setupEsymLsym_fig       import *
+from nucleardatapy.fig.corr_setupEsymDen_fig        import *
+#
+from nucleardatapy.fig.nuc_setupBEExp_fig           import *
+#
+from nucleardatapy.fig.hnuc_setupRE1LExp_fig        import *
+from nucleardatapy.fig.hnuc_setupChart_fig          import *
+#
+from nucleardatapy.fig.eos_setupAM_fig              import *
+from nucleardatapy.fig.eos_setupAMBeq_fig           import *
+from nucleardatapy.fig.eos_setupAMLeq_fig           import *
+#
+from nucleardatapy.fig.crust_setupCrust_fig         import *
+#
+from nucleardatapy.fig.astro_setupGW_fig            import *
+from nucleardatapy.fig.astro_setupMasses_fig        import *
+from nucleardatapy.fig.astro_setupMR_fig            import *
+from nucleardatapy.fig.astro_setupMtov_fig          import *
+from nucleardatapy.fig.astro_setupMup_fig           import *

nucleardatapy/fig/astro_setupGW_fig.py

@@ -0,0 +1,67 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def astro_setupGW_fig( pname, sources ):
+    """
+    Plot Tidal deformabilities for each sources.\
+    The plot is 1x1 with:\
+    [0]: Tidal deformabilities versus sources.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param sources: array of sources names.
+    :type sources: array of str.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    # plot Lambda versus sources
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.14, bottom=0.12, right=None, top=0.8, wspace=0.3, hspace=0.3)
+    #
+    axs.set_xlabel(r'sources',fontsize='12')
+    axs.set_ylabel(r'$\tilde{\Lambda}_{90\%}$',fontsize='12')
+    axs.set_xlim([0.8, 2.5])
+    axs.set_ylim([0, 1200])
+    #
+    isource = 1
+    xlabel = []
+    ilabel = []
+    for source in sources:
+        xlabel.append( source )
+        ilabel.append( isource )
+        #
+        # get the mass associated to `source` and `obs`
+        #
+        hyps = nuda.astro.gw_hyps( source = source )
+        print('obss:',hyps)
+        #
+        ihyp = 0
+        for hyp in hyps:
+            gw = nuda.astro.setupGW( source = source, hyp = hyp )
+            if nuda.env.verb_output: gw.print_output( )
+            if nuda.env.verb_latex: gw.print_latex( )
+            axs.errorbar( isource+ihyp/10, gw.lam, yerr=np.array([(gw.lam_sig_do,gw.lam_sig_up)]).T, label=gw.label, color=nuda.param.col[isource], marker=gw.marker, linestyle = 'solid', linewidth = 1 )
+            ihyp += 1
+            #
+        gwav = nuda.astro.setupGWAverage( source = source )
+        if nuda.env.verb_output: gwav.print_output( )
+        if nuda.env.verb_latex: gwav.print_latex( )
+        axs.errorbar( isource+ihyp/10, gwav.lam_cen, yerr=gwav.lam_sig_std, ms=12, label=gwav.label, color=nuda.param.col[isource], marker='^', linestyle = 'solid', linewidth = 3 )
+        #
+        isource += 1
+    #
+    axs.set_xticks( ilabel )
+    axs.set_xticklabels( xlabel )
+    axs.legend(loc='lower center',bbox_to_anchor=(0.5,1.01),columnspacing=2,fontsize='8', ncol=2,frameon=False)
+    #
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+    #

nucleardatapy/fig/astro_setupMR_fig.py

@@ -0,0 +1,88 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def astro_setupMR_fig( pname, sources, sources_av ):
+    """
+    Plot M-R constraints from NICER analyses.\
+    The plot is 1x1 with:\
+    [0]: Masses versus radii.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param sources: array of sources.
+    :type sources: array of str.
+    :param sources_av: .
+    :type sources_av: .
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.85, wspace=0.3, hspace=0.3)
+    #
+    axs.set_xlabel(r'$R$ (km)',fontsize='12')
+    axs.set_ylabel(r'M (M$_\odot$)',fontsize='12')
+    axs.set_xlim([10.5, 16.5])
+    axs.set_ylim([1.15, 2.2])
+    #
+    isource = 1
+    xlabel = []
+    ilabel = []
+    for source in sources:
+        xlabel.append( source )
+        ilabel.append( isource )
+        #
+        # get the mass associated to `source` and `obs`
+        #
+        obss = nuda.astro.mr_obss( source = source )
+        print(f'source: {source}, obss: {obss}')
+        #
+        iobs = 0
+        for obs in obss:
+            mr = nuda.astro.setupMR( source = source, obs = obs )
+            if nuda.env.verb_output: mr.print_output( )
+            if nuda.env.verb_latex: mr.print_latex( )
+            if source.lower() == 'j0030+0451' and obs == 3:
+                axs.errorbar( mr.rad, mr.mass, xerr=np.array([(mr.rad_sig_lo,mr.rad_sig_up)]).T, yerr=np.array([(mr.mass_sig_lo,mr.mass_sig_up)]).T, label=mr.label, color=nuda.param.col[isource], marker=mr.marker, linewidth = 1 )
+            elif source.lower() == 'j0030+0451' and obs == 4:
+                axs.errorbar( mr.rad, mr.mass, xerr=np.array([(mr.rad_sig_lo,mr.rad_sig_up)]).T, yerr=np.array([(mr.mass_sig_lo,mr.mass_sig_up)]).T, label=mr.label, color=nuda.param.col[isource], marker=mr.marker, linewidth = 1 )
+            else:
+                axs.errorbar( mr.rad, mr.mass, xerr=np.array([(mr.rad_sig_lo,mr.rad_sig_up)]).T, yerr=np.array([(mr.mass_sig_lo,mr.mass_sig_up)]).T, label=mr.label, color=nuda.param.col[isource], marker=mr.marker, linewidth = 1 )
+            iobs += 1
+            #
+        isource += 1
+        #
+    isource = 1
+    for source in sources_av:
+        if source.lower() == 'j0030+0451': 
+            obss = [ 1, 2 ]
+        if source.lower() == 'j0740+6620': 
+            obss = [ 1, 2 ]
+        mrav = nuda.astro.setupMRAverage( source = source, obss = obss )
+        if nuda.env.verb_output: mrav.print_output( )
+        if nuda.env.verb_latex: mrav.print_latex( )
+        axs.errorbar( mrav.rad_cen, mrav.mass_cen, xerr=mrav.rad_sig_std, yerr=mrav.mass_sig_std, ms=12, label=mrav.label, color=nuda.param.col[isource], marker='^', linestyle = 'solid', linewidth = 3 )
+        isource += 1
+    #
+    # write source name in plot:
+    #
+    for source in sources:
+        print('source:',source)
+        if source.lower() == 'j0030+0451':
+            axs.text(13.3,1.25,'J0030+0451')
+        elif source.lower() == 'j0740+6620':
+            axs.text(13,1.95,'J0740+6620')
+        elif source.lower() == 'j0437-4715':
+            axs.text(10.7,1.5,'J0437-4715')
+    #
+    #axs.legend(loc='upper left',fontsize='8', ncol=2)
+    axs.legend(loc='lower center',bbox_to_anchor=(0.48,1.01),columnspacing=2,fontsize='8',ncol=3,frameon=False)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+    #

nucleardatapy/fig/astro_setupMasses_fig.py

@@ -0,0 +1,64 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def astro_setupMasses_fig( pname, sources ):
+    """
+    Plot Masses for massives neutron stars as a function of sources.\
+    The plot is 1x1 with:\
+    [0]: masses versus sources.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param sources: array of sources.
+    :type sources: array of str.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.7, wspace=0.3, hspace=0.3)
+    #
+    axs.set_ylabel(r'M (M$_\odot$)',fontsize='12')
+    axs.set_xlabel(r'sources',fontsize='12')
+    axs.set_ylim([1.7, 3.5])
+    axs.set_xlim([0.8, 5.5])
+    #
+    isource = 1
+    xlabel = []
+    ilabel = []
+    for source in sources:
+        xlabel.append( source )
+        ilabel.append( isource )
+        #
+        # get the mass associated to `source` and `obs`
+        #
+        obss = nuda.astro.masses_obss( source = source )
+        print(f'source: {source}, obss: {obss}')
+        #
+        iobs = 0
+        for obs in obss:
+            m = nuda.astro.setupMasses( source = source, obs = obs )
+            if nuda.env.verb_output: m.print_output( )
+            if nuda.env.verb_latex: m.print_latex( )
+            axs.errorbar( isource+iobs/10, m.mass, yerr=np.array([(m.sig_lo,m.sig_up)]).T, label=m.label, color=nuda.param.col[isource], marker='s', linestyle = 'solid', linewidth = 1 )
+            iobs += 1
+            #
+        mav = nuda.astro.setupMassesAverage( source = source )
+        if nuda.env.verb_output: mav.print_output( )
+        if nuda.env.verb_latex: mav.print_latex( )
+        axs.errorbar( isource+iobs/10, mav.mass_cen, yerr=mav.sig_std, label=mav.label, color=nuda.param.col[isource], marker='o', linestyle = 'solid', linewidth = 3 )
+        isource += 1
+    #
+    axs.set_xticks( ilabel )
+    axs.set_xticklabels( xlabel )
+    #axs.legend(loc='upper left',fontsize='8', ncol=2)
+    axs.legend(loc='lower center',bbox_to_anchor=(0.5,1.01),columnspacing=2,fontsize='8', ncol=2,frameon=False)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+    #

nucleardatapy/fig/astro_setupMtov_fig.py

@@ -0,0 +1,84 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def astro_setupMtov_fig( pname, sources_lo_all, sources_lo_dist, sources_up1, sources_up2, sources_up3 ):
+    """
+    Plot the distribution of maximum masses.\
+    The plot is 1x1 with:\
+    [0]: distribution of maximum mass versus maximum mass.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param sources_lo_all: array of all sources (low boundaries).
+    :type sources_lo_all: array of str.
+    :param sources_lo_dist: array of sources considered for the distribution (low boundaries).
+    :type sources_lo_dist: array of str.
+    :param sources_up1: array of sources (up boundaries).
+    :type sources_up1: array of str.
+    :param sources_up2: array of sources (up boundaries).
+    :type sources_up2: array of str.
+    :param sources_up3: array of sources (up boundaries).
+    :type sources_up3: array of str.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.85, wspace=0.3, hspace=0.3)
+    #
+    axs.set_xlabel(r'M$_\mathrm{tov}$ (M$_\odot$)',fontsize='12')
+    axs.set_ylabel(r'probability (non-normalised)',fontsize='12')
+    axs.set_xlim([1.6, 3.4])
+    axs.set_ylim([0.0, 1.03])
+    #
+    sources_lo_nondist = []
+    for source in sources_lo_all:
+        if source not in sources_lo_dist:
+            sources_lo_nondist.append(source)
+    print('sources removed from distribution:',sources_lo_nondist)
+    #
+    prob0 = nuda.astro.setupMtov( sources_lo_nondist, sources_up1 )
+    prob1 = nuda.astro.setupMtov( sources_lo_dist, sources_up1 )
+    prob2 = nuda.astro.setupMtov( sources_lo_dist, sources_up2 )
+    prob3 = nuda.astro.setupMtov( sources_lo_dist, sources_up3 )
+    #
+    for ind,source in enumerate(sources_lo_nondist):
+        #
+        print(f'source: {source}')
+        axs.plot(prob0.mass, prob0.proba_lo[ind], label=prob0.label_lo[ind], color=nuda.param.col[ind], linestyle='dotted', linewidth = 1 )
+        #
+    for ind,source in enumerate(sources_lo_dist):
+        #
+        print(f'source: {source}')
+        axs.plot(prob1.mass, prob1.proba_lo[ind], label=prob1.label_lo[ind], color=nuda.param.col[ind], linestyle='dashed', linewidth = 3 )
+        #
+    #axs.plot(prob1.mass, prob1.proba_lo_tot, label=prob1.label_lo_tot, marker = 'o', markevery=4, color=nuda.param.col[0], linestyle='None' )
+    #
+    for ind,source in enumerate(sources_up1):
+        #
+        print(f'source: {source}')
+        axs.plot(prob1.mass, prob1.proba_up[ind], label=prob1.label_up[ind], color=nuda.param.col[ind], linestyle='dotted', linewidth = 3 )
+        #
+    #axs.plot(prob1.mass, prob1.proba_up_tot, label=prob1.label_up_tot, marker = 's', markevery=4, color=nuda.param.col[4], linestyle='None' )
+    #axs.plot(prob2.mass, prob2.proba_up_tot, label=prob2.label_up_tot, marker = 'v', markevery=7, color=nuda.param.col[5], linestyle='None' )
+    #axs.plot(prob3.mass, prob3.proba_up_tot, label=prob3.label_up_tot, marker = '^', markevery=10, color=nuda.param.col[6], linestyle='None' )
+    #
+    axs.plot(prob1.mass, prob1.proba_tot, label=prob1.label_tot+' case 1', marker = 's', markevery=4, color=nuda.param.col[4], linestyle='None' )
+    axs.plot(prob2.mass, prob2.proba_tot, label=prob2.label_tot+' case 2', marker = 'v', markevery=5, color=nuda.param.col[5], linestyle='None' )
+    axs.plot(prob3.mass, prob3.proba_tot, label=prob3.label_tot+' case 3', marker = '^', markevery=7, color=nuda.param.col[6], linestyle='None' )
+    #
+    #axs.plot(prob1.mass, prob1.proba_tot, label=prob1.label_tot, color=nuda.param.col[4], linestyle=(1, (2, 10)), linewidth = 3 )
+    #axs.plot(prob2.mass, prob2.proba_tot, label=prob2.label_tot, color=nuda.param.col[4], linestyle=(1, (3, 9)), linewidth = 3 )
+    #axs.plot(prob3.mass, prob3.proba_tot, label=prob3.label_tot, color=nuda.param.col[4], linestyle=(1, (4, 8)), linewidth = 3 )
+    #
+    #axs.legend( loc='upper left',fontsize='8', ncol=1 )
+    axs.legend(loc='lower center',bbox_to_anchor=(0.5,1.01),columnspacing=2,fontsize='8',ncol=4,frameon=False)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+    #

nucleardatapy/fig/astro_setupMup_fig.py

@@ -0,0 +1,68 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def astro_setupMup_fig( pname, sources ):
+    """
+    Plot upper boundaries for the tov mass.\
+    The plot is 1x1 with:\
+    [0]: upper boundary for the mass versus sources.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param sources: array of sources.
+    :type sources: array of str.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.8, wspace=0.3, hspace=0.3)
+    #
+    axs.set_xlabel(r'sources',fontsize='12')
+    axs.set_ylabel(r'M (M$_\odot$)',fontsize='12')
+    axs.set_xlim([0.8, 2.5])
+    axs.set_ylim([2.4, 3.4])
+    #
+    isource = 1
+    xlabel = []
+    ilabel = []
+    for source in sources:
+        xlabel.append( source )
+        ilabel.append( isource )
+        #
+        # get mup associated to `source` and `hyp`
+        #
+        hyps = nuda.astro.mup_hyps( source = source )
+        print('source:',source)
+        print('hyps:',hyps)
+        #
+        ihyp = 0
+        for hyp in hyps:
+            mup = nuda.astro.setupMup( source = source, hyp = hyp )
+            if nuda.env.verb_output: mup.print_output( )
+            if nuda.env.verb_latex: mup.print_latex( )
+            axs.errorbar( isource+ihyp/10, mup.mup, yerr=np.array([(mup.sig_do,mup.sig_up)]).T, label=mup.label, color=nuda.param.col[isource], marker=mup.marker, linestyle = 'solid', linewidth = 1 )
+            ihyp += 1
+            #
+        if source=='GW170817': hyps = [ 3, 4 ]
+        mupav = nuda.astro.setupMupAverage( source = source, hyps = hyps )
+        if nuda.env.verb_output: mupav.print_output( )
+        if nuda.env.verb_latex: mupav.print_latex( )
+        axs.errorbar( isource+ihyp/10, mupav.mup_cen, yerr=mupav.sig_std, label=mupav.label, color=nuda.param.col[isource], marker='o', linestyle = 'solid', linewidth = 3 )
+        #
+        isource += 1
+    #
+    axs.set_xticks( ilabel )
+    axs.set_xticklabels( xlabel )
+    #
+    #axs.legend(loc='upper left',fontsize='8', ncol=2)
+    axs.legend(loc='lower center',bbox_to_anchor=(0.5,1.01),columnspacing=2,fontsize='8',ncol=2,frameon=False)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+    #

nucleardatapy/fig/corr_setupEsymDen_fig.py

@@ -0,0 +1,47 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def corr_setupEsymDen_fig( pname, constraints, Ksym ):
+    """
+    Plot upper boundaries for the tov mass.\
+    The plot is 1x1 with:\
+    [0]: upper boundary for the mass versus sources.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param sources: array of sources.
+    :type sources: array of str.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.98, wspace=0.3, hspace=0.3)
+    #
+    axs.set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)',fontsize='12')
+    axs.set_ylabel(r'$E_\text{sym}$ (MeV)',fontsize='12')
+    axs.set_xlim([0.09, 0.27])
+    axs.set_ylim([10, 60])
+    #
+    for constraint in constraints:
+        #
+        esym = nuda.corr.setupEsymDen( constraint = constraint , Ksym=Ksym )
+        #
+        #print('Den:',esym.esym_den)
+        #print('Esym_max:',esym.esym_e2a_max)
+        #print('Esym_min:',esym.esym_e2a_min)
+        #
+        if esym.plot:
+            axs.fill_between( esym.esym_den, y1=esym.esym_e2a_min, y2=esym.esym_e2a_max, label=esym.label, alpha=esym.alpha )
+    #
+    axs.text(0.15,12,r'$K_\text{sym}$='+str(int(Ksym))+' MeV',fontsize='12')
+    axs.legend(loc='lower right',fontsize='9')
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+    #

nucleardatapy/fig/corr_setupEsymLsym_fig.py

@@ -0,0 +1,56 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def corr_setupEsymLsym_fig( pname, constraints ):
+    """
+    Plot the correlation between Esym and Lsym.\
+    The plot is 1x1 with:\
+    [0]: Esym - Lsym correlation plot
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param constraints: list of constraints to run on.
+    :type constraints: array of str.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.98, wspace=0.3, hspace=0.3)
+    #
+    axs.set_xlabel(r'$E_{\mathrm{sym},2}$ (MeV)')
+    axs.set_ylabel(r'$L_{\mathrm{sym},2}$ (MeV)')
+    axs.set_xlim([23, 44])
+    axs.set_ylim([10, 120])
+    #
+    for constraint in constraints:
+        #
+        print('constraint:',constraint)
+        el = nuda.corr.setupEsymLsym( constraint = constraint )
+        if nuda.env.verb: print('Esym:',el.Esym,'+-',el.Esym_err)
+        if nuda.env.verb: print('Lsym:',el.Lsym,'+-',el.Lsym_err)
+        if nuda.env.verb: print('len(Esym):',el.Esym.size)
+        #
+        if el.plot == 'point_err_xy':
+            axs.errorbar( el.Esym, el.Lsym, xerr=el.Esym_err, yerr=el.Lsym_err, linestyle='solid', label=el.label )
+        elif el.plot == 'curve':
+            axs.plot( el.Esym, el.Lsym, linestyle='solid', linewidth=3, label=el.label )
+        elif el.plot == 'contour':
+            axs.plot( el.Esym, el.Lsym, linestyle='solid', label=el.label )
+        elif el.plot == 'band_y':
+            axs.fill_between( el.Esym, y1=el.Lsym-el.Lsym_err, y2=el.Lsym+el.Lsym_err, label=el.label, alpha=el.alpha )
+            #axs.errorbar( el.Esym, el.Lsym, xerr=el.Esym_err, linestyle='solid', label=el.label )
+        elif el.plot == 'band_x':
+            axs.fill_betweenx( el.Lsym, x1=el.Esym-el.Esym_err, x2=el.Esym+el.Esym_err, label=el.label, alpha=el.alpha )
+            #axs.errorbar( el.Esym, el.Lsym, yerr=el.Lsym_err, linestyle='solid', label=el.label )
+        if nuda.env.verb: el.print_outputs( )
+    #
+    axs.legend(loc='lower right',fontsize='9')
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=300)
+    	plt.close()

nucleardatapy/fig/corr_setupKsatQsat_fig.py

@@ -0,0 +1,64 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def corr_setupKsatQsat_fig( pname, constraints ):
+    """
+    Plot the correlation between Ksat and Qsat.\
+    The plot is 1x1 with:\
+    [0]: Ksat - Qsat correlation plot
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param constraints: list of constraints to run on.
+    :type constraints: array of str.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.15, bottom=0.12, right=None, top=0.98, wspace=0.3, hspace=0.3)
+    #
+    axs.set_xlabel(r'$K_\mathrm{sat}$ (MeV)')
+    axs.set_ylabel(r'$Q_\mathrm{sat}$ (MeV)')
+    axs.set_xlim([190, 360])
+    axs.set_ylim([-1000, 1500])
+    #axs.tick_params(labelbottom=True, labeltop=False, labelleft=True, labelright=False,
+    #                 bottom=True, top=True, left=True, right=True)
+#    axs.xaxis.set_major_locator(MultipleLocator(5))
+    #axs.xaxis.set_major_formatter(FormatStrFormatter('%d'))
+#    axs.xaxis.set_minor_locator(MultipleLocator(1))
+#    axs.yaxis.set_major_locator(MultipleLocator(20))
+#    axs.yaxis.set_minor_locator(MultipleLocator(5))
+#    axs.tick_params(axis = 'both', which='major', length=10, width=1, direction='inout', right = True, left = True, bottom = True, top = True)
+#    axs.tick_params(axis = 'both', which='minor', length=5,  width=1, direction='in', right = True, left = True, bottom = True, top = True )
+    #
+    for k,constraint in enumerate(constraints):
+        #
+        print('constraint:',constraint)
+        kq = nuda.corr.setupKsatQsat( constraint = constraint )
+        if nuda.env.verb: print('Ksat:',kq.Ksat)
+        if nuda.env.verb: print('Qsat:',kq.Qsat)
+        if nuda.env.verb: print('len(Ksat):',kq.Ksat.size)
+        #
+        if k == 2:
+            kk = 0
+        else:
+            kk = k
+        axs.scatter( kq.Ksat, kq.Qsat, label=kq.label, color=nuda.param.col[kk], marker=kq.marker )
+        x = np.linspace(min(kq.Ksat),max(kq.Ksat),10)
+        if k == 3 or k == 4 or k == 5:
+            axs.plot( x, nuda.corr.flinear(x,kq.m,kq.c), color=nuda.param.col[kk], linestyle='dashed' )
+        else:
+            axs.plot( x, nuda.corr.flinear(x,kq.m,kq.c), color=nuda.param.col[kk], linestyle='solid' )
+        #
+        if nuda.env.verb: kq.print_outputs( )
+    #
+    axs.legend(loc='upper left',ncol=3, fontsize='9')
+    #
+    if pname is not None: 
+        plt.savefig(pname, dpi=200)
+        plt.close()

nucleardatapy/fig/crust_setupCrust_fig.py

@@ -0,0 +1,52 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def crust_setupCrust_fig( pname, models ):
+    """
+    Plot hyper-nuclear chart (N versus Z).\
+    The plot is 1x2 with:\
+    [0]: nuclear chart.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param models: table.
+    :type models: 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.8, wspace=0.3, hspace=0.3 )
+    #
+    axs[0].set_xlabel(r'n (fm$^{-3}$)')
+    axs[0].set_ylabel(r'$e_{int}(n)$')
+    axs[0].set_xlim([1e-4, 1e-1])
+    axs[0].set_ylim([-2, 10])
+    axs[0].set_xscale('log')
+    #
+    axs[1].set_xlabel(r'n (fm$^{-3}$)')
+    axs[1].set_ylabel(r'$Z$')
+    axs[1].set_xlim([1e-4, 1e-1])
+    axs[1].set_ylim([10, 100])
+    axs[1].set_xscale('log')
+    #
+    for model in models:
+        #
+        print('model:',model)
+        crust = nuda.crust.setupCrust( model = model )
+        if crust.e2a_int is not None: 
+            axs[0].plot( crust.den, crust.e2a_int, label=crust.label, linestyle=crust.linestyle )
+        if crust.Z is not None: 
+            axs[1].plot( crust.den, crust.Z, linestyle=crust.linestyle )
+    #axs[0].legend(loc='upper left',fontsize='8', ncol=1)
+    #axs[1].legend(loc='upper left',fontsize='8', ncol=1)
+    fig.legend(loc='upper left',bbox_to_anchor=(0.01,1.01),columnspacing=2,fontsize='8',ncol=4,frameon=False)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+    #