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

nucleardatapy/__init__.py

@@ -1,9 +1,11 @@
 """
-This module provides microscopic, phenomenological and experimental data constraints.
+This toolkit is a community-driven tool providing data for nuclear and astro-nuclear research community. It contains microscopic, phenomenological, experimental and observational data for the understanding of neutron stars.
+To get some help, please contact: nucleardatapy@gmail.com.
 """
 from nucleardatapy.cst     import *
 from nucleardatapy.env     import *
 from nucleardatapy.param   import *
+from nucleardatapy.hello   import *
 from nucleardatapy.create_folder import *
 #
 from nucleardatapy.matter  import *

nucleardatapy/astro/setup_gw.py

@@ -95,7 +95,7 @@ class setupGW():
         # fix `file_in` and some properties of the object
         #
         if source.lower()=='gw170817':
-            file_in = nuda.param.path_data+'astro/gw/GW170817.dat'
+            file_in = nuda.param.path_data+'astro/GW/GW170817.dat'
             if hyp==1:
                 #: Attribute providing the full reference to the paper to be citted.
                 self.ref='B.P. Abbott, R. Abbott, T.D. Abbott, et al., PRL 119, 161101 (2017)'
@@ -137,7 +137,7 @@ class setupGW():
                 self.note = "write here notes about this observation."
                 self.marker = 's'
         elif source.lower()=='gw190425':
-            file_in = nuda.param.path_data+'astro/gw/GW190425.dat'
+            file_in = nuda.param.path_data+'astro/GW/GW190425.dat'
             if hyp==1:
                 #: Attribute providing the full reference to the paper to be citted.
                 self.ref='B.P. Abbott, R. Abbott, T.D. Abbott, S. Abraham, et al., ApJL 892, L3 (2020)'
@@ -165,12 +165,12 @@ class setupGW():
                 if int(ele[0]) == hyp:
                     self.mchirp = float(ele[1])
                     self.mchirp_sig_up = float(ele[2])
-                    self.mchirp_sig_do = float(ele[3])
-                    self.q_do = float(ele[4])
+                    self.mchirp_sig_lo = float(ele[3])
+                    self.q_lo = float(ele[4])
                     self.q_up = float(ele[5])
                     self.lam = float(ele[6])
                     self.lam_sig_up = float(ele[7])
-                    self.lam_sig_do = float(ele[8])
+                    self.lam_sig_lo = float(ele[8])
                     self.latexCite = ele[9].replace('\n','').replace(' ','')
         #
         if nuda.env.verb: print("Exit setupGW()")
@@ -188,10 +188,10 @@ class setupGW():
             print("   source:       ",self.source)
             print("   hyp:          ",self.hyp)
             print("   m_chirp:      ",self.mchirp)
-            print("   sig(m_chirp): ",self.mchirp_sig_up,self.mchirp_sig_do)
-            print("   q_do, q_up:   ",self.q_do, self.q_up)
+            print("   sig(m_chirp): ",self.mchirp_sig_up,self.mchirp_sig_lo)
+            print("   q_lo, q_up:   ",self.q_lo, self.q_up)
             print("   lambda:       ",self.lam)
-            print("   sig(lambda):  ",self.lam_sig_up,self.lam_sig_do)
+            print("   sig(lambda):  ",self.lam_sig_up,self.lam_sig_lo)
             print("   latexCite:    ",self.latexCite)
             print("   ref:          ",self.ref)
             print("   label:        ",self.label)
@@ -210,7 +210,7 @@ class setupGW():
         if nuda.env.verb: print("Enter print_latex()")
         #
         if nuda.env.verb_latex:
-            print(rf"- table: {self.source} & {self.hyp} & ${self.mchirp:.4f}^{{{+self.mchirp_sig_up}}}_{{{-self.mchirp_sig_do}}}$ & $[{self.q_do}:{self.q_up}]$ & ${{{self.lam:.2f}}}^{{{+self.lam_sig_up}}}_{{{-self.lam_sig_do}}}$ & \cite{{{self.latexCite}}} \\\\")
+            print(rf"- table: {self.source} & {self.hyp} & ${self.mchirp:.4f}^{{{+self.mchirp_sig_up}}}_{{{-self.mchirp_sig_lo}}}$ & $[{self.q_lo}:{self.q_up}]$ & ${{{self.lam:.2f}}}^{{{+self.lam_sig_up}}}_{{{-self.lam_sig_lo}}}$ & \cite{{{self.latexCite}}} \\\\")
         else:
             print(rf"- No  table for source {self.source}. To get  table, write  'verb_table = True' in env.py.")
         #
@@ -236,14 +236,14 @@ class setupGW():
         #
         self.mchirp = None
         self.mchirp_sig_up = None
-        self.mchirp_sig_do = None
-        self.q_do = None
+        self.mchirp_sig_lo = None
+        self.q_lo = None
         self.q_up = None
         self.lam  = None
         #: Attribute the lower bound of the tidal deformability for the source.
         self.lambda_sig_up = None
         #: Attribute the upper bound of the tidal deformability for the source.
-        self.lambda_sig_do = None
+        self.lambda_sig_lo = None
         #
         if nuda.env.verb: print("Exit init_self()")
         #
@@ -279,9 +279,9 @@ class setupGWAverage():
         for hyp in hyps:
             gw = nuda.astro.setupGW( source = source, hyp = hyp )
             #gw.print_outputs( )
-            lamdo = gw.lam - 3*gw.lam_sig_do
+            lamlo = gw.lam - 3*gw.lam_sig_lo
             lamup = gw.lam + 3*gw.lam_sig_up
-            if lamdo < lammin: lammin = lamdo
+            if lamlo < lammin: lammin = lamlo
             if lamup > lammax: lammax = lamup
         # construct the distribution of observations in ay
         ax = np.linspace(lammin,lammax,300)
@@ -290,7 +290,7 @@ class setupGWAverage():
             #print('hyp:',hyp)
             gw = nuda.astro.setupGW( source = source, hyp = hyp )
             #gw.print_outputs( )
-            ay += gauss(ax,gw.lam,gw.lam_sig_up,gw.lam_sig_do)
+            ay += gauss(ax,gw.lam,gw.lam_sig_up,gw.lam_sig_lo)
         # determine the centroid and standard deviation from the distribution of obs. 
         nor = sum( ay )
         cen = sum( ay*ax )
@@ -355,13 +355,13 @@ class setupGWAverage():
         #
         return self        
 
-def gauss( ax, lam, sig_up, sig_do ):
+def gauss( ax, lam, sig_up, sig_lo ):
     fac = math.sqrt( 2*math.pi )
     gauss = []
     for x in ax:
         if x < lam: 
-            z = ( x - lam ) / sig_do
-            norm = sig_do * fac
+            z = ( x - lam ) / sig_lo
+            norm = sig_lo * fac
         else:
             z = ( x - lam ) / sig_up
             norm = sig_up * fac

nucleardatapy/astro/setup_mr.py

@@ -13,7 +13,7 @@ def mr_sources():
     #
     if nuda.env.verb: print("\nEnter astro_mr()")
     #
-    sources = [ 'J0030+0451', 'J0740+6620', 'J0437-4715' ]
+    sources = [ 'J0030+0451', 'J0740+6620', 'J0437-4715', 'J0614-3329', 'J1731-347' ]
     #
     #print('sources available in the toolkit:',sources)
     sources_lower = [ item.lower() for item in sources ]
@@ -39,9 +39,13 @@ def mr_obss( source ):
     if source.lower()=='j0030+0451':
         obss = [ 1, 2, 3, 4 ]
     elif source.lower()=='j0740+6620':
-        obss = [ 1, 2 ]
+        obss = [ 1, 2, 3 ]
     elif source.lower()=='j0437-4715':
         obss = [ 1 ]
+    elif source.lower()=='j0614-3329':
+        obss = [ 1 ]
+    elif source.lower()=='j1731-347':
+        obss = [ 1 ]
     #
     #print('Observations available in the toolkit:',obss)
     #
@@ -96,7 +100,7 @@ class setupMR():
             file_in = nuda.param.path_data+'astro/NICER/J0030+0451.dat'
             if obs==1:
                 #: Attribute providing the full reference to the paper to be citted.
-                self.ref='M.C. Miller, F.K. Lamb, A.J. Dittmann, aet al., ApJL 887, L24 (2019).'
+                self.ref='M.C. Miller, F.K. Lamb, A.J. Dittmann, aet al., ApJL 887, L24 (2019)'
                 #: Attribute providing the label the data is references for figures.
                 self.label = 'J0030 Miller 2019'
                 #: Attribute providing additional notes about the observation.
@@ -104,7 +108,7 @@ class setupMR():
                 self.marker = 'o'
             elif obs==2:
                 #: Attribute providing the full reference to the paper to be citted.
-                self.ref='T.E. Riley, A.L. Watts, S. Bogdanov, P.S. Ray, et al., ApJ 887, L21 (2019).'
+                self.ref='T.E. Riley, A.L. Watts, S. Bogdanov, P.S. Ray, et al., ApJ 887, L21 (2019)'
                 #: Attribute providing the label the data is references for figures.
                 self.label = 'J0030 Riley 2019'
                 #: Attribute providing additional notes about the observation.
@@ -112,7 +116,7 @@ class setupMR():
                 self.marker = 's'
             elif obs==3:
                 #: Attribute providing the full reference to the paper to be citted.
-                self.ref='S. Vinciguerra, T. Salmi, A.L. Watts, D. Choudhury, et al., ApJ 961, 62 (2024).'
+                self.ref='S. Vinciguerra, T. Salmi, A.L. Watts, D. Choudhury, et al., ApJ 961, 62 (2024)'
                 #: Attribute providing the label the data is references for figures.
                 self.label = 'J0030 Vinciguerra 2024a'
                 #: Attribute providing additional notes about the observation.
@@ -120,7 +124,7 @@ class setupMR():
                 self.marker = 'x'
             elif obs==4:
                 #: Attribute providing the full reference to the paper to be citted.
-                self.ref='S. Vinciguerra, T. Salmi, A.L. Watts, D. Choudhury, et al., ApJ 961, 62 (2024).'
+                self.ref='S. Vinciguerra, T. Salmi, A.L. Watts, D. Choudhury, et al., ApJ 961, 62 (2024)'
                 #: Attribute providing the label the data is references for figures.
                 self.label = 'J0030 Vinciguerra 2024b'
                 #: Attribute providing additional notes about the observation.
@@ -130,7 +134,7 @@ class setupMR():
             file_in = nuda.param.path_data+'astro/NICER/J0740+6620.dat'
             if obs==1:
                 #: Attribute providing the full reference to the paper to be citted.
-                self.ref='M.C. Miller, F.K. Lamb, A.J. Dittmann, S. Bogdanov, et al., ApJL 918, L28 (2021).'
+                self.ref='M.C. Miller, F.K. Lamb, A.J. Dittmann, S. Bogdanov, et al., ApJL 918, L28 (2021)'
                 #: Attribute providing the label the data is references for figures.
                 self.label = 'J0740 Miller 2021'
                 #: Attribute providing additional notes about the observation.
@@ -138,22 +142,50 @@ class setupMR():
                 self.marker = 'o'
             elif obs==2:
                 #: Attribute providing the full reference to the paper to be citted.
-                self.ref='T.E. Riley, A.L. Watts, P.S. Ray, S. Bogdanov, et al., ApJL 918, L27 (2021).'
+                self.ref='T.E. Riley, A.L. Watts, P.S. Ray, S. Bogdanov, et al., ApJL 918, L27 (2021)'
                 #: Attribute providing the label the data is references for figures.
                 self.label = 'J0740 Riley 2021'
                 #: Attribute providing additional notes about the observation.
                 self.note = "write notes about this observation."
                 self.marker = 's'
+            elif obs==3:
+                #: Attribute providing the full reference to the paper to be citted.
+                self.ref='T. Salmi, D. Choudhury, Y. Kini, T.E. Riley et al., ApJ 974, 294 (2024)'
+                #: Attribute providing the label the data is references for figures.
+                self.label = 'J0740 Salmi 2024'
+                #: Attribute providing additional notes about the observation.
+                self.note = "write notes about this observation."
+                self.marker = 's'
         elif source.lower()=='j0437-4715':
             file_in = nuda.param.path_data+'astro/NICER/J0437-4715.dat'
             if obs==1:
                 #: Attribute providing the full reference to the paper to be citted.
-                self.ref='D. Choudhury, T. Salmi, S. Vinciguerra, T.E. Riley, et al., ApJL 971, L20 (2024).'
+                self.ref='D. Choudhury, T. Salmi, S. Vinciguerra, T.E. Riley, et al., ApJL 971, L20 (2024)'
                 #: Attribute providing the label the data is references for figures.
                 self.label = 'J0437 Choudhury 2024'
                 #: Attribute providing additional notes about the observation.
                 self.note = "write notes about this observation."
                 self.marker = 'o'
+        elif source.lower()=='j0614-3329':
+            file_in = nuda.param.path_data+'astro/NICER/J0614-3329.dat'
+            if obs==1:
+                #: Attribute providing the full reference to the paper to be citted.
+                self.ref='L. Mauviard, S. Guillot, T. Salmi, D. Choudhury, et al., arXiv:2506.14883'
+                #: Attribute providing the label the data is references for figures.
+                self.label = 'J0614 Mauviard 2025'
+                #: Attribute providing additional notes about the observation.
+                self.note = "write notes about this observation."
+                self.marker = 'o'
+        elif source.lower()=='j1731-347':
+            file_in = nuda.param.path_data+'astro/HESS/J1731-347.dat'
+            if obs==1:
+                #: Attribute providing the full reference to the paper to be citted.
+                self.ref='V. Doroshenko, V. Suleimanov, G. P\"uhlhofer and A. Santangelo, Nature Astronomy volume 6, pages 1444–1451 (2022)'
+                #: Attribute providing the label the data is references for figures.
+                self.label = 'J1731 HESS 2022'
+                #: Attribute providing additional notes about the observation.
+                self.note = "write notes about this observation."
+                self.marker = 's'
         #
         #: Attribute the observational mass of the source.
         self.mass = None
@@ -176,24 +208,49 @@ class setupMR():
             for line in file:
                 if '#' in line: continue
                 ele = line.split(',')
-                #print('ele[0]:',ele[0],' obs:',obs,' ele[:]:',ele[:])
+                #print('ele[0]:',ele[0],' obs:',obs)
+                #print('ele[:]:',ele[:])
                 if int(ele[0]) == obs:
-                    self.rad = float(ele[1])
-                    self.rad_sig_up = float(ele[2])
-                    self.rad_sig_lo = float(ele[3])
-                    self.mass = float(ele[4])
-                    self.mass_sig_up = float(ele[5])
-                    self.mass_sig_lo = float(ele[6])
-                    #print('ele[7]:',ele[7],' len:',len(ele[7]))
-                    if len(ele[7])>1:
+                    try:
+                        self.rad = float(ele[1])
+                        self.rad_sig_up = float(ele[2])
+                        self.rad_sig_lo = float(ele[3])
+                    except ValueError:
+                        self.rad = 0.0
+                        self.rad_sig_up = 0.0
+                        self.rad_sig_lo = 0.0
+                    try:
+                        self.mass = float(ele[4])
+                        self.mass_sig_up = float(ele[5])
+                        self.mass_sig_lo = float(ele[6])
+                    except ValueError:
+                        self.mass = 0.0
+                        self.mass_sig_up = 0.0
+                        self.mass_sig_lo = 0.0
+                    #print('ele[7]:',ele[7],isinstance(ele[7].replace(' ',''), (int, float)),' len:',len(ele[7]))
+                    #print('ele[10]:',ele[10],isinstance(ele[10].replace(' ',''), (int, float)))
+                    try:
                         self.comp = float(ele[7])
                         self.comp_sig_up = float(ele[8])
                         self.comp_sig_lo = float(ele[9])
-                    else:
+                    except ValueError:
                         self.comp = 0.0
-                        self.comp_sig_up = 1.0
-                        self.comp_sig_lo = 1.0
-                    self.latexCite = ele[10].replace('\n','').replace(' ','')
+                        self.comp_sig_up = 0.0
+                        self.comp_sig_lo = 0.0
+                    try:
+                        self.spin = float(ele[10])
+                    except ValueError:
+                        self.spin = 0.0
+#                    if len(ele[7])>1:
+#                        self.comp = float(ele[7])
+#                        self.comp_sig_up = float(ele[8])
+#                        self.comp_sig_lo = float(ele[9])
+#                    else:
+#                        self.comp = 0.0
+#                        self.comp_sig_up = 1.0
+#                        self.comp_sig_lo = 1.0
+                    self.latexCite = ele[11].replace('\n','').replace(' ','')
+                    #print('outputs:',self.rad,self.mass,self.comp,self.spin,self.latexCite)
         #
         # compute compactness
         #
@@ -310,8 +367,8 @@ class setupMRAverage():
         # construct the distribution of observations in ay
         ar = np.linspace(rmin,rmax,300); ar=np.array( ar )
         am = np.linspace(mmin,mmax,300); am=np.array( am )
-        print('cmax:',cmax)
-        print('cmin:',cmin)
+        #print('cmax:',cmax)
+        #print('cmin:',cmin)
         ac = np.linspace(cmin,cmax,300); ac=np.array( ac )
         #ac = 0.5 * nuda.cst.rshsol_si / 1.e3 * am / ar
         ayr = np.zeros(300); ayr=np.array( ayr )
@@ -333,12 +390,36 @@ class setupMRAverage():
         stdr = sum ( ayr*ar**2 )
         stdm = sum ( aym*am**2 )
         stdc = sum ( ayc*ac**2 )
-        self.rad_cen = cenr / nor
-        self.mass_cen = cenm / nom
-        self.comp_cen = cenc / noc
-        self.rad_sig_std = round( math.sqrt( stdr/nor - self.rad_cen**2 ), 3 )
-        self.mass_sig_std = round( math.sqrt( stdm/nom - self.mass_cen**2 ), 3 )
-        self.comp_sig_std = round( math.sqrt( stdc/noc - self.comp_cen**2 ), 3 )
+        try:
+            if nor != 0.0:            
+                self.rad_cen = cenr / nor
+                self.rad_sig_std = round( math.sqrt( stdr/nor - self.rad_cen**2 ), 3 )
+            else:
+                self.rad_cen = 0.0
+                self.rad_sig_std = 0.0
+        except ZeroDivisionError:
+            self.rad_cen = 0.0
+            self.rad_sig_std = 0.0
+        try:
+            if nom != 0.0:
+                self.mass_cen = cenm / nom
+                self.mass_sig_std = round( math.sqrt( stdm/nom - self.mass_cen**2 ), 3 )
+            else:
+                self.mass_cen = 0.0
+                self.mass_sig_std = 0.0
+        except ZeroDivisionError:
+            self.mass_cen = 0.0
+            self.mass_sig_std = 0.0
+        try:
+            if noc != 0.0:
+                self.comp_cen = cenc / noc
+                self.comp_sig_std = round( math.sqrt( stdc/noc - self.comp_cen**2 ), 3 )
+            else:
+                self.comp_cen = 0.0
+                self.comp_sig_std = 0.0
+        except ZeroDivisionError:
+            self.comp_cen = 0.0
+            self.comp_sig_std = 0.0
         self.rad_cen = round( self.rad_cen, 3)
         self.mass_cen = round( self.mass_cen, 3)
         self.comp_cen = round( self.comp_cen, 3)
@@ -389,11 +470,20 @@ def gauss( ax, cent, sig_up, sig_lo ):
     gauss = []
     for x in ax:
         if x < cent: 
-            z = ( x - cent ) / sig_lo
+            if sig_lo != 0.0:
+                z = ( x - cent ) / sig_lo
+            else:
+                z= x * 0.0
             norm = sig_lo * fac
         else:
-            z = ( x - cent ) / sig_up
+            if sig_up != 0.0:
+                z = ( x - cent ) / sig_up
+            else:
+                z= x * 0.0
             norm = sig_up * fac
-        gauss.append( math.exp( -0.5*z**2 ) / norm )
+        if norm != 0.0:
+            gauss.append( math.exp( -0.5*z**2 ) / norm )
+        else:
+            gauss.append( 0.0 ) 
     return gauss
 

nucleardatapy/astro/setup_mup.py

@@ -141,7 +141,7 @@ class setupMup():
         #: Attribute the positive uncertainty.
         self.sig_up = None
         #: Attribute the negative uncertainty.
-        self.sig_do = None
+        self.sig_lo = None
         #: Attribute latexCite.
         self.latexCite = None
         #
@@ -156,7 +156,7 @@ class setupMup():
                 if int(ele[0]) == hyp:
                     self.mup = float(ele[2])
                     self.sig_up = float(ele[3])
-                    self.sig_do = float(ele[4])
+                    self.sig_lo = float(ele[4])
                     self.latexCite = ele[5].replace('\n','').replace(' ','')
         #
         if nuda.env.verb: print("Exit setupMup()")
@@ -173,7 +173,7 @@ class setupMup():
         print("   source:  ",self.source)
         print("   hyp:",self.hyp)
         print("   mup:",self.mup)
-        print("   sigma(mup):",self.sig_up,self.sig_do)
+        print("   sigma(mup):",self.sig_up,self.sig_lo)
         print("   latexCite:",self.latexCite)
         print("   ref:    ",self.ref)
         print("   label:  ",self.label)
@@ -190,7 +190,7 @@ class setupMup():
         if nuda.env.verb: print("Enter print_latex()")
         #
         if nuda.env.verb_latex:
-            print(rf"- table: {self.source} & {self.hyp} & ${self.mup:.2f}^{{{+self.sig_up}}}_{{{-self.sig_do}}}$ & \cite{{{self.latexCite}}} \\\\")
+            print(rf"- table: {self.source} & {self.hyp} & ${self.mup:.2f}^{{{+self.sig_up}}}_{{{-self.sig_lo}}}$ & \cite{{{self.latexCite}}} \\\\")
         else:
             print(rf"- No table for source {self.source}. To get table, write 'verb_latex = True' in env.py.")
         #
@@ -228,9 +228,9 @@ class setupMupAverage():
         for hyp in hyps:
             mup = nuda.setupMup( source = source, hyp = hyp )
             #mup.print_outputs( )
-            mupdo = mup.mup - 3*mup.sig_do
+            muplo = mup.mup - 3*mup.sig_lo
             mupup = mup.mup + 3*mup.sig_up
-            if mupdo < mupmin: mupmin = mupdo
+            if muplo < mupmin: mupmin = muplo
             if mupup > mupmax: mupmax = mupup
         # construct the distribution of observations in ay
         ax = np.linspace(mupmin,mupmax,300)
@@ -240,7 +240,7 @@ class setupMupAverage():
             #print('hyp:',hyp)
             mup = nuda.setupMup( source = source, hyp = hyp )
             #mup.print_outputs( )
-            ay += gauss(ax,mup.mup,mup.sig_up,mup.sig_do)
+            ay += gauss(ax,mup.mup,mup.sig_up,mup.sig_lo)
         # determine the centroid and standard deviation from the distribution of obs. 
         nor = sum( ay )
         cen = sum( ay*ax )
@@ -288,13 +288,13 @@ class setupMupAverage():
         if nuda.env.verb: print("Exit print_latex()")
         #
 
-def gauss( ax, mass, sig_up, sig_do ):
+def gauss( ax, mass, sig_up, sig_lo ):
     fac = math.sqrt( 2*math.pi )
     gauss = []
     for x in ax:
         if x < mass: 
-            z = ( x - mass ) / sig_do
-            norm = sig_do * fac
+            z = ( x - mass ) / sig_lo
+            norm = sig_lo * fac
         else:
             z = ( x - mass ) / sig_up
             norm = sig_up * fac

nucleardatapy/corr/setup_EsymDen.py

@@ -3,13 +3,8 @@ import sys
 import numpy as np  # 1.15.0
 import random
 
-#nucleardatapy_tk = os.getenv('NUCLEARDATAPY_TK')
-#sys.path.insert(0, nucleardatapy_tk)
-
 import nucleardatapy as nuda
 
-#nsat = 0.16
-
 class setupEsymDen():
     """
     Instantiate the values of Esym and Lsym from the constraint.

nucleardatapy/corr/setup_EsymLsym.py

@@ -12,17 +12,19 @@ def EsymLsym_constraints():
     """
     Return a list of constraints available in this toolkit in the \
     following list: '2009-HIC', '2010-RNP', '2012-FRDM', '2013-NS', \
-    '2014-IAS', '2014-IAS+RNP', '2015-POL-208PB', '2015-POL-120SN', \
-    '2015-POL-68NI', '2017-UG', '2021-PREXII-Reed', \
-    '2021-PREXII-Reinhard', '2023-PREXII+CREX-Zhang'; and \
-    print them all on the prompt.
+    '2014-IAS', '2014-IAS+RNP', '2015-POL-208Pb', '2015-POL-120Sn', \
+    '2015-POL-68Ni', '2017-UG', '2021-PREXII-Reed', \
+    '2021-PREXII-Reinhard', '2023-PREXII+CREX-Zhang', \
+    '2023-EDF-D4', '2023-EDF-D4-IAS', '2023-EDF-D4-IAS-Rnp'; \
+    and print them all on the prompt.
 
     :return: The list of constraints.
     :rtype: list[str].
     """
     constraints = [ '2009-HIC', '2010-RNP', '2012-FRDM', '2013-NS', '2014-IAS', '2014-IAS+RNP', \
-             '2015-POL-208PB', '2015-POL-120SN', '2015-POL-68NI', '2017-UG', \
-              '2021-PREXII-Reed', '2021-PREXII-Reinhard', '2023-PREXII+CREX-Zhang' ]
+             '2015-POL-208Pb', '2015-POL-120Sn', '2015-POL-68Ni', '2017-UG', \
+              '2021-PREXII-Reed', '2021-PREXII-Reinhard', '2023-PREXII+CREX-Zhang',\
+              '2023-EDF-D4', '2023-EDF-D4-IAS', '2023-EDF-D4-IAS-Rnp' ]
     #print('Constraints available in the toolkit:',constraints)
     constraints_lower = [ item.lower() for item in constraints ]
     return constraints, constraints_lower
@@ -47,9 +49,10 @@ class setupEsymLsym():
 
     The name of the constraint to be chosen in the \
     following list: '2009-HIC', '2010-RNP', '2012-FRDM', '2013-NS', \
-    '2014-IAS', '2014-IAS+RNP', '2015-POL-208PB', '2015-POL-120SN', \
-    '2015-POL-68NI', '2017-UG', '2021-PREXII-Reed', \
-    '2021-PREXII-Reinhard', '2021-PREXII+CREX-Zhang'.
+    '2014-IAS', '2014-IAS+RNP', '2015-POL-208Pb', '2015-POL-120Sn', \
+    '2015-POL-68Ni', '2017-UG', '2021-PREXII-Reed', \
+    '2021-PREXII-Reinhard', '2021-PREXII+CREX-Zhang',\
+    '2023-EDF-D4', '2023-EDF-D4-IAS', '2023-EDF-D4-IAS-Rnp'.
 
     :param constraint: Fix the name of `constraint`. Default value: '2014-IAS'.
     :type constraint: str, optional.
@@ -69,9 +72,9 @@ class setupEsymLsym():
         constraints, constraints_lower = EsymLsym_constraints()
         #
         if constraint.lower() not in constraints_lower:
-            print('The constraint ',constraint,' is not in the list of constraints.')
-            print('list of constraints:',constraints)
-            print('-- Exit the code --')
+            print('setup_EsymLsym.py: The constraint ',constraint,' is not in the list of constraints.')
+            print('setup_EsymLsym.py: list of constraints:',constraints)
+            print('setup_EsymLsym.py: -- Exit the code --')
             exit()
         #
         if constraint.lower() == '2009-hic':
@@ -372,7 +375,7 @@ class setupEsymLsym():
             file_in = os.path.join(nuda.param.path_data,'corr/EsymLsym/2021-PREXII-Reed.dat')
             if nuda.env.verb: print('Reads file:',file_in)
             self.ref = 'Reed et al., PRL 126, 172503 (2021)'
-            self.label = 'PREXII-Reed'
+            self.label = 'PREXII-Reed-2021'
             self.note = "."
             self.Esym, self.Esym_err, self.Lsym, self.Lsym_err = \
                 np.loadtxt( file_in, usecols=(0,1,2,3), unpack = True )
@@ -385,7 +388,7 @@ class setupEsymLsym():
             file_in = os.path.join(nuda.param.path_data,'corr/EsymLsym/2021-PREXII-Reinhard.dat')
             if nuda.env.verb: print('Reads file:',file_in)
             self.ref = 'Reinhard et al., PRL 127, 232501 (2021)'
-            self.label = 'PREXII-Reinhard'
+            self.label = 'PREXII-Reinhard-2021'
             self.note = "."
             self.Esym, self.Esym_err, self.Lsym, self.Lsym_err = \
                 np.loadtxt( file_in, usecols=(0,1,2,3), unpack = True )
@@ -398,7 +401,7 @@ class setupEsymLsym():
             file_in = os.path.join(nuda.param.path_data,'corr/EsymLsym/2023-PREXII-Zhang.dat')
             if nuda.env.verb: print('Reads file:',file_in)
             self.ref = 'Z. Zhang, L.W. Chen, Phys. Rev. C 108, 024317 (2023).'
-            self.label = 'PREXII+CREX-Zhang'
+            self.label = 'PREXII+CREX-Zhang-2023'
             self.note = "."
             self.Esym, self.Esym_err, self.Lsym, self.Lsym_err = \
                 np.loadtxt( file_in, usecols=(0,1,2,3), unpack = True )
@@ -406,10 +409,40 @@ class setupEsymLsym():
             self.Lsym_max = self.Lsym + self.Lsym_err
             self.plot = 'point_err_xy'
             #
+        elif constraint.lower() == '2023-edf-d4':
+            #
+            self.ref = 'B.V. Carlson, M. Dutra, O. Lourenco, K. Margueron, Phys. Rev. C 107, 0353805 (2022).'
+            self.label = 'EDF(D4)-2023'
+            self.note = "."
+            # D4:
+            self.Esym = [ 29.85, 31.97, 33.06, 37.26, 39.42, 38.58, 36.44, 30.97, 29.85 ]
+            self.Lsym = [ 50.3, 45.36, 55.45, 99.14, 126.6, 124.57, 115.71, 61.79, 50.3 ]
+            self.plot = 'contour'
+            #
+        elif constraint.lower() == '2023-edf-d4-ias':
+            #
+            self.ref = 'B.V. Carlson, M. Dutra, O. Lourenco, K. Margueron, Phys. Rev. C 107, 0353805 (2022).'
+            self.label = 'EDF(D4+IAS)-2023'
+            self.note = "."
+            # D4+IAS:
+            self.Esym = [ 29.85, 31.97, 33.06, 33.96, 34.54, 30.97, 29.85 ]
+            self.Lsym = [ 50.3, 45.36, 55.45, 71.55, 88.03, 61.79, 50.3 ]
+            self.plot = 'contour'
+            #
+        elif constraint.lower() == '2023-edf-d4-ias-rnp':
+            #
+            self.ref = 'B.V. Carlson, M. Dutra, O. Lourenco, K. Margueron, Phys. Rev. C 107, 0353805 (2022).'
+            self.label = 'EDF(D4+IAS+Rnp)-2023'
+            self.note = "."
+            # D4sym (D4+IAS+\Delta r_{np}):
+            self.Esym = [ 29.85, 32.01, 32.54, 32.74, 31.98, 30.97, 29.85 ]
+            self.Lsym = [ 50.3, 48.15, 61.52, 70.45, 67.44, 61.79, 50.3 ]
+            self.plot = 'contour'
+            #
         else:
             #
-            print('The variable constraint:',constraint)
-            print('does not fit with the options in the code')
+            print('setup_EsymLsym.py: The variable constraint:',constraint)
+            print('setup_EsymLsym.py: does not fit with the options in the code')
 
         if nuda.env.verb: print("Exit setupEsymLsym()")
     #