nucleardatapy 0.2.0__py3-none-any.whl → 1.0.0__py3-none-any.whl

nucleardatapy/matter/setup_ffg.py

@@ -104,10 +104,30 @@ def esymffg_nr( kf ):
     """
     return effg_nr( kf ) * ( nuda.cst.two**nuda.cst.twoThird - 1.0 )
 
+def pre_nr( kf ):
+    """
+    Free Fermi gas pressure as a function of the Fermi momentum.
+
+    :param kf: Fermi momentum.
+    :type kf: float or numpy vector of real numbers.
+
+    """
+    return nuda.cst.twoThird * effg_nr( kf ) * den( kf )
+
+def cs2_nr( kf ):
+    """
+    Free Fermi gas sound speed as a function of the Fermi momentum.
+
+    :param kf: Fermi momentum.
+    :type kf: float or numpy vector of real numbers.
+
+    """
+    return nuda.cst.twoThird
+
 # FFG energy
 def feden(gam, kf, mc2):
     den = gam * kf**3 / ( 6 * nuda.cst.pi2 )
-    e2v = []
+    eps = []
     e2a = []
     for ind, val_kf in enumerate(kf):
         if val_kf > 1e-12:
@@ -120,14 +140,14 @@ def feden(gam, kf, mc2):
                 r = ( pf + ef ) / mc2
         #term = 2.0 * pf * ef**3 - pf * ef * mc2**2 - mc2**4 * np.log(r) - 8.0 / 3.0 * pf**3 * mc2
                 term = 2.0 * pf * ef**3 - pf * ef * mc2**2 - mc2**4 * np.log(r)
-            e2v.append( gam * term / (16.0 * nuda.cst.pi2 * nuda.cst.hbc**3 ) )
-            e2a.append( e2v[-1] / den[ind] )
+            eps.append( gam * term / (16.0 * nuda.cst.pi2 * nuda.cst.hbc**3 ) )
+            e2a.append( eps[-1] / den[ind] )
         else:
-            e2v.append( 0.0 )
+            eps.append( 0.0 )
             e2a.append( 0.0 )
-    e2v = np.array( e2v, dtype = float )
+    eps = np.array( eps, dtype = float )
     e2a = np.array( e2a, dtype = float )
-    return e2v, e2a
+    return eps, e2a
 
 # FFG pressure
 def fpres(gam, kf, mc2):
@@ -210,7 +230,7 @@ class setupFFGNuc():
         #: Attribute proton density
         self.den_p = self.x_p * den
         #: Attribute Fermi momentum for a Fermi system with degeneracy = 4
-        self.kf = (1.5 * nuda.cst.pi2 * self.den)**nuda.cst.third
+        self.kf_nuc = (1.5 * nuda.cst.pi2 * self.den)**nuda.cst.third
         #: Attribute neutron Fermi momentum (degeneracy = 2)
         self.kf_n = (nuda.cst.three * nuda.cst.pi2 * self.den_n)**nuda.cst.third
         #: Attribute proton Fermi momentum (degeneracy = 2)
@@ -225,11 +245,11 @@ class setupFFGNuc():
         self.eF_p_int_nr = nuda.cst.half * nuda.cst.h2m / ms * self.kf_p**nuda.cst.two
         #: Attribute rest mass energy per particle (degeneracy = 2)
         self.e2a_rm = self.x_n * ms * nuda.cst.mnc2 + self.x_p * ms * nuda.cst.mpc2
-        self.e2v_rm = self.e2a_rm * self.den
+        self.eps_rm = self.e2a_rm * self.den
         #: Attribute FFG energy per particle (degeneracy = 2)
-        e2v_n, e2a_n = feden( 2.0, self.kf_n, ms*nuda.cst.mnc2)
-        e2v_p, e2a_p = feden( 2.0, self.kf_p, ms*nuda.cst.mpc2)
-        self.e2v = e2v_n + e2v_p
+        eps_n, e2a_n = feden( 2.0, self.kf_n, ms*nuda.cst.mnc2)
+        eps_p, e2a_p = feden( 2.0, self.kf_p, ms*nuda.cst.mpc2)
+        self.eps = eps_n + eps_p
         self.e2a = self.x_n * e2a_n + self.x_p * e2a_p
         self.e2a_int = self.e2a - self.e2a_rm
         self.e2a_int_nr = nuda.cst.threeFifth * nuda.cst.half * nuda.cst.h2m / ms * \
@@ -239,9 +259,9 @@ class setupFFGNuc():
              (nuda.cst.one-delta)**nuda.cst.fiveThird )
         self.e2a_nr = self.e2a_rm + self.e2a_int_nr 
         #: Attribute FFG energy per unit volum (degeneracy = 2)
-        self.e2v_int = self.e2a_int * self.den
-        self.e2v_int_nr = self.e2a_int_nr * self.den
-        self.e2v_nr = self.e2a_nr * self.den
+        self.eps_int = self.e2a_int * self.den
+        self.eps_int_nr = self.e2a_int_nr * self.den
+        self.eps_nr = self.e2a_nr * self.den
         #: Attribute FFG symmetry energy (degeneracy = 2)
         self.esym_nr = nuda.cst.threeFifth * nuda.cst.half * nuda.cst.h2m / ms * \
            (3*nuda.cst.pi2*nuda.cst.half*den)**nuda.cst.twoThird * \
@@ -257,14 +277,14 @@ class setupFFGNuc():
         #: Attribute FFG pressure (degeneracy = 2)
         self.pre = fpres( 2.0, self.kf_n, ms*nuda.cst.mnc2 ) + fpres( 2.0, self.kf_p, ms*nuda.cst.mpc2 )
         #print('pre:',self.pre)
-        self.pre_nr = nuda.cst.twoThird * self.e2v_int_nr
+        self.pre_nr = nuda.cst.twoThird * self.eps_int_nr
         #print('pre_nr:',self.pre_nr)
         # spline the pressure p(n) to extract dp/dn:
         #cs_pre = CubicSpline( self.den, self.pre )
         #: Attribute enthalpy
-        self.h2a = ( self.e2v + self.pre ) / self.den
+        self.h2a = ( self.eps + self.pre ) / self.den
         #print('h2a:',self.h2a)
-        self.h2a_nr = ( self.e2v_nr + self.pre_nr ) / self.den
+        self.h2a_nr = ( self.eps_nr + self.pre_nr ) / self.den
         #print('h2a_nr:',self.h2a_nr)
         #: Attribute sound speed squared
         dp_dn = self.x_n * f_dp_dn( self.kf_n, ms*nuda.cst.mnc2 ) + self.x_p * f_dp_dn( self.kf_p, ms*nuda.cst.mpc2 )
@@ -273,7 +293,7 @@ class setupFFGNuc():
         #print('dp_dn_num:',dp_dn_num)
         self.cs2 = dp_dn / self.h2a 
         #self.cs2_num = dp_dn_num / self.h2a
-        #self.cs2_nr = 10.0 / 9.0 * self.e2v_nr / self.den / self.h2a_nr
+        #self.cs2_nr = 10.0 / 9.0 * self.eps_nr / self.den / self.h2a_nr
         dp_dn_nr = 10.0 / 9.0 * self.e2a_int_nr
         #print('dp_dn_nr:',dp_dn_nr)
         #print('dp_dn_nr/dp_dn:',dp_dn_nr/dp_dn)
@@ -282,6 +302,7 @@ class setupFFGNuc():
         self.den_unit = 'fm$^{-3}$'
         self.kf_unit = 'fm$^{-1}$'
         self.e2a_unit = 'MeV'
+        self.eps_unit = 'MeV fm$^{-3}$'
         self.pre_unit = 'MeV fm$^{-3}$'
         self.cs2_unit = 'c$^{2}$'
         #
@@ -358,31 +379,44 @@ class setupFFGLep():
         #: Attribute muon Fermi energy (degeneracy = 2)
         self.eF_mu = np.sqrt( nuda.cst.mmuc2**2 + (nuda.cst.hbc*self.kf_mu)**2 )
         #: Attribute FFG energy per particle (degeneracy = 2)
-        self.e2v_el, self.e2a_el = feden( 2.0, self.kf_el, nuda.cst.mec2 )
-        self.e2v_mu, self.e2a_mu = feden( 2.0, self.kf_mu, nuda.cst.mmuc2 )
-        self.e2v_lep = self.e2v_el + self.e2v_mu
-        #self.e2a_el = self.e2v_el / self.den_e
-        self.e2a_el_int = self.e2a_el - nuda.cst.mec2*self.den_el / ( self.den_lep )
-        #self.e2a_mu = self.e2v_mu / self.den_mu
-        self.e2a_mu_int = self.e2a_mu - nuda.cst.mmuc2*self.den_mu / ( self.den_lep )
-        self.e2a_lep = self.e2a_el + self.e2a_mu
-        self.e2a_int_lep = self.e2a_el_int + self.e2a_mu_int
+        # energy
+        self.eps_el, self.e2n_el = feden( 2.0, self.kf_el, nuda.cst.mec2 )
+        self.eps_mu, self.e2n_mu = feden( 2.0, self.kf_mu, nuda.cst.mmuc2 )
+        self.eps_lep = self.eps_el + self.eps_mu
+        self.e2n_lep = self.eps_lep / self.den_lep
+        #self.e2a_el = self.eps_el / self.den_e
+        # internal energy
+        self.eps_el_int = self.eps_el - nuda.cst.mec2*self.den_el
+        self.e2n_el_int = self.eps_el_int / self.den_el
+        self.eps_mu_int = self.eps_mu - nuda.cst.mmuc2*self.den_mu
+        self.e2n_mu_int = np.zeros( np.size(self.den_mu) )
+        for k,n_mu in enumerate(self.den_mu):
+            if n_mu > 0.0:
+                self.e2n_mu_int[k] = self.eps_mu_int[k] / self.den_mu[k]
+        self.eps_lep_int = self.eps_el_int + self.eps_mu_int
+        self.e2n_lep_int = self.eps_lep_int / self.den_lep
         #: Attribute FFG pressure (degeneracy = 2)
         self.pre_el = fpres( 2.0, self.kf_el, nuda.cst.mec2)
         self.pre_mu = fpres( 2.0, self.kf_mu, nuda.cst.mmuc2)
-        self.pre_lep = self.pre_el + self.pre_mu
+        self.pre_lep = self.x_el * self.pre_el + self.x_mu * self.pre_mu
         #: Attribute enthalpy
-        self.h2v_el = self.e2v_el + self.pre_el
-        self.h2v_mu = self.e2v_mu + self.pre_mu
-        self.h2v_lep = self.h2v_el + self.h2v_mu
+        self.h2n_el = self.e2n_el + self.pre_el / self.den_el
+        self.h2n_mu = []
+        for ind,den_mu in enumerate(self.den_mu):
+            if den_mu != 0.0:
+                self.h2n_mu.append( self.e2n_mu[ind] + self.pre_mu[ind] / self.den_mu[ind] )
+            else:
+                self.h2n_mu.append( 0.0 )
+        self.h2n_mu = np.array( self.h2n_mu, dtype = float )
+        self.h2n_lep = self.x_el * self.h2n_el + self.x_mu * self.h2n_mu
         #: Attribute sound speed squared
         dp_dn = self.den_el * f_dp_dn( self.kf_el, nuda.cst.mec2 ) + self.den_mu * f_dp_dn( self.kf_mu, nuda.cst.mmuc2 )
-        self.cs2_lep = dp_dn / self.h2v_lep
+        self.cs2_lep = dp_dn / self.h2n_lep
         #
         self.den_unit = 'fm$^{-3}$'
         self.kf_unit = 'fm$^{-1}$'
-        self.e2a_unit = 'MeV'
-        self.e2v_unit = 'MeV fm$^{-3}$'
+        self.e2n_unit = 'MeV'
+        self.eps_unit = 'MeV fm$^{-3}$'
         self.pre_unit = 'MeV fm$^{-3}$'
         #
         if nuda.env.verb: print("Exit setupFFGLep()")
@@ -396,12 +430,12 @@ class setupFFGLep():
         if nuda.env.verb: print("Enter print_outputs()")
         #
         print("- Print output:")
-        if self.den_e is not None: print(f"   den_e: {np.round(self.den_e,2)} in {self.den_unit}")
+        if self.den_el is not None: print(f"   den_el: {np.round(self.den_el,2)} in {self.den_unit}")
         if self.den_mu is not None: print(f"   den_mu: {np.round(self.den_mu,2)}")
-        if self.kf_e is not None: print(f"   kf_e: {np.round(self.kf_e,2)} in {self.kf_unit}")
-        if self.e2a_int_lep is not None: print(f"   e2a_int_lep: {np.round(self.e2a_int_lep,2)} in {self.e2a_unit}")
+        if self.kf_el is not None: print(f"   kf_el: {np.round(self.kf_el,2)} in {self.kf_unit}")
+        if self.e2n_lep_int is not None: print(f"   e2n_lep_int: {np.round(self.e2n_lep_int,2)} in {self.e2n_unit}")
         if self.pre_lep is not None: print(f"   pre_lep: {np.round(self.pre_lep,2)} in {self.pre_unit}")
-        if self.cs2_lep is not None: print(f"   cs2_lep: {np.round(self.cs2_lep,3)} in {self.cs2_unit}")
+        if self.cs2_lep is not None: print(f"   cs2_lep: {np.round(self.cs2_lep,3)}")
         #
         if nuda.env.verb: print("Exit print_outputs()")
         #

nucleardatapy/matter/setup_hic.py

@@ -9,11 +9,11 @@ sys.path.insert(0, nucleardatapy_tk)
 import nucleardatapy as nuda
 
 
-def hic_constraints():
+def hic_inferences():
     """
-    Return a list of the HIC constraints available in this toolkit 
+    Return a list of the HIC inferences available in this toolkit 
     for the equation of state in SM and NM and print them all on 
-    the prompt. These constraints are the following
+    the prompt. These inferences are the following
     ones: [ '2002-DLL', '2002-KAON', '2016-FOPI', '2011-FOPI-LAND', '2016-ASY-EOS'
     , '2021-SPIRIT','2019-NP-RATIO','2009-ISO-DIFF' ].
 
@@ -21,53 +21,53 @@ def hic_constraints():
     :rtype: list[str].
     """
     #
-    if nuda.env.verb: print("\nEnter hic_constraints()")
+    if nuda.env.verb: print("\nEnter hic_inferences()")
     #
-    constraints = [ '2002-DLL', '2002-KAON', '2016-FOPI', '2009-ISO-DIFF' ,'2011-FOPI-LAND'
+    inferences = [ '2002-DLL', '2002-KAON', '2016-FOPI', '2009-ISO-DIFF' ,'2011-FOPI-LAND'
                    , '2016-ASY-EOS','2019-NP-RATIO', '2021-SPIRIT']
     #
-    constraints_lower = [ item.lower() for item in constraints ]
-    if nuda.env.verb: print('HIC constraints available in the toolkit:',constraints_lower)
+    inferences_lower = [ item.lower() for item in inferences ]
+    if nuda.env.verb: print('HIC inferences available in the toolkit:',inferences_lower)
     #
-    if nuda.env.verb: print("Exit hic_constraints()")
+    if nuda.env.verb: print("Exit hic_inferences()")
     #
-    return constraints, constraints_lower
+    return inferences, inferences_lower
 
 
 class setupHIC():
     """
-    Instantiate the constraints on the EOS from HIC.
+    Instantiate the inferences on the EOS from HIC.
 
-    This choice is defined in the variable `constraint`.
+    This choice is defined in the variable `inference`.
 
-    `constraint` can chosen among the following ones: 
+    `inference` can chosen among the following ones: 
     [ '2002-DLL', '2016-FOPI', '2002-KAON','2009-ISO-DIFF', '2011-FOPI_LAND'
     , '2016-ASY_EOS','2019-NP-RATIO', '2021-SPIRIT' ].
 
-    :param constraint: Default value: '2002-DLL'.
-    :type constraint: str, optional. 
+    :param inference: Default value: '2002-DLL'.
+    :type inferences: str, optional. 
 
     **Attributes:**
     """
-    def __init__(self, constraint = '2002-DLL'):
+    def __init__(self, inference = '2002-DLL'):
         #
         if nuda.env.verb: print("Enter setupHIC()")
         #
-        self.constraint = constraint
-        
-        if nuda.env.verb: print("constraint:",constraint)
+        self.inference = inference
+        #
+        if nuda.env.verb: print("inference:",inference)
         #
         self = setupHIC.init_self( self )
         #
-        constraints, constraints_lower = hic_constraints()
-        if constraint.lower() not in constraints_lower:
-            print('The constraint ',constraint,' is not in the list of EOS HIC constraints.')
-            print('list of EOS HIC constraints:',constraints)
+        inferences, inferences_lower = hic_inferences()
+        if inference.lower() not in inferences_lower:
+            print('The inference ',inference,' is not in the list of EOS HIC inferences.')
+            print('list of EOS HIC inferences:',inferences)
             print('-- Exit the code --')
             exit()
         #
         #
-        if constraint.lower()=='2002-dll':
+        if inference.lower()=='2002-dll':
             #
             file_in1 = nuda.param.path_data+'matter/hic/2002-DLL-SM.dat'
             file_in2 = nuda.param.path_data+'matter/hic/2002-DLL-NM-soft.dat'
@@ -81,10 +81,11 @@ class setupHIC():
             self.label_so = 'DLL-2002-Asy_soft'
             self.label_st = 'DLL-2002-Asy_stiff'
             self.color= 'blue'
+            self.den_err = False
             den2densat, self.sm_pre, self.sm_pre_err = np.loadtxt( file_in1, usecols=(0,1,2), unpack = True )
             den2densat, self.nm_pre_so, self.nm_pre_so_err = np.loadtxt( file_in2, usecols=(0,1,2), unpack = True )
             den2densat, self.nm_pre_st, self.nm_pre_st_err = np.loadtxt( file_in3, usecols=(0,1,2), unpack = True )
-            self.den = den2densat * nuda.cst.nsat # in fm-3
+            self.den_pre = den2densat * nuda.cst.nsat # in fm-3
             self.sm_pre_up = self.sm_pre + self.sm_pre_err
             self.sm_pre_lo = self.sm_pre - self.sm_pre_err
             #self.sm_pre = nuda.cst.half * ( self.sm_pre_up + self.sm_pre_lo )
@@ -92,28 +93,29 @@ class setupHIC():
             #self.nm_pre_so = nuda.cst.half * ( self.nm_pre_so_up + self.nm_pre_so_lo )
             #self.nm_pre_so_err = nuda.cst.half * ( self.nm_pre_so_up - self.nm_pre_so_lo )
             # decide that the NM pressure should be the asy-soft one by default
-            self.nm_pre = self.nm_pre_so
-            self.nm_pre_err = self.nm_pre_so_err
-            self.nm_pre_up = self.nm_pre_so + self.nm_pre_so_err
-            self.nm_pre_lo = self.nm_pre_so - self.nm_pre_so_err
-
+            self.nm_pre_so_up = self.nm_pre_so + self.nm_pre_so_err
+            self.nm_pre_so_lo = self.nm_pre_so - self.nm_pre_so_err
             self.nm_pre_st_up = self.nm_pre_st + self.nm_pre_st_err
             self.nm_pre_st_lo = self.nm_pre_st - self.nm_pre_st_err
+            self.nm_pre = self.nm_pre_so
+            self.nm_pre_err = self.nm_pre_so_err
+            self.nm_pre_up = self.nm_pre_so_up
+            self.nm_pre_lo = self.nm_pre_so_lo
             #
-        elif constraint.lower()=='2002-kaon':
+        elif inference.lower()=='2002-kaon':
             file_in = nuda.param.path_data+'matter/hic/2002-KAON.dat'
             if nuda.env.verb: print('Reads file:',file_in)
             self.ref = 'C. Fuchs, PPNP 56,1 (2006); W. Lynch et al., PPNP 62, 427 (2009).'
             self.note = "Kaon yield ratios studied to contraint symmetric EOS."
             self.label = 'KAON-2002'
             self.color = 'cyan'
+            self.den_err = False
             den2densat, self.sm_pre_lo, self.sm_pre_up = np.loadtxt( file_in, usecols=(0,1,2), unpack = True )
-            self.den = den2densat * nuda.cst.nsat # in fm-3
+            self.den_pre = den2densat * nuda.cst.nsat # in fm-3
             self.sm_pre = nuda.cst.half * ( self.sm_pre_up + self.sm_pre_lo )
             self.sm_pre_err = nuda.cst.half * ( self.sm_pre_up - self.sm_pre_lo )
-            
-
-        elif constraint.lower()=='2016-fopi':
+            #
+        elif inference.lower()=='2016-fopi':
             #
             file_in1 = nuda.param.path_data+'matter/hic/2016-FOPI-E2A.dat'
             file_in2 = nuda.param.path_data+'matter/hic/2016-FOPI-SM.dat'
@@ -122,72 +124,83 @@ class setupHIC():
             self.note = "Elliptical flow data is used to constraint symmetric matter EOS"
             self.label = 'FOPI-2016'
             self.color = 'magenta'
-            den2densat_e2a, self.sm_e2a, self.sm_e2a_err = np.loadtxt( file_in1, usecols=(0,1,2), unpack = True )
+            self.den_err = False
+            den2densat_e2a, self.sm_e2a_int, self.sm_e2a_err = np.loadtxt( file_in1, usecols=(0,1,2), unpack = True )
             den2densat, self.sm_pre_lo, self.sm_pre_up = np.loadtxt( file_in2, usecols=(0,1,2), unpack = True )
             self.den_e2a = den2densat_e2a * nuda.cst.nsat # in fm-3
-            self.den = den2densat * nuda.cst.nsat # in fm-3
-            self.sm_e2a_up = self.sm_e2a + self.sm_e2a_err
-            self.sm_e2a_lo = self.sm_e2a - self.sm_e2a_err
+            self.den_pre = den2densat * nuda.cst.nsat # in fm-3
+            self.sm_e2a_int_up = self.sm_e2a_int + self.sm_e2a_err
+            self.sm_e2a_int_lo = self.sm_e2a_int - self.sm_e2a_err
             self.sm_pre = nuda.cst.half * ( self.sm_pre_up + self.sm_pre_lo )
             self.sm_pre_err = nuda.cst.half * ( self.sm_pre_up - self.sm_pre_lo )
             #
-        elif constraint.lower()=='2009-iso-diff':
+        elif inference.lower()=='2009-iso-diff':
+            #
             file_in = nuda.param.path_data+'matter/hic/2009-iso-diff.dat'
             if nuda.env.verb: print('Reads file:',file_in)
             self.ref = 'M. B. Tsang et al., Phys. Rev. Lett. 102, 122701 (2009); W. Lynch, M. B. Tsang, Phys. Lett. B 830, 137098 (2022).'
             self.note = "Isospin diffusion data studied to constraint symmetry energy at sub-saturation."
             self.label = 'Iso Diff-2009'
             self.color = 'k'
-            den2densat, den2densat_err, self.sym_enr_isodiff, self.sym_enr_isodiff_err = np.loadtxt( file_in, usecols=(0,1,2,3), unpack = True )
-            self.den_isodiff = den2densat * nuda.cst.nsat # in fm-3
-            self.den_isodiff_err = den2densat_err * nuda.cst.nsat # in fm-3
+            self.den_err = True
+            den2densat, den2densat_err, self.esym, self.esym_err = np.loadtxt( file_in, usecols=(0,1,2,3), unpack = True )
+            self.den_esym = den2densat * nuda.cst.nsat # in fm-3
+            self.den_esym_err = den2densat_err * nuda.cst.nsat # in fm-3
+            #
+        elif inference.lower()=='2011-fopi-land':
             #
-        elif constraint.lower()=='2011-fopi-land':
             file_in = nuda.param.path_data+'matter/hic/2011-FOPI-LAND.dat'
             if nuda.env.verb: print('Reads file:',file_in)
             self.ref = ' P. Russotto et al., Physics Letters B 697, 471 (2011).'
             self.note = "Sura-saturation information on symmtery energy usinf n/p elliptical flows"
             self.label = 'FOPI-LAND-2011'
             self.color = 'yellow'
-            den2densat, self.sym_enr_lo, self.sym_enr_up = np.loadtxt( file_in, usecols=(0,1,2), unpack = True )
-            self.den = den2densat * nuda.cst.nsat # in fm-3
-            self.sym_enr = nuda.cst.half * ( self.sym_enr_up + self.sym_enr_lo )
-            self.sym_enr_err = nuda.cst.half * ( self.sym_enr_up - self.sym_enr_lo )
-            
-        elif constraint.lower()=='2016-asy-eos':
+            self.den_err = False
+            den2densat, self.esym_lo, self.esym_up = np.loadtxt( file_in, usecols=(0,1,2), unpack = True )
+            self.den_esym = den2densat * nuda.cst.nsat # in fm-3
+            self.esym = nuda.cst.half * ( self.esym_up + self.esym_lo )
+            self.esym_err = nuda.cst.half * ( self.esym_up - self.esym_lo )
+            #
+        elif inference.lower()=='2016-asy-eos':
+            #
             file_in = nuda.param.path_data+'matter/hic/2016-ASY-EOS.dat'
             if nuda.env.verb: print('Reads file:',file_in)
             self.ref = 'P. Russotto et al., Phys. Rev. C 94, 034608 (2016).'
             self.note = "Sura-saturation information on symmtery energy usinf n/p elliptical flows."
             self.label = 'ASY-EOS-2016'
             self.color = 'red'
-            den2densat, self.sym_enr_lo, self.sym_enr_up = np.loadtxt( file_in, usecols=(0,1,2), unpack = True )
-            self.den = den2densat * nuda.cst.nsat # in fm-3
-            self.sym_enr = nuda.cst.half * ( self.sym_enr_up + self.sym_enr_lo )
-            self.sym_enr_err = nuda.cst.half * ( self.sym_enr_up - self.sym_enr_lo )
-        #   
-        elif constraint.lower()=='2019-np-ratio':
+            self.den_err = False
+            den2densat, self.esym_lo, self.esym_up = np.loadtxt( file_in, usecols=(0,1,2), unpack = True )
+            self.den_esym = den2densat * nuda.cst.nsat # in fm-3
+            self.esym = nuda.cst.half * ( self.esym_up + self.esym_lo )
+            self.esym_err = nuda.cst.half * ( self.esym_up - self.esym_lo )
+            #
+        elif inference.lower()=='2019-np-ratio':
+            #
             file_in = nuda.param.path_data+'matter/hic/2019-n2p-ratio.dat'
             if nuda.env.verb: print('Reads file:',file_in)
             self.ref = 'P. Morfouace et al., Phys. Lett. B 799, 135045 (2019).'
             self.note = "n/p spectral ratios to constraint symmetry energy at sub-saturation densities"
             self.label = 'n/p ratio-2019'
             self.color = 'green'
-            den2densat, den2densat_err, self.sym_enr_np, self.sym_enr_np_err = np.loadtxt( file_in, usecols=(0,1,2,3), unpack = True )
-            self.den_np = den2densat * nuda.cst.nsat # in fm-3
-            self.den_np_err = den2densat_err * nuda.cst.nsat # in fm-3
-        #    
-        elif constraint.lower()=='2021-spirit':
+            self.den_err = True
+            den2densat, den2densat_err, self.esym, self.esym_err = np.loadtxt( file_in, usecols=(0,1,2,3), unpack = True )
+            self.den_esym = den2densat * nuda.cst.nsat # in fm-3
+            self.den_esym_err = den2densat_err * nuda.cst.nsat # in fm-3
+            #
+        elif inference.lower()=='2021-spirit':
+            #
             file_in = nuda.param.path_data+'matter/hic/2021-SPIRIT.dat'
             if nuda.env.verb: print('Reads file:',file_in)
             self.ref = 'J. Estee et al., Phys. Rev. Lett. 126, 162701 (2021).'
             self.note = "Pion double ratios is studied in neutron rich and poor colliding nuclei."
             self.label = 'SPIRIT-2021'
             self.color = 'blue'
-            den2densat, den2densat_err, self.sym_enr_spirit, self.sym_enr_spirit_err, self.sym_pre_spirit, self.sym_pre_spirit_err = np.loadtxt( file_in, usecols=(0,1,2,3,4,5), unpack = True )
-            self.den_spirit = den2densat * nuda.cst.nsat # in fm-3
-            self.den_spirit_err = den2densat_err * nuda.cst.nsat # in fm-3
-    #
+            self.den_err = True
+            den2densat, den2densat_err, self.esym, self.esym_err, self.psym, self.psym_err = np.loadtxt( file_in, usecols=(0,1,2,3,4,5), unpack = True )
+            self.den_esym = den2densat * nuda.cst.nsat # in fm-3
+            self.den_esym_err = den2densat_err * nuda.cst.nsat # in fm-3
+            #
     def print_outputs( self ):
        """
        Method which print outputs on terminal's screen.
@@ -197,7 +210,7 @@ class setupHIC():
        if nuda.env.verb: print("Enter print_outputs()")
        #
        print("- Print output:")
-       print("   constraint:",self.constraint)
+       print("   inference:",self.inference)
        print("   ref:     ",self.ref)
        print("   label:   ",self.label)
        print("   note:    ",self.note)
@@ -235,11 +248,11 @@ class setupHIC():
         #: Attribute the density of the system (in fm^-3).
         self.den = None
         #: Attribute the upper limit of the energy per particle in SM (in MeV).
-        self.sm_e2a_up = None
+        self.sm_e2a_int_up = None
         #: Attribute the lower limit of the energy per particle in SM (in MeV).
-        self.sm_e2a_lo = None
+        self.sm_e2a_int_lo = None
         #: Attribute the energy per particle in SM (in MeV).
-        self.sm_e2a = None
+        self.sm_e2a_int = None
         #: Attribute the uncertainty in the energy per particle in SM (in MeV fm-3).
         self.sm_e2a_err = None
         #: Attribute the upper limit of the pressure in SM (in MeV fm-3).
@@ -277,9 +290,9 @@ class setupHIC():
         #: Attribute the lower limit of the pressure in NM for asy-stiff EOS (in MeV fm-3).
         self.nm_pre_st_lo = None
         #: Attribute the upper edge of the symmetry energy for FOPI/ASY-EOS constraints (in MeV).
-        self.sym_enr_lo = None
+        self.esym_lo = None
         #: Attribute the upper edge of the symmetry energy for FOPI/ASY-EOS constraints (in MeV).
-        self.sym_enr_up = None
+        self.esym_up = None
         #: Attribute the centroids of the symmetry energy for FOPI/ASY-EOS constraints (in MeV).
         self.sym_enr = None
         #: Attribute the uncertainty of the symmetry energy for FOPI/ASY-EOS constraints (in MeV).
@@ -288,10 +301,14 @@ class setupHIC():
         self.sym_enr_spirit = None
         #: Attribute the uncertainty of the symmetry energy for SPIRIT constraints (in MeV).
         self.sym_enr_spirit_err = None
-        #: Attribute the centroid of the symmetry energy for isospin difusion (in MeV).
-        self.sym_enr_isodiff = None
-        #: Attribute the uncertainty of the symmetry energy for isospin difusion (in MeV).
-        self.sym_enr_isodiff_err = None
+        #: Attribute the upper boundary of the symmetry energy.
+        self.esym_up = None
+        #: Attribute the lower boundary of the symmetry energy.
+        self.esym_lo = None
+        #: Attribute the centroid of the symmetry energy.
+        self.esym = None
+        #: Attribute the uncertainty of the symmetry energy.
+        self.esym_err = None
         #: Attribute the centroid of the symmetry energy for isospin difusion (in MeV).
         self.sym_enr_np = None
         #: Attribute the uncertainty of the symmetry energy for isospin difusion (in MeV).