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

nucleardatapy/fig/nuc_setupBEExp_chart_fig.py

@@ -5,8 +5,10 @@ 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:\
+    Plot nuclear chart (N versus Z).
+
+    The plot is 1x1 with:
+
     [0]: nuclear chart.
 
     :param pname: name of the figure (*.png)
@@ -17,7 +19,6 @@ def nuc_setupBEExp_chart_lt_fig( pname, table, version, theo_tables ):
     :type version: str.
     :param theo_tables: object instantiated on the reference band.
     :type theo_tables: object.
-
     """
     #
     print(f'Plot name: {pname}')
@@ -31,8 +32,8 @@ def nuc_setupBEExp_chart_lt_fig( pname, table, version, theo_tables ):
     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_ylabel(r'Z',fontsize='14')
+    axs.set_xlabel(r'N',fontsize='14')
     axs.set_xlim([0, 200])
     axs.set_ylim([0, 132])
     axs.text(10,120,'Number of nuclei:')
@@ -41,21 +42,28 @@ def nuc_setupBEExp_chart_lt_fig( pname, table, version, theo_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.scatter( ustbl.sel_nucN, ustbl.sel_nucZ, marker='s', s = 3, linewidth=0, color = 'grey', label='long-lived ('+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.scatter( ustbl.sel_nucN, ustbl.sel_nucZ, marker='s', s = 3, linewidth=0, color = 'r', label='short-lived ('+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.scatter( ustbl.sel_nucN, ustbl.sel_nucZ, marker='s', s = 3, linewidth=0, color = 'b', label='very-short-lived ('+str(ustbl.sel_nbNucSel)+')' )
+    #axs.text(10,80,'very short live: '+str(ustbl.sel_nbNucSel))
+    #
+    # hypershortlive nuclei
+    #
+    mas = nuda.nuc.setupBEExp( table = table, version = version )
+    ustbl = mas.select( state= 'gs', interp = 'n', nucleus = 'hypershortlive' )
+    axs.scatter( ustbl.sel_nucN, ustbl.sel_nucZ, marker='s', s = 3, linewidth=0, color = 'g', label='hyper-short-lived ('+str(ustbl.sel_nbNucSel)+')' )
     #axs.text(10,80,'very short live: '+str(ustbl.sel_nbNucSel))
     #
     # unstable nuclei:
@@ -70,16 +78,20 @@ def nuc_setupBEExp_chart_lt_fig( pname, table, version, theo_tables ):
     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' )
+        for zref in range(1,96):
+            itp = theo.isotopes( Zref= zref)
+            s2n = itp.S2n( Zref= zref )
+            drip_S2n = itp.drip_S2n(Zref = zref)
+            if legend == 0:
+                axs.scatter( drip_S2n.drip_S2n_N, drip_S2n.drip_S2n_Z, marker='o', s = 3, linewidth=0, color = 'purple', 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 = 'purple' )
+        # for nref in range(1,170):
+        #     itn = theo.isotones( Nref= nref)
+        #     s2p = itn.S2p( Nref= nref)
+        # drip_S2p = theo.drip_S2p( Nmin=1, Nmax=169)
+        # axs.scatter( drip_S2p.drip_S2p_N, drip_S2p.drip_S2p_Z, marker='o', s = 3, linewidth=0, color = 'purple' )
     #
     # First and last isotopes
     #
@@ -94,7 +106,7 @@ def nuc_setupBEExp_chart_lt_fig( pname, table, version, theo_tables ):
     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))
+    axs.text(49,120,str(ustbl.sel_nbNucSel+stbl.sel_nbNucSel))
     #
     # plot N=Z dotted line
     #
@@ -103,8 +115,10 @@ def nuc_setupBEExp_chart_lt_fig( pname, table, version, theo_tables ):
     #
     # plot stable_fit
     #
-    N, Z = nuda.stable_fit( Zmin = 1, Zmax = 120 )
+    N, Z = nuda.stable_fit_Z( Zmin = 1, Zmax = 120 )
     axs.plot( N, Z, linestyle='dashed', linewidth=1, color='k')
+    #N, Z = nuda.stable_fit_N( Nmin = 1, Nmax = 200 )
+    #axs.plot( N, Z, linestyle='dashed', linewidth=1, color='r')
     #
     # plot shells for isotopes and isotones
     #
@@ -121,26 +135,27 @@ def nuc_setupBEExp_chart_lt_fig( pname, table, version, theo_tables ):
         plt.close()
     #
 
-def nuc_setupBEExp_chart_Rch_fig( pname, table, version, Rch_table ):
+def nuc_setupBEExp_chart_year_fig( pname, sYear, year_min, year_max ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
+    Plot nuclear chart (N versus Z) for a range of discovery years.
+
+    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.
-
+    :param sYear: select nuclei for given discovery years.
+    :type sYear: object.
+    :param year_min: lower range of the discovery years.
+    :type year_min: real.
+    :param year_max: upper range of the discovery years.
+    :type year_max: real.
     """
     #
     print(f'Plot name: {pname}')
     #
-    print('Table:',table)
+    #print('Table:',table)
     #
     # plot
     #
@@ -148,35 +163,23 @@ def nuc_setupBEExp_chart_Rch_fig( pname, table, version, Rch_table ):
     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_title(r''+table+' mass table version '+version)
     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' )
+    #axs.text(10,120,'Number of nuclei:')
+    axs.text(10,100,'Year: '+str(year_min)+' to '+str(year_max))
     #
-    # stable nuclei:
+    # nuclei given in argument: sYear
     #
-    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))
+    axs.scatter( sYear.sel_N, sYear.sel_Z, marker='s', s = 3, linewidth=0, color = 'grey' )
     #
-    # Charge radii:
+    # First and last isotopes
     #
-    rch = nuda.nuc.setupRchExp( table = Rch_table )
-    axs.scatter( rch.nucN, rch.nucZ, marker='s', s = 3, linewidth=0, color = 'blue', label = 'charge radii' )
+    #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
     #
@@ -185,8 +188,10 @@ def nuc_setupBEExp_chart_Rch_fig( pname, table, version, Rch_table ):
     #
     # plot stable_fit
     #
-    N, Z = nuda.stable_fit( Zmin = 1, Zmax = 120 )
+    N, Z = nuda.stable_fit_Z( Zmin = 1, Zmax = 120 )
     axs.plot( N, Z, linestyle='dashed', linewidth=1, color='k')
+    N, Z = nuda.stable_fit_N( Nmin = 1, Nmax = 170 )
+    axs.plot( N, Z, linestyle='dashed', linewidth=1, color='r')
     #
     # plot shells for isotopes and isotones
     #
@@ -194,7 +199,7 @@ def nuc_setupBEExp_chart_Rch_fig( pname, table, version, Rch_table ):
     #
     # set legend
     #
-    axs.legend(loc='lower right',fontsize='10')
+    #axs.legend(loc='lower right',fontsize='10')
     #
     # set plot name and close
     #
@@ -203,11 +208,12 @@ def nuc_setupBEExp_chart_Rch_fig( pname, table, version, Rch_table ):
         plt.close()
     #
 
-
-def nuc_setupBEExp_chart_year_fig( pname, sYear, year_min, year_max ):
+def nuc_setupBEExp_chart_Rch_fig( pname, table, version, Rch_table ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
+    Plot nuclear chart (N versus Z).
+
+    The plot is 1x1 with:
+
     [0]: nuclear chart.
 
     :param pname: name of the figure (*.png)
@@ -216,14 +222,14 @@ def nuc_setupBEExp_chart_year_fig( pname, sYear, year_min, year_max ):
     :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.
+    :param Rch_table: table for Rch.
+    :type Rch_table: str.
 
     """
     #
     print(f'Plot name: {pname}')
     #
-    #print('Table:',table)
+    print('Table:',table)
     #
     # plot
     #
@@ -231,23 +237,36 @@ def nuc_setupBEExp_chart_year_fig( pname, sYear, year_min, year_max ):
     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_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.text(10,120,'Number of nuclei:')
-    axs.text(10,100,'Year: '+str(year_min)+' to '+str(year_max))
+    axs.set_yticks([0,20,40,60,80,100,120])
+    axs.text(10,120,'Number of nuclei:')
     #
-    # nuclei given in argument: sYear
+    # First and last isotopes
     #
-    axs.scatter( sYear.sel_N, sYear.sel_Z, marker='s', s = 3, linewidth=0, color = 'grey' )
+    mas = nuda.nuc.setupBEExp( table = table, version = version )
+    ustbl = mas.select( state= 'gs', interp = 'n', nucleus = 'unstable' )
+    for Zref in range(1,95):
+        iso = ustbl.isotopes( Zref = Zref )
+        axs.scatter( iso.itp_nucNmin, iso.itp_nucZ, marker='s', s = 3, linewidth=0, color = 'green' )
+        axs.scatter( iso.itp_nucNmax, iso.itp_nucZ, marker='s', s = 3, linewidth=0, color = 'green' )
     #
-    # First and last isotopes
+    # stable nuclei:
     #
-    #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' )
+    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
     #
@@ -256,8 +275,8 @@ def nuc_setupBEExp_chart_year_fig( pname, sYear, year_min, year_max ):
     #
     # plot stable_fit
     #
-    N, Z = nuda.stable_fit( Zmin = 1, Zmax = 120 )
-    axs.plot( N, Z, linestyle='dashed', linewidth=1, color='k')
+    #N, Z = nuda.stable_fit_N( Zmin = 1, Zmax = 120 )
+    #axs.plot( N, Z, linestyle='dashed', linewidth=1, color='k')
     #
     # plot shells for isotopes and isotones
     #
@@ -265,7 +284,7 @@ def nuc_setupBEExp_chart_year_fig( pname, sYear, year_min, year_max ):
     #
     # set legend
     #
-    #axs.legend(loc='lower right',fontsize='10')
+    axs.legend(loc='lower right',fontsize='10')
     #
     # set plot name and close
     #
@@ -273,3 +292,4 @@ def nuc_setupBEExp_chart_year_fig( pname, sYear, year_min, year_max ):
         plt.savefig(pname, dpi=300)
         plt.close()
     #
+

nucleardatapy/fig/nuc_setupBEExp_fig.py

@@ -5,9 +5,11 @@ import nucleardatapy as nuda
 
 def nuc_setupBEExp_year_fig( pname, table, version ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot the histogram for the discovery year.
+
+    The plot is 1x2 with:
+
+    [0]: full range of years. [1]: last two decades.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
@@ -15,9 +17,6 @@ def nuc_setupBEExp_year_fig( pname, table, version ):
     :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(50*'-')
@@ -34,11 +33,11 @@ def nuc_setupBEExp_year_fig( pname, table, 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)
+    fig.subplots_adjust(left=0.12, bottom=0.15, right=0.95, 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_ylabel(r'number of discovered nuclei',fontsize='14')
+    axs[0].set_xlabel(r'year',fontsize='14')
     axs[0].set_xlim([1890, 2020])
     #axs.set_yscale('log')
     axs[0].set_ylim([0, 250])
@@ -49,7 +48,7 @@ def nuc_setupBEExp_year_fig( pname, table, version ):
     #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_xlabel(r'year',fontsize='14')
     axs[1].set_xlim([2000, 2020])
     axs[1].set_ylim([0, 100])
     axs[1].hist( mas.nucYear, bins=100 )
@@ -67,19 +66,20 @@ def nuc_setupBEExp_year_fig( pname, table, version ):
 
 def nuc_setupBEExp_S2n_fig( pname, tables, versions, Zref = 50 ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot S2n from the data extracted from the nuclear chart.
+
+    The plot is 1x1 with:
+
+    [0]: show S2n as a function of N.
 
     :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.
-
+    :param tables: tables.
+    :type tables: array of str.
+    :param versions: versions of the tables.
+    :type versions: array of str.
+    :param Zref: Reference value for Z.
+    :type Zref: integer.
     """
     #
     print(f'Plot name: {pname}')
@@ -89,26 +89,27 @@ def nuc_setupBEExp_S2n_fig( pname, tables, versions, Zref = 50 ):
     #
     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)
+    fig.subplots_adjust(left=0.14, bottom=0.15, right=0.95, top=0.85, wspace=0.3, hspace=0.3)
     #
-    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_title(r'Experimental mass tables',fontsize='14')
+    axs.set_ylabel(r'$S_{2n}$ (MeV)',fontsize='14')
+    axs.set_xlabel(r'N',fontsize='14')
     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')
+    axs.text(int(Zref),10,'For Z='+str(Zref),fontsize='14')
     #
     # 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 )
-        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 )
+        mas_exp2 = mas_exp.select( state = 'gs', interp = 'n' )
+        mas_exp3 = mas_exp2.isotopes( Zref = Zref )
+        s2n_exp = mas_exp3.S2n( Zref = Zref )
+        axs.errorbar( s2n_exp.S2n_N, s2n_exp.S2n_E, yerr=s2n_exp.S2n_E_err, fmt='o', label=table+' '+version )
+        #axs.scatter( s2n_exp.S2n_N, s2n_exp.S2n_E, label=exp_table+' '+exp_version )
     #
     axs.legend(loc='upper right',fontsize='10', ncol=1)
     #
@@ -119,19 +120,20 @@ def nuc_setupBEExp_S2n_fig( pname, tables, versions, Zref = 50 ):
 
 def nuc_setupBEExp_S2p_fig( pname, tables, versions, Nref = 50 ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot S2p from the data extracted from the nuclear chart.
+
+    The plot is 1x1 with:
+
+    [0]: show S2p as a function of Z.
 
     :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.
-
+    :param tables: tables.
+    :type tables: array of str.
+    :param versions: versions of the tables.
+    :type versions: array of str.
+    :param Nref: Reference value for N.
+    :type Nref: integer.
     """
     #
     print(f'Plot name: {pname}')
@@ -139,31 +141,29 @@ def nuc_setupBEExp_S2p_fig( pname, tables, versions, Nref = 50 ):
     print('Tables:',tables)
     print('Nref:',Nref)
     #
-    print(f'Plot name: {pname}')
-    #
-    # 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)
+    fig.subplots_adjust(left=0.14, bottom=0.15, right=0.95, top=0.85, wspace=0.3, hspace=0.3)
     #
-    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_title(r'Experimental mass tables',fontsize='14')
+    axs.set_ylabel(r'$S_{2p}$ (MeV)',fontsize='14')
+    axs.set_xlabel(r'Z',fontsize='14')
     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')
+    axs.text(int(0.7*Nref),10,'For N='+str(Nref),fontsize='14')
     #
     # 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 )
-        s2p_exp = mas_exp.S2p( Nmin = Nref, Nmax = Nref )
-        axs.scatter( s2p_exp.S2p_Z, s2p_exp.S2p, label=table+' '+version )
+        mas_exp2 = mas_exp.select( state = 'gs', interp = 'n' )
+        mas_exp3 = mas_exp2.isotones( Nref = Nref )
+        s2p_exp = mas_exp3.S2p( Nref = Nref )
+        axs.errorbar( s2p_exp.S2p_Z, s2p_exp.S2p_E, yerr=s2p_exp.S2p_E_err, fmt='o', label=table+' '+version )
+        #axs.scatter( s2p_exp.S2p_Z, s2p_exp.S2p_E, label=table+' '+version )
     #
     axs.legend(loc='upper right',fontsize='10', ncol=1)
     #
@@ -174,19 +174,20 @@ def nuc_setupBEExp_S2p_fig( pname, tables, versions, Nref = 50 ):
 
 def nuc_setupBEExp_D3n_fig( pname, tables, versions, Zref = 50 ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot D3n (3-point formula for the odd-even mass staggering) from the data extracted from the nuclear chart.
+
+    The plot is 1x1 with:
+
+    [0]: show D3n as a function of N.
 
     :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.
-
+    :param tables: tables.
+    :type tables: array of str.
+    :param versions: versions of the tables.
+    :type versions: array of str.
+    :param Zref: Reference value for Z.
+    :type Zref: integer.
     """
     #
     print(f'Plot name: {pname}')
@@ -196,15 +197,15 @@ def nuc_setupBEExp_D3n_fig( pname, tables, versions, Zref = 50 ):
     #
     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)
+    fig.subplots_adjust(left=0.14, bottom=0.15, right=0.95, 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'N',fontsize='12')
+    axs.set_title(r'Experimental mass tables',fontsize='14')
+    axs.set_ylabel(r'$\Delta_{3n}$ (MeV)',fontsize='14')
+    axs.set_xlabel(r'N',fontsize='14')
     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')
+    axs.text(int(Zref),1.0,'For Z='+str(Zref),fontsize='14')
     #
     # loop over the tables
     #
@@ -213,9 +214,13 @@ def nuc_setupBEExp_D3n_fig( pname, tables, versions, Zref = 50 ):
         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)' )
+        mas_exp2 = mas_exp.select( state = 'gs', interp = 'n' )
+        mas_exp3 = mas_exp2.isotopes( Zref = Zref )
+        D3n_exp = mas_exp3.D3n( Zref = Zref )
+        axs.errorbar( D3n_exp.D3n_N_even, D3n_exp.D3n_E_even, yerr=D3n_exp.D3n_E_err_even, fmt='o', label=table+' '+version )
+        axs.errorbar( D3n_exp.D3n_N_odd,  D3n_exp.D3n_E_odd, yerr=D3n_exp.D3n_E_err_odd, fmt='o', label=table+' '+version )
+        #axs.scatter( d3n_exp.D3n_N_even, d3n_exp.D3n_E_even, label=table+' '+version+'(even)' )
+        #axs.scatter( d3n_exp.D3n_N_odd,  d3n_exp.D3n_E_odd,  label=table+' '+version+'(odd)' )
     #
     axs.legend(loc='upper right',fontsize='10', ncol=1)
     #
@@ -226,19 +231,20 @@ def nuc_setupBEExp_D3n_fig( pname, tables, versions, Zref = 50 ):
 
 def nuc_setupBEExp_D3p_fig( pname, tables, versions, Nref = 50 ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot D3p (3-point formula for the odd-even mass staggering) from the data extracted from the nuclear chart.
+
+    The plot is 1x1 with:
+
+    [0]: show D3p as a function of Z.
 
     :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.
-
+    :param tables: tables.
+    :type tables: array of str.
+    :param versions: versions of the tables.
+    :type versions: array of str.
+    :param Nref: Reference value for N.
+    :type Nref: integer.
     """
     #
     print(f'Plot name: {pname}')
@@ -248,15 +254,15 @@ def nuc_setupBEExp_D3p_fig( pname, tables, versions, Nref = 50 ):
     #
     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)
+    fig.subplots_adjust(left=0.14, bottom=0.15, right=0.95, 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_title(r'Experimental mass tables',fontsize='14')
+    axs.set_ylabel(r'$\Delta_{3p}$ (MeV)',fontsize='14')
+    axs.set_xlabel(r'Z',fontsize='14')
     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')
+    axs.text(int(0.7*Nref),1.4,'For N='+str(Nref),fontsize='14')
     #
     # loop over the tables
     #
@@ -265,9 +271,13 @@ def nuc_setupBEExp_D3p_fig( pname, tables, versions, Nref = 50 ):
         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)' )
+        mas_exp2 = mas_exp.select( state = 'gs', interp = 'n' )
+        mas_exp3 = mas_exp2.isotones( Nref = Nref )
+        D3p_exp = mas_exp3.D3p( Nref = Nref )
+        axs.errorbar( D3p_exp.D3p_Z_even, D3p_exp.D3p_E_even, yerr=D3p_exp.D3p_E_err_even, fmt='o', label=table+' '+version )
+        axs.errorbar( D3p_exp.D3p_Z_odd,  D3p_exp.D3p_E_odd,  yerr=D3p_exp.D3p_E_err_odd,  fmt='o', label=table+' '+version )
+        #axs.scatter( d3p_exp.D3p_Z_even, d3p_exp.D3p_E_even, label=table+' '+version+'(even)' )
+        #axs.scatter( d3p_exp.D3p_Z_odd,  d3p_exp.D3p_E_odd,  label=table+' '+version+'(odd)' )
     #
     axs.legend(loc='upper right',fontsize='10', ncol=1)
     #