nucleardatapy 0.2.0__py3-none-any.whl → 1.0.0__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

@@ -39,7 +39,7 @@ 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 ]
     #
@@ -144,6 +144,14 @@ class setupMR():
                 #: 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:

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()")
     #