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

nucleardatapy/fig/matter_setupNEPModelDist_fig.py

@@ -1,68 +0,0 @@
-import numpy as np
-import matplotlib.pyplot as plt
-
-import nucleardatapy as nuda
-
-def matter_setupNEPModelDist_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).\
-
-    :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=None, 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.setupNEPModelDist( model )
-        #
-        axs[0,0].hist( dist.Esat, bins=10, label=model )
-        axs[1,0].hist( dist.nsat, bins=10 )
-        axs[2,0].hist( dist.Ksat, bins=10 )
-        axs[3,0].hist( dist.Qsat, bins=10 )
-        axs[4,0].hist( dist.msat, bins=10 )
-        axs[0,1].hist( dist.Esym, bins=10 )
-        axs[1,1].hist( dist.Lsym, bins=10 )
-        axs[2,1].hist( dist.Ksym, bins=10 )
-        axs[3,1].hist( dist.Qsym, bins=10 )
-        axs[4,1].hist( dist.Dmsat, bins=10 )
-    #
-    #axs[0,0].legend(loc='lower right',fontsize='10',ncol=2)
-    fig.legend(loc='upper left',bbox_to_anchor=(0.2,0.99),columnspacing=2,fontsize='8',ncol=4,frameon=False)
-    #
-    if pname is not None:
-        plt.savefig(pname, dpi=300)
-        plt.close()