nucleardatapy 0.2.0__py3-none-any.whl

nucleardatapy/nuc/setup_nskin_theo.py

@@ -0,0 +1,263 @@
+import os
+import sys
+import numpy as np  # 1.15.0
+from scipy.interpolate import CubicSpline
+
+# nucleardatapy_tk = os.getenv('NUCLEARDATAPY_TK')
+# sys.path.insert(0, nucleardatapy_tk)
+
+import nucleardatapy as nuda
+
+def nskin_theo_models():
+    """
+    Return a the neutron skin values of the models available in this toolkit and print them all on the prompt.
+
+    :return: The list of models with can be 'Skyrme', 'NLRH', 'DDRH'.
+    :rtype: list[str].
+    """
+    #
+    if nuda.env.verb: print("\nEnter nskin_theo_models()")
+    #
+    models = [ 'Skyrme', 'NLRH', 'DDRH' ]
+    #print('Phenomenological models available in the toolkit:',models)
+    models_lower = [ item.lower() for item in models ]
+    #
+    if nuda.env.verb: print("Exit nskin_theo_models()")
+    #
+    return models, models_lower
+
+# def nskin_theo_nucleus():
+#     #
+#     if nuda.env.verb: print("\nEnter nskin_theo_nucleus()")
+#     #
+#     nucleus = [ '48Ca', '208Pb']
+#     nucleus_lower = [ item.lower() for item in nucleus ]
+#     #
+#     if nuda.env.verb: print("Exit nskin_theo_nucleus()")
+#     #
+#     return nucleus, nucleus_lower
+# 
+
+def nskin_theo_params( model ):
+    """
+    Return a list with the parameterizations available in 
+    this toolkit for a given model and print them all on the prompt.
+
+    :param model: The type of model for which there are parametrizations. \
+    They should be chosen among the following options: 'Skyrme', 'NLRH', \
+    'DDRH', 'DDRHF'.
+    :type model: str.
+    :return: The list of parametrizations. \
+    If `models` == 'skyrme': 'BSK14', \
+    'BSK16', 'BSK17', 'F-', 'F+', 'F0', 'LNS1', 'LNS5', \
+    'NRAPR', 'RATP', 'SAMI', 'SGII', 'SIII', 'SKGSIGMA', 'SKI2', 'SKI4', 'SKMP', \
+    'SKMS', 'SKO', 'SKOP', 'SKRSIGMA', 'SKX', 'Skz2', 'SLY4', 'SLY5', \
+    'SLY230A', 'SLY230B', 'T6', 'UNEDF0', 'UNEDF1'. \
+    If `models` == 'NLRH': 'NL-SH', 'NL3', 'NL3II', 'PK1', 'TM1'. \
+    If `models` == 'DDRH': 'DDME1', 'DDME2', 'DDMEd', 'PKDD', 'TW99'. \
+    :rtype: list[str].
+    """
+    #
+    if nuda.env.verb: print("\nEnter nskin_theo_params()")
+    #
+    #print('For model:',model)
+    if model.lower() == 'skyrme':
+        params = [ 'BSK14', 'BSK16', 'BSK17', 'F-', 'F+', 'F0',\
+            'LNS1', 'LNS5', 'NRAPR', 'RATP', 'SAMI', 'SGII', 'SIII', \
+            'SKGSIGMA', 'SKI2', 'SKI4', 'SKMP','SKMS', 'SKO', 'SKOP', \
+            'SKRSIGMA', 'SKX', 'Skz2', 'SLY4', 'SLY5', 'SLY230A', \
+            'SLY230B', 'T6', 'UNEDF0', 'UNEDF1']
+        nucleus = ['48Ca', '208Pb']    
+    elif model.lower() == 'nlrh':
+        params = [ 'NL-SH', 'NL3', 'NL3II', 'PK1', 'TM1' ]
+        nucleus = ['48Ca', '208Pb']
+    elif model.lower() == 'ddrh':
+        params = [ 'DDME1', 'DDME2', 'DDMEd', 'PKDD', 'TW99' ]
+        nucleus = ['48Ca', '208Pb']
+    #print('Parameters available in the toolkit:',params)
+    params_lower = [ item.lower() for item in params ]
+    #
+    if nuda.env.verb: print("Exit nskin_theo_params()")
+    #
+    return params, params_lower
+
+class SetupNeutronSkinTheo():
+    """
+    Instantiate the object with results based on phenomenological\
+    interactions and choosen by the toolkit practitioner. \
+    This choice is defined in the variables `model` and `param`.
+
+    If `models` == 'skyrme', `param` can be: 'BSK14', \
+    'BSK16', 'BSK17', 'F-', 'F+', 'F0', 'LNS1', 'LNS5', \
+    'NRAPR', 'RATP', 'SAMI', 'SGII', 'SIII', 'SKGSIGMA', 'SKI2', 'SKI4', 'SKMP', \
+    'SKMS', 'SKO', 'SKOP', 'SKRSIGMA', 'SKX', 'Skz2', 'SLY4', 'SLY5', \
+    'SLY230A', 'SLY230B', 'T6', 'UNEDF0', 'UNEDF1'. 
+
+    If `models` == 'NLRH', `param` can be: 'NL-SH', 'NL3', 'NL3II', 'PK1', 'TM1'. 
+
+    If `models` == 'DDRH', `param` can be: 'DDME1', 'DDME2', 'DDMEd', 'PKDD', 'TW99'. 
+
+    
+    :param model: Fix the name of model: 'Skyrme', 'NLRH', \
+    'DDRH', 'DDRHF'. Default value: 'Skyrme'.
+    :type model: str, optional. 
+    :param param: Fix the parameterization associated to model. \
+    Default value: 'SLY5'.
+    :type param: str, optional. 
+
+    **Attributes:**
+    """
+    #
+    def __init__( self, model = 'Skyrme', param = 'SLY5', nucleus = '208Pb' ):
+        #
+        if nuda.env.verb: print("\nEnter SetupNeutronSkinTheo()")
+        #
+        #: Attribute model.
+        self.model = model
+        if nuda.env.verb: print("model:",model)
+        print("-> model:",model)
+        #: Attribute param.
+        self.param = param
+        if nuda.env.verb: print("param:",param)
+        print("-> param:",param)
+        #
+        #: Attribute nucleus.
+        self.nucleus = nucleus
+        if nuda.env.verb: print("nucleus:",nucleus)
+        print("-> nucleus:",nucleus)
+        # 
+        self = SetupNeutronSkinTheo.init_self( self )
+        #
+        models, models_lower = nskin_theo_models( )
+        #
+        if model.lower() not in models_lower:
+            print('The model name ',model,' is not in the list of models.')
+            print('list of models:',models)
+            print('-- Exit the code --')
+            exit()
+        #
+        params, params_lower = nskin_theo_params( model = model )
+        #
+        if param.lower() not in params_lower:
+            print('The param set ',param,' is not in the list of param.')
+            print('list of param:',params)
+            print('-- Exit the code --')
+            exit()
+        #
+        # nucleus, nucleus_lower = nskin_theo_nucleus( )
+        # #
+        # if nucleus.lower() not in nucleus_lower:
+        #     print('The param set ',nucleus,' is not in the list of param.')
+        #     print('list of param:',nucleus)
+        #     print('-- Exit the code --')
+        #     exit()
+        # #            
+        if model.lower() == 'skyrme':
+            #
+             file_in1 = os.path.join(nuda.param.path_data,'NeutronSkin/skyrmeNskin-'+nucleus+'.dat')
+             if nuda.env.verb: print('Reads file1:',file_in1)
+             name = np.loadtxt( file_in1, usecols=(0), comments='#', unpack = True, dtype=str )
+             Rn, Rp, Rskin = np.loadtxt( file_in1, usecols=(1,2,3), comments='#', unpack = True )
+             #: Attribute providing the label the data is references for figures.
+             self.label = 'SKY-'+param
+             #: Attribute providing additional notes about the data.
+             self.note = "write here notes about this EOS."
+             #
+             if param in name:
+                 self.nskin = True
+                 ind = np.where(name == param )
+                 self.Rn = Rn[ind][0]; 
+                 self.Rp = Rp[ind][0]; 
+                 self.Rskin = Rskin[ind][0];
+             else:
+                 self.nskin = False
+            #
+        elif model.lower() == 'nlrh':
+            #
+             file_in1 = os.path.join(nuda.param.path_data,'NeutronSkin/nlrhNskin-'+nucleus+'.dat')
+             if nuda.env.verb: print('Reads file1:',file_in1)
+             name = np.loadtxt( file_in1, usecols=(0), comments='#', unpack = True, dtype=str )
+             Rn, Rp, Rskin = np.loadtxt( file_in1, usecols=(1,2,3), comments='#', unpack = True )
+             #: Attribute providing the label the data is references for figures.
+             self.label = 'NLRH-'+param
+             #: Attribute providing additional notes about the data.
+             self.note = "write here notes about this EOS."
+             #
+             if param in name:
+                 self.nskin = True
+                 ind = np.where(name == param )
+                 self.Rn = Rn[ind][0]; 
+                 self.Rp = Rp[ind][0]; 
+                 self.Rskin = Rskin[ind][0];
+             else:
+                 self.nskin = False
+            #
+        elif model.lower() == 'ddrh':
+            #
+             file_in1 = os.path.join(nuda.param.path_data,'NeutronSkin/ddrhNskin-'+nucleus+'.dat')
+             if nuda.env.verb: print('Reads file1:',file_in1)
+             name = np.loadtxt( file_in1, usecols=(0), comments='#', unpack = True, dtype=str )
+             Rn, Rp, Rskin = np.loadtxt( file_in1, usecols=(1,2,3), comments='#', unpack = True )
+             #: Attribute providing the label the data is references for figures.
+             self.label = 'DDRH-'+param
+             #: Attribute providing additional notes about the data.
+             self.note = "write here notes about this EOS."
+             #
+             if param in name:
+                 self.nskin = True
+                 ind = np.where(name == param )
+                 self.Rn = Rn[ind][0]; 
+                 self.Rp = Rp[ind][0]; 
+                 self.Rskin = Rskin[ind][0];
+             else:
+                 self.nskin = False
+            #
+        # 
+        self.Rn_unit = 'fm'
+        self.Rp_unit = 'fm'
+        self.Rskin_unit = 'fm'
+        #
+        if nuda.env.verb: print("Exit SetupNeutronSkinThe()")
+    #
+    def print_outputs( self ):
+        """
+        Method which print outputs on terminal's screen.
+        """
+        print("")
+        #
+        if nuda.env.verb: print("Enter print_outputs()")
+        #
+        print("- Print output:")
+        print("   nucleus:",self.nucleus)
+        print("   model:",self.model)
+        print("   param:",self.param)
+        #
+        print(f"   Rn: {np.round(self.Rn,4)} in {self.Rn_unit}")
+        print(f"   Rp: {np.round(self.Rn,4)} in {self.Rp_unit}")
+        print(f"   Rskin: {np.round(self.Rskin,4)} in {self.Rskin_unit}")
+        #
+        if nuda.env.verb: print("Exit print_outputs()")
+        #
+    def init_self( self ):
+        """
+        Initialize variables in self.
+        """
+        #
+        if nuda.env.verb: print("Enter init_self()")
+        #
+        #: Attribute the neutron matter density.
+        self.Rn = []
+        #: Attribute the symmetric matter density.
+        self.Rp = []
+        #: Attribute the neutron matter neutron Fermi momentum.
+        self.Rskin = []
+        #: Attribute the plot linestyle.
+        self.linestyle = None
+        #: Attribute the plot marker.
+        self.marker = 'o'
+        #: Attribute the plot every data.
+        self.every = 1
+        #
+        if nuda.env.verb: print("Exit init_self()")
+        #
+        return self        

nucleardatapy/nuc/setup_rch_exp.py

@@ -0,0 +1,159 @@
+import os
+import sys
+import numpy as np  # 1.15.0
+
+#nucleardatapy_tk = os.getenv('NUCLEARDATAPY_TK')
+#sys.path.insert(0, nucleardatapy_tk)
+
+import nucleardatapy as nuda
+
+def rch_exp_tables():
+    """
+    Return a list of the tables available in this toolkit for the charge radiuus and
+    print them all on the prompt.  These tables are the following
+    ones: '2013-Angeli'.
+
+    :return: The list of tables.
+    :rtype: list[str].
+    """
+    #
+    if nuda.env.verb: print("\nEnter rch_exp_tables()")
+    #
+    tables = [ '2013-Angeli' ]
+    #
+    print('tables available in the toolkit:',tables)
+    tables_lower = [ item.lower() for item in tables ]
+    print('tables available in the toolkit:',tables_lower)
+    #
+    if nuda.env.verb: print("Exit rch_exp_tables()")
+    #
+    return tables, tables_lower
+
+class setupRchExp():
+   """
+   Instantiate the object with charge radii choosen \
+   from a table.
+
+   This choice is defined in the variable `table`.
+
+   The tables can chosen among the following ones: \
+   '2013-Angeli'.
+
+   :param table: Fix the name of `table`. Default value: '2013-Angeli'.
+   :type table: str, optional. 
+
+   **Attributes:**
+   """
+   #
+   def __init__( self, table = '2013-Angeli' ):
+      """
+      Parameters
+      ----------
+      model : str, optional
+      The model to consider. Choose between: 1998-VAR-AM-APR (default), 2008-AFDMC-NM, ...
+      """
+      #
+      if nuda.env.verb: print("\nEnter setupRchExp()")
+      #
+      self.table = table
+      if nuda.env.verb: print("table:",table)
+      #
+      #: Attribute Z (charge of the nucleus).
+      self.nucZ = []
+      #: Attribute symb (symbol) of the element, e.g., Fe.
+      self.nucSymb = []
+      #: Attribute N (number of neutrons of the nucleus).
+      self.nucN = []
+      #: Attribute A (mass of the nucleus).
+      self.nucA = []
+      #: Attribue R_ch (charge radius) in fm.
+      self.nucRch = []
+      #: Attribue uncertainty in R_ch (charge radius) in fm.
+      self.nucRch_err = []
+      #
+      tables, tables_lower = rch_exp_tables()
+      #
+      if table.lower() not in tables_lower:
+         print('Table ',table,' is not in the list of tables.')
+         print('list of tables:',tables)
+         print('-- Exit the code --')
+         exit()
+      #
+      if table.lower() == '2013-angeli':
+         #
+         file_in = os.path.join(nuda.param.path_data,'nuclei/radch/2013-Angeli.csv')
+         if nuda.env.verb: print('Reads file:',file_in)
+         #: Attribute providing the full reference to the paper to be citted.
+         self.ref = 'I. Angeli and K.P. Marinova, Table of experimental nuclear ground state charge radii: An update, Atomic Data and Nuclear Data Tables 69, 69 (2013)'
+         #: Attribute providing the label the data is references for figures.
+         self.label = 'Angeli-Marinova-2013'
+         #: 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:
+               #print('line:',line)
+               if '#' in line:
+                  continue
+               linesplit = line.split(',')
+               #print('line.split:',linesplit)
+               self.nucZ.append(linesplit[0])
+               self.nucSymb.append(linesplit[1])
+               self.nucN.append(linesplit[2])
+               self.nucA.append(linesplit[3])
+               self.nucRch.append(linesplit[4])
+               self.nucRch_err.append(linesplit[5])
+               #
+         #: Attribute radius unit.
+         self.R_unit = 'fm'
+         #
+         if nuda.env.verb: print("Exit setupRchExp()")
+   #
+   def Rch_isotopes(self, Zref = 50 ):
+      """
+      This method provide a list if radii for an isotopic chain defined by Zref.
+
+      """
+      #
+      if nuda.env.verb: print("Enter Rch_isotopes()")
+      #
+      Nref = []
+      Aref = []
+      Rchref = []
+      Rchref_err = []
+      for k in range(len(self.nucZ)):
+         if int( self.nucZ[k] ) == Zref:
+            Nref.append( self.nucN[k] )
+            Aref.append( self.nucA[k] )
+            Rchref.append( self.nucRch[k] )
+            Rchref_err.append( self.nucRch_err[k] )
+      Nref = np.array( Nref, dtype = int )
+      Aref = np.array( Aref, dtype = int )
+      Rchref = np.array( Rchref, dtype = float )
+      Rchref_err = np.array( Rchref_err, dtype = float )
+      #
+      return Nref, Aref, Rchref, Rchref_err
+      #
+      if nuda.env.verb: print("Exit Rch_isotopes()")
+   #
+   def print_outputs( self ):
+      """
+      Method which print outputs on terminal's screen.
+      """
+      print("")
+      #
+      if nuda.env.verb: print("Enter print_outputs()")
+      #
+      print("- Print output:")
+      print("   table:",self.table)
+      print("   ref:",self.ref)
+      print("   label:",self.label)
+      print("   note:",self.note)
+      if any(self.nucZ): print(f"   Z: {self.nucZ}")
+      if any(self.nucA): print(f"   A: {self.nucA}")
+      if any(self.nucRch): print(f"   Rch: {self.nucRch}")
+      if any(self.nucRch_err): print(f"   Rch_err: {self.nucRch_err}")
+      #
+      if nuda.env.verb: print("Exit print_outputs()")
+      #
+

nucleardatapy/nuc/setup_rch_theo.py

@@ -0,0 +1,198 @@
+import os
+import sys
+import numpy as np  # 1.15.0
+
+#nucleardatapy_tk = os.getenv('NUCLEARDATAPY_TK')
+#sys.path.insert(0, nucleardatapy_tk)
+
+import nucleardatapy as nuda
+
+def rch_theo_tables():
+    """
+    Return a list of the tables available in this toolkit for the charge radiuus and
+    print them all on the prompt.  These tables are the following
+    ones: '2013-Angeli'.
+
+    :return: The list of tables.
+    :rtype: list[str].
+    """
+    #
+    if nuda.env.verb: print("\nEnter rch_theo_tables()")
+    #
+    tables = [ '2021-BSkG1', '2022-BSkG2', '2023-BSkG3', '2025-BSkG4'  ]
+    #
+    print('tables available in the toolkit:',tables)
+    tables_lower = [ item.lower() for item in tables ]
+    print('tables available in the toolkit:',tables_lower)
+    #
+    if nuda.env.verb: print("Exit rch_theo_tables()")
+    #
+    return tables, tables_lower
+
+class setupRchTheo():
+   """
+   Instantiate the object with charge radii choosen \
+   from a table.
+
+   This choice is defined in the variable `table`.
+
+   The tables can chosen among the following ones: \
+   '2013-Angeli'.
+
+   :param table: Fix the name of `table`. Default value: '2013-Angeli'.
+   :type table: str, optional. 
+
+   **Attributes:**
+   """
+   #
+   def __init__( self, table = '2021-BSkG1' ):
+      """
+      Parameters
+      ----------
+      table : str, optional
+      The theoretical table to consider. Choose between: 2021-BSkG1 (default), 2022-BSkG2, 2023-BSkG3, 2025-BSkG4 ...
+      """
+      #
+      if nuda.env.verb: print("\nEnter setupRchTheo()")
+      #
+      self.table = table
+      if nuda.env.verb: print("table:",table)
+      #
+      #: Attribute Z (charge of the nucleus).
+      self.nucZ = []
+      #: Attribute symb (symbol) of the element, e.g., Fe.
+      self.nucSymb = []
+      #: Attribute N (number of neutrons of the nucleus).
+      self.nucN = []
+      #: Attribute A (mass of the nucleus).
+      self.nucA = []
+      #: Attribue R_ch (charge radius) in fm.
+      self.nucRch = []
+      #
+      tables, tables_lower = rch_theo_tables()
+      #
+      if table.lower() not in tables_lower:
+         print('Table ',table,' is not in the list of tables.')
+         print('list of tables:',tables)
+         print('-- Exit the code --')
+         exit()
+      #
+      if table.lower() == '2021-bskg1':
+         #
+         file_in = os.path.join(nuda.param.path_data,'nuclei/masses/Theory/2021-BSkG1.txt')
+         if nuda.env.verb: print('Reads file:',file_in)
+         #: Attribute providing the full reference to the paper to be citted.
+         self.ref = 'G. Scamps, S. Goriely, E. Olsen, M. Bender, and W. Ryssens, EPJA 57, 333 (2021).'
+         #: Attribute providing additional notes about the data.
+         self.note = "write here notes about this EOS."
+         #: Attribute providing the label the data is references for figures.
+         self.label = 'BSkG1-2021'
+         #self.nucZr, self.nucNr, self.nucBE2A  = np.loadtxt( file_in, usecols=(0,1,2), unpack = True )
+         self.nucZr, self.nucNr, self.nucMass, self.Ebind, self.beta20, self.beta22, self.beta2, self.Erot, \
+         self.gap_n, self.gap_p, self.nucRch, self.moi = \
+         np.loadtxt( file_in, usecols=(0,1,3,5,6,7,8,9,10,11,12,15), unpack = True )
+         self.nucZ = np.array( [ int(ele) for ele in self.nucZr ] )
+         self.nucN = np.array( [ int(ele) for ele in self.nucNr ] )
+         self.nucA = self.nucZ + self.nucN
+         #
+      elif table.lower() == '2022-bskg2':
+         #
+         file_in = os.path.join(nuda.param.path_data,'nuclei/masses/Theory/2022-BSkG2.txt')
+         if nuda.env.verb: print('Reads file:',file_in)
+         #: Attribute providing the full reference to the paper to be citted.
+         self.ref = 'W. Ryssens, G. Scamps, S. Goriely, and M. Bender, EPJA 58, 246 (2022).'
+         #: Attribute providing additional notes about the data.
+         self.note = "write here notes about this EOS."
+         #: Attribute providing the label the data is references for figures.
+         self.label = 'BSkG2-2022'
+         #self.nucZr, self.nucNr, self.nucBE2A  = np.loadtxt( file_in, usecols=(0,1,2), unpack = True )
+         self.nucZr, self.nucNr, self.nucMass, self.Ebind, self.beta20, self.beta22, self.beta2, self.Erot, \
+         self.gap_n, self.gap_p, self.nucRch, self.moi = \
+         np.loadtxt( file_in, usecols=(0,1,3,5,6,7,8,9,10,11,12,15), unpack = True )
+         self.nucZ = np.array( [ int(ele) for ele in self.nucZr ] )
+         self.nucN = np.array( [ int(ele) for ele in self.nucNr ] )
+         self.nucA = self.nucZ + self.nucN
+         #
+      elif table.lower() == '2023-bskg3':
+         #
+         file_in = os.path.join(nuda.param.path_data,'nuclei/masses/Theory/2023-BSkG3.txt')
+         if nuda.env.verb: print('Reads file:',file_in)
+         #: Attribute providing the full reference to the paper to be citted.
+         self.ref = 'G. Grams, W. Ryssens, G. Scamps, S. Goriely, and N. Chamel, EPJA 59, 270 (2023).'
+         #: Attribute providing additional notes about the data.
+         self.note = "write here notes about this EOS."
+         #: Attribute providing the label the data is references for figures.
+         self.label = 'BSkG3-2023'
+         self.nucZr, self.nucNr, self.nucMass, self.Ebind, self.beta20, self.beta22, self.beta2, \
+         self.beta30, self.beta32, self.Erot, self.gap_n, self.gap_p, self.nucRch, self.moi = \
+         np.loadtxt( file_in, usecols=(0,1,3,5,6,7,8,9,10,11,12,13,14,17), unpack = True )
+         self.nucZ = np.array( [ int(ele) for ele in self.nucZr ] )
+         self.nucN = np.array( [ int(ele) for ele in self.nucNr ] )
+         self.nucA = self.nucZ + self.nucN
+         #
+      elif table.lower() == '2025-bskg4':
+         #
+         file_in = os.path.join(nuda.param.path_data,'nuclei/masses/Theory/2025-BSkG4.txt')
+         if nuda.env.verb: print('Reads file:',file_in)
+         #: Attribute providing the full reference to the paper to be citted.
+         self.ref = 'G. Grams, W. Ryssens, N. Shchechilin, A. Sanchez-Fernandez, N. Chamel, and S. Goriely,  EPJA 61, 35 (2025).'
+         #: Attribute providing additional notes about the data.
+         self.note = "write here notes about this EOS."
+         #: Attribute providing the label the data is references for figures.
+         self.label = 'BSkG4-2025'
+         self.nucZr, self.nucNr, self.nucMass, self.Ebind, self.beta20, self.beta22, self.beta2, \
+         self.beta30, self.beta32, self.Erot, self.gap_n, self.gap_p, self.nucRch, self.moi = \
+         np.loadtxt( file_in, usecols=(0,1,3,5,6,7,8,9,10,11,12,13,14,17), unpack = True )
+         self.nucZ = np.array( [ int(ele) for ele in self.nucZr ] )
+         self.nucN = np.array( [ int(ele) for ele in self.nucNr ] )
+         self.nucA = self.nucZ + self.nucN
+         #
+      #: Attribute radius unit.
+      self.R_unit = 'fm'
+      #
+      if nuda.env.verb: print("Exit setupChTheo()")
+   #
+   def Rch_isotopes(self, Zref = 50 ):
+      """
+      This method provide a list if radii for an isotopic chain defined by Zref.
+
+      """
+      #
+      if nuda.env.verb: print("Enter Rch_isotopes()")
+      #
+      Nref = []
+      Aref = []
+      Rchref = []
+      for k in range(len(self.nucZ)):
+         if int( self.nucZ[k] ) == Zref:
+            Nref.append( self.nucN[k] )
+            Aref.append( self.nucA[k] )
+            Rchref.append( self.nucRch[k] )
+      Nref = np.array( Nref, dtype = int )
+      Aref = np.array( Aref, dtype = int )
+      Rchref = np.array( Rchref, dtype = float )
+      #
+      return Nref, Aref, Rchref
+      #
+      if nuda.env.verb: print("Exit Rch_isotopes()")
+   #
+   def print_outputs( self ):
+      """
+      Method which print outputs on terminal's screen.
+      """
+      print("")
+      #
+      if nuda.env.verb: print("Enter print_outputs()")
+      #
+      print("- Print output:")
+      print("   table:",self.table)
+      print("   ref:",self.ref)
+      print("   label:",self.label)
+      print("   note:",self.note)
+      if any(self.nucZ): print(f"   Z: {self.nucZ}")
+      if any(self.nucA): print(f"   A: {self.nucA}")
+      if any(self.nucRch): print(f"   Rch: {self.nucRch}")
+      #
+      if nuda.env.verb: print("Exit print_outputs()")
+      #
+