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

nucleardatapy/fig/nuc_setupRnpExp_fig.py

@@ -0,0 +1,88 @@
+import os
+import numpy as np
+import matplotlib.pyplot as plt
+from matplotlib.ticker import AutoMinorLocator
+
+import nucleardatapy as nuda
+
+SOURCE_LABELS_LATEX = {
+    "48Ca": r"$^{48}\mathrm{Ca}$",
+    "208Pb": r"$^{208}\mathrm{Pb}$"
+}
+
+def nuc_setupRnpExp_fig(pname=None, source=None):
+    print(f"Using source: {source}")
+    
+    # subplot_label = "(a)"
+    subplot_label = " "
+
+    if source is None:
+        print("Erro: nenhum source fornecido.")
+        return
+
+    labels = []
+    rskin_values = []
+    error_lower = []
+    error_upper = []
+    markers = []
+
+    cals = nuda.nuc.rnp_exp_source(source)
+
+    for cal in cals:
+        neutron_skin_calc = nuda.nuc.setupRnpExp(source=source, cal=cal)
+
+        if neutron_skin_calc.rnp is not None:
+            labels.append(neutron_skin_calc.label)
+            rskin_values.append(neutron_skin_calc.rnp)
+
+            err_down = neutron_skin_calc.rnp_sig_lo if neutron_skin_calc.rnp_sig_lo is not None else 0.0
+            err_up = neutron_skin_calc.rnp_sig_up if neutron_skin_calc.rnp_sig_up is not None else 0.0
+            error_lower.append(err_down)
+            error_upper.append(err_up)
+
+            marker = getattr(neutron_skin_calc, "marker", 'o')
+            markers.append(marker)
+
+    if not rskin_values:
+        print(f"Nenhum dado disponível para {source}.")
+        return
+
+    fig, ax = plt.subplots(figsize=(10, 8))
+    x_positions = range(len(labels) + 1)
+
+    for i, (x, y, err_down, err_up, marker) in enumerate(zip(x_positions, rskin_values, error_lower, error_upper, markers)):
+        adjusted_err_down = min(err_down, 0.2)
+        adjusted_err_up = min(err_up, 0.2)
+
+        ax.errorbar(x, y, yerr=[[adjusted_err_down], [adjusted_err_up]], fmt=marker, markersize=8, capsize=0, label=labels[i])
+
+        if err_down >= 1000:
+            ax.plot([x], [y - adjusted_err_down], marker="v", color="grey", markersize=8)
+        if err_up >= 1000:
+            ax.plot([x], [y + adjusted_err_up], marker="^", color="grey", markersize=8)
+
+    nsav = nuda.nuc.setupRnpAverage(source=source)
+    if nsav.rnp_cen is not None:
+        ax.errorbar(len(labels), nsav.rnp_cen, yerr=nsav.sig_std, label=nsav.label,
+                    color='red', marker='o', markersize=10, linestyle='solid', linewidth=3)
+        labels.append(nsav.label)
+
+    ax.set_ylim([0, 0.5])
+    ax.set_xticks(x_positions)
+    ax.set_xticklabels(labels, rotation=45, ha="right", fontsize=15)
+    ax.tick_params(axis='y', labelsize=15)
+    ax.set_ylabel(rf"$R_{{\rm{{skin}}}}$ {SOURCE_LABELS_LATEX[source]} (fm)", fontsize=15)
+    ax.text(0.95, 0.95, subplot_label, transform=ax.transAxes, fontsize=15,
+            verticalalignment='top', horizontalalignment='right')
+    ax.yaxis.set_minor_locator(AutoMinorLocator())
+    ax.tick_params(axis='y', which='minor', length=4, color='gray')
+
+    if pname is None:
+        output_dir = os.path.abspath("figs/")
+        os.makedirs(output_dir, exist_ok=True)
+        pname = os.path.join(output_dir, f"plot_nuc_setupRnp_Exp_{source.replace(' ', '_')}.png")
+
+    plt.tight_layout()
+    plt.savefig(pname, dpi=200)
+    plt.close()
+    print(f"Plot saved: {pname}")

nucleardatapy/fig/nuc_setupRnpTheo_fig.py

@@ -0,0 +1,133 @@
+import numpy as np
+import matplotlib.pyplot as plt
+from matplotlib.ticker import AutoMinorLocator  # Import para minor ticks
+import os
+
+import nucleardatapy as nuda
+
+# Dictionary to map sources to LaTeX names
+SOURCE_LABELS_LATEX = {
+    "48Ca": r"$^{48}\mathrm{Ca}$",
+    "208Pb": r"$^{208}\mathrm{Pb}$"
+}
+
+# Directory containing the model data tables
+MODEL_TABLES_DIR = nuda.param.path_data + '/rnp/'
+
+# Define markers and colors for each model
+MODEL_STYLES = {
+    "skyrme": {"color": "blue", "marker": "s", "label": "Skyrme"},
+    "nlrh": {"color": "red", "marker": "^", "label": "NLRH"},
+    "ddrh": {"color": "green", "marker": "D", "label": "DDRH"}
+}
+
+def read_model_data(directory, source):
+    model_data = {}
+    for filename in os.listdir(directory):
+        if filename.endswith(".dat") and source in filename:
+            model_name = filename.split("rnp")[0].lower()
+            filepath = os.path.join(directory, filename)
+            data = []
+            with open(filepath, 'r') as file:
+                for line in file:
+                    if line.startswith("#") or not line.strip():
+                        continue
+                    parts = line.split()
+                    param = parts[0]
+                    rn, rp, rskin = map(float, parts[1:4])
+                    data.append((param, rn, rp, rskin))
+            model_data[model_name] = np.array(data, dtype=object)
+            if nuda.env.verb:
+                print(f"Loaded model data for {model_name}: {model_data[model_name]}")
+    return model_data
+
+def nuc_setupRnpTheo_fig(pname, source):
+    print(f'Plot name: {pname}')
+    print(f'Using source: {source}')
+
+    # subplot_label = "(a)"  # remove this if you don't want any labels
+    subplot_label = " "
+
+    labels = []
+    rskin_values = []
+    error_lower = []
+    error_upper = []
+    xexp = []
+
+    cals = nuda.nuc.rnp_exp_source(source)
+    for i, cal in enumerate(cals):
+        neutron_skin_calc = nuda.nuc.setupRnpExp(source=source, cal=cal)
+        if neutron_skin_calc.rnp is not None:
+            labels.append(neutron_skin_calc.label)
+            rskin_values.append(neutron_skin_calc.rnp)
+            xexp.append(i)
+            err_down = neutron_skin_calc.rnp_sig_lo if neutron_skin_calc.rnp_sig_lo is not None else 0.0
+            err_up = neutron_skin_calc.rnp_sig_up if neutron_skin_calc.rnp_sig_up is not None else 0.0
+            error_lower.append(err_down)
+            error_upper.append(err_up)
+
+    model_data = read_model_data(MODEL_TABLES_DIR, source)
+    combined_rskin = []
+    combined_errors = []
+    combined_markers = []
+    combined_colors = []
+    xtheo = []
+
+    for model_name, data in model_data.items():
+        for j, (_, _, _, rskin) in enumerate(data):
+            x_position = xexp[j % len(xexp)] + 0.5  # Cycle through xexp and offset by 0.5
+            xtheo.append(x_position)
+            combined_rskin.append(rskin)
+            combined_errors.append((0.0, 0.0))
+            combined_markers.append(MODEL_STYLES[model_name]["marker"])
+            combined_colors.append(MODEL_STYLES[model_name]["color"])
+
+    if nuda.env.verb:
+        print(f"Experimental positions for {source}: {xexp}")
+        print(f"Theoretical positions for {source}: {xtheo}")
+
+    fig, ax = plt.subplots(figsize=(10, 8))
+    for i, (x, y, err_down, err_up) in enumerate(zip(xexp, rskin_values, error_lower, error_upper)):
+        adjusted_err_down = min(err_down, 0.2)
+        adjusted_err_up = min(err_up, 0.2)
+        ax.errorbar(x, y, yerr=[[adjusted_err_down], [adjusted_err_up]], fmt='o', markersize=8, capsize=0, color='black', markerfacecolor='none')
+        if err_down >= 1000:
+            ax.plot([x], [y - adjusted_err_down], marker="v", color="black", markersize=8)
+        if err_up >= 1000:
+            ax.plot([x], [y + adjusted_err_up], marker="^", color="black", markersize=8)
+
+    nsav = nuda.nuc.setupRnpAverage(source=source)
+    if nsav.rnp_cen is not None:
+        ax.errorbar(len(labels), nsav.rnp_cen, yerr=nsav.sig_std, label=nsav.label, 
+                   color='k', marker='o', markersize=10, linestyle='solid', linewidth=3)           
+
+    for i, (x, y, marker, color) in enumerate(zip(xtheo, combined_rskin, combined_markers, combined_colors)):
+        ax.plot(x, y, marker=marker, markersize=8, color=color)
+
+    # Add legend for experimental points
+    ax.scatter([], [], color='black', marker='o', facecolors='none', label='Experimental/Analysis')
+
+    # Add legend for theoretical models
+    for model_name, style in MODEL_STYLES.items():
+        ax.scatter([], [], color=style["color"], marker=style["marker"], label=style["label"])
+
+    ax.set_ylim([0, 0.5])
+    ax.set_xticks(xexp)
+    ax.set_xticklabels(labels, rotation=45, ha="right", fontsize=15)
+    ax.set_ylabel(rf"$R_{{\rm{{skin}}}}$ {SOURCE_LABELS_LATEX[source]} (fm)", fontsize=15)
+    ax.text(0.95, 0.95, subplot_label, transform=ax.transAxes, fontsize=15, verticalalignment='top', horizontalalignment='right')
+
+    # Add minor ticks on y-axis
+    ax.yaxis.set_minor_locator(AutoMinorLocator())
+    ax.tick_params(axis='y', which='minor', length=4, color='gray')
+
+    ax.legend(loc="upper right", bbox_to_anchor=(0.5, 1), fontsize=12)
+
+    output_dir = "figs/"
+    os.makedirs(output_dir, exist_ok=True)
+    fig_name = os.path.join(output_dir, f"plot_nuc_setupRnpTheo_source{source.replace(' ', '_')}.png")    
+    plt.tight_layout()
+    plt.savefig(fig_name, dpi=200)
+    plt.close()
+
+    print(f"Plot saved: {fig_name}")

nucleardatapy/hello.py

@@ -0,0 +1,6 @@
+
+def hello():
+	print('hello world!')
+
+if __name__ == "__main__":
+    hello()

nucleardatapy/hnuc/__init__.py

@@ -2,7 +2,7 @@
 This module provides microscopic, phenomenological and experimental data constraints.
 """
 #
-from nucleardatapy.hnuc.setup_be1L_exp     import *
-from nucleardatapy.hnuc.setup_be2L_exp     import *
-from nucleardatapy.hnuc.setup_be1Xi_exp    import *
+from nucleardatapy.hnuc.setup_re1L_exp     import *
+from nucleardatapy.hnuc.setup_re2L_exp     import *
+from nucleardatapy.hnuc.setup_re1Xi_exp    import *
 #

nucleardatapy/hnuc/{setup_be1L_exp.py → setup_re1L_exp.py}

@@ -150,16 +150,16 @@ class setupRE1LExp():
       #: Attribute A (mass of the nucleus).
       self.A = self.Z + self.N + np.ones(len(self.N),dtype=int)
       #: charge of the hypernuclei (=Z, since Lamnda is charged 0)
-      self.ch = self.Z
+      self.Q = self.Z
       #: Strangness number
-      self.S = -2*np.ones(len(self.N),dtype=int)
+      self.S = -1*np.ones(len(self.N),dtype=int)
       #: symbol representing the nucleus
       self.symb = nucSymb
-      #: Attribute the angular momentum of the state.
+      #: Attribute the s.p. state.
       self.sps = nucsps
       #: Attribute the angular momentum of the state.
       self.ell = np.array( nucell, dtype = int )
-      #: Attribute 1L binding energy in MeV.
+      #: Attribute 1L removal energy in MeV.
       self.lre = np.array( nuclre, dtype = float )
       #: Attribute 1L binding energy error in MeV.
       self.lre_err = np.array( nuclre_err, dtype = float )
@@ -204,7 +204,7 @@ class setupRE1LExp():
       if any(self.Z): print(f"   Z: {self.Z}")
       if any(self.N): print(f"   N: {self.N}")
       if any(self.S): print(f"   S: {self.S}")
-      if any(self.ch): print(f"  ch: {self.ch}")
+      if any(self.Q): print(f"   Q: {self.Q}")
       if any(self.symb): print(f" symb: {self.symb}")
       if any(self.ell): print(f" ell: {self.ell}")
       if any(self.lre): print(f"  re: {self.lre}")
@@ -225,7 +225,7 @@ class setupRE1LExp():
          print(rf" index & Z & N & S & ch & symb & $RE$  & Ref. \\\\")
          print(rf"       &   &   &   &    &      & (MeV) &      \\\\")
          for ind,A in enumerate(self.A):
-            print(rf" {ind} & {self.Z[ind]} & {self.N[ind]} & {self.S[ind]} & {self.ch[ind]} & {self.symb[ind]} & ${self.lre[ind]:.3f}\pm {self.lre_err[ind]:.3f}$ & \cite{{"+self.keyref+"} \\\\")
+            print(rf" {ind} & {self.Z[ind]} & {self.N[ind]} & {self.S[ind]} & {self.Q[ind]} & {self.symb[ind]} & ${self.lre[ind]:.3f}\pm {self.lre_err[ind]:.3f}$ & \cite{{"+self.keyref+"} \\\\")
       else:
          print(f"- No  table for source {self.table} (average). To get table, write 'verb_latex = True' in env.py.")
       #

nucleardatapy/hnuc/{setup_be1Xi_exp.py → setup_re1Xi_exp.py}

@@ -117,12 +117,12 @@ class setupRE1XiExp():
       #: Strangness number (to be checked)
       self.S = -2*np.ones(len(self.N),dtype=int)
       #: charge of the hypernuclei (=Z-1, since Xi is charged -1)
-      self.ch = self.Z - np.ones(len(self.N),dtype=int)
+      self.Q = self.Z - np.ones(len(self.N),dtype=int)
       #: symbol representing the nucleus
       self.symb = nucSymb
-      #: Attribute 1L binding energy in MeV.
+      #: Attribute 1Xi removal energy in MeV.
       self.xire = np.array( nucxire, dtype = float )
-      #: Attribute 1L binding energy error in MeV.
+      #: Attribute 1Xi removal energy error in MeV.
       self.xire_err = np.array( nucxire_err, dtype = float )
       #: Attribute the probe.
       self.probe = probe
@@ -162,7 +162,7 @@ class setupRE1XiExp():
       if any(self.Z): print(f"   Z: {self.Z}")
       if any(self.N): print(f"   N: {self.N}")
       if any(self.S): print(f"   S: {self.S}")
-      if any(self.ch): print(f"  ch: {self.ch}")
+      if any(self.Q): print(f"   Q: {self.Q}")
       if any(self.symb): print(f" symb: {self.symb}")
       if any(self.xire): print(f" xibe: {self.xire}")
       if any(self.xire_err): print(f" xibe_err: {self.xire_err}")
@@ -183,7 +183,7 @@ class setupRE1XiExp():
          print(rf" index & Z & N & S & ch & symb & RE    & & Ref. \\\\")
          print(rf"       &   &   &   &    &      & (MeV) & & \\\\")
          for ind,A in enumerate(self.A):
-            print(rf" {ind} & {self.Z[ind]} & {self.N[ind]} & {self.S[ind]} & {self.ch[ind]} & {self.symb[ind]} & ${self.xire[ind]:.3f}\pm {self.xire_err[ind]:.3f}$ & & \\cite{{"+self.keyref+"} \\\\")
+            print(rf" {ind} & {self.Z[ind]} & {self.N[ind]} & {self.S[ind]} & {self.Q[ind]} & {self.symb[ind]} & ${self.xire[ind]:.3f}\pm {self.xire_err[ind]:.3f}$ & & \\cite{{"+self.keyref+"} \\\\")
       else:
          print(f"- No  table for source {self.table} (average). To get table, write 'verb_latex = True' in env.py.")
       #

nucleardatapy/hnuc/{setup_be2L_exp.py → setup_re2L_exp.py}

@@ -119,14 +119,14 @@ class setupRE2LExp():
       #: Attribute A (mass of the nucleus).
       self.A = self.Z + self.N + 2*np.ones(len(self.N),dtype=int)
       #: charge of the hypernuclei (=Z, since Lamnda is charged 0)
-      self.ch = self.Z
+      self.Q = self.Z
       #: Strangness number
-      self.S = -2*2*np.ones(len(self.N),dtype=int)
+      self.S = -2*np.ones(len(self.N),dtype=int)
       #: symbol representing the nucleus
       self.symb = nucSymb
-      #: Attribute 2L binding energy in MeV.
+      #: Attribute 2L removal energy in MeV.
       self.llre = np.array( nuclre, dtype = float )
-      #: Attribute 2L binding energy error in MeV.
+      #: Attribute 2L removal energy error in MeV.
       self.llre_err = np.array( nuclre_err, dtype = float )
       #: Attribute 2L bond energy in MeV.
       self.lldre = np.array( nucldre, dtype = float )
@@ -171,7 +171,7 @@ class setupRE2LExp():
       if any(self.Z): print(f"   Z: {self.Z}")
       if any(self.N): print(f"   N: {self.N}")
       if any(self.S): print(f"   S: {self.S}")
-      if any(self.ch): print(f"  ch: {self.ch}")
+      if any(self.Q): print(f"   Q: {self.Q}")
       if any(self.symb): print(f" symb: {self.symb}")
       if any(self.llre): print(f" be: {self.llre}")
       if any(self.llre_err): print(f" be_err: {self.llre_err}")
@@ -193,7 +193,7 @@ class setupRE2LExp():
          print(rf" index & Z & N & S & ch & symb & $RE$  & $\Delta RE$ & Ref. \\\\")
          print(rf"       &   &   &   &    &      & (MeV) & (MeV)       & \\\\")
          for ind,A in enumerate(self.A):
-            print(rf" {ind} & {self.Z[ind]} & {self.N[ind]} & {self.S[ind]} & {self.ch[ind]} & {self.symb[ind]} & ${self.llre[ind]:.3f}\pm {self.llre_err[ind]:.3f}$ & ${self.lldre[ind]:.3f}\pm {self.lldre_err[ind]:.3f}$ & \\cite{{"+self.keyref+"}  \\\\")
+            print(rf" {ind} & {self.Z[ind]} & {self.N[ind]} & {self.S[ind]} & {self.Q[ind]} & {self.symb[ind]} & ${self.llre[ind]:.3f}\pm {self.llre_err[ind]:.3f}$ & ${self.lldre[ind]:.3f}\pm {self.lldre_err[ind]:.3f}$ & \\cite{{"+self.keyref+"}  \\\\")
       else:
          print(f"- No  table for source {self.table} (average). To get table, write 'verb_latex = True' in env.py.")
       #

nucleardatapy/matter/__init__.py

@@ -2,16 +2,17 @@
 This module provides microscopic, phenomenological and experimental data constraints.
 """
 #
-from nucleardatapy.matter.setup_ffg            import *
-from nucleardatapy.matter.setup_micro          import *
-from nucleardatapy.matter.setup_micro_effmass  import *
-from nucleardatapy.matter.setup_micro_gap      import *
-from nucleardatapy.matter.setup_micro_band     import *
-from nucleardatapy.matter.setup_micro_esym     import *
-from nucleardatapy.matter.setup_micro_lp       import *
-from nucleardatapy.matter.setup_pheno          import *
-from nucleardatapy.matter.setup_pheno_esym     import *
-from nucleardatapy.matter.setup_hic            import *
-from nucleardatapy.matter.setup_nep            import *
-from nucleardatapy.matter.setup_nep_dist       import *
-from nucleardatapy.matter.setup_nep_model_dist import *
+from nucleardatapy.matter.setup_ffg             import *
+from nucleardatapy.matter.setup_micro           import *
+from nucleardatapy.matter.setup_micro_effmass   import *
+from nucleardatapy.matter.setup_micro_gap       import *
+from nucleardatapy.matter.setup_micro_band      import *
+from nucleardatapy.matter.setup_micro_esym      import *
+from nucleardatapy.matter.setup_micro_lp        import *
+from nucleardatapy.matter.setup_pheno           import *
+from nucleardatapy.matter.setup_pheno_esym      import *
+from nucleardatapy.matter.setup_hic             import *
+from nucleardatapy.matter.setup_nep             import *
+from nucleardatapy.matter.setup_nep_stat_model  import *
+from nucleardatapy.matter.setup_nep_stat_models import *
+from nucleardatapy.matter.setup_check           import *

nucleardatapy/matter/setup_check.py

@@ -0,0 +1,97 @@
+
+import numpy as np  # 1.15.0
+from scipy.interpolate import CubicSpline
+
+import nucleardatapy as nuda
+
+class setupCheck():
+    """
+    Instantiate a flag reflecting if e2a from `eos` passes through the reference `band` or not.
+
+    :param eos: object containing the eos to check.
+    :type eos: object. 
+    :param band: object containing the band employed to check the eos.
+    :type band: object. 
+    
+    **Attributes:**
+    
+    """
+    #
+    def __init__( self, eos, band ):
+        """
+        Parameters
+        ----------
+        eos : object.
+        Object containing the eos to check.
+        band: object.
+        Object containing the band employed to check the eos.
+        matter: string which can be: 'nm' (default), 'sm', or 'esym'.
+        """
+        #
+        if nuda.env.verb: print("Enter setupCheck()")
+        #
+        #: Attribute the object `eos`.
+        self.eos = eos
+        #: Attribute the object `band`.
+        self.band = band
+        #: Attribute the value for the variable `matter`.
+        self.matter = band.matter
+        #
+        if 'fit' in eos.model:
+            if self.matter.lower() == 'nm' or self.matter.lower() == 'sm':
+                self.x = np.insert( eos.den, 0, 0.0 )
+                self.y = np.insert( eos.e2a_int, 0, 0.0 )
+            elif self.matter.lower() == 'esym':
+                self.x = np.insert( eos.den, 0, 0.0 )
+                self.y = np.insert( eos.esym, 0, 0.0 )
+        else:
+            if self.matter.lower() == 'nm':
+                self.x = np.insert( eos.nm_den, 0, 0.0 )
+                self.y = np.insert( eos.nm_e2a_int, 0, 0.0 )
+            elif self.matter.lower() == 'sm':
+                self.x = np.insert( eos.sm_den, 0, 0.0 )
+                self.y = np.insert( eos.sm_e2a_int, 0, 0.0 )
+            elif self.matter.lower() == 'esym':
+                self.x = np.insert( eos.den, 0, 0.0 )
+                self.y = np.insert( eos.esym, 0, 0.0 )
+            else:
+                print('setup_check: issue with matter:',self.matter)
+                exit()
+        cs_e2a = CubicSpline( self.x, self.y )
+        self.eos_e2a_int = cs_e2a(band.den)
+        flag = True
+        for ind,den in enumerate(band.den):
+            #if abs(cs_e2a(den)-band.e2a[ind]) > band.e2a_std[ind]:
+            if abs(self.eos_e2a_int[ind]-band.e2a_int[ind]) > band.e2a_std[ind]:
+                flag = False
+        #: Attribute is eos is inside the band.
+        self.isInside = flag
+        #: Attribute is eos is outside the band.
+        self.isOutside = not flag
+        #
+        self.den_unit = 'fm$^{-3}$'
+        self.e2a_unit = 'MeV'
+        #
+        if nuda.env.verb: print("Exit setupCheck()")
+    #
+    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('For matter:',self.matter)
+        print('EOS:')
+        if self.x is not None: print(f"   den: {np.round(self.x,2)} in {self.den_unit}")
+        if self.y is not None: print(f"   e2a: {np.round(self.y,2)} in {self.e2a_unit}")
+        print('BAND:')
+        if self.band.den is not None: print(f"   den: {np.round(self.band.den,2)} in {self.den_unit}")
+        if self.band.e2a_int is not None: print(f"   e2a_int: {np.round(self.band.e2a_int,2)} in {self.e2a_unit}")
+        if self.band.e2a_std is not None: print(f"   e2a_std: {np.round(self.band.e2a_std,2)} in {self.e2a_unit}")
+        #
+        if nuda.env.verb: print("Exit print_outputs()")
+        #
+