running code to produce results

Author: ja-vazquez

simplemc/likelihoods/BBNLikelihood.py

@@ -0,0 +1,31 @@
+
+from simplemc.likelihoods.BaseLikelihood import BaseLikelihood
+
+
+class BBNLikelihood (BaseLikelihood):
+    def __init__(self, mean=0.02218, sigma=0.00055):
+        """
+        A.-K. Burns, T.M.P. Tait and M. Valli, PRyMordial: the first
+        three minutes, within and beyond the standard model,
+        European Physical Journal C 84 (2024) 86 [2307.07061].
+        .
+        Parameters
+        ----------
+        mean
+        sigma
+
+        Returns
+        -------
+
+        """
+        BaseLikelihood.__init__(self, "BBN_%g_pm_%g" % (mean, sigma))
+        self.mean = mean
+        self.pr = 1/(2*sigma**2)
+
+    def loglike(self):
+        return -self.pr*(self.mean-self.theory_.Obh2)**2
+
+
+class BBN(BBNLikelihood):
+    def __init__(self):
+        BBNLikelihood.__init__(self, 0.02218, 0.00055)

simplemc/models/HolographicCosmology.py

@@ -51,7 +51,7 @@ def __init__(self, varyOk=False, varychde=True, mean=-2,
 
         # = 0 for CC, = -2 for Barrow
         mean = mean
-        priors = (-3.0 , 0.0)
+        priors = (-3.5 , -0.01)
         sigma = 0.1
 
         self.pname = 'amp_'

simplemc/runbase.py

@@ -50,6 +50,7 @@
 from .likelihoods.CompressedHDLikelihood import HubbleDiagram, HD23
 from .likelihoods.Compressedfs8Likelihood import fs8Diagram
 from .likelihoods.HubbleParameterLikelihood import RiessH0, RiessH0_21
+from .likelihoods.BBNLikelihood import BBN
 
 from .likelihoods.StrongLensingLikelihood import StrongLensing
 
@@ -174,8 +175,8 @@ def ParseModel(model, **kwargs):
         T = RotationCurves()
     elif model == "Chaplygin_s":
         T = ChaplyginCosmology(usesigmoid=True)
-    elif model == "HDE":
-        T = HolographicCosmology()
+    elif model == "HDE_node1":
+        T = HolographicCosmology(nodes=1)
     elif model == 'simple':
         T = SimpleModel(custom_parameters, custom_function)
     elif model == 'simple_cosmo':
@@ -284,7 +285,7 @@ def ParseDataset(datasets, **kwargs):
             L.addLikelihood(PLK15())
         elif name == 'PLK18':
             L.addLikelihood(PLK18())
-        elif name == 'PLKw':
+        elif name == 'PLKW':
             from .likelihoods.WangWangCMB import PlanckLikelihood
             L.addLikelihood(PlanckLikelihood())
         elif name == 'WMAP9':
@@ -324,6 +325,8 @@ def ParseDataset(datasets, **kwargs):
             L.addLikelihood(HubbleDiagram())
         elif name == 'HD23':
             L.addLikelihood(HD23())
+        elif name == 'BBN':
+            L.addLikelihood(BBN())
         elif name == 'fs8':
             L.addLikelihood(fs8Diagram())
         elif name == 'dline':