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

nucleardatapy/eos/setup_am_Leq.py

@@ -87,11 +87,11 @@ class setupAMLeq():
         self.marker = esym.marker
         #print('type esym:',type(esym.esym))
         if esym.esym is not None:
-            self.nm_den = esym.nm_den
-            self.nm_e2a = esym.nm_e2a
-            self.sm_den = esym.sm_den
-            self.sm_e2a = esym.sm_e2a
             self.den = esym.den
+            self.nm_den = esym.den
+            self.sm_den = esym.den
+            self.nm_e2a_int = esym.esym_nm_e2a_int
+            self.sm_e2a_int = esym.esym_sm_e2a_int
             self.esym = esym.esym
             #print('esym:',self.esym)
             self.x_n = ( 1.0 + self.asy ) / 2.0
@@ -120,29 +120,54 @@ class setupAMLeq():
             self.n_mu = self.x_mu * self.den
             #
             # Thermodynamical variables
-            self.e2a_nuc = esym.esym_sm_e2a + esym.esym * self.asy**2
-            self.e2v_nuc = self.e2a_nuc * self.den
+            # nucleons
+            self.rmass = self.x_n * nuda.cst.mnc2 + self.x_p * nuda.cst.mpc2
+            self.e2a_int_nuc = esym.esym_sm_e2a_int + esym.esym * self.asy**2
+            self.e2a_nuc = self.rmass + self.e2a_int_nuc
+            self.eps_int_nuc = self.e2a_int_nuc * self.den
+            self.eps_nuc = self.e2a_nuc * self.den
             self.pre_nuc = esym.esym_sm_pre + esym.esym_sym_pre * self.asy**2
+            # leptons
             lep = nuda.matter.setupFFGLep( den_el = self.n_el, den_mu = self.n_mu )
             self.e2a_el = self.x_el * lep.e2n_el
+            self.e2a_int_el = self.e2a_el - self.x_el * nuda.cst.mec2
             self.e2a_mu = self.x_mu * lep.e2n_mu
+            self.e2a_int_mu = self.e2a_mu - self.x_mu * nuda.cst.mmuc2
             self.e2a_lep = self.x_lep * lep.e2n_lep
-            self.e2v_lep = self.e2a_lep / self.den
+            self.eps_lep = self.e2a_lep / self.den
             self.pre_el = lep.pre_el
             self.pre_mu = lep.pre_mu
             self.pre_lep = lep.pre_el + lep.pre_mu
-            # self.h2a
-            # self.h2v
-            # self.cs2
             # total
+            self.e2a_int_tot = self.e2a_int_nuc + self.e2a_lep
             self.e2a_tot = self.e2a_nuc + self.e2a_lep
-            self.e2v_tot = self.e2v_nuc + self.e2v_lep
+            self.eps_int_tot = self.eps_int_nuc + self.eps_lep
+            self.eps_tot = self.eps_nuc + self.eps_lep
             self.pre_tot = self.pre_nuc + self.pre_lep
-            #
+            # enthalpy self.h2a
+            self.h2a_lep = self.e2a_lep + self.pre_lep / self.den
+            self.h2a_nuc = self.e2a_nuc + self.pre_nuc / self.den
+            self.h2a_tot = self.e2a_tot + self.pre_tot / self.den
+            # enthaply density self.h2v
+            self.h2v_lep = self.h2a_lep * self.den
+            self.h2v_nuc = self.h2a_nuc * self.den
+            self.h2v_tot = self.h2a_tot * self.den
+            # sound speed self.cs2
+            x = np.insert(self.den, 0, 0.0)
+            y = np.insert(self.pre_lep, 0, 0.0)
+            cs_pre = CubicSpline(x, y)
+            self.cs2_lep = cs_pre(self.den, 1) / self.h2a_lep
+            y = np.insert(self.pre_nuc, 0, 0.0)
+            cs_pre = CubicSpline(x, y)
+            self.cs2_nuc = cs_pre(self.den, 1) / self.h2a_nuc
+            y = np.insert(self.pre_tot, 0, 0.0)
+            cs_pre = CubicSpline(x, y)
+            self.cs2_tot = cs_pre(self.den, 1) / self.h2a_tot
+            #            #
         self.den_unit = 'fm$^{-3}$'
         self.kf_unit = 'fm$^{-1}$'
         self.e2a_unit = 'MeV'
-        self.e2v_unit = 'MeV fm$^{-3}$'
+        self.eps_unit = 'MeV fm$^{-3}$'
         self.pre_unit = 'MeV fm$^{-3}$'
         self.gap_unit = 'MeV'
         #
@@ -165,7 +190,7 @@ class setupAMLeq():
         if self.kfn is not None: print(f"   kfn: {np.round(self.den,3)} in {self.kf_unit}")
         if self.asy is not None: print(f"   asy: {np.round(self.asy,3)}")
         if self.e2a is not None: print(f"   e2a: {np.round(self.e2a,3)} in {self.e2a_unit}")
-        if self.e2v is not None: print(f"   e2v: {np.round(self.e2v,3)} in {self.e2v_unit}")
+        if self.eps is not None: print(f"   eps: {np.round(self.eps,3)} in {self.eps_unit}")
         if self.pre is not None: print(f"   pre: {np.round(self.pre,3)} in {self.pre_unit}")
         if self.cs2 is not None: print(f"   cs2: {np.round(self.cs2,2)}")
         #
@@ -202,7 +227,7 @@ class setupAMLeq():
         #: Attribute the energy per particle.
         self.e2a = None
         #: Attribute the energy per unit volume.
-        self.e2v = None
+        self.eps = None
         #: Attribute the enthalpy per particle.
         self.h2a = None
         #: Attribute the enthalpy per unit volume.

nucleardatapy/fig/__init__.py

@@ -2,7 +2,8 @@
 This module provides microscopic, phenomenological and experimental data constraints.
 """
 from nucleardatapy.fig.matter_setupFFGNuc_fig        import *
-from nucleardatapy.fig.matter_setupMicro_fig import *
+from nucleardatapy.fig.matter_setupFFGLep_fig        import *
+from nucleardatapy.fig.matter_setupMicro_fig         import *
 from nucleardatapy.fig.matter_setupMicro_err_NM_fig  import *
 from nucleardatapy.fig.matter_setupMicro_band_fig    import *
 from nucleardatapy.fig.matter_setupMicro_effmass_fig import *
@@ -10,19 +11,26 @@ from nucleardatapy.fig.matter_setupMicro_gap_fig     import *
 from nucleardatapy.fig.matter_setupMicro_LP_fig      import *
 from nucleardatapy.fig.matter_setupMicroEsym_fig     import *
 #
-from nucleardatapy.fig.matter_setupPheno_fig        import *
-from nucleardatapy.fig.matter_setupPhenoEsym_fig    import *
+from nucleardatapy.fig.matter_setupPheno_fig         import *
+from nucleardatapy.fig.matter_setupPhenoEsym_fig     import *
 #
-from nucleardatapy.fig.matter_setupHIC_fig          import *
+from nucleardatapy.fig.matter_setupHIC_fig           import *
 #
-from nucleardatapy.fig.matter_setupNEPModelDist_fig import *
+from nucleardatapy.fig.matter_setupNEPStats_fig      import *
 #
-from nucleardatapy.fig.matter_Esym_fig              import *
-from nucleardatapy.fig.matter_ESM_fig               import *
-from nucleardatapy.fig.matter_ENM_fig               import *
-from nucleardatapy.fig.matter_cs2_fig               import *
-from nucleardatapy.fig.matter_preNM_fig             import *
-from nucleardatapy.fig.matter_preSM_fig             import *
+from nucleardatapy.fig.matter_setupCheck_fig         import *
+#
+from nucleardatapy.fig.matter_all_fig                import *
+#
+#from nucleardatapy.fig.matter_Esym_fig              import *
+#from nucleardatapy.fig.matter_ESM_fig               import *
+#from nucleardatapy.fig.matter_ENM_fig               import *
+#from nucleardatapy.fig.matter_cs2_fig               import *
+#from nucleardatapy.fig.matter_preNM_fig             import *
+#from nucleardatapy.fig.matter_preSM_fig             import *
+#from nucleardatapy.fig.matter_eosNM_fig             import *
+#from nucleardatapy.fig.matter_eosSM_fig             import *
+#from nucleardatapy.fig.matter_cs2NM_fig             import *
 #
 from nucleardatapy.fig.corr_setupKsatQsat_fig       import *
 from nucleardatapy.fig.corr_setupEsymLsym_fig       import *
@@ -40,12 +48,13 @@ from nucleardatapy.fig.nuc_setupRnpTheo_fig         import *
 from nucleardatapy.fig.hnuc_setupRE1LExp_fig        import *
 from nucleardatapy.fig.hnuc_setupChart_fig          import *
 #
-from nucleardatapy.fig.eos_setupAM_e_fig            import *
-from nucleardatapy.fig.eos_setupAM_e_asy_nuc_fig    import *
-from nucleardatapy.fig.eos_setupAM_e_asy_lep_fig    import *
-from nucleardatapy.fig.eos_setupAM_e_asy_tot_fig    import *
+from nucleardatapy.fig.eos_setupAM_fig              import *
+from nucleardatapy.fig.eos_setupAM_asy_nuc_fig      import *
+from nucleardatapy.fig.eos_setupAM_asy_lep_fig      import *
+from nucleardatapy.fig.eos_setupAM_asy_tot_fig      import *
 from nucleardatapy.fig.eos_setupAMBeq_fig           import *
 from nucleardatapy.fig.eos_setupAMLeq_fig           import *
+from nucleardatapy.fig.eos_setupCC_fig              import *
 #
 from nucleardatapy.fig.crust_setupCrust_fig         import *
 #

nucleardatapy/fig/astro_setupGW_fig.py

@@ -5,8 +5,10 @@ import nucleardatapy as nuda
 
 def astro_setupGW_fig( pname, sources ):
     """
-    Plot Tidal deformabilities for each sources.\
-    The plot is 1x1 with:\
+    Plot Tidal deformabilities as a function of sources.
+
+    The plot is 1x1 with:
+    
     [0]: Tidal deformabilities versus sources.
 
     :param pname: name of the figure (*.png)
@@ -22,10 +24,10 @@ def astro_setupGW_fig( pname, sources ):
     #
     fig, axs = plt.subplots(1,1)
     fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.14, bottom=0.12, right=None, top=0.8, wspace=0.3, hspace=0.3)
+    fig.subplots_adjust(left=0.14, bottom=0.12, right=0.95, top=0.85, wspace=0.3, hspace=0.3)
     #
-    axs.set_xlabel(r'sources',fontsize='12')
-    axs.set_ylabel(r'$\tilde{\Lambda}_{90\%}$',fontsize='12')
+    axs.set_xlabel(r'sources',fontsize='14')
+    axs.set_ylabel(r'$\tilde{\Lambda}_{90\%}$',fontsize='14')
     axs.set_xlim([0.8, 2.5])
     axs.set_ylim([0, 1200])
     #
@@ -46,7 +48,7 @@ def astro_setupGW_fig( pname, sources ):
             gw = nuda.astro.setupGW( source = source, hyp = hyp )
             if nuda.env.verb_output: gw.print_output( )
             if nuda.env.verb_latex: gw.print_latex( )
-            axs.errorbar( isource+ihyp/10, gw.lam, yerr=np.array([(gw.lam_sig_do,gw.lam_sig_up)]).T, label=gw.label, color=nuda.param.col[isource], marker=gw.marker, linestyle = 'solid', linewidth = 1 )
+            axs.errorbar( isource+ihyp/10, gw.lam, yerr=np.array([(gw.lam_sig_lo,gw.lam_sig_up)]).T, label=gw.label, color=nuda.param.col[isource], marker=gw.marker, linestyle = 'solid', linewidth = 1 )
             ihyp += 1
             #
         gwav = nuda.astro.setupGWAverage( source = source )
@@ -58,7 +60,7 @@ def astro_setupGW_fig( pname, sources ):
     #
     axs.set_xticks( ilabel )
     axs.set_xticklabels( xlabel )
-    axs.legend(loc='lower center',bbox_to_anchor=(0.5,1.01),columnspacing=2,fontsize='8', ncol=2,frameon=False)
+    axs.legend(loc='lower center',bbox_to_anchor=(0.5,1.01),columnspacing=2,fontsize='8', ncol=3,frameon=False)
     #
     #
     if pname is not None:

nucleardatapy/fig/astro_setupMR_fig.py

@@ -5,16 +5,18 @@ import nucleardatapy as nuda
 
 def astro_setupMR_fig( pname, sources, sources_av ):
     """
-    Plot M-R constraints from NICER analyses.\
-    The plot is 1x1 with:\
+    Plot M-R constraints from observation measurements.
+
+    The plot is 1x1 with:
+    
     [0]: Masses versus radii.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
     :param sources: array of sources.
     :type sources: array of str.
-    :param sources_av: .
-    :type sources_av: .
+    :param sources_av: array of averaged sources.
+    :type sources_av: array of str.
 
     """
     #
@@ -22,12 +24,12 @@ def astro_setupMR_fig( pname, sources, sources_av ):
     #
     fig, axs = plt.subplots(1,1)
     fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.85, wspace=0.3, hspace=0.3)
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=0.95, top=0.85, wspace=0.3, hspace=0.3)
     #
-    axs.set_xlabel(r'$R$ (km)',fontsize='12')
-    axs.set_ylabel(r'M (M$_\odot$)',fontsize='12')
-    axs.set_xlim([10.5, 16.5])
-    axs.set_ylim([1.15, 2.2])
+    axs.set_xlabel(r'$R$ (km)',fontsize='14')
+    axs.set_ylabel(r'M (M$_\odot$)',fontsize='14')
+    axs.set_xlim([9.2, 16.5])
+    axs.set_ylim([0.58, 2.2])
     #
     isource = 1
     xlabel = []
@@ -56,31 +58,37 @@ def astro_setupMR_fig( pname, sources, sources_av ):
             #
         isource += 1
         #
-    isource = 1
-    for source in sources_av:
+    #isource = 1
+    for isource,source in enumerate(sources_av):
+        #if source.lower() == 'j0030+0451': 
+        #    obss = [ 1, 2 ]
+        #if source.lower() == 'j0740+6620': 
+        #    obss = [ 1, 2, 3 ]
+        obss = nuda.astro.mr_obss( source = source )
         if source.lower() == 'j0030+0451': 
             obss = [ 1, 2 ]
-        if source.lower() == 'j0740+6620': 
-            obss = [ 1, 2 ]
         mrav = nuda.astro.setupMRAverage( source = source, obss = obss )
         if nuda.env.verb_output: mrav.print_output( )
         if nuda.env.verb_latex: mrav.print_latex( )
-        axs.errorbar( mrav.rad_cen, mrav.mass_cen, xerr=mrav.rad_sig_std, yerr=mrav.mass_sig_std, ms=12, label=mrav.label, color=nuda.param.col[isource], marker='^', linestyle = 'solid', linewidth = 3 )
-        isource += 1
+        axs.errorbar( mrav.rad_cen, mrav.mass_cen, xerr=mrav.rad_sig_std, yerr=mrav.mass_sig_std, ms=12, label=mrav.label, color=nuda.param.col[isource+1], marker='^', linewidth = 3 )
     #
     # write source name in plot:
     #
     for source in sources:
         print('source:',source)
         if source.lower() == 'j0030+0451':
-            axs.text(13.3,1.25,'J0030+0451')
+            axs.text(13.3,1.55,'J0030+0451')
         elif source.lower() == 'j0740+6620':
             axs.text(13,1.95,'J0740+6620')
         elif source.lower() == 'j0437-4715':
-            axs.text(10.7,1.5,'J0437-4715')
+            axs.text(11.0,1.48,'J0437-4715')
+        elif source.lower() == 'j0614-3329':
+            axs.text(9.8,1.55,'J0614-3329')
+        elif source.lower() == 'j1731-347':
+            axs.text(10.5,0.8,'J1731-347')
     #
     #axs.legend(loc='upper left',fontsize='8', ncol=2)
-    axs.legend(loc='lower center',bbox_to_anchor=(0.48,1.01),columnspacing=2,fontsize='8',ncol=3,frameon=False)
+    axs.legend(loc='lower center',bbox_to_anchor=(0.48,1.01),columnspacing=2,fontsize='8',ncol=4,frameon=False)
     #
     if pname is not None:
     	plt.savefig(pname, dpi=200)

nucleardatapy/fig/astro_setupMasses_fig.py

@@ -5,8 +5,10 @@ import nucleardatapy as nuda
 
 def astro_setupMasses_fig( pname, sources ):
     """
-    Plot Masses for massives neutron stars as a function of sources.\
-    The plot is 1x1 with:\
+    Plot Masses for massives neutron stars as a function of sources.
+
+    The plot is 1x1 with:
+    
     [0]: masses versus sources.
 
     :param pname: name of the figure (*.png)
@@ -20,10 +22,10 @@ def astro_setupMasses_fig( pname, sources ):
     #
     fig, axs = plt.subplots(1,1)
     fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.7, wspace=0.3, hspace=0.3)
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=0.95, top=0.75, wspace=0.3, hspace=0.3)
     #
-    axs.set_ylabel(r'M (M$_\odot$)',fontsize='12')
-    axs.set_xlabel(r'sources',fontsize='12')
+    axs.set_ylabel(r'M (M$_\odot$)',fontsize='14')
+    axs.set_xlabel(r'sources',fontsize='14')
     axs.set_ylim([1.7, 3.5])
     axs.set_xlim([0.8, 5.5])
     #
@@ -56,7 +58,7 @@ def astro_setupMasses_fig( pname, sources ):
     axs.set_xticks( ilabel )
     axs.set_xticklabels( xlabel )
     #axs.legend(loc='upper left',fontsize='8', ncol=2)
-    axs.legend(loc='lower center',bbox_to_anchor=(0.5,1.01),columnspacing=2,fontsize='8', ncol=2,frameon=False)
+    axs.legend(loc='lower center',bbox_to_anchor=(0.5,1.01),columnspacing=2,fontsize='8', ncol=3,frameon=False)
     #
     if pname is not None:
     	plt.savefig(pname, dpi=200)

nucleardatapy/fig/astro_setupMtov_fig.py

@@ -5,8 +5,12 @@ import nucleardatapy as nuda
 
 def astro_setupMtov_fig( pname, sources_lo_all, sources_lo_dist, sources_up1, sources_up2, sources_up3 ):
     """
-    Plot the distribution of maximum masses.\
-    The plot is 1x1 with:\
+    Plot the probability distribution functions associated to maximum masses.
+
+    Includes the observation uncertainty in the PDF by using error-function.
+
+    The plot is 1x1 with:
+    
     [0]: distribution of maximum mass versus maximum mass.
 
     :param pname: name of the figure (*.png)
@@ -28,10 +32,10 @@ def astro_setupMtov_fig( pname, sources_lo_all, sources_lo_dist, sources_up1, so
     #
     fig, axs = plt.subplots(1,1)
     fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.85, wspace=0.3, hspace=0.3)
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=0.95, top=0.9, wspace=0.3, hspace=0.3)
     #
-    axs.set_xlabel(r'M$_\mathrm{tov}$ (M$_\odot$)',fontsize='12')
-    axs.set_ylabel(r'probability (non-normalised)',fontsize='12')
+    axs.set_xlabel(r'M$_\mathrm{tov}$ (M$_\odot$)',fontsize='14')
+    axs.set_ylabel(r'probability (non-normalised)',fontsize='14')
     axs.set_xlim([1.6, 3.4])
     axs.set_ylim([0.0, 1.03])
     #
@@ -76,7 +80,7 @@ def astro_setupMtov_fig( pname, sources_lo_all, sources_lo_dist, sources_up1, so
     #axs.plot(prob3.mass, prob3.proba_tot, label=prob3.label_tot, color=nuda.param.col[4], linestyle=(1, (4, 8)), linewidth = 3 )
     #
     #axs.legend( loc='upper left',fontsize='8', ncol=1 )
-    axs.legend(loc='lower center',bbox_to_anchor=(0.5,1.01),columnspacing=2,fontsize='8',ncol=4,frameon=False)
+    axs.legend(loc='lower center',bbox_to_anchor=(0.5,1.01),columnspacing=2,fontsize='8',ncol=5,frameon=False)
     #
     if pname is not None:
     	plt.savefig(pname, dpi=200)

nucleardatapy/fig/astro_setupMup_fig.py

@@ -5,8 +5,10 @@ import nucleardatapy as nuda
 
 def astro_setupMup_fig( pname, sources ):
     """
-    Plot upper boundaries for the tov mass.\
-    The plot is 1x1 with:\
+    Plot mass upper boundaries from GW measurements as a function of GW sources.
+
+    The plot is 1x1 with:
+    
     [0]: upper boundary for the mass versus sources.
 
     :param pname: name of the figure (*.png)
@@ -20,10 +22,10 @@ def astro_setupMup_fig( pname, sources ):
     #
     fig, axs = plt.subplots(1,1)
     fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.8, wspace=0.3, hspace=0.3)
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=0.95, top=0.85, wspace=0.3, hspace=0.3)
     #
-    axs.set_xlabel(r'sources',fontsize='12')
-    axs.set_ylabel(r'M (M$_\odot$)',fontsize='12')
+    axs.set_xlabel(r'sources',fontsize='14')
+    axs.set_ylabel(r'M (M$_\odot$)',fontsize='14')
     axs.set_xlim([0.8, 2.5])
     axs.set_ylim([2.4, 3.4])
     #
@@ -45,7 +47,7 @@ def astro_setupMup_fig( pname, sources ):
             mup = nuda.astro.setupMup( source = source, hyp = hyp )
             if nuda.env.verb_output: mup.print_output( )
             if nuda.env.verb_latex: mup.print_latex( )
-            axs.errorbar( isource+ihyp/10, mup.mup, yerr=np.array([(mup.sig_do,mup.sig_up)]).T, label=mup.label, color=nuda.param.col[isource], marker=mup.marker, linestyle = 'solid', linewidth = 1 )
+            axs.errorbar( isource+ihyp/10, mup.mup, yerr=np.array([(mup.sig_lo,mup.sig_up)]).T, label=mup.label, color=nuda.param.col[isource], marker=mup.marker, linestyle = 'solid', linewidth = 1 )
             ihyp += 1
             #
         if source=='GW170817': hyps = [ 3, 4 ]
@@ -60,7 +62,7 @@ def astro_setupMup_fig( pname, sources ):
     axs.set_xticklabels( xlabel )
     #
     #axs.legend(loc='upper left',fontsize='8', ncol=2)
-    axs.legend(loc='lower center',bbox_to_anchor=(0.5,1.01),columnspacing=2,fontsize='8',ncol=2,frameon=False)
+    axs.legend(loc='lower center',bbox_to_anchor=(0.5,1.01),columnspacing=2,fontsize='8',ncol=3,frameon=False)
     #
     if pname is not None:
     	plt.savefig(pname, dpi=200)

nucleardatapy/fig/corr_setupEsymDen_fig.py

@@ -3,16 +3,22 @@ import matplotlib.pyplot as plt
 
 import nucleardatapy as nuda
 
-def corr_setupEsymDen_fig( pname, constraints, Ksym ):
+def corr_setupEsymDen_fig( pname, constraints, Ksym, origine ):
     """
-    Plot upper boundaries for the tov mass.\
-    The plot is 1x1 with:\
+    Plot Esym or Esym,2 (depending on `origine`) as a function of the density.
+
+    The plot is 1x1 with:
+    
     [0]: upper boundary for the mass versus sources.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
-    :param sources: array of sources.
-    :type sources: array of str.
+    :param constraints: array of constraints to run on.
+    :type constraints: array of str.
+    :param Ksym: Value (in MeV) of Ksym.
+    :type Ksym: real number.
+    :param origine: can be 'finiteNuclei' or 'neutronStar'.
+    :type origine: str.
 
     """
     #
@@ -20,10 +26,17 @@ def corr_setupEsymDen_fig( pname, constraints, Ksym ):
     #
     fig, axs = plt.subplots(1,1)
     fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.98, wspace=0.3, hspace=0.3)
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=0.95, top=0.98, wspace=0.3, hspace=0.3)
     #
-    axs.set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)',fontsize='12')
-    axs.set_ylabel(r'$E_\text{sym}$ (MeV)',fontsize='12')
+    axs.set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)',fontsize='14')
+    if origine == 'finiteNuclei':
+        axs.set_ylabel(r'$e_{\text{sym},2}(n_\text{nuc})$ (MeV)',fontsize='14')
+    elif origine == 'neutronStar':
+        axs.set_ylabel(r'$e_\text{sym}(n_\text{nuc})$ (MeV)',fontsize='14')
+    else:
+        print('corr_setupEsymDen_fig.py: origine is not well documented, origine=',origine)
+        print('corr_setupEsymDen_fig.py: exit()')
+        exit()
     axs.set_xlim([0.09, 0.27])
     axs.set_ylim([10, 60])
     #
@@ -38,7 +51,7 @@ def corr_setupEsymDen_fig( pname, constraints, Ksym ):
         if esym.plot:
             axs.fill_between( esym.esym_den, y1=esym.esym_e2a_min, y2=esym.esym_e2a_max, label=esym.label, alpha=esym.alpha )
     #
-    axs.text(0.15,12,r'$K_\text{sym}$='+str(int(Ksym))+' MeV',fontsize='12')
+    axs.text(0.15,12,r'$K_\text{sym}$='+str(int(Ksym))+' MeV',fontsize='14')
     axs.legend(loc='lower right',fontsize='9')
     #
     if pname is not None:

nucleardatapy/fig/corr_setupEsymLsym_fig.py

@@ -3,27 +3,41 @@ import matplotlib.pyplot as plt
 
 import nucleardatapy as nuda
 
-def corr_setupEsymLsym_fig( pname, constraints ):
+def corr_setupEsymLsym_fig( pname, constraints, origine ):
     """
-    Plot the correlation between Esym and Lsym.\
-    The plot is 1x1 with:\
-    [0]: Esym - Lsym correlation plot
+    Plot the correlation between Esym and Lsym.
+
+    The plot is 1x1 with:
+
+    [0]: Esym - Lsym correlation plot.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
     :param constraints: list of constraints to run on.
     :type constraints: array of str.
+    :param origine: can be 'finiteNuclei' or 'neutronStar'.
+    :type origine: str.
 
     """
     #
     print(f'Plot name: {pname}')
+    print('list of constraints:',constraints)
     #
     fig, axs = plt.subplots(1,1)
     fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.98, wspace=0.3, hspace=0.3)
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=0.95, top=0.98, wspace=0.3, hspace=0.3)
+    #
+    if origine == 'finiteNuclei':
+        axs.set_xlabel(r'$E_{\mathrm{sym},2}$ (MeV)',fontsize='14')
+        axs.set_ylabel(r'$L_{\mathrm{sym},2}$ (MeV)',fontsize='14')
+    elif origine == 'neutronStar':
+        axs.set_xlabel(r'$E_\mathrm{sym}$ (MeV)',fontsize='14')
+        axs.set_ylabel(r'$L_\mathrm{sym}$ (MeV)',fontsize='14')
+    else:
+        print('corr_setupEsymLsym_fig.py: origine is not well documented, origine=',origine)
+        print('corr_setupEsymLsym_fig.py: exit()')
+        exit()
     #
-    axs.set_xlabel(r'$E_{\mathrm{sym},2}$ (MeV)')
-    axs.set_ylabel(r'$L_{\mathrm{sym},2}$ (MeV)')
     axs.set_xlim([23, 44])
     axs.set_ylim([10, 120])
     #
@@ -40,6 +54,9 @@ def corr_setupEsymLsym_fig( pname, constraints ):
         elif el.plot == 'curve':
             axs.plot( el.Esym, el.Lsym, linestyle='solid', linewidth=3, label=el.label )
         elif el.plot == 'contour':
+            print('plot:',el.plot)
+            print('Esym:',el.Esym)
+            print('Lsym:',el.Lsym)
             axs.plot( el.Esym, el.Lsym, linestyle='solid', label=el.label )
         elif el.plot == 'band_y':
             axs.fill_between( el.Esym, y1=el.Lsym-el.Lsym_err, y2=el.Lsym+el.Lsym_err, label=el.label, alpha=el.alpha )
@@ -49,7 +66,7 @@ def corr_setupEsymLsym_fig( pname, constraints ):
             #axs.errorbar( el.Esym, el.Lsym, yerr=el.Lsym_err, linestyle='solid', label=el.label )
         if nuda.env.verb: el.print_outputs( )
     #
-    axs.legend(loc='lower right',fontsize='9')
+    axs.legend(loc='lower right',fontsize='10')
     #
     if pname is not None:
     	plt.savefig(pname, dpi=300)

nucleardatapy/fig/corr_setupKsatQsat_fig.py

@@ -5,8 +5,10 @@ import nucleardatapy as nuda
 
 def corr_setupKsatQsat_fig( pname, constraints ):
     """
-    Plot the correlation between Ksat and Qsat.\
-    The plot is 1x1 with:\
+    Plot the correlation between Ksat and Qsat.
+
+    The plot is 1x1 with:
+    
     [0]: Ksat - Qsat correlation plot
 
     :param pname: name of the figure (*.png)
@@ -20,12 +22,12 @@ def corr_setupKsatQsat_fig( pname, constraints ):
     #
     fig, axs = plt.subplots(1,1)
     fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.15, bottom=0.12, right=None, top=0.98, wspace=0.3, hspace=0.3)
+    fig.subplots_adjust(left=0.15, bottom=0.12, right=0.95, top=0.98, wspace=0.3, hspace=0.3)
     #
-    axs.set_xlabel(r'$K_\mathrm{sat}$ (MeV)')
-    axs.set_ylabel(r'$Q_\mathrm{sat}$ (MeV)')
-    axs.set_xlim([190, 360])
-    axs.set_ylim([-1000, 1500])
+    axs.set_xlabel(r'$K_\mathrm{sat}$ (MeV)',fontsize='14')
+    axs.set_ylabel(r'$Q_\mathrm{sat}$ (MeV)',fontsize='14')
+    axs.set_xlim([100, 360])
+    axs.set_ylim([-2500, 2200])
     #axs.tick_params(labelbottom=True, labeltop=False, labelleft=True, labelright=False,
     #                 bottom=True, top=True, left=True, right=True)
 #    axs.xaxis.set_major_locator(MultipleLocator(5))
@@ -44,20 +46,24 @@ def corr_setupKsatQsat_fig( pname, constraints ):
         if nuda.env.verb: print('Qsat:',kq.Qsat)
         if nuda.env.verb: print('len(Ksat):',kq.Ksat.size)
         #
-        if k == 2:
-            kk = 0
-        else:
-            kk = k
-        axs.scatter( kq.Ksat, kq.Qsat, label=kq.label, color=nuda.param.col[kk], marker=kq.marker )
-        x = np.linspace(min(kq.Ksat),max(kq.Ksat),10)
-        if k == 3 or k == 4 or k == 5:
-            axs.plot( x, nuda.corr.flinear(x,kq.m,kq.c), color=nuda.param.col[kk], linestyle='dashed' )
+        #if k == 2:
+        #    kk = 0
+        #else:
+        #    kk = k
+        if k == 5 or k == 9 or k == 10 or k == 11: 
+            lstyle = 'dashed'
         else:
-            axs.plot( x, nuda.corr.flinear(x,kq.m,kq.c), color=nuda.param.col[kk], linestyle='solid' )
+            lstyle = 'solid'
+        if kq.Ksat is not None:
+            axs.scatter( kq.Ksat, kq.Qsat, label=kq.label, color = nuda.param.col[k], marker = kq.marker )
+        if kq.Ksat_lin is not None:
+            axs.plot( kq.Ksat_lin, kq.Qsat_lin, color = nuda.param.col[k], linestyle = lstyle )
+        if kq.Ksat_band is not None:
+            axs.fill_between( kq.Ksat_band, kq.Qsat_lo, kq.Qsat_up, label=kq.label, color = nuda.param.col[k], alpha = 0.2 )
         #
         if nuda.env.verb: kq.print_outputs( )
     #
-    axs.legend(loc='upper left',ncol=3, fontsize='9')
+    axs.legend(loc='upper left',ncol=4, fontsize='10')
     #
     if pname is not None: 
         plt.savefig(pname, dpi=200)

nucleardatapy/fig/crust_setupCrust_fig.py

@@ -5,13 +5,15 @@ import nucleardatapy as nuda
 
 def crust_setupCrust_fig( pname, models ):
     """
-    Plot hyper-nuclear chart (N versus Z).\
-    The plot is 1x2 with:\
-    [0]: nuclear chart.
+    Plot crust predictions for the models given in `models`.
+
+    The plot is 1x2 with:
+
+    [0]: internal energy per nucleon as a function of the density n. [1]: Z as a function of the density n.
 
     :param pname: name of the figure (*.png)
     :type pname: str.
-    :param models: table.
+    :param models: list of different models.
     :type models: str.
 
     """
@@ -22,14 +24,14 @@ def crust_setupCrust_fig( pname, models ):
     fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
     fig.subplots_adjust(left=0.10, bottom=0.12, right=None, top=0.78, wspace=0.3, hspace=0.3 )
     #
-    axs[0].set_xlabel(r'n (fm$^{-3}$)')
-    axs[0].set_ylabel(r'$e_{int}(n)$')
+    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)',fontsize='14')
+    axs[0].set_ylabel(r'$e_\text{int}(n_\text{nuc})$ (MeV)',fontsize='14')
     axs[0].set_xlim([1e-4, 1e-1])
     axs[0].set_ylim([-2, 10])
     axs[0].set_xscale('log')
     #
-    axs[1].set_xlabel(r'n (fm$^{-3}$)')
-    axs[1].set_ylabel(r'$Z$')
+    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)',fontsize='14')
+    axs[1].set_ylabel(r'$Z$',fontsize='14')
     axs[1].set_xlim([1e-4, 1e-1])
     axs[1].set_ylim([10, 100])
     axs[1].set_xscale('log')
@@ -44,7 +46,7 @@ def crust_setupCrust_fig( pname, models ):
             axs[1].plot( crust.den, crust.Z, linestyle=crust.linestyle )
     #axs[0].legend(loc='upper left',fontsize='8', ncol=1)
     #axs[1].legend(loc='upper left',fontsize='8', ncol=1)
-    fig.legend(loc='upper left',bbox_to_anchor=(0.01,1.01),columnspacing=2,fontsize='7.5',ncol=4,frameon=False)
+    fig.legend(loc='upper left',bbox_to_anchor=(0.08,1.01),columnspacing=2,fontsize='7.5',ncol=4,frameon=False)
     #
     if pname is not None:
     	plt.savefig(pname, dpi=200)