MAP model from Qin+2025 included

Author: sambit-giri

src/tools21cm/compare_constraints.py

@@ -283,8 +283,66 @@ def _load_models(self):
                                     "method": method, "name": "Evolving IMF"},
         ]
 
+        method = "21cmFAST"
+        MAPmod = np.array([
+                ##  z  MAP  16th 84th 
+                [15.20,0.99,0.99,0.99],
+                [14.88,0.99,0.99,0.99],
+                [14.57,0.99,0.99,0.99],
+                [14.26,0.98,0.98,0.99],
+                [13.96,0.98,0.98,0.98],
+                [13.67,0.98,0.98,0.98],
+                [13.38,0.97,0.97,0.98],
+                [13.10,0.97,0.97,0.97],
+                [12.82,0.96,0.96,0.97],
+                [12.55,0.96,0.96,0.96],
+                [12.29,0.95,0.95,0.96],
+                [12.03,0.94,0.94,0.95],
+                [11.77,0.93,0.93,0.94],
+                [11.52,0.92,0.92,0.93],
+                [11.28,0.91,0.91,0.92],
+                [11.04,0.90,0.90,0.91],
+                [10.80,0.89,0.89,0.90],
+                [10.57,0.87,0.87,0.89],
+                [10.34,0.86,0.86,0.87],
+                [10.12,0.84,0.84,0.85],
+                [9.90,0.82,0.82,0.84],
+                [9.69,0.80,0.80,0.82],
+                [9.48,0.77,0.77,0.79],
+                [9.27,0.75,0.75,0.77],
+                [9.07,0.72,0.72,0.75],
+                [8.87,0.70,0.70,0.72],
+                [8.68,0.67,0.67,0.69],
+                [8.49,0.64,0.64,0.66],
+                [8.30,0.61,0.61,0.63],
+                [8.12,0.57,0.58,0.60],
+                [7.94,0.54,0.54,0.57],
+                [7.77,0.51,0.51,0.54],
+                [7.60,0.47,0.48,0.50],
+                [7.43,0.44,0.44,0.47],
+                [7.26,0.40,0.41,0.44],
+                [7.10,0.37,0.37,0.40],
+                [6.94,0.33,0.34,0.36],
+                [6.78,0.29,0.30,0.33],
+                [6.63,0.26,0.26,0.29],
+                [6.48,0.22,0.23,0.25],
+                [6.34,0.19,0.19,0.21],
+                [6.19,0.15,0.15,0.17],
+                [6.05,0.11,0.12,0.14],
+                [5.91,0.08,0.08,0.10],
+                [5.78,0.05,0.05,0.07],
+                [5.64,0.03,0.03,0.04],
+                [5.51,0.01,0.01,0.02],
+                [5.39,0.00,0.00,0.00],
+            ])
+        self.neutral_fraction_models["Qin+2025"] = [
+                {"redshift": MAPmod[:,0], 
+                 "value": MAPmod[:,1], "error": [MAPmod[:,2]-MAPmod[:,1],MAPmod[:,3]-MAPmod[:,1]],
+                                    "method": method, "name": "MAP model"},
+        ]
+
         if self.verbose:
-            print(f'Theory:\n{list(self.neutral_fraction_constraints.keys())}')
+            print(f'Theory:\n{list(self.neutral_fraction_models.keys())}')
 
 class EscapeFractionLymanAlpha:
     def __init__(self, verbose=True):