Merge pull request #364 from SNEWS2/JostMigenda/updateOscillationPara…

…meters Add NuFIT 6.0 and PDG 2024 oscillation parameter presets
SNEWS2 · Nov 25, 2024 · 1ad71b8 · 1ad71b8
2 parents 2856f12 + 49b3dc5
commit 1ad71b8
Showing 1 changed file with 56 additions and 12 deletions.
diff --git a/python/snewpy/neutrino.py b/python/snewpy/neutrino.py
@@ -163,49 +163,71 @@ def __post_init__(self):
         assert np.isclose(dm2_sum,0), f'dm41_2 - dm43_2 - dm31_2 = {dm2_sum} !=0'
 
 parameter_presets = {
-    # Values from JHEP 09 (2020) 178 [arXiv:2007.14792] and www.nu-fit.org.
     'NuFIT5.0': {
+        # Values from http://www.nu-fit.org/?q=node/228; cite as JHEP 09 (2020) 178 [arXiv:2007.14792]
         MassHierarchy.NORMAL:
         MixingParameters3Flavor(
-            theta12 = 33.44<< u.deg,
+            theta12 = 33.44 << u.deg,
             theta13 = 8.57 << u.deg,
-            theta23 = 49.20 << u.deg,
-            deltaCP = 197  << u.deg,
-            dm21_2 =  7.42e-5  << u.eV**2,
-            dm31_2 =  2.517e-3 << u.eV**2
+            theta23 = 49.2 << u.deg,
+            deltaCP = 197 << u.deg,
+            dm21_2 = 7.42e-5 << u.eV**2,
+            dm31_2 = 2.517e-3 << u.eV**2
         ),
         MassHierarchy.INVERTED:
-        MixingParameters3Flavor(    
+        MixingParameters3Flavor(
             theta12 = 33.45 << u.deg,
             theta13 = 8.60 << u.deg,
-            theta23 = 49.30 << u.deg,
+            theta23 = 49.3 << u.deg,
             deltaCP = 282 << u.deg,
             dm21_2 = 7.42e-5 << u.eV**2,
             dm32_2 = -2.498e-3 << u.eV**2
         )
     },
     'NuFIT5.2': {
+        # Values from http://www.nu-fit.org/?q=node/256; cite as JHEP 09 (2020) 178 [arXiv:2007.14792]
         MassHierarchy.NORMAL:
         MixingParameters3Flavor(
             theta12 = 33.41 << u.deg,
             theta13 = 8.58 << u.deg,
-            theta23 = 42.20 << u.deg,
+            theta23 = 42.2 << u.deg,
             deltaCP = 232 << u.deg,
             dm21_2 = 7.41e-5 << u.eV**2,
             dm31_2 = 2.507e-3 << u.eV**2
         ),
         MassHierarchy.INVERTED:
-        MixingParameters3Flavor(        
+        MixingParameters3Flavor(
             theta12 = 33.41 << u.deg,
             theta13 = 8.57 << u.deg,
-            theta23 = 49.00 << u.deg,
+            theta23 = 49.0 << u.deg,
             deltaCP = 276 << u.deg,
             dm21_2 = 7.41e-5 << u.eV**2,
             dm32_2 = -2.486e-3 << u.eV**2
         )
     },
+    'NuFIT6.0': {
+        # Values from http://www.nu-fit.org/?q=node/294; cite as arXiv:2410.05380
+        MassHierarchy.NORMAL:
+        MixingParameters3Flavor(
+            theta12 = 33.68 << u.deg,
+            theta13 = 8.56 << u.deg,
+            theta23 = 43.3 << u.deg,
+            deltaCP = 212 << u.deg,
+            dm21_2 = 7.49e-5 << u.eV**2,
+            dm31_2 = 2.513e-3 << u.eV**2
+        ),
+        MassHierarchy.INVERTED:
+        MixingParameters3Flavor(
+            theta12 = 33.68 << u.deg,
+            theta13 = 8.59 << u.deg,
+            theta23 = 47.9 << u.deg,
+            deltaCP = 274 << u.deg,
+            dm21_2 = 7.49e-5 << u.eV**2,
+            dm32_2 = -2.484e-3 << u.eV**2
+        )
+    },
     'PDG2022':{
-    # Values from https://pdg.lbl.gov
+        # Cite as R.L. Workman et al. (Particle Data Group), Prog. Theor. Exp. Phys. 2022, 083C01 (2022)
         MassHierarchy.NORMAL:
         MixingParameters3Flavor(
             theta12 = 33.65 << u.deg,
@@ -224,6 +246,28 @@ def __post_init__(self):
             dm21_2 = 7.53e-5 << u.eV**2,
             dm32_2 = -2.536e-3 << u.eV**2
         )
+    },
+    'PDG2024':{
+        # Values from https://pdglive.lbl.gov/Particle.action?node=S067&init=0
+        # Cite as S. Navas et al. (Particle Data Group), Phys. Rev. D 110, 030001 (2024)
+        MassHierarchy.NORMAL:
+        MixingParameters3Flavor(
+            theta12 = 33.65 << u.deg,
+            theta13 = 8.51 << u.deg,
+            theta23 = 48.33 << u.deg,
+            deltaCP = 214 << u.deg,
+            dm21_2 = 7.53e-5 << u.eV**2,
+            dm32_2 = 2.455e-3 << u.eV**2
+        ),
+        MassHierarchy.INVERTED:
+        MixingParameters3Flavor(
+            theta12 = 33.65 << u.deg,
+            theta13 = 8.51 << u.deg,
+            theta23 = 48.04 << u.deg,
+            deltaCP = 214 << u.deg,
+            dm21_2 = 7.53e-5 << u.eV**2,
+            dm32_2 = -2.529e-3 << u.eV**2
+        )
     }
 }