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

nucleardatapy/fig/matter_all_fig.py

@@ -5,9 +5,11 @@ import nucleardatapy as nuda
 
 def matter_all_e2a_fig( pname, micro_mbs, pheno_models, band_check, band_plot, matter ):
     """
-    Plot nucleonic energy per particle E/A in matter.\
-    The plot is 1x2 with:\
-    [0,0]: E/A versus den (micro). [0,1]: E/A versus den (pheno).\
+    Plot nucleonic energy per particle E/A in matter.
+
+    The plot is 1x2 with:
+
+    [0]: E/A versus den (micro). [1]: E/A versus den (pheno).
 
     :param pname: name of the figure (*.png)
     :type pname: str.
@@ -15,9 +17,12 @@ def matter_all_e2a_fig( pname, micro_mbs, pheno_models, band_check, band_plot, m
     :type micro_mbs: str.
     :param pheno_models: models to run on.
     :type pheno_models: array of str.
-    :param band: object instantiated on the reference band.
-    :type band: object.
-
+    :param band_check: object instantiated on the reference band.
+    :type band_check: object.
+    :param band_plot: object instantiated on the reference band.
+    :type band_plot: object.
+    :param matter: can be 'SM' or 'NM'.
+    :type matter: str.
     """
     #
     print(f'Plot name: {pname}')
@@ -51,6 +56,7 @@ def matter_all_e2a_fig( pname, micro_mbs, pheno_models, band_check, band_plot, m
     for kmb,mb in enumerate(micro_mbs):
         #
         models, models_lower = nuda.matter.micro_models_mb( mb )
+        print('models in matter code:',models)
         #
         for model in models:
             #
@@ -202,27 +208,27 @@ def matter_all_e2a_fig( pname, micro_mbs, pheno_models, band_check, band_plot, m
     if pname is not None:
     	plt.savefig(pname, dpi=200)
     #
-    	plt.close()
+    plt.close()
     #
 
 def matter_all_Esym_fig( pname, micro_mbs, pheno_models, band_check, band_plot ):
     """
-    Plot nucleonic symmetry energy.\
-    The plot is 2x2 with:\
-    [0,0]: Esym versus den.            [0,1]: Esym versus kfn.\
-    [1,0]: Esym/Esym_NRFFG versus den. [1,1]: Esym/Esym_NRFFG versus kfn.\
+    Plot nucleonic symmetry energy.
+
+    The plot is 1x2 with:
+
+    [0]: Esym versus den (micro).  [1]: Esym versus den (pheno).
 
     :param pname: name of the figure (*.png)
     :type pname: str.
-    :param mb: many-body (mb) approach considered.
-    :type mb: str.
-    :param models: models to run on.
-    :type models: array of str.
-    :param band: object instantiated on the reference band.
-    :type band: object.
-    :param matter: can be 'SM' or 'NM'.
-    :type matter: str.
-
+    :param micro_mbs: many-body (mb) approach considered.
+    :type micro_mbs: str.
+    :param pheno_models: models to run on.
+    :type pheno_models: array of str.
+    :param band_check: object instantiated on the reference band.
+    :type band_check: object.
+    :param band_plot: object instantiated on the reference band.
+    :type band_plot: object.
     """
     #
     print(f'Plot name: {pname}')
@@ -340,9 +346,11 @@ def matter_all_Esym_fig( pname, micro_mbs, pheno_models, band_check, band_plot )
 
 def matter_all_pre_fig( pname, micro_mbs, pheno_models, band_check, matter ):
     """
-    Plot nucleonic pressure in matter.\
-    The plot is 1x2 with:\
-    [0]: pre versus den (micro). [1]: pre versus den (pheno).\
+    Plot nucleonic pressure in matter.
+
+    The plot is 1x2 with:
+
+    [0]: pre versus den (micro). [1]: pre versus den (pheno).
 
     :param pname: name of the figure (*.png)
     :type pname: str.
@@ -350,9 +358,12 @@ def matter_all_pre_fig( pname, micro_mbs, pheno_models, band_check, matter ):
     :type micro_mbs: str.
     :param pheno_models: models to run on.
     :type pheno_models: array of str.
-    :param band: object instantiated on the reference band.
-    :type band: object.
-
+    :param band_check: object instantiated on the reference band.
+    :type band_check: object.
+    :param band_plot: object instantiated on the reference band.
+    :type band_plot: object.
+    :param matter: can be 'SM' or 'NM'.
+    :type matter: str.
     """
     #
     print(f'Plot name: {pname}')
@@ -551,9 +562,11 @@ def matter_all_pre_fig( pname, micro_mbs, pheno_models, band_check, matter ):
 
 def matter_all_eos_fig( pname, micro_mbs, pheno_models, band_check, matter ):
     """
-    Plot EoS in matter.\
-    The plot is 1x2 with:\
-    [0]: pre versus energy-density (micro). [1]: pre versus energy-density (pheno).\
+    Plot EoS in matter.
+
+    The plot is 1x2 with:
+
+    [0]: pre versus eps (micro). [1]: pre versus eps (pheno).
 
     :param pname: name of the figure (*.png)
     :type pname: str.
@@ -561,9 +574,10 @@ def matter_all_eos_fig( pname, micro_mbs, pheno_models, band_check, matter ):
     :type micro_mbs: str.
     :param pheno_models: models to run on.
     :type pheno_models: array of str.
-    :param band: object instantiated on the reference band.
-    :type band: object.
-
+    :param band_check: object instantiated on the reference band.
+    :type band_check: object.
+    :param matter: can be 'SM' or 'NM'.
+    :type matter: str.
     """
     #
     print(f'Plot name: {pname}')
@@ -762,9 +776,11 @@ def matter_all_eos_fig( pname, micro_mbs, pheno_models, band_check, matter ):
 
 def matter_all_cs2_fig( pname, micro_mbs, pheno_models, band_check, matter ):
     """
-    Plot nucleonic sound speed in matter.\
-    The plot is 1x2 with:\
-    [0]: cs2 versus den (micro). [1]: cs2 versus den (pheno).\
+    Plot nucleonic sound speed in matter.
+
+    The plot is 1x2 with:
+
+    [0]: cs2 versus den (micro). [1]: cs2 versus den (pheno).
 
     :param pname: name of the figure (*.png)
     :type pname: str.
@@ -772,9 +788,10 @@ def matter_all_cs2_fig( pname, micro_mbs, pheno_models, band_check, matter ):
     :type micro_mbs: str.
     :param pheno_models: models to run on.
     :type pheno_models: array of str.
-    :param band: object instantiated on the reference band.
-    :type band: object.
-
+    :param band_check: object instantiated on the reference band.
+    :type band_check: object.
+    :param matter: can be 'SM' or 'NM'.
+    :type matter: str.
     """
     #
     print(f'Plot name: {pname}')

nucleardatapy/fig/matter_setupCheck_fig.py

@@ -3,10 +3,13 @@ import matplotlib.pyplot as plt
 
 import nucleardatapy as nuda
 
-def matter_setupCheck_fig( pname, mb, models, band ):
+def matter_setupCheck_fig( pname, mb, models, band, matter ):
     """
-    Plot nucleonic energy per particle E/A in matter.\
-    The plot is 1x1 with: E/A versus den.
+    Check E/A from models in `models` and show the reference band.
+
+    The plot is 1x1 with: 
+    
+    [0]: E/A versus den.
 
     :param pname: name of the figure (*.png)
     :type pname: str.

nucleardatapy/fig/matter_setupFFGLep_fig.py

@@ -7,18 +7,22 @@ mpl.use("Agg")  # Use a non-interactive backend
 
 def matter_setupFFGLep_fig( pname, den_el=None, den_mu1=None, den_mu2=None, den_mu3=None ):
     """
-    Plot leptonic FFG energy per particle E/A and pressure in NM and SM.\
-    The plot is 2x2 with:\
-    [0,0]: E/A versus den.     [0,1]: E/A versus kfn.\
-    [1,0]: pre versus den.     [1,1]: pre versus kfn.\
+    Plot leptonic FFG energy per particle E/A and pressure in NM and SM.
+
+    The plot is 2x1 with:
+
+    [0]: E/A versus den.    [1]: pre versus den.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
-    :param den: density.
-    :type den: float or numpy vector of real numbers.
-    :param kfn: neutron Fermi momentum.
-    :type kfn: float or numpy vector of real numbers.
-
+    :param den_el: electron densities.
+    :type den_el: float or numpy vector of real numbers.
+    :param den_mu1: muon densities (set 1).
+    :type den_mu1: float or numpy vector of real numbers.
+    :param den_mu2: muon densities (set 2).
+    :type den_mu2: float or numpy vector of real numbers.
+    :param den_mu3: muon densities (set 3).
+    :type den_mu3: float or numpy vector of real numbers.
     """
     #
     print(f"Plot name: {pname}")

nucleardatapy/fig/matter_setupFFGNuc_fig.py

@@ -5,21 +5,24 @@ import nucleardatapy as nuda
 
 mpl.use("Agg")  # Use a non-interactive backend
 
-def matter_setupFFGNuc_EP_fig(
-    pname, mss=[1.0], den=np.linspace(0.01, 0.35, 10), kf=np.linspace(0.5, 2.0, 10)
-):
+def matter_setupFFGNuc_EP_fig( pname, mss=[1.0], den=np.linspace(0.01, 0.35, 10), kf=np.linspace(0.5, 2.0, 10) ):
     """
-    Plot nucleonic FFG energy per particle E/A and pressure in NM and SM.\
-    The plot is 2x2 with:\
-    [0,0]: E/A versus den.     [0,1]: E/A versus kfn.\
-    [1,0]: pre versus den.     [1,1]: pre versus kfn.\
+    Plot nucleonic FFG energy per particle and pressure in NM and SM.
+
+    The plot is 2x2 with:
+
+    [0,0]: E/A versus den.     [0,1]: E/A versus kfn.
+
+    [1,0]: pre versus den.     [1,1]: pre versus kfn.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
-    :param den: density.
+    :param mss: effective mass/bare mass of the nucleons.
+    :type mss: float or numpy vector of real numbers.
+    :param den: densities.
     :type den: float or numpy vector of real numbers.
-    :param kfn: neutron Fermi momentum.
-    :type kfn: float or numpy vector of real numbers.
+    :param kf: Fermi momenta.
+    :type kf: float or numpy vector of real numbers.
 
     """
     #
@@ -204,13 +207,18 @@ def matter_setupFFGNuc_EP_fig(
 
 def matter_setupFFGNuc_EOS_fig(pname, mss=[1.0], den=np.linspace(0.01, 0.35, 10)):
     """
-    Plot nucleonic FFG EOS in NM and SM.\
-    The plot is 1x2 with:\
-    [0]: EOS (pre) versus energy density rho.\
-    [1]: Sound speed c_s^2 versus energy density rho.\
+    Plot nucleonic FFG EOS in NM and SM.
+
+    The plot is 2x1 with:
+
+    [0]: EOS (pre) versus energy density rho.
+
+    [1]: Sound speed c_s^2 versus energy density rho.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
+    :param mss: effective mass/bare mass of the nucleons.
+    :type mss: float or numpy vector of real numbers.
     :param den: density.
     :type den: float or numpy vector of real numbers.
 

nucleardatapy/fig/matter_setupHIC_fig.py

@@ -14,18 +14,18 @@ import nucleardatapy as nuda
 
 def matter_setupHIC_fig(pname, inferences):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x2 with:\
-    [0]: nuclear chart.
+    Plot the inferences from HIC.
+
+    The plot is 2x2 with:
+
+    [0,0]: pressure in SM versus den. [0,1]: E/A in SM versus den.
+    
+    [1,0]: pressure in NM versus den. [1,1]: Esym versus den.
 
     :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 inferences: inferences for HIC.
+    :type inferences: array of str.
 
     """
     #

nucleardatapy/fig/matter_setupMicroEsym_fig.py

@@ -6,19 +6,20 @@ import nucleardatapy as nuda
 
 def matter_setupMicroEsym_fig(pname, mbs, band):
     """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x2 with:\
-    [0]: nuclear chart.
+    Plot the symmetry energy esym for microscopic 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 mbs: list of many-body approaches.
+    :type mbs: array of str.
+    :param band: object instantiated on the reference band.
+    :type band: object.
     """
     #
     print(f"Plot name: {pname}")
@@ -40,18 +41,18 @@ def matter_setupMicroEsym_fig(pname, mbs, band):
     axs[1, 0].set_ylabel(r"$E_\mathrm{sym}/E_\mathrm{sym, FFG, NR}$", fontsize="14")
     axs[1, 0].set_xlabel(r"$n_\mathrm{nuc}$ (fm$^{-3}$)", fontsize="14")
     axs[1, 0].set_xlim([0,   0.33])
-    axs[1, 0].set_ylim([1.5, 3.6])
+    axs[1, 0].set_ylim([1.8, 2.9])
     #
     axs[1, 1].set_xlabel(r"$k_F$ (fm$^{-1}$)", fontsize="14")
     axs[1, 1].set_xlim([0.5, 2.0])
-    axs[1, 1].set_ylim([1.5, 3.6])
+    axs[1, 1].set_ylim([1.8, 2.9])
     axs[1, 1].tick_params("y", labelleft=False)
     #
     mb_check = []
     #
     for kmb, mb in enumerate(mbs):
         #
-        models, models_lower = nuda.matter.micro_models_mb(mb)
+        models, models_lower = nuda.matter.micro_esym_models_mb(mb)
         #
         for model in models:
             #
@@ -234,55 +235,55 @@ def matter_setupMicroEsym_fig(pname, mbs, 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), color="k", linestyle="dashed", visible=True, zorder = 100 )
+        band.den, (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), color="k", linestyle="dashed", visible=True, zorder = 100 )
+        band.den, (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), color="k", linestyle="dashed", visible=True, zorder = 100 )
+        band.kfn, (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), color="k", linestyle="dashed", visible=True, zorder = 100 )
+        band.kfn, (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.kf),
-        y2=(band.e2a + band.e2a_std) / nuda.esymffg_nr(band.kf),
+        y1=(band.e2a_int - band.e2a_std) / nuda.esymffg_nr(band.kf),
+        y2=(band.e2a_int + band.e2a_std) / nuda.esymffg_nr(band.kf),
         color=band.color,
         alpha=band.alpha,
         visible=True,
     )
     axs[1, 0].plot(
-        band.den, (band.e2a - band.e2a_std) / nuda.esymffg_nr(band.kf),
+        band.den, (band.e2a_int - band.e2a_std) / nuda.esymffg_nr(band.kf),
         color="k", linestyle="dashed", visible=True, zorder = 100 )
     axs[1, 0].plot(
-        band.den, (band.e2a + band.e2a_std) / nuda.esymffg_nr(band.kf),
+        band.den, (band.e2a_int + band.e2a_std) / nuda.esymffg_nr(band.kf),
         color="k", linestyle="dashed", visible=True, zorder = 100 )
     axs[1, 1].fill_between(
         band.kfn,
-        y1=(band.e2a - band.e2a_std) / nuda.esymffg_nr(band.kf),
-        y2=(band.e2a + band.e2a_std) / nuda.esymffg_nr(band.kf),
+        y1=(band.e2a_int - band.e2a_std) / nuda.esymffg_nr(band.kf),
+        y2=(band.e2a_int + band.e2a_std) / nuda.esymffg_nr(band.kf),
         color=band.color,
         alpha=band.alpha,
         visible=True,
     )
     axs[1, 1].plot(
-        band.kfn, (band.e2a - band.e2a_std) / nuda.esymffg_nr(band.kf),
+        band.kfn, (band.e2a_int - band.e2a_std) / nuda.esymffg_nr(band.kf),
         color="k", linestyle="dashed", visible=True, zorder = 100 )
     axs[1, 1].plot(
-        band.kfn, (band.e2a + band.e2a_std) / nuda.esymffg_nr(band.kf),
+        band.kfn, (band.e2a_int + band.e2a_std) / nuda.esymffg_nr(band.kf),
         color="k", linestyle="dashed", visible=True, zorder = 100 )
 
     # axs[1,0].legend(loc='upper right',fontsize='8')
@@ -291,7 +292,7 @@ def matter_setupMicroEsym_fig(pname, mbs, band):
         bbox_to_anchor=(0.1, 1.0),
         columnspacing=2,
         fontsize="8",
-        ncol=4,
+        ncol=6,
         frameon=False,
     )
     #

nucleardatapy/fig/matter_setupMicro_LP_fig.py

@@ -3,13 +3,19 @@ import matplotlib.pyplot as plt
 
 import nucleardatapy as nuda
 
-
 def matter_setupMicro_LP_fig( pname, models, matter="SM", ell=0 ):
     """
-    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 energy per particle E/A in matter.
+
+    The plot is 2x2 in SM with:
+
+    [0,0]: F_ell versus den.    [0,1]: G_ell versus den.
+
+    [1,0]: F'_ell' versus den.  [1,1]: G'_ell versus den.
+
+    Or the plot is 2x1 in NM with:
+
+    [0]: F_ell versus den.    [1]: G_ell versus den.
 
     :param pname: name of the figure (*.png)
     :type pname: str.

nucleardatapy/fig/matter_setupMicro_band_fig.py

@@ -6,15 +6,20 @@ import nucleardatapy as nuda
 
 def matter_setupMicro_band_fig(pname, models, den, matter):
     """
-    Plot the correlation between Esym and Lsym.\
-    The plot is 1x1 with:\
-    [0]: Esym - Lsym correlation plot
+    Plot the reference band together with the models on wich it is constructed.
+
+    The plot is 1x2 with:
+
+    [0]: E versus den.  [1]: E/EFFG versus den.
 
     :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: list of models to run on.
+    :type models: array of str.
+    :param den: numpy array with densities to define the reference band.
+    :type den: array of reals.
+    :param matter: chose between 'SM', 'NM' or 'Esym'.
+    :type matter: str.
     """
     #
     print(f"Plot name: {pname}")

nucleardatapy/fig/matter_setupMicro_effmass_fig.py

@@ -6,10 +6,11 @@ import nucleardatapy as nuda
 
 def matter_setupMicro_effmass_fig(pname, models, matter="NM"):
     """
-    Plot the effective mass as function of the density and Fermi momentum.\
-    The plot is 1x2 with:\
-    [0]: effmass(den)
-    [1]: effmass(kF)
+    Plot the effective mass as function of the density and Fermi momentum.
+
+    The plot is 1x2 with:
+
+    [0]: effmass(den).  [1]: effmass(kF).
 
     :param pname: name of the figure (*.png)
     :type pname: str.

nucleardatapy/fig/matter_setupMicro_err_NM_fig.py

@@ -6,8 +6,10 @@ import nucleardatapy as nuda
 
 def matter_setupMicro_err_NM_fig(pname, models):
     """
-    Plot uncertainties (err) estimated by different authors in NM.\
-    The plot is 1x1 with:\
+    Plot uncertainties (err) estimated by different authors in NM.
+
+    The plot is 1x1 with:
+
     [0]: uncertainty for E/A versus den.
 
     :param pname: name of the figure (*.png)

nucleardatapy/fig/matter_setupMicro_fig.py

@@ -3,12 +3,15 @@ import matplotlib.pyplot as plt
 
 import nucleardatapy as nuda
 
-def matter_setupMicro_e2a_fig( pname, mb, models, band ):
+def matter_setupMicro_e2a_fig( pname, mb, models, 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 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.
@@ -36,10 +39,10 @@ def matter_setupMicro_e2a_fig( pname, mb, models, band ):
     axs[1,1].tick_params('y', labelleft=False)
     axs[0,0].tick_params('x', labelbottom=False)
     axs[0,1].tick_params('x', labelbottom=False)
-    axs[0,0].set_xlim([0, 0.33])
-    axs[1,0].set_xlim([0, 0.33])
-    axs[0,1].set_xlim([0, 1.9])
-    axs[1,1].set_xlim([0, 1.9])
+    axs[0,0].set_xlim([0,0.33]) #([0.06, 0.28]) #[0, 0.33]
+    axs[1,0].set_xlim([0,0.33]) #([0.06, 0.28]) #[0, 0.33]
+    axs[0,1].set_xlim([0,1.9]) #([0.8, 1.6]) #[0, 1.9]
+    axs[1,1].set_xlim([0,1.9]) #([0.8, 1.6]) #[0, 1.9]
     if matter.lower() == 'nm':
         axs[0,0].set_ylabel(r'$E_\text{int,NM}/A$ (MeV)', fontsize = '14' )
         axs[1,0].set_ylabel(r'$E_\text{int,NM}/E_\text{int,NM}^\text{NRFFG}$', fontsize = '14' )
@@ -51,10 +54,10 @@ def matter_setupMicro_e2a_fig( pname, mb, models, band ):
     elif matter.lower() == 'sm':
         axs[0,0].set_ylabel(r'$E_\text{int,SM}/A$ (MeV)', fontsize = '14' )
         axs[1,0].set_ylabel(r'$E_\text{int,SM}/E_\text{int,SM}^\text{NRFFG}$', fontsize = '14' )
-        axs[0,0].set_ylim([-20, 10])
-        axs[0,1].set_ylim([-20, 10])
-        axs[1,0].set_ylim([-1.2, 0.5])
-        axs[1,1].set_ylim([-1.2, 0.5])
+        axs[0,0].set_ylim([-20,10]) #([-18, -8]) #[-20, 10]
+        axs[0,1].set_ylim([-20,10]) #([-18, -8]) #[-20, 10]
+        axs[1,0].set_ylim([-1.2, 0.5]) #([-1.0, -0.5]) #[-1.2, 0.5]
+        axs[1,1].set_ylim([-1.2, 0.5]) #([-1.0, -0.5]) #[-1.2, 0.5]
         delta = 0.0
     #
     for model in models:
@@ -166,12 +169,15 @@ def matter_setupMicro_e2a_fig( pname, mb, models, band ):
         plt.savefig(pname, dpi=300)
         plt.close()
 
-def matter_setupMicro_pre_fig( pname, mb, models, band ):
+def matter_setupMicro_pre_fig( pname, mb, models, 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.
@@ -295,11 +301,13 @@ def matter_setupMicro_pre_fig( pname, mb, models, band ):
         plt.savefig(pname, dpi=300)
         plt.close()
 
-def matter_setupMicro_cs2_fig( pname, mb, models, band ):
+def matter_setupMicro_cs2_fig( pname, mb, models, band, matter ):
     """
-    Plot nucleonic pressure in matter.\
-    The plot is 1x2 with:\
-    [0]: cs2 versus den.   [1]: cs2 versus kfn.\
+    Plot nucleonic pressure 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.