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

nucleardatapy/fig/nuc_setupBETheo_fig.py

@@ -5,19 +5,20 @@ import nucleardatapy as nuda
 
 def nuc_setupBETheo_diff_fig( pname, tables, table_ref = '1995-DZ', Zref = 50 ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot the energy difference between theoretical mass tables with respect to `table_ref`.
+
+    The plot is 1x1 with:
+
+    [0]: show the energy difference 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 table_ref: reference table.
+    :type table_ref: str.
+    :param Zref: Reference value for Z.
+    :type Zref: integer.
     """
     #
     print(f'Plot name: {pname}')
@@ -34,16 +35,16 @@ def nuc_setupBETheo_diff_fig( pname, tables, table_ref = '1995-DZ', Zref = 50 ):
     fig.subplots_adjust(left=0.15, bottom=0.13, right=None, top=0.8)#, wspace=0.3, hspace=0.3)
     #
     #axs.set_title(r'Comparison of theoretical mass models',fontsize='12')
-    axs.set_ylabel(r'$E-E_{DZ}$ (MeV)',fontsize='12')
-    axs.set_xlabel(r'N',fontsize='12')
+    axs.set_ylabel(r'$E-E_{DZ}$ (MeV)',fontsize='14')
+    axs.set_xlabel(r'N',fontsize='14')
     axs.set_ylim([-5, 5])
     #axs.text(int(Zref)+5,-7,'For Z='+str(Zref),fontsize='12')
     if Zref == 50:
         axs.set_xlim( [ 40, 100 ] )
-        axs.text(55,4,'For Z='+str(Zref),fontsize='12')
+        axs.text(55,4,'For Z='+str(Zref),fontsize='14')
     elif Zref == 82:
         axs.set_xlim( [ 90, 150 ] )
-        axs.text(110,4,'For Z='+str(Zref),fontsize='12')
+        axs.text(110,4,'For Z='+str(Zref),fontsize='14')
     #
     # loop over the tables
     #
@@ -68,19 +69,20 @@ def nuc_setupBETheo_diff_fig( pname, tables, table_ref = '1995-DZ', Zref = 50 ):
 
 def nuc_setupBETheo_S2n_fig( pname, tables, Zref = 50 ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot S2n from the data extracted from the theoretical mass table.
+
+    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}')
@@ -95,27 +97,33 @@ def nuc_setupBETheo_S2n_fig( pname, tables, Zref = 50 ):
     fig.subplots_adjust(left=0.15, bottom=0.13, right=None, top=0.8)#, wspace=0.3, hspace=0.3)
     #
     #axs.set_title(r'Comparison of theoretical mass models',fontsize='12')
-    axs.set_ylabel(r'$S_{2n}$ (MeV)',fontsize='12')
-    axs.set_xlabel(r'N',fontsize='12')
+    axs.set_ylabel(r'$S_{2n}$ (MeV)',fontsize='14')
+    axs.set_xlabel(r'N',fontsize='14')
     axs.set_xlim([Zref-5, int(2.3*Zref)])
     axs.set_xticks(np.arange(start=Zref-5,stop=2.3*Zref,step=5))
-    axs.set_ylim([0, 40])
-    axs.text(int(Zref),10,'For Z='+str(Zref),fontsize='12')
+    axs.set_ylim([-40, 0])
+    axs.text(int(Zref),-10,'For Z='+str(Zref),fontsize='14')
     #
     # loop over the tables
     #
     for i,table in enumerate( tables ):
         #
         mas = nuda.nuc.setupBETheo( table = table )
-        s2n = mas.S2n( Zmin = Zref, Zmax = Zref )
+        mas2 = mas.isotopes( Zref = Zref )
+        s2n = mas2.S2n( Zref = Zref )
         #
-        axs.plot( s2n.S2n_N, s2n.S2n, linestyle='solid', linewidth=1, label=table )
+        axs.plot( s2n.S2n_N, s2n.S2n_E, linestyle='solid', linewidth=1, label=table )
     #
     exp_table = 'AME'
     exp_version = '2020'
+    #
     mas_exp = nuda.nuc.setupBEExp( table = exp_table, version = exp_version )
-    s2n_exp = mas_exp.S2n( Zmin = Zref, Zmax = Zref )
-    axs.scatter( s2n_exp.S2n_N, s2n_exp.S2n, 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=exp_table+' '+exp_version )
+    #axs.scatter( s2n_exp.S2n_N, s2n_exp.S2n_E, label=exp_table+' '+exp_version )
     #axs.plot( s2n_exp.S2n_N, s2n_exp.S2n, linestyle='solid', linewidth=1, label=exp_table+' '+exp_version )
     #N_diff, A_diff, BE_diff, BE_diff = mas.diff_exp( table_exp = 'AME', version_exp = '2020', Zref = Zref )
     #axs.scatter( N_diff, BE_diff, label='AME2020' )
@@ -130,19 +138,20 @@ def nuc_setupBETheo_S2n_fig( pname, tables, Zref = 50 ):
 
 def nuc_setupBETheo_S2p_fig( pname, tables, Nref = 50 ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot S2p from the data extracted from the theoretical mass table.
+
+    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}')
@@ -156,29 +165,34 @@ def nuc_setupBETheo_S2p_fig( pname, tables, Nref = 50 ):
     fig.subplots_adjust(left=0.15, bottom=0.13, right=None, top=0.8)#, wspace=0.3, hspace=0.3)
     #
     #axs.set_title(r'Comparison of theoretical mass models',fontsize='12')
-    axs.set_ylabel(r'$S_{2p}$ (MeV)',fontsize='12')
-    axs.set_xlabel(r'Z',fontsize='12')
+    axs.set_ylabel(r'$S_{2p}$ (MeV)',fontsize='14')
+    axs.set_xlabel(r'Z',fontsize='14')
     #axs.set_xlim([0.4*Nref, 1.3*Nref])
     axs.set_xlim([0.5*Nref, 1.05*Nref])
     #axs.set_xticks(np.arange(start=int(0.4*Nref),stop=1.3*Nref,step=5))
     axs.set_xticks(np.arange(start=int(0.5*Nref),stop=1.05*Nref,step=5))
-    axs.set_ylim([0, 46])
-    axs.text(int(0.7*Nref),35,'For N='+str(Nref),fontsize='12')
+    axs.set_ylim([-46, 0])
+    axs.text(int(0.7*Nref),-35,'For N='+str(Nref),fontsize='14')
     #
     # loop over the tables
     #
     for i,table in enumerate( tables ):
         #
         mas = nuda.nuc.setupBETheo( table = table )
-        s2p = mas.S2p( Nmin = Nref, Nmax = Nref )
+        mas2 = mas.isotones( Nref = Nref )
+        s2p = mas2.S2p( Nref = Nref )
         #
-        axs.plot( s2p.S2p_Z, s2p.S2p, linestyle='solid', linewidth=1, label=table )
+        axs.plot( s2p.S2p_Z, s2p.S2p_E, linestyle='solid', linewidth=1, label=table )
     #
     exp_table = 'AME'
     exp_version = '2020'
     mas_exp = nuda.nuc.setupBEExp( table = exp_table, version = exp_version )
-    s2p_exp = mas_exp.S2p( Nmin = Nref, Nmax = Nref )
-    axs.scatter( s2p_exp.S2p_Z, s2p_exp.S2p, label=exp_table+' '+exp_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=exp_table+' '+exp_version )
+    #axs.scatter( s2p_exp.S2p_Z, s2p_exp.S2p_E, label=exp_table+' '+exp_version )
     #
     #axs.legend(loc='upper right',fontsize='10', ncol=4)
     fig.legend(loc='upper left',bbox_to_anchor=(0.15,1.0),columnspacing=2,fontsize='8',ncol=4,frameon=False)
@@ -190,19 +204,20 @@ def nuc_setupBETheo_S2p_fig( pname, tables, Nref = 50 ):
 
 def nuc_setupBETheo_D3n_fig( pname, tables, 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 theoretical mass table.
+
+    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}')
@@ -220,7 +235,7 @@ def nuc_setupBETheo_D3n_fig( pname, tables, Zref = 50 ):
     fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.8)#, wspace=0.3, hspace=0.3)
     #
     #axs.set_title(r'Comparison of theoretical mass models',fontsize='12')
-    axs.set_ylabel(r'$\Delta_{3,n}$ (MeV)',fontsize='12')
+    axs.set_ylabel(r'$\Delta_{3n}$ (MeV)',fontsize='12')
     axs.set_xlabel(r'N',fontsize='12')
     axs.set_xlim([Zref-5, int(2.3*Zref)])
     axs.set_xticks(np.arange(start=Zref-5,stop=2.3*Zref,step=5))
@@ -232,15 +247,24 @@ def nuc_setupBETheo_D3n_fig( pname, tables, Zref = 50 ):
     for i,table in enumerate( tables ):
         #
         mas = nuda.nuc.setupBETheo( table = table )
-        d3p = mas.D3n( Zmin = Zref, Zmax = Zref )
+        mas2 = mas.isotopes( Zref = Zref )
+        D3n = mas2.D3n( Zref = Zref )
         #
-        axs.plot( d3p.D3n_N_even, d3p.D3n_even, linestyle='solid', linewidth=1, label=table+'(even)' )
+        axs.plot( D3n.D3n_N_even, D3n.D3n_E_even, linestyle='solid', linewidth=1, label=table+'(even)' )
     #
     exp_table = 'AME'
     exp_version = '2020'
+    #
     mas_exp = nuda.nuc.setupBEExp( table = exp_table, version = exp_version )
-    d3n_exp = mas_exp.D3n( Zmin = Zref, Zmax = Zref )
-    axs.scatter( d3n_exp.D3n_N_even, d3n_exp.D3n_even, label=exp_table+' '+exp_version+'(even)' )
+    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=exp_table+' '+exp_version+'(even)' )
+    #axs.scatter( D3n_exp.D3n_N_even, D3n_exp.D3n_E_even, label=exp_table+' '+exp_version+'(even)' )
+    #
+    # empirical relations:
+    axs.plot( D3n.D3n_N_even, nuda.nuc.delta_emp(D3n.D3n_N_even,Zref,formula='classic'), linestyle='solid', linewidth=3, label='classic' )
+    axs.plot( D3n.D3n_N_even, nuda.nuc.delta_emp(D3n.D3n_N_even,Zref,formula='Vogel'), linestyle='solid', linewidth=3, label='Vogel' )
     #
     #axs.legend(loc='upper right',fontsize='10', ncol=4)
     fig.legend(loc='upper left',bbox_to_anchor=(0.1,1.0),columnspacing=2,fontsize='7',ncol=4,frameon=False)
@@ -252,19 +276,20 @@ def nuc_setupBETheo_D3n_fig( pname, tables, Zref = 50 ):
 
 def nuc_setupBETheo_D3p_fig( pname, tables, 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 theoretical mass table.
+
+    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}')
@@ -282,7 +307,7 @@ def nuc_setupBETheo_D3p_fig( pname, tables, Nref = 50 ):
     fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.8)#, wspace=0.3, hspace=0.3)
     #
     #axs.set_title(r'Comparison of theoretical mass models',fontsize='12')
-    axs.set_ylabel(r'$\Delta_{3,p}$ (MeV)',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.1*Nref])
     axs.set_xticks(np.arange(start=int(0.4*Nref),stop=1.2*Nref,step=5))
@@ -294,20 +319,28 @@ def nuc_setupBETheo_D3p_fig( pname, tables, Nref = 50 ):
     for i,table in enumerate( tables ):
         #
         mas = nuda.nuc.setupBETheo( table = table )
-        d3p = mas.D3p( Nmin = Nref, Nmax = Nref )
+        mas2 = mas.isotones( Nref = Nref )
+        D3p = mas2.D3p( Nref = Nref )
         #
-        axs.plot( d3p.D3p_Z_even, d3p.D3p_even, linestyle='solid', linewidth=1, label=table+'(even)' )
+        axs.plot( D3p.D3p_Z_even, D3p.D3p_E_even, linestyle='solid', linewidth=1, label=table+'(even)' )
     #
     exp_table = 'AME'
     exp_version = '2020'
+    #
     mas_exp = nuda.nuc.setupBEExp( table = exp_table, version = exp_version )
-    d3p_exp = mas_exp.D3p( Nmin = Nref, Nmax = Nref )
-    axs.scatter( d3p_exp.D3p_Z_even, d3p_exp.D3p_even, label=exp_table+' '+exp_version+'(even)' )
+    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=exp_table+' '+exp_version+'(even)' )
+    #axs.scatter( D3p_exp.D3p_Z_even, D3p_exp.D3p_E_even, label=exp_table+' '+exp_version+'(even)' )
+    #
+    # empirical relations:
+    axs.plot( D3p.D3p_Z_even, nuda.nuc.delta_emp(Nref,D3p.D3p_Z_even,formula='classic'), linestyle='solid', linewidth=3, label='classic' )
+    axs.plot( D3p.D3p_Z_even, nuda.nuc.delta_emp(Nref,D3p.D3p_Z_even,formula='Vogel'), linestyle='solid', linewidth=3, label='Vogel' )
     #
     #axs.legend(loc='upper right',fontsize='8', ncol=4)
     fig.legend(loc='upper left',bbox_to_anchor=(0.1,1.0),columnspacing=2,fontsize='7',ncol=4,frameon=False)
     #
-    #
     if pname is not None: 
         plt.savefig(pname, dpi=200)
         plt.close()

nucleardatapy/fig/nuc_setupISGMRExp_fig.py

@@ -5,19 +5,16 @@ import nucleardatapy as nuda
 
 def nuc_setupISGMRExp_fig( pname, tables ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot the experimental ISGMR energy from the tables.
+
+    The plot is 1x3 with:
+
+    [0]: for Zr. [1]: for Sn. [2]: for Pb.
 
     :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.
     """
     #
     print(f'Plot name: {pname}')
@@ -28,14 +25,14 @@ def nuc_setupISGMRExp_fig( pname, tables ):
     #
     fig, axs = plt.subplots(1,3)
     fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.12, bottom=0.15, right=None, top=0.85, wspace=0.1, hspace=0.3)
+    fig.subplots_adjust(left=0.12, bottom=0.15, right=0.95, top=0.85, wspace=0.1, hspace=0.3)
     #
     if obs == 'M12M0':
-        axs[0].set_ylabel(r'$E_{ISGMR}$ ($m_1/m_0$)')
+        axs[0].set_ylabel(r'$E_{ISGMR}$ from $m_1/m_0$ (MeV)',fontsize='14')
     elif obs == 'M12Mm1':
-        axs[0].set_ylabel(r'$E_{ISGMR}$ ($\sqrt{m_1/m_{-1}}$)')
+        axs[0].set_ylabel(r'$E_{ISGMR}$ from $\sqrt{m_1/m_{-1}}$ (MeV)',fontsize='14')
     elif obs == 'M12Mm1':
-        axs[0].set_ylabel(r'$E_{ISGMR}$ ($\sqrt{m_3/m_1}$)')
+        axs[0].set_ylabel(r'$E_{ISGMR}$ from $\sqrt{m_3/m_1}$ (MeV)',fontsize='14')
     #
     for k,table in enumerate( tables ):
         #
@@ -44,7 +41,7 @@ def nuc_setupISGMRExp_fig( pname, tables ):
         for i in [0,1,2]:
             print('For Z = ',nucZ[i])
             axs[i].set_title(nuda.param.elements[nucZ[i]-1])
-            axs[i].set_xlabel(r'A')
+            axs[i].set_xlabel(r'A',fontsize='14')
             axs[i].set_ylim([13, 18])
             if i>0: axs[i].tick_params('y', labelleft=False)
             gmrs = gmr.select( Zref = nucZ[i], obs = obs )

nucleardatapy/fig/nuc_setupRchExp_fig.py

@@ -5,19 +5,16 @@ import nucleardatapy as nuda
 
 def nuc_setupRchExp_fig( pname, tables ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot the experimental charge radii from the tables.
+
+    The plot is 1x1 with:
+
+    [0]: Rch 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.
     """
     #
     print(f'Plot name: {pname}')
@@ -26,8 +23,8 @@ def nuc_setupRchExp_fig( pname, tables ):
     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_ylabel(r'$R_{ch}$')
-    axs.set_xlabel(r'N')
+    axs.set_ylabel(r'$R_{ch}$ (fm)',fontsize='14')
+    axs.set_xlabel(r'N',fontsize='14')
     axs.set_xlim([10, 140])
     axs.set_ylim([3.2, 5.8])
     #
@@ -53,19 +50,14 @@ def nuc_setupRchExp_fig( pname, tables ):
 
 def nuc_setupRchExp_3Zref_fig( pname, tables ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot the experimental charge radii from the tables.\\
+    The plot is 1x1 with:\\
+    [0]: Rch 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.
     """
     #
     print(f'Plot name: {pname}')

nucleardatapy/fig/nuc_setupRchTheo_fig.py

@@ -5,19 +5,18 @@ import nucleardatapy as nuda
 
 def nuc_setupRchTheo_fig( pname, tables, table_exp ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot the theoretical charge radii from the tables and compare to the experimental value.
+
+    The plot is 1x1 with:
+
+    [0]: Rch 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: theoretical tables.
+    :type tables: array of str.
+    :param table_exp: experimental table.
+    :type table_exp: str.
     """
     #
     print(f'Plot name: {pname}')
@@ -28,28 +27,13 @@ def nuc_setupRchTheo_fig( pname, tables, table_exp ):
     fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
     fig.subplots_adjust(left=0.12, bottom=0.15, right=None, top=0.9, wspace=0.35, hspace=0.3)
     #
-    axs.set_ylabel(r'$R_{ch}$')
-    axs.set_xlabel(r'N')
+    axs.set_ylabel(r'$R_{ch}$ (fm)',fontsize='14')
+    axs.set_xlabel(r'N',fontsize='14')
     axs.set_xlim([10, 160])
     axs.set_ylim([3.0, 6.5])
     #
     Zrefs = [ 20, 28, 40, 50, 60, 70, 82, 90 ]
     #
-    for indZ,Zref in enumerate(Zrefs):
-        rchExpIsot = nuda.nuc.setupRchExpIsotopes( rch_exp, Zref = Zref )
-        if indZ == 0:
-            axs.errorbar( rchExpIsot.N, rchExpIsot.Rch, yerr=rchExpIsot.Rch_err, color='k', fmt='s', markersize=3, label=table_exp )
-            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, 'classic' ), linestyle='dashed', color='k', label='empirical(classic)' )
-            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, '1994-NPP' ), linestyle='dashed', color='red', label='empirical(NPP-1994)' )
-            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, '2013-BAKS-1' ), linestyle='dotted', color='k', label='empirical(BAKS-1-2013)' )
-            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, '2013-BAKS-3' ), linestyle='dotted', color='red', label='empirical(BAKS-3-2013)' )
-        else:
-            axs.errorbar( rchExpIsot.N, rchExpIsot.Rch, yerr=rchExpIsot.Rch_err, color='k', fmt='s', markersize=3 )
-            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, 'classic' ), linestyle='dashed', color='k' )
-            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, '1994-NPP' ), linestyle='dashed', color='red' )
-            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, '2013-BAKS-1' ), linestyle='dotted', color='k' )
-            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, '2013-BAKS-3' ), linestyle='dotted', color='red' )
-    #
     for indT,table in enumerate(tables):
         #
         rch = nuda.nuc.setupRchTheo( table = table )
@@ -66,7 +50,23 @@ def nuc_setupRchTheo_fig( pname, tables, table_exp ):
             #else:
             #    axs.plot( rchIsot.N, nuda.nuc.rch_emp( rchIsot.A, rchIsot.Z, 'classic' ), linestyle='dashed', color='k' )
             #
-    axs.legend(loc='upper left',fontsize='7',frameon=False)
+    #
+    for indZ,Zref in enumerate(Zrefs):
+        rchExpIsot = nuda.nuc.setupRchExpIsotopes( rch_exp, Zref = Zref )
+        if indZ == 0:
+            axs.errorbar( rchExpIsot.N, rchExpIsot.Rch, yerr=rchExpIsot.Rch_err, color='k', fmt='s', markersize=3, label=table_exp )
+            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, 'classic' ), linestyle='dashed', color='k', label='empirical(classic)' )
+            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, '1994-NPP' ), linestyle='dashed', color='red', label='empirical(NPP-1994)' )
+            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, '2013-BAKS-1' ), linestyle='dotted', color='k', label='empirical(BAKS-1-2013)' )
+            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, '2013-BAKS-3' ), linestyle='dotted', color='red', label='empirical(BAKS-3-2013)' )
+        else:
+            axs.errorbar( rchExpIsot.N, rchExpIsot.Rch, yerr=rchExpIsot.Rch_err, color='k', fmt='s', markersize=3 )
+            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, 'classic' ), linestyle='dashed', color='k' )
+            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, '1994-NPP' ), linestyle='dashed', color='red' )
+            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, '2013-BAKS-1' ), linestyle='dotted', color='k' )
+            axs.plot( rchExpIsot.N, nuda.nuc.rch_emp( rchExpIsot.A, rchExpIsot.Z, '2013-BAKS-3' ), linestyle='dotted', color='red' )
+    #
+    axs.legend(loc='lower right',fontsize='10',ncol=2,frameon=False)
     #
     axs.text( 50,3.7,'Ca')
     axs.text( 70,4.2,'Ni')
@@ -84,19 +84,16 @@ def nuc_setupRchTheo_fig( pname, tables, table_exp ):
 
 def nuc_setupRchTheo_3Zref_fig( pname, tables, table_exp ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot the theoretical charge radii from the tables and compare to the experimental value.\\
+    The plot is 1x1 with:\\
+    [0]: Rch 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: theoretical tables.
+    :type tables: array of str.
+    :param table_exp: experimental table.
+    :type table_exp: str.
     """
     #
     print(f'Plot name: {pname}')
@@ -108,8 +105,8 @@ def nuc_setupRchTheo_3Zref_fig( pname, tables, table_exp ):
     fig.subplots_adjust(left=0.12, bottom=0.15, right=None, top=0.9, wspace=0.35, hspace=0.3)
     #
     axs[0].set_title(r'Zr')
-    axs[0].set_ylabel(r'$R_{ch}$')
-    axs[0].set_xlabel(r'A')
+    axs[0].set_ylabel(r'$R_{ch}$ (fm)',fontsize='14')
+    axs[0].set_xlabel(r'A',fontsize='14')
     axs[0].set_xlim([88, 98])
     axs[0].set_ylim([4.25, 4.45])
     #