nucleardatapy 1.0.0__py3-none-any.whl → 1.0.2__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}')
@@ -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}')
@@ -136,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}')
@@ -201,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}')
@@ -272,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}')

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}')

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}')
@@ -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}')
@@ -85,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}')

nucleardatapy/fig/nuc_setupRnpExp_fig.py

@@ -11,6 +11,18 @@ SOURCE_LABELS_LATEX = {
 }
 
 def nuc_setupRnpExp_fig(pname=None, source=None):
+    """
+    Plot the experimental np radii (neutron skin).
+
+    The plot is 1x1 with:
+
+    [0]: Rch as a function of N. 
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param source: experimental table.
+    :type source: str.
+    """
     print(f"Using source: {source}")
     
     # subplot_label = "(a)"

nucleardatapy/fig/nuc_setupRnpTheo_fig.py

@@ -7,8 +7,8 @@ import nucleardatapy as nuda
 
 # Dictionary to map sources to LaTeX names
 SOURCE_LABELS_LATEX = {
-    "48Ca": r"$^{48}\mathrm{Ca}$",
-    "208Pb": r"$^{208}\mathrm{Pb}$"
+    "48Ca": r"$^{48}$Ca",
+    "208Pb": r"$^{208}$Pb"
 }
 
 # Directory containing the model data tables
@@ -42,6 +42,18 @@ def read_model_data(directory, source):
     return model_data
 
 def nuc_setupRnpTheo_fig(pname, source):
+    """
+    Plot the theoretical np radii (neutron skin).
+
+    The plot is 1x1 with:
+
+    [0]: Rch as a function of N. 
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param source: experimental table.
+    :type source: str.
+    """
     print(f'Plot name: {pname}')
     print(f'Using source: {source}')
 

nucleardatapy/hnuc/setup_re2L_exp.py

@@ -16,7 +16,8 @@ def re2L_exp_tables():
     #
     if nuda.env.verb: print("\nEnter re2L_exp_tables()")
     #
-    tables = [ '2013-2L-Ahn' ]
+    #tables = [ '2013-2L-Ahn' ]
+    tables = [ '1991-2L-Yamamoto', '2013-2L-Ahn', '2019-2L-Ekawa' ]
     #
     #print('tables available in the toolkit:',tables)
     tables_lower = [ item.lower() for item in tables ]
@@ -78,7 +79,18 @@ class setupRE2LExp():
       color = []
       mark = []
       #
-      if table.lower() == '2013-2l-ahn':
+      if table.lower() == '1991-2l-yamamoto':
+         #
+         file_in = os.path.join(nuda.param.path_data,'hnuclei/1991-2L-Yamamoto.csv')
+         if nuda.env.verb: print('Reads file:',file_in)
+         #: Attribute providing the full reference to the paper to be citted.
+         self.ref = 'Y. Yamamoto, H. Takaki, K. Ikeda, Prog. Theo. Phys. 86, 867 (1991).'
+         self.keyref = 'YYamamoto:1991'
+         #: Attribute providing additional notes about the data.
+         self.note = "https://doi.org/10.1143/ptp/86.4.867"
+         #
+
+      elif table.lower() == '2013-2l-ahn':
          #
          file_in = os.path.join(nuda.param.path_data,'hnuclei/2013-2L-Ahn.csv')
          if nuda.env.verb: print('Reads file:',file_in)
@@ -88,6 +100,16 @@ class setupRE2LExp():
          #: Attribute providing additional notes about the data.
          self.note = "write here notes about this table."
          #
+      if table.lower() == '2019-2l-ekawa':
+         #
+         file_in = os.path.join(nuda.param.path_data,'hnuclei/2019-2L-Ekawa.csv')
+         if nuda.env.verb: print('Reads file:',file_in)
+         #: Attribute providing the full reference to the paper to be citted.
+         self.ref = 'H. Ekawa et al., PTEP 2019, 021D02.'
+         self.keyref = 'HEkawa:2019'
+         #: Attribute providing additional notes about the data.
+         self.note = "write here notes about this table."
+         #
       #
       with open(file_in,'r') as file:
          for line in file:
@@ -102,11 +124,12 @@ class setupRE2LExp():
                nuclre_err.append(linesplit[4].strip())
                nucldre.append(linesplit[5].strip())
                nucldre_err.append(linesplit[6].strip())
-               probe.append(linesplit[6].strip().strip('\n'))
-               if probe[-1] == 'emul':
-                  label.append("Ahn-2013 Emul")
-                  color.append('blue')
-                  mark.append('s')
+               label.append(linesplit[7].strip('\n'))
+               probe.append(linesplit[7].strip().strip('\n'))
+               color.append('blue')
+               mark.append('s')
+               print('label:',label[-1])
+               print('probe:',probe[-1])
             else:
                break
       #