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

nucleardatapy/fig/nuc_setupBEExp_chart_fig.py

@@ -0,0 +1,275 @@
+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()
+    #
+
+def nuc_setupBEExp_chart_Rch_fig( pname, table, version, Rch_table ):
+    """
+    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+' + '+Rch_table+' charge radius table')
+    axs.set_ylabel(r'Z')
+    axs.set_xlabel(r'N')
+    axs.set_xlim([0, 200])
+    axs.set_xticks([0,25,50,75,100,125,150,175,200])
+    axs.set_ylim([0, 132])
+    axs.set_yticks([0,20,40,60,80,100,120])
+    axs.text(10,120,'Number of nuclei:')
+    #
+    # First and last isotopes
+    #
+    mas = nuda.nuc.setupBEExp( table = table, version = version )
+    ustbl = mas.select( state= 'gs', interp = 'n', nucleus = 'unstable' )
+    iso = ustbl.isotopes( Zmin=1, Zmax = 95 )
+    axs.scatter( iso.isotopes_Nmin, iso.isotopes_Z, marker='s', s = 3, linewidth=0, color = 'green', label='AME boundaries' )
+    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))
+    #
+    # Charge radii:
+    #
+    rch = nuda.nuc.setupRchExp( table = Rch_table )
+    axs.scatter( rch.nucN, rch.nucZ, marker='s', s = 3, linewidth=0, color = 'blue', label = 'charge radii' )
+    #
+    # 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()
+    #
+
+
+def nuc_setupBEExp_chart_year_fig( pname, sYear, year_min, year_max ):
+    """
+    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:')
+    axs.text(10,100,'Year: '+str(year_min)+' to '+str(year_max))
+    #
+    # nuclei given in argument: sYear
+    #
+    axs.scatter( sYear.sel_N, sYear.sel_Z, marker='s', s = 3, linewidth=0, color = 'grey' )
+    #
+    # 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' )
+    #
+    # 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/fig/nuc_setupBEExp_fig.py

@@ -3,7 +3,7 @@ import matplotlib.pyplot as plt
 
 import nucleardatapy as nuda
 
-def nuc_setupBEExp_chart_lt_fig( pname, table, version, theo_tables ):
+def nuc_setupBEExp_year_fig( pname, table, version ):
     """
     Plot nuclear chart (N versus Z).\
     The plot is 1x1 with:\
@@ -20,103 +20,258 @@ def nuc_setupBEExp_chart_lt_fig( pname, table, version, theo_tables ):
 
     """
     #
+    print(50*'-')
+    print("Enter nuc_setupBEExp_year_fig.py:")
+    print(50*'-')
+    #
     print(f'Plot name: {pname}')
     #
     print('Table:',table)
     #
-    # plot
+    # read all the mass table:
+    #
+    mas = nuda.nuc.setupBEExp( table = table, version = version )
+    #
+    fig, axs = plt.subplots(1,2)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.15, right=None, top=0.85, wspace=0.2, hspace=0.2)
+    #
+    axs[0].set_title(r''+table+' mass table version '+version)
+    axs[0].set_ylabel(r'number of discovered nuclei')
+    axs[0].set_xlabel(r'year')
+    axs[0].set_xlim([1890, 2020])
+    #axs.set_yscale('log')
+    axs[0].set_ylim([0, 250])
+    #axs.text(10,120,'Number of nuclei:')
+    #
+    axs[0].hist( mas.nucYear, bins=100 )
+    #axs.hist( mas.year, bins=100, linestyle='solid', linewidth=1, color='k')
+    #axs.plot( mas.dist_year*10, mas.dist_nbNuc, linestyle='solid', linewidth=1, color='k')
+    #
+    axs[1].set_title(r''+table+' mass table version '+version)
+    axs[1].set_xlabel(r'year')
+    axs[1].set_xlim([2000, 2020])
+    axs[1].set_ylim([0, 100])
+    axs[1].hist( mas.nucYear, bins=100 )
+    #
+    #axs.legend(loc='lower right',fontsize='10')
+    #
+    if pname is not None: 
+        plt.savefig(pname, dpi=200)
+        plt.close()
+    #
+    print(50*'-')
+    print("Exit nuc_setupBEExp_year_fig.py:")
+    print(50*'-')
+    #
+
+def nuc_setupBEExp_S2n_fig( pname, tables, versions, Zref = 50 ):
+    """
+    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('Tables:',tables)
+    print('Zref:',Zref)
     #
     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)
+    fig.subplots_adjust(left=0.14, 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:')
+    axs.set_title(r'Experimental mass tables',fontsize='12')
+    axs.set_ylabel(r'$S_{2n}$ (MeV)',fontsize='12')
+    axs.set_xlabel(r'N',fontsize='12')
+    axs.set_xlim([Zref-5, int(1.85*Zref)])
+    axs.set_xticks(np.arange(start=Zref-5,stop=2*Zref,step=5))
+    #axs.set_ylim([-10, 10])
+    axs.text(int(Zref),10,'For Z='+str(Zref),fontsize='12')
     #
-    # longlive nuclei
+    # loop over the tables
     #
-    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))
+    for i,table in enumerate( tables ):
+        #
+        version = versions[i]
+        # plot nuclear chart:
+        mas_exp = nuda.nuc.setupBEExp( table = table, version = version )
+        s2n_exp = mas_exp.S2n( Zmin = Zref, Zmax = Zref )
+        axs.scatter( s2n_exp.S2n_N, s2n_exp.S2n, label=table+' '+version )
+        #axs.plot( s2n_exp.S2n_N, s2n_exp.S2n, linestyle='solid', linewidth=1, label=exp_table+' '+exp_version )
     #
-    # shortlive nuclei
+    axs.legend(loc='upper right',fontsize='10', ncol=1)
     #
-    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))
+    if pname is not None: 
+        plt.savefig(pname, dpi=200)
+        plt.close()
     #
-    # veryshortlive nuclei
+
+def nuc_setupBEExp_S2p_fig( pname, tables, versions, Nref = 50 ):
+    """
+    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.
+
+    """
     #
-    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))
+    print(f'Plot name: {pname}')
     #
-    # unstable nuclei:
+    print('Tables:',tables)
+    print('Nref:',Nref)
     #
-    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:
+    print(f'Plot name: {pname}')
     #
-    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))
+    # plot
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.14, bottom=0.15, right=None, top=0.85, wspace=0.3, hspace=0.3)
     #
-    axs.text(60,120,str(ustbl.sel_nbNucSel+stbl.sel_nbNucSel))
+    axs.set_title(r'Experimental mass tables',fontsize='12')
+    axs.set_ylabel(r'$S_{2p}$ (MeV)',fontsize='12')
+    axs.set_xlabel(r'Z',fontsize='12')
+    axs.set_xlim([0.4*Nref, 1.2*Nref])
+    axs.set_xticks(np.arange(start=int(0.4*Nref),stop=1.2*Nref,step=5))
+    #axs.set_ylim([-10, 10])
+    axs.text(int(0.7*Nref),10,'For N='+str(Nref),fontsize='12')
     #
-    # plot N=Z dotted line
+    # loop over the tables
     #
-    axs.plot( [0, 130], [0, 130], linestyle='dotted', linewidth=1, color='k')
-    axs.text(105,120,'N=Z')
+    for i,table in enumerate( tables ):
+        #
+        version = versions[i]
+        # plot nuclear chart:
+        mas_exp = nuda.nuc.setupBEExp( table = table, version = version )
+        s2p_exp = mas_exp.S2p( Nmin = Nref, Nmax = Nref )
+        axs.scatter( s2p_exp.S2p_Z, s2p_exp.S2p, label=table+' '+version )
     #
-    # plot stable_fit
+    axs.legend(loc='upper right',fontsize='10', ncol=1)
     #
-    N, Z = nuda.stable_fit( Zmin = 1, Zmax = 120 )
-    axs.plot( N, Z, linestyle='dashed', linewidth=1, color='k')
+    if pname is not None: 
+        plt.savefig(pname, dpi=200)
+        plt.close()
     #
-    # plot shells for isotopes and isotones
+
+def nuc_setupBEExp_D3n_fig( pname, tables, versions, Zref = 50 ):
+    """
+    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('Tables:',tables)
+    print('Zref:',Zref)
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.14, bottom=0.15, right=None, top=0.85, wspace=0.3, hspace=0.3)
     #
-    axs = nuda.nuc.plot_shells(axs)
+    axs.set_title(r'Experimental mass tables',fontsize='12')
+    axs.set_ylabel(r'$\Delta_{3p}$ (MeV)',fontsize='12')
+    axs.set_xlabel(r'N',fontsize='12')
+    axs.set_xlim([Zref-5, int(1.85*Zref)])
+    axs.set_xticks(np.arange(start=Zref-5,stop=2*Zref,step=5))
+    #axs.set_ylim([-10, 10])
+    axs.text(int(Zref),1.0,'For Z='+str(Zref),fontsize='12')
     #
-    # set legend
+    # loop over the tables
     #
-    axs.legend(loc='lower right',fontsize='10')
+    for i,table in enumerate( tables ):
+        #
+        version = versions[i]
+        # plot nuclear chart:
+        mas_exp = nuda.nuc.setupBEExp( table = table, version = version )
+        d3n_exp = mas_exp.D3n( Zmin = Zref, Zmax = Zref )
+        axs.scatter( d3n_exp.D3n_N_even, d3n_exp.D3n_even, label=table+' '+version+'(even)' )
+        axs.scatter( d3n_exp.D3n_N_odd,  d3n_exp.D3n_odd,  label=table+' '+version+'(odd)' )
     #
-    # set plot name and close
+    axs.legend(loc='upper right',fontsize='10', ncol=1)
     #
     if pname is not None: 
-        plt.savefig(pname, dpi=300)
+        plt.savefig(pname, dpi=200)
         plt.close()
     #
+
+def nuc_setupBEExp_D3p_fig( pname, tables, versions, Nref = 50 ):
+    """
+    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('Tables:',tables)
+    print('Nref:',Nref)
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.14, bottom=0.15, right=None, top=0.85, wspace=0.3, hspace=0.3)
+    #
+    axs.set_title(r'Experimental mass tables',fontsize='12')
+    axs.set_ylabel(r'$\Delta_{3p}$ (MeV)',fontsize='12')
+    axs.set_xlabel(r'Z',fontsize='12')
+    axs.set_xlim([0.4*Nref, 1.2*Nref])
+    axs.set_xticks(np.arange(start=int(0.4*Nref),stop=1.2*Nref,step=5))
+    #axs.set_ylim([-10, 10])
+    axs.text(int(0.7*Nref),1.4,'For N='+str(Nref),fontsize='12')
+    #
+    # loop over the tables
+    #
+    for i,table in enumerate( tables ):
+        #
+        version = versions[i]
+        # plot nuclear chart:
+        mas_exp = nuda.nuc.setupBEExp( table = table, version = version )
+        d3p_exp = mas_exp.D3p( Nmin = Nref, Nmax = Nref )
+        axs.scatter( d3p_exp.D3p_Z_even, d3p_exp.D3p_even, label=table+' '+version+'(even)' )
+        axs.scatter( d3p_exp.D3p_Z_odd,  d3p_exp.D3p_odd,  label=table+' '+version+'(odd)' )
+    #
+    axs.legend(loc='upper right',fontsize='10', ncol=1)
+    #
+    if pname is not None: 
+        plt.savefig(pname, dpi=200)
+        plt.close()
+    #