nucleardatapy 0.2.0__py3-none-any.whl

nucleardatapy/fig/matter_setupMicro_gap_fig.py

@@ -0,0 +1,161 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_setupMicro_gap_1s0_fig( pname, models, matter = 'NM' ):
+    """
+    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}')
+    #
+    print('models:',models)
+    #
+    fig, axs = plt.subplots(2,2)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.8, wspace=0.05, hspace=0.05 )
+    #
+    axs[0,0].set_ylabel(r'$\Delta_{1S0}$ (MeV)')
+    axs[0,0].set_xlim([0, 0.048])
+    axs[0,0].set_ylim([0, 3.0])
+    axs[0,0].tick_params('x', labelbottom=False)
+    #
+    axs[0,1].set_xlim([0, 1.6])
+    axs[0,1].set_ylim([0, 3.0])
+    axs[0,1].tick_params('x', labelbottom=False)
+    axs[0,1].tick_params('y', labelleft=False)
+    #
+    axs[1,0].set_ylabel(r'$\Delta_{1S0}/E_F$')
+    axs[1,0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs[1,0].set_xlim([0, 0.048])
+    axs[1,0].set_ylim([0, 0.78])
+    #
+    axs[1,1].set_xlabel(r'$k_{F_n}$ (fm$^{-1}$)')
+    axs[1,1].set_xlim([0, 1.6])
+    axs[1,1].set_ylim([0, 0.78])
+    axs[1,1].tick_params('y', labelleft=False)
+    #
+    for model in models:
+        #
+        gap = nuda.matter.setupMicroGap( model = model, matter = matter )
+        #
+        if matter.lower() == 'nm':
+            if gap.nm_gap_1s0 is not None:
+                if gap.nm_gap_1s0_err is not None:
+                    axs[0,0].errorbar( gap.nm_den_1s0, gap.nm_gap_1s0, yerr=gap.nm_gap_1s0_err, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
+                    axs[0,1].errorbar( gap.nm_kfn_1s0, gap.nm_gap_1s0, yerr=gap.nm_gap_1s0_err, marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,0].errorbar( gap.nm_den_1s0, gap.nm_gap_1s0/nuda.eF_n(gap.nm_kfn_1s0), yerr=gap.nm_gap_1s0_err/nuda.eF_n(gap.nm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,1].errorbar( gap.nm_kfn_1s0, gap.nm_gap_1s0/nuda.eF_n(gap.nm_kfn_1s0), yerr=gap.nm_gap_1s0_err/nuda.eF_n(gap.nm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                else:
+                    axs[0,0].plot( gap.nm_den_1s0, gap.nm_gap_1s0, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
+                    axs[0,1].plot( gap.nm_kfn_1s0, gap.nm_gap_1s0, marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,0].plot( gap.nm_den_1s0, gap.nm_gap_1s0/nuda.eF_n(gap.nm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,1].plot( gap.nm_kfn_1s0, gap.nm_gap_1s0/nuda.eF_n(gap.nm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
+        elif matter.lower() == 'sm':
+            if gap.sm_gap_1s0 is not None:
+                if gap.sm_gap_1s0_err is not None:
+                    axs[0,0].errorbar( gap.sm_den_1s0, gap.sm_gap_1s0, yerr=gap.sm_gap_1s0_err, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
+                    axs[0,1].errorbar( gap.sm_kfn_1s0, gap.sm_gap_1s0, yerr=gap.sm_gap_1s0_err, marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,0].errorbar( gap.sm_den_1s0, gap.sm_gap_1s0/nuda.eF_n(gap.sm_kfn_1s0), yerr=gap.sm_gap_1s0_err/nuda.eF_n(gap.sm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,1].errorbar( gap.sm_kfn_1s0, gap.sm_gap_1s0/nuda.eF_n(gap.sm_kfn_1s0), yerr=gap.sm_gap_1s0_err/nuda.eF_n(gap.sm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                else:
+                    axs[0,0].plot( gap.sm_den_1s0, gap.sm_gap_1s0, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
+                    axs[0,1].plot( gap.sm_kfn_1s0, gap.sm_gap_1s0, marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,0].plot( gap.sm_den_1s0, gap.sm_gap_1s0/nuda.eF_n(gap.sm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,1].plot( gap.sm_kfn_1s0, gap.sm_gap_1s0/nuda.eF_n(gap.sm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
+        if nuda.env.verb_output: gap.print_outputs( )
+    #
+    #axs[1,0].legend(loc='upper right',fontsize='8')
+    fig.legend(loc='upper left',bbox_to_anchor=(0.08,1.0),columnspacing=2,fontsize='8',ncol=3,frameon=False)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=300)
+    	plt.close()
+
+def matter_setupMicro_gap_3pf2_fig( pname, models, matter = 'NM' ):
+    """
+    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}')
+    #
+    # plot 3PF2 pairing gap in NM
+    print('models:',models)
+    #
+    fig, axs = plt.subplots(2,2)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.8, wspace=0.05, hspace=0.05 )
+    #
+    axs[0,0].set_ylabel(r'$\Delta_{3PF2}$ (MeV)')
+    axs[0,0].set_xlim([0, 0.38])
+    axs[0,0].set_ylim([0, 0.6])
+    axs[0,0].tick_params('x', labelbottom=False)
+    #
+    axs[0,1].set_xlim([0.6, 2.3])
+    axs[0,1].set_ylim([0, 0.6])
+    axs[0,1].tick_params('x', labelbottom=False)
+    axs[0,1].tick_params('y', labelleft=False)
+    #
+    axs[1,0].set_ylabel(r'$100\times \Delta_{3PF2}/E_F$')
+    axs[1,0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs[1,0].set_xlim([0, 0.38])
+    axs[1,0].set_ylim([0, 0.78])
+    #
+    axs[1,1].set_xlabel(r'$k_{F_n}$ (fm$^{-1}$)')
+    axs[1,1].set_xlim([0.6, 2.3])
+    axs[1,1].set_ylim([0, 0.78])
+    axs[1,1].tick_params('y', labelleft=False)
+    #
+    for model in models:
+        #
+        gap = nuda.matter.setupMicroGap( model = model, matter = matter )
+        #
+        if matter.lower() == 'nm':
+            if gap.nm_gap_3pf2 is not None:
+                if gap.nm_gap_3pf2_err is not None:
+                    axs[0,0].errorbar( gap.nm_den_3pf2, gap.nm_gap_3pf2, yerr=gap.nm_gap_3pf2_err, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
+                    axs[0,1].errorbar( gap.nm_kfn_3pf2, gap.nm_gap_3pf2, yerr=gap.nm_gap_3pf2_err, marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,0].errorbar( gap.nm_den_3pf2, 100*gap.nm_gap_3pf2/nuda.eF_n(gap.nm_kfn_3pf2), yerr=gap.nm_gap_3pf2_err/nuda.eF_n(gap.nm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,1].errorbar( gap.nm_kfn_3pf2, 100*gap.nm_gap_3pf2/nuda.eF_n(gap.nm_kfn_3pf2), yerr=gap.nm_gap_3pf2_err/nuda.eF_n(gap.nm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                else:
+                    axs[0,0].plot( gap.nm_den_3pf2, gap.nm_gap_3pf2, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
+                    axs[0,1].plot( gap.nm_kfn_3pf2, gap.nm_gap_3pf2, marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,0].plot( gap.nm_den_3pf2, 100*gap.nm_gap_3pf2/nuda.eF_n(gap.nm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,1].plot( gap.nm_kfn_3pf2, 100*gap.nm_gap_3pf2/nuda.eF_n(gap.nm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
+        elif matter.lower() == 'sm':
+            if gap.sm_gap_3pf2 is not None:
+                if gap.sm_gap_3pf2_err is not None:
+                    axs[0,0].errorbar( gap.sm_den_3pf2, gap.sm_gap_3pf2, yerr=gap.sm_gap_3pf2_err, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
+                    axs[0,1].errorbar( gap.sm_kfn_3pf2, gap.sm_gap_3pf2, yerr=gap.sm_gap_3pf2_err, marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,0].errorbar( gap.sm_den_3pf2, 100*gap.sm_gap_3pf2/nuda.eF_n(gap.sm_kfn_3pf2), yerr=gap.sm_gap_3pf2_err/nuda.eF_n(gap.sm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,1].errorbar( gap.sm_kfn_3pf2, 100*gap.sm_gap_3pf2/nuda.eF_n(gap.sm_kfn_3pf2), yerr=gap.sm_gap_3pf2_err/nuda.eF_n(gap.sm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                else:
+                    axs[0,0].plot( gap.sm_den_3pf2, gap.sm_gap_3pf2, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
+                    axs[0,1].plot( gap.sm_kfn_3pf2, gap.sm_gap_3pf2, marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,0].plot( gap.sm_den_3pf2, 100*gap.sm_gap_3pf2/nuda.eF_n(gap.sm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[1,1].plot( gap.sm_kfn_3pf2, 100*gap.sm_gap_3pf2/nuda.eF_n(gap.sm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
+        if nuda.env.verb_output: gap.print_outputs( )
+    #
+    #axs[1,0].legend(loc='upper right',fontsize='8')
+    fig.legend(loc='upper left',bbox_to_anchor=(0.04,1.0),columnspacing=2,fontsize='8',ncol=3,frameon=False)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=300)
+    	plt.close()

nucleardatapy/fig/matter_setupNEPModelDist_fig.py

@@ -0,0 +1,68 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_setupNEPModelDist_fig( pname, models ):
+    """
+    Plot nucleonic energy per particle E/A in matter.\
+    The plot is 5x2 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 models: models to run on.
+    :type models: array of str.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    fig, axs = plt.subplots(5,2)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.06, right=None, top=0.9, wspace=0.3, hspace=0.4 )
+    #
+    axs[0,0].set_ylabel(r'$E_\text{sat}$')
+    axs[0,0].set_xlim([-16.5, -15])
+    axs[1,0].set_ylabel(r'$n_\text{sat}$')
+    axs[1,0].set_xlim([0.14, 0.18])
+    axs[2,0].set_ylabel(r'$K_\text{sat}$')
+    axs[2,0].set_xlim([180, 360])
+    axs[3,0].set_ylabel(r'$Q_\text{sat}$')
+    axs[3,0].set_xlim([-1000, 1000])
+    axs[4,0].set_ylabel(r'$m_\text{sat}^{*}/m$')
+    axs[4,0].set_xlim([0.35, 1.2])
+    axs[0,1].set_ylabel(r'$E_\text{sym}$')
+    axs[0,1].set_xlim([26, 40])
+    axs[1,1].set_ylabel(r'$L_\text{sym}$')
+    axs[1,1].set_xlim([0, 120])
+    axs[2,1].set_ylabel(r'$K_\text{sym}$')
+    axs[2,1].set_xlim([-400, 220])
+    axs[3,1].set_ylabel(r'$Q_\text{sym}$')
+    axs[3,1].set_xlim([-50, 900])
+    axs[4,1].set_ylabel(r'$\Delta m_\text{sat}^{*}/m$')
+    axs[4,1].set_xlim([-0.5, 1.1])
+    #
+    # Built distribution of NEP
+    #
+    for model in models:
+        #
+        dist = nuda.matter.setupNEPModelDist( model )
+        #
+        axs[0,0].hist( dist.Esat, bins=10, label=model )
+        axs[1,0].hist( dist.nsat, bins=10 )
+        axs[2,0].hist( dist.Ksat, bins=10 )
+        axs[3,0].hist( dist.Qsat, bins=10 )
+        axs[4,0].hist( dist.msat, bins=10 )
+        axs[0,1].hist( dist.Esym, bins=10 )
+        axs[1,1].hist( dist.Lsym, bins=10 )
+        axs[2,1].hist( dist.Ksym, bins=10 )
+        axs[3,1].hist( dist.Qsym, bins=10 )
+        axs[4,1].hist( dist.Dmsat, bins=10 )
+    #
+    #axs[0,0].legend(loc='lower right',fontsize='10',ncol=2)
+    fig.legend(loc='upper left',bbox_to_anchor=(0.2,0.99),columnspacing=2,fontsize='8',ncol=4,frameon=False)
+    #
+    if pname is not None:
+        plt.savefig(pname, dpi=300)
+        plt.close()

nucleardatapy/fig/matter_setupPhenoEsym_fig.py

@@ -0,0 +1,100 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_setupPhenoEsym_fig( pname, models, band ):
+    """
+    Plot nuclear chart (N versus Z).\
+    The plot is 1x2 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(2,2)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.8, wspace=0.2, hspace=0.2 )
+    #
+    axs[0,0].set_ylabel(r'$E_\mathrm{sym}$ (MeV)')
+    axs[0,0].set_xlim([0, 0.3])
+    axs[0,0].set_ylim([0, 50])
+    #
+    axs[0,1].set_xlim([0.5, 2.0])
+    axs[0,1].set_ylim([0, 50])
+    #
+    axs[1,0].set_ylabel(r'$E_\mathrm{sym}/E_\mathrm{sym, FFG}$')
+    axs[1,0].set_xlabel(r'n (fm$^{-3}$)')
+    axs[1,0].set_xlim([0, 0.3])
+    axs[1,0].set_ylim([1, 4])
+    #
+    axs[1,1].set_xlabel(r'$k_F$ (fm$^{-1}$)')
+    axs[1,1].set_xlim([0.5, 2.0])
+    axs[1,1].set_ylim([1, 4])
+    #
+    for model in models:
+        #
+        params, params_lower = nuda.matter.pheno_params( model = model )
+        #
+        for param in params:
+            #
+            print('in Sample: model, param',model,param)
+            esym = nuda.matter.setupPhenoEsym( model = model, param = param )
+            #
+            #print("esym:",esym.esym)
+            #print("den:",esym.den)
+            if esym.esym is not None:
+                print("esym_err:",esym.esym_err)
+                if esym.esym_err is not None:
+                    if not nuda.matter.checkPheno(esym,band,'Esym'):
+                        axs[0,0].errorbar( esym.den, esym.esym, yerr=esym.esym_err, marker=esym.marker, linestyle='none', errorevery=esym.every )
+                    else:
+                        axs[0,0].errorbar( esym.den, esym.esym, yerr=esym.esym_err, marker=esym.marker, linestyle='none', errorevery=esym.every, label=esym.label )
+                    axs[0,1].errorbar( esym.kf,  esym.esym, yerr=esym.esym_err, marker=esym.marker, 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, 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, linestyle='none', errorevery=esym.every )
+                else:
+                    if not nuda.matter.checkPheno(esym,band,'Esym'):
+                        axs[0,0].plot( esym.den, esym.esym, marker=esym.marker, linestyle='dotted' )
+                        axs[0,1].plot( esym.kf,  esym.esym, marker=esym.marker, linestyle='dotted' )
+                        axs[1,0].plot( esym.den, esym.esym/nuda.esymffg_nr(esym.kf), marker=esym.marker, linestyle='dotted' )
+                        axs[1,1].plot( esym.kf,  esym.esym/nuda.esymffg_nr(esym.kf), marker=esym.marker, linestyle='dotted' )
+                    else:
+                        axs[0,0].plot( esym.den, esym.esym, marker=esym.marker, linestyle=esym.linestyle, label=esym.label )
+                        axs[0,1].plot( esym.kf,  esym.esym, marker=esym.marker, linestyle=esym.linestyle )
+                        axs[1,0].plot( esym.den, esym.esym/nuda.esymffg_nr(esym.kf), marker=esym.marker, linestyle=esym.linestyle )
+                        axs[1,1].plot( esym.kf,  esym.esym/nuda.esymffg_nr(esym.kf), marker=esym.marker, linestyle=esym.linestyle )
+            if nuda.env.verb_output: esym.print_outputs( )
+            #
+        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  )
+        axs[0,0].plot( band.den, (band.e2a+band.e2a_std), color='k', linestyle='dashed', visible=True  )
+        axs[0,1].fill_between( band.kfn, y1=(band.e2a-band.e2a_std), y2=(band.e2a+band.e2a_std), color=band.color, alpha=band.alpha, visible=True  )
+        axs[0,1].plot( band.kfn, (band.e2a-band.e2a_std), color='k', linestyle='dashed', visible=True  )
+        axs[0,1].plot( band.kfn, (band.e2a+band.e2a_std), color='k', linestyle='dashed', visible=True  )
+        axs[1,0].fill_between( band.den, y1=(band.e2a-band.e2a_std)/nuda.esymffg_nr(band.kfn), y2=(band.e2a+band.e2a_std)/nuda.esymffg_nr(band.kfn), color=band.color, alpha=band.alpha, visible=True  )
+        axs[1,0].plot( band.den, (band.e2a-band.e2a_std)/nuda.esymffg_nr(band.kfn), color='k', linestyle='dashed', visible=True  )
+        axs[1,0].plot( band.den, (band.e2a+band.e2a_std)/nuda.esymffg_nr(band.kfn), color='k', linestyle='dashed', visible=True  )
+        axs[1,1].fill_between( band.kf, y1=(band.e2a-band.e2a_std)/nuda.esymffg_nr(band.kfn), y2=(band.e2a+band.e2a_std)/nuda.esymffg_nr(band.kfn), color=band.color, alpha=band.alpha, visible=True  )
+        axs[1,1].plot( band.kfn, (band.e2a-band.e2a_std)/nuda.esymffg_nr(band.kfn), color='k', linestyle='dashed', visible=True  )
+        axs[1,1].plot( band.kfn, (band.e2a+band.e2a_std)/nuda.esymffg_nr(band.kfn), color='k', linestyle='dashed', visible=True  )
+    #
+    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[1,0].legend(loc='upper right',fontsize='8')
+    fig.legend(loc='upper left',bbox_to_anchor=(0.1,1.0),columnspacing=2,fontsize='8',ncol=4,frameon=False)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=300)
+    	plt.close()

nucleardatapy/fig/matter_setupPheno_fig.py

@@ -0,0 +1,121 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def matter_setupPheno_E_fig( pname, model, band, matter ):
+    """
+    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 model: class of model considered.
+    :type model: 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}')
+    #
+    # plot
+    #
+    fig, axs = plt.subplots(2,2)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.15, bottom=0.12, right=None, top=0.88, wspace=0.05, hspace=0.05)
+    #
+    axs[0,0].set_xlim([0, 0.33])
+    axs[1,0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)',fontsize='12')
+    axs[1,0].set_xlim([0, 0.33])
+    axs[0,1].set_xlim([0.5, 1.5])
+    axs[1,1].set_xlabel(r'$k_{F_n}$ (fm$^{-1}$)',fontsize='12')
+    axs[1,1].set_xlim([0.5, 1.5])
+    #
+    axs[0,0].tick_params('x', labelbottom=False)
+    axs[0,1].tick_params('x', labelbottom=False)
+    axs[0,1].tick_params('y', labelleft=False)
+    axs[1,1].tick_params('y', labelleft=False)
+    #
+    if matter.lower() == 'nm':
+        axs[0,0].set_ylabel(r'$E_\text{NM}/A$ (MeV)',fontsize='12')
+        axs[0,0].set_ylim([0, 35])
+        axs[0,1].set_ylim([0, 15])
+        axs[1,0].set_ylabel(r'$E_\text{NM}/E_\text{NRFFG,NM}$',fontsize='12')
+        axs[1,0].set_ylim([0.2, 0.84])
+        axs[1,1].set_ylim([0.2, 0.84])
+    elif matter.lower() == 'sm':
+        axs[0,0].set_ylabel(r'$E_\text{SM}/A$ (MeV)',fontsize='12')
+        axs[0,0].set_ylim([-20, 10])
+        axs[0,1].set_ylim([-20, 10])
+        axs[1,0].set_ylabel(r'$E_\text{SM}/E_\text{NRFFG,SM}$',fontsize='12')
+        axs[1,0].set_ylim([-2.0, 0.1])
+        axs[1,1].set_ylim([-2.0, 0.1])
+    #
+    params, params_lower = nuda.matter.pheno_params( model = model )
+    #
+    for param in params:
+        #
+        pheno = nuda.matter.setupPheno( model = model, param = param )
+        if matter.lower() == 'nm':
+            if any(pheno.nm_e2a): 
+                if model == 'Skyrme' and not nuda.matter.checkPheno(pheno,band,matter):
+                    axs[0,0].plot( pheno.nm_den, pheno.nm_e2a, linestyle='solid' )
+                else:
+                    axs[0,0].plot( pheno.nm_den, pheno.nm_e2a, linestyle='solid', label=pheno.label )
+                    #print('model:',model,' param:',param,' inside band')
+                axs[0,1].plot( pheno.nm_kfn, pheno.nm_e2a, linestyle='solid' )
+                axs[1,0].plot( pheno.nm_den, pheno.nm_e2a/nuda.effg_nr(pheno.nm_kfn), linestyle='solid' )
+                axs[1,1].plot( pheno.nm_kfn, pheno.nm_e2a/nuda.effg_nr(pheno.nm_kfn), linestyle='solid' )
+        elif matter.lower() == 'sm':
+            if any(pheno.sm_e2a): 
+                if model == 'Skyrme' and not nuda.matter.checkPheno(pheno,band,matter):
+                    axs[0,0].plot( pheno.sm_den, pheno.sm_e2a, linestyle='solid' )
+                else:
+                    axs[0,0].plot( pheno.sm_den, pheno.sm_e2a, linestyle='solid', label=pheno.label )
+                    #print('model:',model,' param:',param,' inside band')
+                #axs[0,0].plot( pheno.sm_den, pheno.sm_e2a, linestyle='solid', label=pheno.label )
+                axs[0,1].plot( pheno.sm_kf, pheno.sm_e2a, linestyle='solid' )
+                axs[1,0].plot( pheno.sm_den, pheno.sm_e2a/nuda.effg_nr(pheno.sm_kf), linestyle='solid' )
+                axs[1,1].plot( pheno.sm_kf, pheno.sm_e2a/nuda.effg_nr(pheno.sm_kf), linestyle='solid' )
+        if nuda.env.verb_output: pheno.print_outputs( )
+    if matter.lower() == 'nm':
+        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 )
+        axs[0,0].plot( band.den, (band.e2a-band.e2a_std), color='k', linestyle='dashed' )
+        axs[0,0].plot( band.den, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
+        axs[0,1].fill_between( band.kfn, y1=(band.e2a-band.e2a_std), y2=(band.e2a+band.e2a_std), color=band.color, alpha=band.alpha )
+        axs[0,1].plot( band.kfn, (band.e2a-band.e2a_std), color='k', linestyle='dashed' )
+        axs[0,1].plot( band.kfn, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
+        axs[1,0].fill_between( band.den, y1=(band.e2a-band.e2a_std)/nuda.effg_nr(band.kfn), y2=(band.e2a+band.e2a_std)/nuda.effg_nr(band.kfn), color=band.color, alpha=band.alpha )
+        axs[1,0].plot( band.den, (band.e2a-band.e2a_std)/nuda.effg_nr(band.kfn), color='k', linestyle='dashed' )
+        axs[1,0].plot( band.den, (band.e2a+band.e2a_std)/nuda.effg_nr(band.kfn), color='k', linestyle='dashed' )
+        axs[1,1].fill_between( band.kfn, y1=(band.e2a-band.e2a_std)/nuda.effg_nr(band.kfn), y2=(band.e2a+band.e2a_std)/nuda.effg_nr(band.kfn), color=band.color, alpha=band.alpha )
+        axs[1,1].plot( band.kfn, (band.e2a-band.e2a_std)/nuda.effg_nr(band.kfn), color='k', linestyle='dashed' )
+        axs[1,1].plot( band.kfn, (band.e2a+band.e2a_std)/nuda.effg_nr(band.kfn), color='k', linestyle='dashed' )
+    elif matter.lower() == 'sm':
+        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' )
+        axs[0,0].plot( band.den, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
+        axs[0,1].fill_between( band.kfn, y1=(band.e2a-band.e2a_std), y2=(band.e2a+band.e2a_std), color=band.color, alpha=band.alpha, visible=True )
+        axs[0,1].plot( band.kfn, (band.e2a-band.e2a_std), color='k', linestyle='dashed' )
+        axs[0,1].plot( band.kfn, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
+        axs[1,0].fill_between( band.den, y1=(band.e2a-band.e2a_std)/nuda.effg_nr(band.kfn), y2=(band.e2a+band.e2a_std)/nuda.effg_nr(band.kfn), color=band.color, alpha=band.alpha, visible=True )
+        axs[1,0].plot( band.den, (band.e2a-band.e2a_std)/nuda.effg_nr(band.kfn), color='k', linestyle='dashed' )
+        axs[1,0].plot( band.den, (band.e2a+band.e2a_std)/nuda.effg_nr(band.kfn), color='k', linestyle='dashed' )
+        axs[1,1].fill_between( band.kfn, y1=(band.e2a-band.e2a_std)/nuda.effg_nr(band.kfn), y2=(band.e2a+band.e2a_std)/nuda.effg_nr(band.kfn), color=band.color, alpha=band.alpha, visible=True )
+        axs[1,1].plot( band.kfn, (band.e2a-band.e2a_std)/nuda.effg_nr(band.kfn), color='k', linestyle='dashed' )
+        axs[1,1].plot( band.kfn, (band.e2a+band.e2a_std)/nuda.effg_nr(band.kfn), color='k', linestyle='dashed' )
+    #
+    #if model != 'Skyrme':
+    #    axs[0,0].legend(loc='upper right',fontsize='8', ncol=2)
+    fig.legend(loc='upper left',bbox_to_anchor=(0.1,1.0),columnspacing=2,fontsize='8',ncol=4,frameon=False)
+    #
+    #
+    if pname is not None:
+        plt.savefig(pname, dpi=300)
+        plt.close()
+    #

nucleardatapy/fig/nuc_setupBEExp_fig.py

@@ -0,0 +1,122 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def nuc_setupBEExp_chart_lt_fig( pname, table, version, theo_tables ):
+    """
+    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}')
+    #
+    print('Table:',table)
+    #
+    # 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.85, wspace=0.3, hspace=0.3)
+    #
+    axs.set_title(r''+table+' mass table version '+version)
+    axs.set_ylabel(r'Z')
+    axs.set_xlabel(r'N')
+    axs.set_xlim([0, 200])
+    axs.set_ylim([0, 132])
+    axs.text(10,120,'Number of nuclei:')
+    #
+    # longlive nuclei
+    #
+    mas = nuda.nuc.setupBEExp( table = table, version = version )
+    ustbl = mas.select( state= 'gs', interp = 'n', nucleus = 'longlive' )
+    axs.scatter( ustbl.sel_nucN, ustbl.sel_nucZ, marker='s', s = 3, linewidth=0, color = 'grey', label='longlive ('+str(ustbl.sel_nbNucSel)+')' )
+    #axs.text(10,96,'long live: '+str(ustbl.sel_nbNucSel))
+    #
+    # shortlive nuclei
+    #
+    mas = nuda.nuc.setupBEExp( table = table, version = version )
+    ustbl = mas.select( state= 'gs', interp = 'n', nucleus = 'shortlive' )
+    axs.scatter( ustbl.sel_nucN, ustbl.sel_nucZ, marker='s', s = 3, linewidth=0, color = 'r', label='shortlive ('+str(ustbl.sel_nbNucSel)+')' )
+    #axs.text(10,88,'short live: '+str(ustbl.sel_nbNucSel))
+    #
+    # veryshortlive nuclei
+    #
+    mas = nuda.nuc.setupBEExp( table = table, version = version )
+    ustbl = mas.select( state= 'gs', interp = 'n', nucleus = 'veryshortlive' )
+    axs.scatter( ustbl.sel_nucN, ustbl.sel_nucZ, marker='s', s = 3, linewidth=0, color = 'b', label='veryshortlive ('+str(ustbl.sel_nbNucSel)+')' )
+    #axs.text(10,80,'very short live: '+str(ustbl.sel_nbNucSel))
+    #
+    # unstable nuclei:
+    #
+    mas = nuda.nuc.setupBEExp( table = table, version = version )
+    ustbl = mas.select( state= 'gs', interp = 'n', nucleus = 'unstable' )
+    #axs.scatter( ustbl.sel_nucN, ustbl.sel_Z, marker='.', s = 1, linewidth=0, color = 'b' )
+    axs.text(10,104,'unstable: '+str(ustbl.sel_nbNucSel))
+    #
+    # drip line nuclei
+    #
+    legend = 0
+    for i,theo_table in enumerate( theo_tables ):
+        theo = nuda.nuc.setupBETheo( table = theo_table )
+        s2n = theo.S2n( Zmin=1, Zmax = 95 )
+        drip_S2n = s2n.drip_S2n( Zmin = 1, Zmax = 95 )
+        if legend == 0:
+            axs.scatter( drip_S2n.drip_S2n_N, drip_S2n.drip_S2n_Z, marker='o', s = 3, linewidth=0, color = 'green', label='Drip Lines' )
+            legend = 1
+        else:
+            axs.scatter( drip_S2n.drip_S2n_N, drip_S2n.drip_S2n_Z, marker='o', s = 3, linewidth=0, color = 'green' )
+        s2p = theo.S2p( Nmin=1, Nmax = 150 )
+        drip_S2p = s2p.drip_S2p( Nmin = 1, Nmax = 150 )
+        axs.scatter( drip_S2p.drip_S2p_N, drip_S2p.drip_S2p_Z, marker='o', s = 3, linewidth=0, color = 'green' )
+    #
+    # First and last isotopes
+    #
+    #iso = ustbl.isotopes( Zmin=1, Zmax = 95 )
+    #axs.scatter( iso.isotopes_Nmin, iso.isotopes_Z, marker='s', s = 3, linewidth=0, color = 'green', label='Isotope bounds' )
+    #axs.scatter( iso.isotopes_Nmax, iso.isotopes_Z, marker='s', s = 3, linewidth=0, color = 'green' )
+    #
+    # stable nuclei:
+    #
+    mas = nuda.nuc.setupBEExp( table = table, version = version )
+    stbl = mas.select( state= 'gs', interp = 'n', nucleus = 'stable' )
+    axs.scatter( stbl.sel_nucN, stbl.sel_nucZ, marker='s', s = 3, linewidth=0, color = 'k' )
+    axs.text(10,112,'stable: '+str(stbl.sel_nbNucSel))
+    #
+    axs.text(60,120,str(ustbl.sel_nbNucSel+stbl.sel_nbNucSel))
+    #
+    # plot N=Z dotted line
+    #
+    axs.plot( [0, 130], [0, 130], linestyle='dotted', linewidth=1, color='k')
+    axs.text(105,120,'N=Z')
+    #
+    # plot stable_fit
+    #
+    N, Z = nuda.stable_fit( Zmin = 1, Zmax = 120 )
+    axs.plot( N, Z, linestyle='dashed', linewidth=1, color='k')
+    #
+    # plot shells for isotopes and isotones
+    #
+    axs = nuda.nuc.plot_shells(axs)
+    #
+    # set legend
+    #
+    axs.legend(loc='lower right',fontsize='10')
+    #
+    # set plot name and close
+    #
+    if pname is not None: 
+        plt.savefig(pname, dpi=300)
+        plt.close()
+    #

nucleardatapy/hnuc/__init__.py

@@ -0,0 +1,8 @@
+"""
+This module provides microscopic, phenomenological and experimental data constraints.
+"""
+#
+from nucleardatapy.hnuc.setup_be1L_exp     import *
+from nucleardatapy.hnuc.setup_be2L_exp     import *
+from nucleardatapy.hnuc.setup_be1Xi_exp    import *
+#