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

nucleardatapy/fig/matter_setupMicro_gap_fig.py

@@ -5,15 +5,20 @@ import nucleardatapy as nuda
 
 def matter_setupMicro_gap_1s0_fig(pname, models, matter="NM"):
     """
-    Plot the correlation between Esym and Lsym.\
-    The plot is 2x2 with:\
-    [0,0]: Esym - Lsym correlation plot
+    Plot the 1S0 pairing gap in `matter`.
+
+    The plot is 2x2 with:
+
+    [0,0]: gap versus den.    [0,1]: gap versus kfn.
+
+    [1,0]: gap/EF versus den. [1,1]: gap/EF versus kfn.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
-    :param constraints: list of constraints to run on.
-    :type constraints: array of str.
-
+    :param models: models to run on.
+    :type models: array of str.
+    :param matter: can be 'SM' or 'NM' (default).
+    :type matter: str.
     """
     #
     print(f"Plot name: {pname}")
@@ -100,15 +105,20 @@ def matter_setupMicro_gap_1s0_fig(pname, models, matter="NM"):
 
 def matter_setupMicro_gap_3pf2_fig(pname, models, matter="NM"):
     """
-    Plot the correlation between Esym and Lsym.\
-    The plot is 2x2 with:\
-    [0,0]: Esym - Lsym correlation plot
+    Plot the 3PF2 pairing gap in `matter`.
+
+    The plot is 2x2 with:
+
+    [0,0]: gap versus den.    [0,1]: gap versus kfn.
+
+    [1,0]: gap/EF versus den. [1,1]: gap/EF versus kfn.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
-    :param constraints: list of constraints to run on.
-    :type constraints: array of str.
-
+    :param models: models to run on.
+    :type models: array of str.
+    :param matter: can be 'SM' or 'NM' (default).
+    :type matter: str.
     """
     #
     print(f"Plot name: {pname}")

nucleardatapy/fig/matter_setupNEPStats_fig.py

@@ -6,9 +6,19 @@ import nucleardatapy as nuda
 
 def matter_setupNEPStats_fig(pname, models):
     """
-    Plot nucleonic energy per particle E/A in matter.\
-    The plot is 5x2 with:\
-    [0,0]: E/A versus den (micro). [0,1]: E/A versus den (pheno).\
+    Plot the PDF for NEPs.
+
+    The plot is 5x2 with:
+
+    [0,0]: Esat. [0,1]: Esym.
+
+    [1,0]: nsat. [1,1]: Lsym.
+
+    [2,0]: Ksat. [2,1]: Ksym.
+
+    [3,0]: Qsat. [3,1]: Qsym.
+
+    [4,0]: m*sat/m. [4,1]: Delta m*sat/m.
 
     :param pname: name of the figure (*.png)
     :type pname: str.

nucleardatapy/fig/matter_setupPhenoEsym_fig.py

@@ -6,26 +6,27 @@ import nucleardatapy as nuda
 
 def matter_setupPhenoEsym_fig(pname, models, band):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x2 with:\
-    [0]: nuclear chart.
+    Plot the symmetry energy esym for phenomenologic models.
+
+    The plot is 2x2 with:
+
+    [0,0]: esym function of the density.          [0,1]: esym function of the Fermi momentum.
+
+    [0,0]: esym/esym,FFG function of the density. [0,1]: esym/esym,FFG function of the Fermi momentum.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
-    :param table: table.
-    :type table: str.
-    :param version: version of table to run on.
-    :type version: str.
-    :param theo_tables: object instantiated on the reference band.
-    :type theo_tables: object.
-
+    :param models: list of phenomenological models.
+    :type models: array of str.
+    :param band: object instantiated on the reference band.
+    :type band: object.
     """
     #
     print(f"Plot name: {pname}")
     #
     fig, axs = plt.subplots(2, 2)
     fig.tight_layout()  # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust( left=0.12, bottom=0.12, right=0.95, top=0.9, wspace=0.05, hspace=0.05 )
+    fig.subplots_adjust( left=0.12, bottom=0.12, right=0.95, top=0.88, wspace=0.05, hspace=0.05 )
     #
     axs[0, 0].set_ylabel(r"$E_\mathrm{sym}$ (MeV)")
     axs[0, 0].set_xlim([0, 0.33])
@@ -134,88 +135,88 @@ def matter_setupPhenoEsym_fig(pname, models, band):
             #
         axs[0, 0].fill_between(
             band.den,
-            y1=(band.e2a - band.e2a_std),
-            y2=(band.e2a + band.e2a_std),
+            y1=(band.e2a_int - band.e2a_std),
+            y2=(band.e2a_int + band.e2a_std),
             color=band.color,
             alpha=band.alpha,
             visible=True,
         )
         axs[0, 0].plot(
             band.den,
-            (band.e2a - band.e2a_std),
+            (band.e2a_int - band.e2a_std),
             color="k",
             linestyle="dashed",
             visible=True, zorder = 100
         )
         axs[0, 0].plot(
             band.den,
-            (band.e2a + band.e2a_std),
+            (band.e2a_int + band.e2a_std),
             color="k",
             linestyle="dashed",
             visible=True, zorder = 100
         )
         axs[0, 1].fill_between(
             band.kfn,
-            y1=(band.e2a - band.e2a_std),
-            y2=(band.e2a + band.e2a_std),
+            y1=(band.e2a_int - band.e2a_std),
+            y2=(band.e2a_int + band.e2a_std),
             color=band.color,
             alpha=band.alpha,
             visible=True,
         )
         axs[0, 1].plot(
             band.kfn,
-            (band.e2a - band.e2a_std),
+            (band.e2a_int - band.e2a_std),
             color="k",
             linestyle="dashed",
             visible=True, zorder = 100
         )
         axs[0, 1].plot(
             band.kfn,
-            (band.e2a + band.e2a_std),
+            (band.e2a_int + band.e2a_std),
             color="k",
             linestyle="dashed",
             visible=True, zorder = 100
         )
         axs[1, 0].fill_between(
             band.den,
-            y1=(band.e2a - band.e2a_std) / nuda.esymffg_nr(band.kfn),
-            y2=(band.e2a + band.e2a_std) / nuda.esymffg_nr(band.kfn),
+            y1=(band.e2a_int - band.e2a_std) / nuda.esymffg_nr(band.kfn),
+            y2=(band.e2a_int + band.e2a_std) / nuda.esymffg_nr(band.kfn),
             color=band.color,
             alpha=band.alpha,
             visible=True,
         )
         axs[1, 0].plot(
             band.den,
-            (band.e2a - band.e2a_std) / nuda.esymffg_nr(band.kfn),
+            (band.e2a_int - band.e2a_std) / nuda.esymffg_nr(band.kfn),
             color="k",
             linestyle="dashed",
             visible=True, zorder = 100
         )
         axs[1, 0].plot(
             band.den,
-            (band.e2a + band.e2a_std) / nuda.esymffg_nr(band.kfn),
+            (band.e2a_int + band.e2a_std) / nuda.esymffg_nr(band.kfn),
             color="k",
             linestyle="dashed",
             visible=True, zorder = 100
         )
         axs[1, 1].fill_between(
             band.kf,
-            y1=(band.e2a - band.e2a_std) / nuda.esymffg_nr(band.kfn),
-            y2=(band.e2a + band.e2a_std) / nuda.esymffg_nr(band.kfn),
+            y1=(band.e2a_int - band.e2a_std) / nuda.esymffg_nr(band.kfn),
+            y2=(band.e2a_int + band.e2a_std) / nuda.esymffg_nr(band.kfn),
             color=band.color,
             alpha=band.alpha,
             visible=True,
         )
         axs[1, 1].plot(
             band.kfn,
-            (band.e2a - band.e2a_std) / nuda.esymffg_nr(band.kfn),
+            (band.e2a_int - band.e2a_std) / nuda.esymffg_nr(band.kfn),
             color="k",
             linestyle="dashed",
             visible=True, zorder = 100
         )
         axs[1, 1].plot(
             band.kfn,
-            (band.e2a + band.e2a_std) / nuda.esymffg_nr(band.kfn),
+            (band.e2a_int + band.e2a_std) / nuda.esymffg_nr(band.kfn),
             color="k",
             linestyle="dashed",
             visible=True, zorder = 100

nucleardatapy/fig/matter_setupPheno_fig.py

@@ -3,13 +3,15 @@ import matplotlib.pyplot as plt
 
 import nucleardatapy as nuda
 
-
-def matter_setupPheno_e2a_fig(pname, model, band):
+def matter_setupPheno_e2a_fig(pname, model, band, matter):
     """
-    Plot nucleonic energy per particle E/A in matter.\
-    The plot is 2x2 with:\
-    [0,0]: E/A versus den.       [0,1]: E/A versus kfn.\
-    [1,0]: E/E_NRFFG versus den. [1,1]: E/E_NRFFG versus kfn.\
+    Plot nucleonic internal energy per particle E/A in `matter`.
+
+    The plot is 2x2 with:
+
+    [0,0]: E/A versus den.       [0,1]: E/A versus kfn.
+
+    [1,0]: E/E_NRFFG versus den. [1,1]: E/E_NRFFG versus kfn.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
@@ -19,7 +21,6 @@ def matter_setupPheno_e2a_fig(pname, model, band):
     :type band: object.
     :param matter: can be 'SM' or 'NM'.
     :type matter: str.
-
     """
     #
     print(f"Plot name: {pname}")
@@ -218,12 +219,15 @@ def matter_setupPheno_e2a_fig(pname, model, band):
         plt.close()
     #
 
-def matter_setupPheno_pre_fig(pname, model, band):
+def matter_setupPheno_pre_fig(pname, model, band, matter):
     """
-    Plot nucleonic pressure in matter.\
-    The plot is 2x2 with:\
-    [0,0]: pre versus den.           [0,1]: pre versus kfn.\
-    [1,0]: pre/pre_NRFFG versus den. [1,1]: pre/pre_NRFFG versus kfn.\
+    Plot nucleonic pressure in matter.
+
+    The plot is 2x2 with:
+
+    [0,0]: pre versus den.           [0,1]: pre versus kfn.
+
+    [1,0]: pre/pre_NRFFG versus den. [1,1]: pre/pre_NRFFG versus kfn.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
@@ -329,11 +333,13 @@ def matter_setupPheno_pre_fig(pname, model, band):
         plt.close()
     #
 
-def matter_setupPheno_cs2_fig(pname, model, band):
+def matter_setupPheno_cs2_fig(pname, model, band, matter):
     """
-    Plot nucleonic sound speed in matter.\
-    The plot is 1x2 with:\
-    [0]: cs2 versus den.    [1]: cs2 versus kfn.\
+    Plot nucleonic sound speed in matter.
+
+    The plot is 1x2 with:
+
+    [0]: cs2 versus den.    [1]: cs2 versus kfn.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
@@ -343,7 +349,6 @@ def matter_setupPheno_cs2_fig(pname, model, band):
     :type band: object.
     :param matter: can be 'SM' or 'NM'.
     :type matter: str.
-
     """
     #
     print(f"Plot name: {pname}")

nucleardatapy/fig/nuc_setupBEExp_chart_fig.py

@@ -5,8 +5,10 @@ import nucleardatapy as nuda
 
 def nuc_setupBEExp_chart_lt_fig( pname, table, version, theo_tables ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
+    Plot nuclear chart (N versus Z).
+
+    The plot is 1x1 with:
+
     [0]: nuclear chart.
 
     :param pname: name of the figure (*.png)
@@ -17,7 +19,6 @@ def nuc_setupBEExp_chart_lt_fig( pname, table, version, theo_tables ):
     :type version: str.
     :param theo_tables: object instantiated on the reference band.
     :type theo_tables: object.
-
     """
     #
     print(f'Plot name: {pname}')
@@ -77,16 +78,20 @@ def nuc_setupBEExp_chart_lt_fig( pname, table, version, theo_tables ):
     legend = 0
     for i,theo_table in enumerate( theo_tables ):
         theo = nuda.nuc.setupBETheo( table = theo_table )
-        s2n = theo.S2n( Zmin=1, Zmax = 95 )
-        drip_S2n = s2n.drip_S2n( Zmin = 1, Zmax = 95 )
-        if legend == 0:
-            axs.scatter( drip_S2n.drip_S2n_N, drip_S2n.drip_S2n_Z, marker='o', s = 3, linewidth=0, color = 'purple', label='Drip Lines' )
-            legend = 1
-        else:
-            axs.scatter( drip_S2n.drip_S2n_N, drip_S2n.drip_S2n_Z, marker='o', s = 3, linewidth=0, color = 'purple' )
-        s2p = theo.S2p( Nmin=1, Nmax = 150 )
-        drip_S2p = s2p.drip_S2p( Nmin = 1, Nmax = 150 )
-        axs.scatter( drip_S2p.drip_S2p_N, drip_S2p.drip_S2p_Z, marker='o', s = 3, linewidth=0, color = 'purple' )
+        for zref in range(1,96):
+            itp = theo.isotopes( Zref= zref)
+            s2n = itp.S2n( Zref= zref )
+            drip_S2n = itp.drip_S2n(Zref = zref)
+            if legend == 0:
+                axs.scatter( drip_S2n.drip_S2n_N, drip_S2n.drip_S2n_Z, marker='o', s = 3, linewidth=0, color = 'purple', label='Drip Lines' )
+                legend = 1
+            else:
+                axs.scatter( drip_S2n.drip_S2n_N, drip_S2n.drip_S2n_Z, marker='o', s = 3, linewidth=0, color = 'purple' )
+        # for nref in range(1,170):
+        #     itn = theo.isotones( Nref= nref)
+        #     s2p = itn.S2p( Nref= nref)
+        # drip_S2p = theo.drip_S2p( Nmin=1, Nmax=169)
+        # axs.scatter( drip_S2p.drip_S2p_N, drip_S2p.drip_S2p_Z, marker='o', s = 3, linewidth=0, color = 'purple' )
     #
     # First and last isotopes
     #
@@ -132,19 +137,20 @@ def nuc_setupBEExp_chart_lt_fig( pname, table, version, theo_tables ):
 
 def nuc_setupBEExp_chart_year_fig( pname, sYear, year_min, year_max ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
+    Plot nuclear chart (N versus Z) for a range of discovery years.
+
+    The plot is 1x1 with:
+
     [0]: nuclear chart.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
-    :param table: table.
-    :type table: str.
-    :param version: version of table to run on.
-    :type version: str.
-    :param theo_tables: object instantiated on the reference band.
-    :type theo_tables: object.
-
+    :param sYear: select nuclei for given discovery years.
+    :type sYear: object.
+    :param year_min: lower range of the discovery years.
+    :type year_min: real.
+    :param year_max: upper range of the discovery years.
+    :type year_max: real.
     """
     #
     print(f'Plot name: {pname}')
@@ -204,8 +210,10 @@ def nuc_setupBEExp_chart_year_fig( pname, sYear, year_min, year_max ):
 
 def nuc_setupBEExp_chart_Rch_fig( pname, table, version, Rch_table ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
+    Plot nuclear chart (N versus Z).
+
+    The plot is 1x1 with:
+
     [0]: nuclear chart.
 
     :param pname: name of the figure (*.png)
@@ -214,8 +222,8 @@ def nuc_setupBEExp_chart_Rch_fig( pname, table, version, Rch_table ):
     :type table: str.
     :param version: version of table to run on.
     :type version: str.
-    :param theo_tables: object instantiated on the reference band.
-    :type theo_tables: object.
+    :param Rch_table: table for Rch.
+    :type Rch_table: str.
 
     """
     #
@@ -242,9 +250,10 @@ def nuc_setupBEExp_chart_Rch_fig( pname, table, version, Rch_table ):
     #
     mas = nuda.nuc.setupBEExp( table = table, version = version )
     ustbl = mas.select( state= 'gs', interp = 'n', nucleus = 'unstable' )
-    iso = ustbl.isotopes( Zmin=1, Zmax = 95 )
-    axs.scatter( iso.isotopes_Nmin, iso.isotopes_Z, marker='s', s = 3, linewidth=0, color = 'green', label='AME boundaries' )
-    axs.scatter( iso.isotopes_Nmax, iso.isotopes_Z, marker='s', s = 3, linewidth=0, color = 'green' )
+    for Zref in range(1,95):
+        iso = ustbl.isotopes( Zref = Zref )
+        axs.scatter( iso.itp_nucNmin, iso.itp_nucZ, marker='s', s = 3, linewidth=0, color = 'green' )
+        axs.scatter( iso.itp_nucNmax, iso.itp_nucZ, marker='s', s = 3, linewidth=0, color = 'green' )
     #
     # stable nuclei:
     #
@@ -266,8 +275,8 @@ def nuc_setupBEExp_chart_Rch_fig( pname, table, version, Rch_table ):
     #
     # plot stable_fit
     #
-    N, Z = nuda.stable_fit( Zmin = 1, Zmax = 120 )
-    axs.plot( N, Z, linestyle='dashed', linewidth=1, color='k')
+    #N, Z = nuda.stable_fit_N( Zmin = 1, Zmax = 120 )
+    #axs.plot( N, Z, linestyle='dashed', linewidth=1, color='k')
     #
     # plot shells for isotopes and isotones
     #

nucleardatapy/fig/nuc_setupBEExp_fig.py

@@ -5,9 +5,11 @@ import nucleardatapy as nuda
 
 def nuc_setupBEExp_year_fig( pname, table, version ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot the histogram for the discovery year.
+
+    The plot is 1x2 with:
+
+    [0]: full range of years. [1]: last two decades.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
@@ -15,9 +17,6 @@ def nuc_setupBEExp_year_fig( pname, table, version ):
     :type table: str.
     :param version: version of table to run on.
     :type version: str.
-    :param theo_tables: object instantiated on the reference band.
-    :type theo_tables: object.
-
     """
     #
     print(50*'-')
@@ -67,19 +66,20 @@ def nuc_setupBEExp_year_fig( pname, table, version ):
 
 def nuc_setupBEExp_S2n_fig( pname, tables, versions, Zref = 50 ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot S2n from the data extracted from the nuclear chart.
+
+    The plot is 1x1 with:
+
+    [0]: show S2n as a function of N.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
-    :param table: table.
-    :type table: str.
-    :param version: version of table to run on.
-    :type version: str.
-    :param theo_tables: object instantiated on the reference band.
-    :type theo_tables: object.
-
+    :param tables: tables.
+    :type tables: array of str.
+    :param versions: versions of the tables.
+    :type versions: array of str.
+    :param Zref: Reference value for Z.
+    :type Zref: integer.
     """
     #
     print(f'Plot name: {pname}')
@@ -120,19 +120,20 @@ def nuc_setupBEExp_S2n_fig( pname, tables, versions, Zref = 50 ):
 
 def nuc_setupBEExp_S2p_fig( pname, tables, versions, Nref = 50 ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot S2p from the data extracted from the nuclear chart.
+
+    The plot is 1x1 with:
+
+    [0]: show S2p as a function of Z.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
-    :param table: table.
-    :type table: str.
-    :param version: version of table to run on.
-    :type version: str.
-    :param theo_tables: object instantiated on the reference band.
-    :type theo_tables: object.
-
+    :param tables: tables.
+    :type tables: array of str.
+    :param versions: versions of the tables.
+    :type versions: array of str.
+    :param Nref: Reference value for N.
+    :type Nref: integer.
     """
     #
     print(f'Plot name: {pname}')
@@ -173,19 +174,20 @@ def nuc_setupBEExp_S2p_fig( pname, tables, versions, Nref = 50 ):
 
 def nuc_setupBEExp_D3n_fig( pname, tables, versions, Zref = 50 ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot D3n (3-point formula for the odd-even mass staggering) from the data extracted from the nuclear chart.
+
+    The plot is 1x1 with:
+
+    [0]: show D3n as a function of N.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
-    :param table: table.
-    :type table: str.
-    :param version: version of table to run on.
-    :type version: str.
-    :param theo_tables: object instantiated on the reference band.
-    :type theo_tables: object.
-
+    :param tables: tables.
+    :type tables: array of str.
+    :param versions: versions of the tables.
+    :type versions: array of str.
+    :param Zref: Reference value for Z.
+    :type Zref: integer.
     """
     #
     print(f'Plot name: {pname}')
@@ -229,19 +231,20 @@ def nuc_setupBEExp_D3n_fig( pname, tables, versions, Zref = 50 ):
 
 def nuc_setupBEExp_D3p_fig( pname, tables, versions, Nref = 50 ):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x1 with:\
-    [0]: nuclear chart.
+    Plot D3p (3-point formula for the odd-even mass staggering) from the data extracted from the nuclear chart.
+
+    The plot is 1x1 with:
+
+    [0]: show D3p as a function of Z.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
-    :param table: table.
-    :type table: str.
-    :param version: version of table to run on.
-    :type version: str.
-    :param theo_tables: object instantiated on the reference band.
-    :type theo_tables: object.
-
+    :param tables: tables.
+    :type tables: array of str.
+    :param versions: versions of the tables.
+    :type versions: array of str.
+    :param Nref: Reference value for N.
+    :type Nref: integer.
     """
     #
     print(f'Plot name: {pname}')