nucleardatapy 0.2.0__py3-none-any.whl

nucleardatapy/nuc/setup_isgmr_exp.py

@@ -0,0 +1,494 @@
+import os
+import sys
+import numpy as np  # 1.15.0
+import math
+
+#nucleardatapy_tk = os.getenv('NUCLEARDATAPY_TK')
+#sys.path.insert(0, nucleardatapy_tk)
+
+import nucleardatapy as nuda
+
+def isgmr_exp_tables():
+    """
+    Return a list of tables available in this toolkit for the ISGMR energy and
+    print them all on the prompt. These tables are the following
+    ones: '2010-ISGMR-LI', '2018-ISGMR-GARG', '2018-ISGMR-GARG-LATEX'.
+
+    :return: The list of tables.
+    :rtype: list[str].    
+    """
+    #
+    if nuda.env.verb: print("\nEnter isgmr_exp_tables()")
+    #
+    tables = [ '2010-ISGMR-LI', '2018-ISGMR-GARG', '2018-ISGMR-GARG-LATEX', '2022-ISGMR-average' ]
+    #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 isgmr_exp_tables()")
+    #
+    return tables, tables_lower
+
+class setupISGMRExp():
+   """
+   Instantiate the object with microscopic results choosen \
+   by the toolkit practitioner. \
+
+   This choice is defined in the variable `table`.
+
+   The `table` can chosen among the following ones: \
+   '2010-ISGMR-LI', '2018-ISGMR-GARG'.
+
+   :param table: Fix the name of `table`. Default value: '2018-ISGMR-GARG', '2018-ISGMR-GARG-LATEX'.
+   :type table: str, optional. 
+
+   **Attributes:**
+   """
+   #
+   def __init__( self, table = '2018-ISGMR-GARG' ):
+      """
+      Parameters
+      ----------
+      table : str, optional
+      The table to consider. \
+      Choose between: 2018-ISGMR-GARG (default) and 2010-ISGMR-LI.
+      """
+      #
+      if nuda.env.verb: print("\nEnter setupISGMRExp()")
+      #: Attribute table.
+      self.table = table
+      if nuda.env.verb: print("table:",table)
+      #
+      tables, tables_lower = isgmr_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()
+      #
+      nucA=[]; nucZ=[]; nucN=[]; nucSymbol=[]; nucEprobe=[]; nucProj=[]; 
+      nucE0=[]; nucE0_errp=[]; nucE0_errm=[];
+      nucG=[]; nucG_errp=[]; nucG_errm=[]; 
+      nucEWSR=[]; nucEWSR_errp=[]; nucEWSR_errm=[]; 
+      nucM12M0=[]; nucM12M0_errp=[]; nucM12M0_errm=[]; 
+      nucM12Mm1=[]; nucM12Mm1_errp=[]; nucM12Mm1_errm=[]; 
+      nucM32M1=[]; nucM32M1_errp=[]; nucM32M1_errm=[];
+      nucRef=[];
+      #
+      if table.lower() == '2010-isgmr-li':
+         #
+         file_in = os.path.join(nuda.param.path_data,'nuclei/isgmr/2010-ISGMR-Li.dat')
+         if nuda.env.verb: print('Reads file:',file_in)
+         #: Attribute providing the full reference to the paper to be citted.
+         self.ref = 'T. Li, U. Garg, Y. Liu et al., Phys. Rev. C 81, 034309 (2010)'
+         #: Attribute providing the label the data is references for figures.
+         self.label = 'Li-Garg-Liu-2010'
+         #: Attribute providing additional notes about the data.
+         self.note = "write here notes about this table."
+         nucZ, nucA, nucM12Mm1, nucM12Mm1_errp, nucM12Mm1_errm = \
+            np.loadtxt( file_in, usecols=(0,1,4,5,6), comments='#', unpack = True )
+         nucN = nucA - nucZ
+         for k,Z in enumerate(nucZ):
+            nucSymbol.append( nuda.param.elements[int(Z)-1] )
+            nucEprobe.append( '386' )
+            nucProj.append( '$\alpha$' )
+            nucE0.append( None ); nucE0_errp.append( None ); nucE0_errm.append( None )
+            nucG.append( None ); nucG_errp.append( None ); nucG_errm.append( None )
+            nucEWSR.append( None ); nucEWSR_errp.append( None ); nucEWSR_errm.append( None )
+            nucM12M0.append( None ); nucM12M0_errp.append( None ); nucM12M0_errm.append( None )
+            nucM32M1.append( None ); nucM32M1_errp.append( None ); nucM32M1_errm.append( None )
+         nuc = len( nucZ ); nbk = nuc
+         #
+      elif table.lower() == '2018-isgmr-garg':
+         #
+         file_in = os.path.join(nuda.param.path_data,'nuclei/isgmr/2018-ISGMR-Garg.dat')
+         if nuda.env.verb: print('Reads file:',file_in)
+         self.ref = 'U. Garg and G. Colo, Prog. Part. Nucl. Phys. 101, 55 (2018)'
+         self.label = 'Garg-Colo-2018'
+         self.note = "write here notes about this table."
+         nucZ, nucA, nucM12Mm1, nucM12Mm1_errp, nucM12Mm1_errm = \
+            np.loadtxt( file_in, usecols=(0,1,2,3,4), comments='#', unpack = True )
+         nucN = nucA - nucZ
+         #print('elements:',nuda.param.elements)
+         for k,Z in enumerate(nucZ):
+            nucSymbol.append( nuda.param.elements[int(Z)-1] )
+            nucEprobe.append( '100' )
+            nucProj.append( '$\alpha$' )
+            nucE0.append( None ); nucE0_errp.append( None ); nucE0_errm.append( None )
+            nucG.append( None ); nucG_errp.append( None ); nucG_errm.append( None )
+            nucEWSR.append( None ); nucEWSR_errp.append( None ); nucEWSR_errm.append( None )
+            nucM12M0.append( None ); nucM12M0_errp.append( None ); nucM12M0_errm.append( None )
+            nucM32M1.append( None ); nucM32M1_errp.append( None ); nucM32M1_errm.append( None )
+         nuc = len( nucZ ); nbk = nuc
+         #
+      elif table.lower() == '2018-isgmr-garg-latex':
+         #
+         file_in = os.path.join(nuda.param.path_data,'nuclei/isgmr/2018-ISGMR-Garg.tex')
+         if nuda.env.verb: print('Reads file:',file_in)
+         self.ref = 'U. Garg and G. Colo, Prog. Part. Nucl. Phys. 101, 55 (2018)'
+         self.label = 'Garg-Colo-2018'
+         self.note = "Parameters of the ISGMR peaks and moment ratios of the ISGMR strength distributions in stable nuclei as reported by the TAMU and RCNP groups. The probes employed in the measurements are listed for each case. Entries marked with $\\star$ indicate that the $\\Gamma$ is an RMS width, not that of a fitted peak. Entries marked with $\\dagger$ indicate a multimodal strength distribution; in those cases the parameters for only the ``main'' ISGMR peak are included. For the TAMU data, the peak parameters correspond to a Gaussian fit, whereas for the RCNP  data, the corresponding parameters are for a Lorentzian fit."
+         #
+         nbk = 0
+         nuc = -1
+         with open(file_in,'r') as file:
+            for line in file:
+               #print('line:',line)
+               if '#' in line[0]: 
+                  continue
+               ele = line.split('&')
+               #print('ele:',ele)
+               # ele[0]: nucleus
+               if ele[0] == '  ' or ele[0] == ' ':
+                  nucSymbol.append( nucSymbol[-1] )
+                  nucA.append( nucA[-1] )
+                  nucZ.append( nucZ[-1] )
+                  nucN.append( nucN[-1] )
+               else:
+                  nuc += 1
+                  symbol = ele[0].split('$')[2].strip()
+                  ZZ, = np.where( nuda.param.elements == symbol )[0] + 1
+                  AA = int( ele[0].split('$')[1].strip('^').strip('{').strip('}') )
+                  NN = AA - ZZ
+                  #ZZ += 1
+                  #print('Z=',ZZ,' symbol:',symbol,' A=',AA,' N=',NN)
+                  nucSymbol.append( symbol )
+                  nucA.append( str( AA ) )
+                  nucZ.append( str( ZZ ) )
+                  nucN.append( str( NN ) )
+                  #print('A=',AA)
+               #print('Z=',nucZ[-1],' symbol:',nucSymbol[-1],' A=',nucA[-1],' N=',nucN[-1])
+               # ele[1]: probe
+               #print('ele[1]:',ele[1])
+               if ele[1] == ' ' or ele[1] == '  ':
+                  nucEprobe.append( nucEprobe[-1] )
+                  nucProj.append( nucProj[-1] )
+               else:
+                  Eprobe = int( ele[1].split('MeV-')[0].strip() )
+                  proj = ele[1].split('MeV-')[1].strip()
+                  nucEprobe.append( Eprobe )
+                  nucProj.append( proj )
+               #print('Z=',nucZ[-1],' symbol:',nucSymbol[-1],' A=',nucA[-1],' N=',nucN[-1],' Eprobe=',nucEprobe[-1],' proj:',nucProj[-1])
+               # ele[3]: E0
+               cent, errp, errm = nuda.param.tex2str( ele[3] )
+               nucE0.append( cent ); nucE0_errp.append( errp ); nucE0_errm.append( errm );
+               #print('Z=',nucZ[-1],' symbol:',nucSymbol[-1],' A=',nucA[-1],' N=',nucN[-1],' Eprobe=',nucEprobe[-1],' proj:',nucProj[-1],' E0:',nucE0[-1],nucE0_errp[-1],nucE0_errm[-1])
+               # ele[3]: Gamma
+               cent, errp, errm = nuda.param.tex2str( ele[4] )
+               nucG.append( cent ); nucG_errp.append( errp ); nucG_errm.append( errm );
+               # ele[4]: EWSR
+               cent, errp, errm = nuda.param.tex2str( ele[5] )
+               nucEWSR.append( cent ); nucEWSR_errp.append( errp ); nucEWSR_errm.append( errm );
+               # ele[5]: nada
+               # ele[6]: M12M0
+               cent, errp, errm = nuda.param.tex2str( ele[7] )
+               nucM12M0.append( cent ); nucM12M0_errp.append( errp ); nucM12M0_errm.append( errm );
+               # ele[7]: M12Mm1
+               cent, errp, errm = nuda.param.tex2str( ele[8] )
+               nucM12Mm1.append( cent ); nucM12Mm1_errp.append( errp ); nucM12Mm1_errm.append( errm );
+               print('nbk:',nbk,' nuc:',nuc,' Z=',nucZ[-1],' symbol:',nucSymbol[-1],' A=',nucA[-1],' N=',nucN[-1],' Eprobe=',nucEprobe[-1],' proj:',nucProj[-1],' E0:',nucE0[-1],nucE0_errp[-1],nucE0_errm[-1])
+               # ele[8]: M32M1
+               cent, errp, errm = nuda.param.tex2str( ele[9] )
+               nucM32M1.append( cent ); nucM32M1_errp.append( errp ); nucM32M1_errm.append( errm );
+               # ele[9]: nada
+               # ele[10]: ref
+               nucRef.append( ele[11] )
+               #print('nuc:',nuc,nucA,nucSymbol,nucProbe,nucTarget,nucG,nucEWSR,nucM12M0,nucM12Mm1,nucM32M1,nucRef)
+               #exit()
+               nbk += 1
+         #
+         nbk -= 1
+         #
+      elif table.lower() == '2022-isgmr-average':
+         #
+         file_in = os.path.join(nuda.param.path_data,'nuclei/isgmr/2022-ISGMR-average.dat')
+         if nuda.env.verb: print('Reads file:',file_in)
+         self.ref = 'U. Garg and G. Colo, Prog. Part. Nucl. Phys. 101, 55 (2018)'
+         self.label = 'Average-2022'
+         self.note = "write here notes about this table."
+         nucZ, nucA, nucM12Mm1, nucM12Mm1_errp, nucM12Mm1_errm = \
+            np.loadtxt( file_in, usecols=(0,1,2,3,4), comments='#', unpack = True )
+         nucN = nucA - nucZ
+         #print('elements:',nuda.param.elements)
+         for k,Z in enumerate(nucZ):
+            nucSymbol.append( nuda.param.elements[int(Z)-1] )
+            nucEprobe.append( '100' )
+            nucProj.append( '$\alpha$' )
+            nucE0.append( None ); nucE0_errp.append( None ); nucE0_errm.append( None )
+            nucG.append( None ); nucG_errp.append( None ); nucG_errm.append( None )
+            nucEWSR.append( None ); nucEWSR_errp.append( None ); nucEWSR_errm.append( None )
+            nucM12M0.append( None ); nucM12M0_errp.append( None ); nucM12M0_errm.append( None )
+            nucM32M1.append( None ); nucM32M1_errp.append( None ); nucM32M1_errm.append( None )
+         nuc = len( nucZ ); nbk = nuc
+      #
+      print('\nnumber of different nuclei:',nuc)
+      print('\nnumber of total entries:   ',nbk)
+      #
+      isgmr = {}
+      isgmr['A'] = nucA; isgmr['Z'] = nucZ; isgmr['N'] = nucN; isgmr['symbol'] = nucSymbol
+      isgmr['Eprobe'] = nucEprobe; isgmr['proj'] = nucProj
+      isgmr['E0'] = nucE0; isgmr['E0_errp'] = nucE0_errp; isgmr['E0_errm'] = nucE0_errm
+      isgmr['G'] = nucG; isgmr['G_errp'] = nucG_errp; isgmr['G_errm'] = nucG_errm
+      isgmr['EWSR'] = nucEWSR; isgmr['EWSR_errp'] = nucEWSR_errp; isgmr['EWSR_errm'] = nucEWSR_errm
+      isgmr['M12M0'] = nucM12M0; isgmr['M12M0_errp'] = nucM12M0_errp; isgmr['M12M0_errm'] = nucM12M0_errm
+      isgmr['M12Mm1'] = nucM12Mm1; isgmr['M12Mm1_errp'] = nucM12Mm1_errp; isgmr['M12Mm1_errm'] = nucM12Mm1_errm
+      isgmr['M32M1'] = nucM32M1; isgmr['M32M1_errp'] = nucM32M1_errp; isgmr['M32M1_errm'] = nucM32M1_errm
+      isgmr['ref'] = nucRef
+      self.isgmr = isgmr
+      #
+      #: Attribute energy unit.
+      self.E_unit = 'MeV'
+      #
+      if nuda.env.verb: print("Exit setupISGMRExp()")
+   #
+   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)
+      print('\nZ:',self.isgmr['Z'])
+      print('\nA:',self.isgmr['A'])
+      for ind,Z in enumerate( self.isgmr['Z'] ):
+         print('For Z:',Z,' A:',self.isgmr['A'][ind])
+         for A in self.isgmr['A'][ind]:
+            print('Centroid energy:',self.isgmr['M12Mm1'][ind])
+            print('   with errp:',self.isgmr['M12Mm1_errp'][ind])
+            print('   with errm:',self.isgmr['M12Mm1_errm'][ind])
+      #
+      if nuda.env.verb: print("Exit print_outputs()")
+      #
+   #
+   def average( self ):
+      """
+      Method to average the data when same target is given.
+
+      **Attributes:**
+      """
+      print("")
+      #
+      if nuda.env.verb: print("Enter average()")
+      #
+      k = 0
+      AAm1 = 0
+      ZZm1 = 0
+      #
+      nA=[]; nZ=[]; nN=[]; nSymbol=[];  
+      nE0=[]; nE0_errp=[]; nE0_errm=[];
+      nG=[]; nG_errp=[]; nG_errm=[]; 
+      nEWSR=[]; nEWSR_errp=[]; nEWSR_errm=[]; 
+      nM12M0=[]; nM12M0_errp=[]; nM12M0_errm=[]; 
+      nM12Mm1=[]; nM12Mm1_errp=[]; nM12Mm1_errm=[]; 
+      nM32M1=[]; nM32M1_errp=[]; nM32M1_errm=[];
+      #
+      while k < nbk:
+         AA   = nucA[k]
+         ZZ   = nucZ[k]
+         if k>0: 
+            AAm1 = nucA[k-1]
+            ZZm1 = nucZ[k-1]
+         if k < nbk-1:
+            AAp1 = nucA[k+1]
+            ZZp1 = nucZ[k+1]
+         else:
+            AAp1 = 0
+            ZZp1 = 0
+         #
+         if AA != AAm1 or ZZ != ZZm1:
+            #
+            # Initialisation
+            #
+            nbE0 = 0
+            if nucE0[k] is not None:
+               nbE0 += 1
+               E0m = float(nucE0[k])
+               E0m_errp = float(nucE0_errp[k])**2
+               E0m_errm = float(nucE0_errm[k])**2
+            nbG = 0
+            if nucG[k] is not None:
+               nbG += 1
+               Gm = float(nucG[k])
+               Gm_errp = float(nucG_errp[k])**2
+               Gm_errm = float(nucG_errm[k])**2
+            nbEWSR = 0
+            if nucEWSR[k] is not None:
+               nbEWSR += 1
+               EWSRm = float(nucEWSR[k])
+               EWSRm_errp = float(nucEWSR_errp[k])**2
+               EWSRm_errm = float(nucEWSR_errm[k])**2
+            nbM12M0 = 0
+            if nucM12M0[k] is not None:
+               nbM12M0 += 1
+               M12M0m = float(nucM12M0[k])
+               M12M0m_errp = float(nucM12M0_errp[k])**2
+               M12M0m_errm = float(nucM12M0_errm[k])**2
+            nbM12Mm1 = 0
+            if nucM12Mm1[k] is not None:
+               nbM12Mm1 += 1
+               M12Mm1m = float(nucM12Mm1[k])
+               M12Mm1m_errp = float(nucM12Mm1_errp[k])**2
+               M12Mm1m_errm = float(nucM12Mm1_errm[k])**2
+            nbM32M1 = 0
+            if nucM32M1[k] is not None:
+               nbM32M1 += 1
+               M32M1m = float(nucM32M1[k])
+               M32M1m_errp = float(nucM32M1_errp[k])**2
+               M32M1m_errm = float(nucM32M1_errm[k])**2
+            #
+         if AA == AAp1 and ZZ == ZZp1:
+            #
+            if nucE0[k+1] is not None:
+               nbE0 += 1
+               E0m += float(nucE0[k+1])
+               E0m_errp += float(nucE0_errp[k+1])**2
+               E0m_errm += float(nucE0_errm[k+1])**2
+            if nucG[k+1] is not None:
+               nbG += 1
+               Gm += float(nucG[k+1])
+               Gm_errp += float(nucG_errp[k+1])**2
+               Gm_errm += float(nucG_errm[k+1])**2
+            if nucEWSR[k+1] is not None:
+               nbEWSR += 1
+               EWSRm += float(nucEWSR[k+1])
+               EWSRm_errp += float(nucEWSR_errp[k+1])**2
+               EWSRm_errm += float(nucEWSR_errm[k+1])**2
+            if nucM12M0[k+1] is not None:
+               nbM12M0 += 1
+               M12M0m += float(nucM12M0[k+1])
+               M12M0m_errp += float(nucM12M0_errp[k+1])**2
+               M12M0m_errm += float(nucM12M0_errm[k+1])**2
+            if nucM12Mm1[k+1] is not None:
+               nbM12Mm1 += 1
+               M12Mm1m += float(nucM12Mm1[k+1])
+               M12Mm1m_errp += float(nucM12Mm1_errp[k+1])**2
+               M12Mm1m_errm += float(nucM12Mm1_errm[k+1])**2
+            if nucM32M1[k+1] is not None:
+               nbM32M1 += 1
+               M32M1m += float(nucM32M1[k+1])
+               M32M1m_errp += float(nucM32M1_errp[k+1])**2
+               M32M1m_errm += float(nucM32M1_errm[k+1])**2
+            #
+         else:
+            #
+            nA.append( nucA[k] )
+            nZ.append( nucZ[k] )
+            nN.append( nucN[k] )
+            nSymbol.append( nucSymbol[k] )
+            if nbE0 == 0:
+               nE0.append( None )
+               nE0_errp.append( None )
+               nE0_errm.append( None )
+            else:
+               nE0.append( E0m / nbE0 )
+               nE0_errp.append( math.sqrt( E0m_errp / nbE0 ) )
+               nE0_errm.append( math.sqrt( E0m_errm / nbE0 ) )
+            if nbG == 0:
+               nG.append( None )
+               nG_errp.append( None )
+               nG_errm.append( None )
+            else:
+               nG.append( Gm / nbG )
+               nG_errp.append( math.sqrt( Gm_errp / nbG ) )
+               nG_errm.append( math.sqrt( Gm_errm / nbG ) )
+            if nbEWSR == 0:
+               nEWSR.append( None )
+               nEWSR_errp.append( None )
+               nEWSR_errm.append( None )
+            else:
+               nEWSR.append( EWSRm / nbEWSR )
+               nEWSR_errp.append( math.sqrt( EWSRm_errp / nbEWSR ) )
+               nEWSR_errm.append( math.sqrt( EWSRm_errm / nbEWSR ) )
+            if nbM12M0 == 0:
+               nM12M0.append( None )
+               nM12M0_errp.append( None )
+               nM12M0_errm.append( None )
+            else:
+               nM12M0.append( M12M0m / nbM12M0 )
+               nM12M0_errp.append( math.sqrt( M12M0m_errp / nbM12M0 ) )
+               nM12M0_errm.append( math.sqrt( M12M0m_errm / nbM12M0 ) )
+            if nbM12Mm1 == 0:
+               nM12Mm1.append( None )
+               nM12Mm1_errp.append( None )
+               nM12Mm1_errm.append( None )
+            else:
+               nM12Mm1.append( M12Mm1m / nbM12Mm1 )
+               nM12Mm1_errp.append( math.sqrt( M12Mm1m_errp / nbM12Mm1 ) )
+               nM12Mm1_errm.append( math.sqrt( M12Mm1m_errm / nbM12Mm1 ) )
+            if nbM32M1 == 0:
+               nM32M1.append( None )
+               nM32M1_errp.append( None )
+               nM32M1_errm.append( None )
+            else:
+               nM32M1.append( M32M1m / nbM32M1 )
+               nM32M1_errp.append( math.sqrt( M32M1m_errp / nbM32M1 ) )
+               nM32M1_errm.append( math.sqrt( M32M1m_errm / nbM32M1 ) )
+         k += 1
+      print('End k:',k)
+      isgmrm = {}
+      isgmrm['A'] = nA; isgmrm['Z'] = nZ; isgmrm['N'] = nN; isgmrm['symbol'] = nSymbol
+      isgmrm['E0'] = nE0; isgmrm['E0_errp'] = nE0_errp; isgmrm['E0_errm'] = nE0_errm
+      isgmrm['G'] = nG; isgmrm['G_errp'] = nG_errp; isgmrm['G_errm'] = nG_errm
+      isgmrm['EWSR'] = nEWSR; isgmrm['EWSR_errp'] = nEWSR_errp; isgmrm['EWSR_errm'] = nEWSR_errm
+      isgmrm['M12M0'] = nM12M0; isgmrm['M12M0_errp'] = nM12M0_errp; isgmrm['M12M0_errm'] = nM12M0_errm
+      isgmrm['M12Mm1'] = nM12Mm1; isgmrm['M12Mm1_errp'] = nM12Mm1_errp; isgmrm['M12Mm1_errm'] = nM12Mm1_errm
+      isgmrm['M32M1'] = nM32M1; isgmrm['M32M1_errp'] = nM32M1_errp; isgmrm['M32M1_errm'] = nM32M1_errm
+      self.isgmrm = isgmrm
+      #
+      for k in range(len(isgmrm['A'])):
+         print('Z=',isgmrm['Z'][k],' symbol:',isgmrm['symbol'][k],' A=',isgmrm['A'][k],' N=',isgmrm['N'][k],' E0:',isgmrm['E0'][k],isgmrm['E0_errp'][k],isgmrm['E0_errm'][k])
+      #
+      return self
+      #
+      if nuda.env.verb: print("Exit average()")
+      #
+   #
+   def select( self, Zref=50, obs = 'M12Mm1' ):
+      """
+      Method to select a subset of data.
+
+      :param Zref: Fix the reference charge for the search of isotopes.
+      :type Zref: int, optional. Default: 1.
+      :param obs: kind of observable to extract: 'M12M0', 'M12Mm1', 'M32M1'.
+      :type obs: str
+      **Attributes:**
+      """
+      print("")
+      #
+      if nuda.env.verb: print("Enter select()")
+      #
+      nucA = []; cent = []; errp = []; errm = [];
+      for ind,A in enumerate(self.isgmr['A']):
+         if obs == 'M12M0' and int( self.isgmr['Z'][ind] ) == Zref and self.isgmr['M12M0'][ind] is not None:
+            nucA.append( int(A) )
+            cent.append( float( self.isgmr['M12M0'][ind] ) )
+            errp.append( float( self.isgmr['M12M0_errp'][ind] ) )
+            errm.append( float( self.isgmr['M12M0_errm'][ind] ) )
+         if obs == 'M12Mm1' and int( self.isgmr['Z'][ind] ) == Zref and self.isgmr['M12Mm1'][ind] is not None:
+            nucA.append( int(A) )
+            cent.append( float( self.isgmr['M12Mm1'][ind] ) )
+            errp.append( float( self.isgmr['M12Mm1_errp'][ind] ) )
+            errm.append( float( self.isgmr['M12Mm1_errm'][ind] ) )
+         if obs == 'M32M1' and int( self.isgmr['Z'][ind] ) == Zref and self.isgmr['M32M1'][ind] is not None:
+            nucA.append( int(A) )
+            cent.append( float( self.isgmr['M32M1'][ind] ) )
+            errp.append( float( self.isgmr['M32M1_errp'][ind] ) )
+            errm.append( float( self.isgmr['M32M1_errm'][ind] ) )
+      erra = 0.5 * np.add( errp, errm )
+      self.nucA = nucA
+      self.cent = cent
+      self.errp = errp
+      self.errm = errm
+      self.erra = erra
+      #
+      return self
+      #
+      if nuda.env.verb: print("Exit select()")
+      #