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

nucleardatapy/fig/matter_setupMicro_gap_fig.py

@@ -3,159 +3,296 @@ import matplotlib.pyplot as plt
 
 import nucleardatapy as nuda
 
-def matter_setupMicro_gap_1s0_fig( pname, models, matter = 'NM' ):
+def matter_setupMicro_gap_1s0_fig(pname, models, matter="NM"):
     """
-    Plot the correlation between Esym and Lsym.\
-    The plot is 1x1 with:\
-    [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}')
+    print(f"Plot name: {pname}")
     #
-    print('models:',models)
+    print("models:", models)
     #
-    fig, axs = plt.subplots(2,2)
-    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.8, wspace=0.05, hspace=0.05 )
+    fig, axs = plt.subplots(2, 2)
+    fig.subplots_adjust( left=0.12, bottom=0.12, right=0.95, top=0.85, wspace=0.05, hspace=0.05 )
     #
-    axs[0,0].set_ylabel(r'$\Delta_{1S0}$ (MeV)')
-    axs[0,0].set_xlim([0, 0.048])
+    axs[0,0].set_ylabel(r"$\Delta_{1S0}$ (MeV)", fontsize="14")
+    axs[0,0].set_xlim([0, 0.09])
     axs[0,0].set_ylim([0, 3.0])
-    axs[0,0].tick_params('x', labelbottom=False)
+    axs[0,0].tick_params("x", labelbottom=False)
     #
     axs[0,1].set_xlim([0, 1.6])
     axs[0,1].set_ylim([0, 3.0])
-    axs[0,1].tick_params('x', labelbottom=False)
-    axs[0,1].tick_params('y', labelleft=False)
+    axs[0,1].tick_params("x", labelbottom=False)
+    axs[0,1].tick_params("y", labelleft=False)
     #
-    axs[1,0].set_ylabel(r'$\Delta_{1S0}/E_F$')
-    axs[1,0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    axs[1,0].set_xlim([0, 0.048])
-    axs[1,0].set_ylim([0, 0.78])
+    axs[1,0].set_ylabel(r"$\Delta_{1S0}/E_F$", fontsize="14")
+    axs[1,0].set_xlabel(r"$n_\text{nuc}$ (fm$^{-3}$)", fontsize="14")
+    axs[1,0].set_xlim([0, 0.09])
+    axs[1,0].set_ylim([0, 0.65])
     #
-    axs[1,1].set_xlabel(r'$k_{F_n}$ (fm$^{-1}$)')
+    axs[1,1].set_xlabel(r"$k_{F}$ (fm$^{-1}$)", fontsize="14")
     axs[1,1].set_xlim([0, 1.6])
-    axs[1,1].set_ylim([0, 0.78])
-    axs[1,1].tick_params('y', labelleft=False)
+    axs[1,1].set_ylim([0, 0.65])
+    axs[1,1].tick_params("y", labelleft=False)
     #
     for model in models:
         #
-        gap = nuda.matter.setupMicroGap( model = model, matter = matter )
+        gap = nuda.matter.setupMicroGap(model=model, matter=matter)
         #
-        if matter.lower() == 'nm':
+        if matter.lower() == "nm":
             if gap.nm_gap_1s0 is not None:
                 if gap.nm_gap_1s0_err is not None:
-                    axs[0,0].errorbar( gap.nm_den_1s0, gap.nm_gap_1s0, yerr=gap.nm_gap_1s0_err, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
-                    axs[0,1].errorbar( gap.nm_kfn_1s0, gap.nm_gap_1s0, yerr=gap.nm_gap_1s0_err, marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,0].errorbar( gap.nm_den_1s0, gap.nm_gap_1s0/nuda.eF_n(gap.nm_kfn_1s0), yerr=gap.nm_gap_1s0_err/nuda.eF_n(gap.nm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,1].errorbar( gap.nm_kfn_1s0, gap.nm_gap_1s0/nuda.eF_n(gap.nm_kfn_1s0), yerr=gap.nm_gap_1s0_err/nuda.eF_n(gap.nm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[0, 0].errorbar(gap.nm_den_1s0, gap.nm_gap_1s0, yerr=gap.nm_gap_1s0_err,
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle, label=gap.label )
+                    axs[0, 1].errorbar(gap.nm_kfn_1s0, gap.nm_gap_1s0, yerr=gap.nm_gap_1s0_err,
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle )
+                    axs[1, 0].errorbar(gap.nm_den_1s0, gap.nm_gap_1s0 / nuda.eF_n(gap.nm_kfn_1s0), yerr=gap.nm_gap_1s0_err / nuda.eF_n(gap.nm_kfn_1s0),
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle )
+                    axs[1, 1].errorbar(gap.nm_kfn_1s0, gap.nm_gap_1s0 / nuda.eF_n(gap.nm_kfn_1s0), yerr=gap.nm_gap_1s0_err / nuda.eF_n(gap.nm_kfn_1s0),
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle )
                 else:
-                    axs[0,0].plot( gap.nm_den_1s0, gap.nm_gap_1s0, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
-                    axs[0,1].plot( gap.nm_kfn_1s0, gap.nm_gap_1s0, marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,0].plot( gap.nm_den_1s0, gap.nm_gap_1s0/nuda.eF_n(gap.nm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,1].plot( gap.nm_kfn_1s0, gap.nm_gap_1s0/nuda.eF_n(gap.nm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
-        elif matter.lower() == 'sm':
+                    axs[0, 0].plot(gap.nm_den_1s0, gap.nm_gap_1s0,
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle, label=gap.label)
+                    axs[0, 1].plot(gap.nm_kfn_1s0, gap.nm_gap_1s0,
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle )
+                    axs[1, 0].plot(gap.nm_den_1s0, gap.nm_gap_1s0 / nuda.eF_n(gap.nm_kfn_1s0),
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle )
+                    axs[1, 1].plot(gap.nm_kfn_1s0, gap.nm_gap_1s0 / nuda.eF_n(gap.nm_kfn_1s0),
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle )
+        elif matter.lower() == "sm":
             if gap.sm_gap_1s0 is not None:
                 if gap.sm_gap_1s0_err is not None:
-                    axs[0,0].errorbar( gap.sm_den_1s0, gap.sm_gap_1s0, yerr=gap.sm_gap_1s0_err, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
-                    axs[0,1].errorbar( gap.sm_kfn_1s0, gap.sm_gap_1s0, yerr=gap.sm_gap_1s0_err, marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,0].errorbar( gap.sm_den_1s0, gap.sm_gap_1s0/nuda.eF_n(gap.sm_kfn_1s0), yerr=gap.sm_gap_1s0_err/nuda.eF_n(gap.sm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,1].errorbar( gap.sm_kfn_1s0, gap.sm_gap_1s0/nuda.eF_n(gap.sm_kfn_1s0), yerr=gap.sm_gap_1s0_err/nuda.eF_n(gap.sm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[0, 0].errorbar(gap.sm_den_1s0, gap.sm_gap_1s0, yerr=gap.sm_gap_1s0_err,
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle, label=gap.label )
+                    axs[0, 1].errorbar(gap.sm_kfn_1s0, gap.sm_gap_1s0, yerr=gap.sm_gap_1s0_err,
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle )
+                    axs[1, 0].errorbar(gap.sm_den_1s0, gap.sm_gap_1s0 / (2*nuda.eF_n(gap.sm_kfn_1s0)), yerr=gap.sm_gap_1s0_err / nuda.eF_n(gap.sm_kfn_1s0),
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle )
+                    axs[1, 1].errorbar(gap.sm_kfn_1s0, gap.sm_gap_1s0 / (2*nuda.eF_n(gap.sm_kfn_1s0)), yerr=gap.sm_gap_1s0_err / nuda.eF_n(gap.sm_kfn_1s0),
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle )
                 else:
-                    axs[0,0].plot( gap.sm_den_1s0, gap.sm_gap_1s0, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
-                    axs[0,1].plot( gap.sm_kfn_1s0, gap.sm_gap_1s0, marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,0].plot( gap.sm_den_1s0, gap.sm_gap_1s0/nuda.eF_n(gap.sm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,1].plot( gap.sm_kfn_1s0, gap.sm_gap_1s0/nuda.eF_n(gap.sm_kfn_1s0), marker=gap.marker, markevery=gap.every, linestyle='none' )
-        if nuda.env.verb_output: gap.print_outputs( )
+                    axs[0, 0].plot(gap.sm_den_1s0, gap.sm_gap_1s0,
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle, label=gap.label )
+                    axs[0, 1].plot(gap.sm_kfn_1s0, gap.sm_gap_1s0,
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle )
+                    axs[1, 0].plot(gap.sm_den_1s0, gap.sm_gap_1s0 / (2*nuda.eF_n(gap.sm_kfn_1s0)),
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle )
+                    axs[1, 1].plot(gap.sm_kfn_1s0, gap.sm_gap_1s0 / (2*nuda.eF_n(gap.sm_kfn_1s0)),
+                        marker=gap.marker, markevery=gap.every, linestyle=gap.lstyle )
+        if nuda.env.verb_output:
+            gap.print_outputs()
     #
-    #axs[1,0].legend(loc='upper right',fontsize='8')
-    fig.legend(loc='upper left',bbox_to_anchor=(0.08,1.0),columnspacing=2,fontsize='8',ncol=3,frameon=False)
+    # axs[1,0].legend(loc='upper right',fontsize='8')
+    fig.legend(loc="upper left", bbox_to_anchor=(0.1, 1.0), columnspacing=2, fontsize="8", ncol=3, frameon=False )
     #
     if pname is not None:
-    	plt.savefig(pname, dpi=300)
-    	plt.close()
+        plt.savefig(pname, dpi=300)
+        plt.close()
+
 
-def matter_setupMicro_gap_3pf2_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 1x1 with:\
-    [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}')
+    print(f"Plot name: {pname}")
     #
     # plot 3PF2 pairing gap in NM
-    print('models:',models)
+    print("models:", models)
     #
-    fig, axs = plt.subplots(2,2)
-    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.8, wspace=0.05, hspace=0.05 )
+    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.85, wspace=0.05, hspace=0.05 )
     #
-    axs[0,0].set_ylabel(r'$\Delta_{3PF2}$ (MeV)')
-    axs[0,0].set_xlim([0, 0.38])
-    axs[0,0].set_ylim([0, 0.6])
-    axs[0,0].tick_params('x', labelbottom=False)
+    axs[0, 0].set_ylabel(r"$\Delta_{3PF2}$ (MeV)", fontsize="14")
+    axs[0, 0].set_xlim([0, 0.38])
+    axs[0, 0].set_ylim([0, 0.6])
+    axs[0, 0].tick_params("x", labelbottom=False)
     #
-    axs[0,1].set_xlim([0.6, 2.3])
-    axs[0,1].set_ylim([0, 0.6])
-    axs[0,1].tick_params('x', labelbottom=False)
-    axs[0,1].tick_params('y', labelleft=False)
+    axs[0, 1].set_xlim([0.6, 2.3])
+    axs[0, 1].set_ylim([0, 0.6])
+    axs[0, 1].tick_params("x", labelbottom=False)
+    axs[0, 1].tick_params("y", labelleft=False)
     #
-    axs[1,0].set_ylabel(r'$100\times \Delta_{3PF2}/E_F$')
-    axs[1,0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    axs[1,0].set_xlim([0, 0.38])
-    axs[1,0].set_ylim([0, 0.78])
+    axs[1, 0].set_ylabel(r"$100\times \Delta_{3PF2}/E_F$", fontsize="14")
+    axs[1, 0].set_xlabel(r"$n_\text{nuc}$ (fm$^{-3}$)", fontsize="14")
+    axs[1, 0].set_xlim([0, 0.38])
+    axs[1, 0].set_ylim([0, 0.78])
     #
-    axs[1,1].set_xlabel(r'$k_{F_n}$ (fm$^{-1}$)')
-    axs[1,1].set_xlim([0.6, 2.3])
-    axs[1,1].set_ylim([0, 0.78])
-    axs[1,1].tick_params('y', labelleft=False)
+    axs[1, 1].set_xlabel(r"$k_{F}$ (fm$^{-1}$)", fontsize="14")
+    axs[1, 1].set_xlim([0.6, 2.3])
+    axs[1, 1].set_ylim([0, 0.78])
+    axs[1, 1].tick_params("y", labelleft=False)
     #
     for model in models:
         #
-        gap = nuda.matter.setupMicroGap( model = model, matter = matter )
+        gap = nuda.matter.setupMicroGap(model=model, matter=matter)
         #
-        if matter.lower() == 'nm':
+        if matter.lower() == "nm":
             if gap.nm_gap_3pf2 is not None:
                 if gap.nm_gap_3pf2_err is not None:
-                    axs[0,0].errorbar( gap.nm_den_3pf2, gap.nm_gap_3pf2, yerr=gap.nm_gap_3pf2_err, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
-                    axs[0,1].errorbar( gap.nm_kfn_3pf2, gap.nm_gap_3pf2, yerr=gap.nm_gap_3pf2_err, marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,0].errorbar( gap.nm_den_3pf2, 100*gap.nm_gap_3pf2/nuda.eF_n(gap.nm_kfn_3pf2), yerr=gap.nm_gap_3pf2_err/nuda.eF_n(gap.nm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,1].errorbar( gap.nm_kfn_3pf2, 100*gap.nm_gap_3pf2/nuda.eF_n(gap.nm_kfn_3pf2), yerr=gap.nm_gap_3pf2_err/nuda.eF_n(gap.nm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[0, 0].errorbar(
+                        gap.nm_den_3pf2,
+                        gap.nm_gap_3pf2,
+                        yerr=gap.nm_gap_3pf2_err,
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                        label=gap.label,
+                    )
+                    axs[0, 1].errorbar(
+                        gap.nm_kfn_3pf2,
+                        gap.nm_gap_3pf2,
+                        yerr=gap.nm_gap_3pf2_err,
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                    )
+                    axs[1, 0].errorbar(
+                        gap.nm_den_3pf2,
+                        100 * gap.nm_gap_3pf2 / nuda.eF_n(gap.nm_kfn_3pf2),
+                        yerr=gap.nm_gap_3pf2_err / nuda.eF_n(gap.nm_kfn_3pf2),
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                    )
+                    axs[1, 1].errorbar(
+                        gap.nm_kfn_3pf2,
+                        100 * gap.nm_gap_3pf2 / nuda.eF_n(gap.nm_kfn_3pf2),
+                        yerr=gap.nm_gap_3pf2_err / nuda.eF_n(gap.nm_kfn_3pf2),
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                    )
                 else:
-                    axs[0,0].plot( gap.nm_den_3pf2, gap.nm_gap_3pf2, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
-                    axs[0,1].plot( gap.nm_kfn_3pf2, gap.nm_gap_3pf2, marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,0].plot( gap.nm_den_3pf2, 100*gap.nm_gap_3pf2/nuda.eF_n(gap.nm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,1].plot( gap.nm_kfn_3pf2, 100*gap.nm_gap_3pf2/nuda.eF_n(gap.nm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
-        elif matter.lower() == 'sm':
+                    axs[0, 0].plot(
+                        gap.nm_den_3pf2,
+                        gap.nm_gap_3pf2,
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                        label=gap.label,
+                    )
+                    axs[0, 1].plot(
+                        gap.nm_kfn_3pf2,
+                        gap.nm_gap_3pf2,
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                    )
+                    axs[1, 0].plot(
+                        gap.nm_den_3pf2,
+                        100 * gap.nm_gap_3pf2 / nuda.eF_n(gap.nm_kfn_3pf2),
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                    )
+                    axs[1, 1].plot(
+                        gap.nm_kfn_3pf2,
+                        100 * gap.nm_gap_3pf2 / nuda.eF_n(gap.nm_kfn_3pf2),
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                    )
+        elif matter.lower() == "sm":
             if gap.sm_gap_3pf2 is not None:
                 if gap.sm_gap_3pf2_err is not None:
-                    axs[0,0].errorbar( gap.sm_den_3pf2, gap.sm_gap_3pf2, yerr=gap.sm_gap_3pf2_err, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
-                    axs[0,1].errorbar( gap.sm_kfn_3pf2, gap.sm_gap_3pf2, yerr=gap.sm_gap_3pf2_err, marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,0].errorbar( gap.sm_den_3pf2, 100*gap.sm_gap_3pf2/nuda.eF_n(gap.sm_kfn_3pf2), yerr=gap.sm_gap_3pf2_err/nuda.eF_n(gap.sm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,1].errorbar( gap.sm_kfn_3pf2, 100*gap.sm_gap_3pf2/nuda.eF_n(gap.sm_kfn_3pf2), yerr=gap.sm_gap_3pf2_err/nuda.eF_n(gap.sm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
+                    axs[0, 0].errorbar(
+                        gap.sm_den_3pf2,
+                        gap.sm_gap_3pf2,
+                        yerr=gap.sm_gap_3pf2_err,
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                        label=gap.label,
+                    )
+                    axs[0, 1].errorbar(
+                        gap.sm_kfn_3pf2,
+                        gap.sm_gap_3pf2,
+                        yerr=gap.sm_gap_3pf2_err,
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                    )
+                    axs[1, 0].errorbar(
+                        gap.sm_den_3pf2,
+                        100 * gap.sm_gap_3pf2 / (2*nuda.eF_n(gap.sm_kfn_3pf2)),
+                        yerr=gap.sm_gap_3pf2_err / (2*nuda.eF_n(gap.sm_kfn_3pf2)),
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                    )
+                    axs[1, 1].errorbar(
+                        gap.sm_kfn_3pf2,
+                        100 * gap.sm_gap_3pf2 / (2*nuda.eF_n(gap.sm_kfn_3pf2)),
+                        yerr=gap.sm_gap_3pf2_err / (2*nuda.eF_n(gap.sm_kfn_3pf2)),
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                    )
                 else:
-                    axs[0,0].plot( gap.sm_den_3pf2, gap.sm_gap_3pf2, marker=gap.marker, markevery=gap.every, linestyle='none', label=gap.label )
-                    axs[0,1].plot( gap.sm_kfn_3pf2, gap.sm_gap_3pf2, marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,0].plot( gap.sm_den_3pf2, 100*gap.sm_gap_3pf2/nuda.eF_n(gap.sm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
-                    axs[1,1].plot( gap.sm_kfn_3pf2, 100*gap.sm_gap_3pf2/nuda.eF_n(gap.sm_kfn_3pf2), marker=gap.marker, markevery=gap.every, linestyle='none' )
-        if nuda.env.verb_output: gap.print_outputs( )
+                    axs[0, 0].plot(
+                        gap.sm_den_3pf2,
+                        gap.sm_gap_3pf2,
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                        label=gap.label,
+                    )
+                    axs[0, 1].plot(
+                        gap.sm_kfn_3pf2,
+                        gap.sm_gap_3pf2,
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                    )
+                    axs[1, 0].plot(
+                        gap.sm_den_3pf2,
+                        100 * gap.sm_gap_3pf2 / (2*nuda.eF_n(gap.sm_kfn_3pf2)),
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                    )
+                    axs[1, 1].plot(
+                        gap.sm_kfn_3pf2,
+                        100 * gap.sm_gap_3pf2 / (2*nuda.eF_n(gap.sm_kfn_3pf2)),
+                        marker=gap.marker,
+                        markevery=gap.every,
+                        linestyle=gap.lstyle,
+                    )
+        if nuda.env.verb_output:
+            gap.print_outputs()
     #
-    #axs[1,0].legend(loc='upper right',fontsize='8')
-    fig.legend(loc='upper left',bbox_to_anchor=(0.04,1.0),columnspacing=2,fontsize='8',ncol=3,frameon=False)
+    # axs[1,0].legend(loc='upper right',fontsize='8')
+    fig.legend(loc="upper left", bbox_to_anchor=(0.1, 0.95), columnspacing=2, fontsize="8", ncol=3, frameon=False )
     #
     if pname is not None:
-    	plt.savefig(pname, dpi=300)
-    	plt.close()
+        plt.savefig(pname, dpi=300)
+        plt.close()
+    #

nucleardatapy/fig/matter_setupNEPStats_fig.py

@@ -0,0 +1,106 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+
+def matter_setupNEPStats_fig(pname, models):
+    """
+    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.
+    :param models: models to run on.
+    :type models: array of str.
+
+    """
+    #
+    print(f"Plot name: {pname}")
+    #
+    fig, axs = plt.subplots(5, 2)
+    fig.tight_layout()  # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust( left=0.12, bottom=0.06, right=0.95, top=0.9, wspace=0.3, hspace=0.4 )
+    #
+    axs[0, 0].set_ylabel(r"$E_\text{sat}$")
+    axs[0, 0].set_xlim([-16.5, -15])
+    axs[1, 0].set_ylabel(r"$n_\text{sat}$")
+    axs[1, 0].set_xlim([0.14, 0.18])
+    axs[2, 0].set_ylabel(r"$K_\text{sat}$")
+    axs[2, 0].set_xlim([180, 360])
+    axs[3, 0].set_ylabel(r"$Q_\text{sat}$")
+    axs[3, 0].set_xlim([-1000, 1000])
+    axs[4, 0].set_ylabel(r"$m_\text{sat}^{*}/m$")
+    axs[4, 0].set_xlim([0.35, 1.2])
+    axs[0, 1].set_ylabel(r"$E_\text{sym}$")
+    axs[0, 1].set_xlim([26, 40])
+    axs[1, 1].set_ylabel(r"$L_\text{sym}$")
+    axs[1, 1].set_xlim([0, 120])
+    axs[2, 1].set_ylabel(r"$K_\text{sym}$")
+    axs[2, 1].set_xlim([-400, 220])
+    axs[3, 1].set_ylabel(r"$Q_\text{sym}$")
+    axs[3, 1].set_xlim([-50, 900])
+    axs[4, 1].set_ylabel(r"$\Delta m_\text{sat}^{*}/m$")
+    axs[4, 1].set_xlim([-0.5, 1.1])
+    #
+    # Built distribution of NEP
+    #
+    for model in models:
+        #
+        dist = nuda.matter.setupNEPStat_model(model)
+        #
+        xbins = np.arange(-16.5, -15.0, 0.15)
+        if len(dist.Esat) != 0:
+            axs[0, 0].hist(dist.Esat, bins=xbins, alpha=0.5, weights=1/len(dist.Esat) * np.ones(len(dist.Esat)) )
+        xbins = np.arange(0.14, 0.18, 0.004)
+        if len(dist.nsat) != 0:
+            axs[1, 0].hist(dist.nsat, bins=xbins, alpha=0.5, weights=1/len(dist.nsat) * np.ones(len(dist.nsat)) )
+        xbins = np.arange(180, 320, 20)
+        if len(dist.Ksat) != 0:
+            axs[2, 0].hist(dist.Ksat, bins=xbins, alpha=0.5, weights=1/len(dist.Ksat) * np.ones(len(dist.Ksat)), label=model )
+        xbins = np.arange(-1000, 1000, 200)
+        if len(dist.Qsat) != 0:
+            axs[3, 0].hist(dist.Qsat, bins=xbins, alpha=0.5, weights=1/len(dist.Qsat) * np.ones(len(dist.Qsat)) )
+        xbins = np.arange(0.5, 1.2, 0.1)
+        if len(dist.msat) != 0:
+            axs[4, 0].hist(dist.msat, bins=xbins, alpha=0.5, weights=1/len(dist.msat) * np.ones(len(dist.msat)) )
+        xbins = np.arange(25, 40, 1.0)
+        if len(dist.Esym) != 0:
+            axs[0, 1].hist(dist.Esym, bins=xbins, alpha=0.5, weights=1/len(dist.Esym) * np.ones(len(dist.Esym)) )
+        xbins = np.arange(0, 120, 10)
+        if len(dist.Lsym) != 0:
+            axs[1, 1].hist(dist.Lsym, bins=xbins, alpha=0.5, weights=1/len(dist.Lsym) * np.ones(len(dist.Lsym)) )
+        xbins = np.arange(-400, 400, 100)
+        if len(dist.Ksym) != 0:
+            axs[2, 1].hist(dist.Ksym, bins=xbins, alpha=0.5, weights=1/len(dist.Ksym) * np.ones(len(dist.Ksym)) )
+        xbins = np.arange(0, 1100, 100)
+        if len(dist.Qsym) != 0:
+            axs[3, 1].hist(dist.Qsym, bins=xbins, alpha=0.5, weights=1/len(dist.Qsym) * np.ones(len(dist.Qsym)) )
+        xbins = np.arange(-0.6, 1.0, 0.1)
+        if len(dist.Dmsat) != 0:
+            axs[4, 1].hist(dist.Dmsat, bins=xbins, alpha=0.5, weights=1/len(dist.Dmsat) * np.ones(len(dist.Dmsat)) )
+    #
+    # axs[0,0].legend(loc='lower right',fontsize='10',ncol=2)
+    fig.legend(
+        loc="upper left",
+        bbox_to_anchor=(0.02, 0.99),
+        columnspacing=2,
+        fontsize="7.5",
+        ncol=7,
+        frameon=False,
+    )
+    #
+    if pname is not None:
+        plt.savefig(pname, dpi=300)
+        plt.close()