def __init__(self):
obh2 = 0.0224
Om = 0.274
h = 0.7
mnu = 0
fidTheory = LCDMCosmology(obh2, Om, h, mnu)
GaussBAODVLikelihood.__init__(self, "DR11LOWZ", 0.32, 1264.0, 25.0,
fidTheory)
def ParseModel(model):
# if model=="LCDM_2":
if model == "LCDM" or "LCDM_" in model:
T = LCDMCosmology()
elif model == "LCDMmasslessnu" or "LCDMmasslessnu_" in model:
T = LCDMCosmology(mnu=0)
elif model == "nuLCDM" or "nuLCDM_" in model:
T = LCDMCosmology()
T.setVaryMnu()
elif model == "NnuLCDM" or "NnuLCDM_" in model:
LCDMCosmology.rd_approx = "tabulated_Nnu"
T = LCDMCosmology()
T.setVaryNnu()
elif model == "noradLCDM" or "noradLCDM_" in model:
T = LCDMCosmology(disable_radiation=True)
elif model == "oLCDM" or "oLCDM_" in model:
T = oLCDMCosmology()
elif model == "wCDM" or "wCDM_" in model:
T = wLCDMCosmology()
elif model == "waCDM" or "waCDM_" in model:
T = owa0CDMCosmology(varyOk=False)
elif model == "owCDM" or "owCDM_" in model:
T = owa0CDMCosmology(varywa=False)
elif model == "owaCDM" or "owaCDM_" in model:
T = owa0CDMCosmology()
elif model == "JordiCDM" or "JordiCDM_" in model:
T = JordiCDMCosmology()
elif model == "WeirdCDM" or "WeirdCDM_" in model:
T = WeirdCDMCosmology()
elif model == "TLight" or "TLight_" in model:
T = TLightCosmology()
elif model == "StepCDM" or "StepCDM_" in model:
T = StepCDMCosmology()
elif model == "Spline" or "Spline_" in model:
T = SplineLCDMCosmology()
elif model == "Decay" or "Decay_" in model:
T = DecayLCDMCosmology()
elif model == "PolyCDM" or "PolyCDM_" in model:
T = PolyCDMCosmology()
elif model == "EarlyDE" or "EarlyDE_" in model:
T = EarlyDECosmology()
else:
print "Cannot recognize model", model
sys.exit(1)
return T
def __init__(self):
obh2 = 0.021547
Om = 0.31
h = 0.67
mnu = 0
fidTheory = LCDMCosmology(obh2, Om, h, mnu)
TabulatedBAODVLikelihood.__init__(
self, "MGS", "data/chidavexi8stavePk5staverec.dat", fidTheory,
0.15)
def __init__(self):
obh2 = 0.0227
Om = 0.27
h = 0.7
mnu = 0 #;# rd=149.77
fidTheory = LCDMCosmology(obh2, Om, h, mnu)
TabulatedBAOLikelihood.__init__(self, "DR11LyaCross",
'data/lyabaocross.scan', 2, fidTheory,
2.36)
def __init__(self):
obh2 = 0.022
Om = 0.31
h = 0.676
mnu = 0.06
fidTheory = LCDMCosmology(obh2, Om, h, mnu)
TabulatedBAOLikelihood.__init__(
self, "ACPNS_1",
'data/aratdatacombrecACPNScombaccarec15st01covchi.dat', 2,
fidTheory, 0.39)
def __init__(self):
## fiducial cosmology for LOWZ/CMASS data.
## see Anderson et al, page 28
obh2 = 0.0224
Om = 0.274
h = 0.7
mnu = 0 #rd=149.28
fidTheory = LCDMCosmology(obh2, Om, h, mnu)
## negative col means the cols is probability and not chi2
TabulatedBAOLikelihood.__init__(self, "DR11CMASS",
'data/sdss_DR11CMASS_consensus.dat',
-2, fidTheory, 0.57)
def __init__(self):
obh2 = 0.02227
Om = 0.27
h = 0.7
mnu = 0
fidTheory = LCDMCosmology(obh2, Om, h, mnu)
GaussBAODVLikelihood.__init__(self,
"SixdFGS",
0.106,
456.0,
27.0,
fidTheory,
maxchi2=4)
def __init__(self):
## fiducial cosmology for Lya data.
## see e.g. Busca's email on 12/3/13
obh2 = 0.0227
Om = 0.27
h = 0.7
mnu = 0.06 #;# rd=149.77
fidTheory = LCDMCosmology(obh2, Om, h, mnu)
# File from 5/16 from Nicolas.
TabulatedBAOLikelihood.__init__(
self, "DR11LyaAuto", 'data/chi2_surface_dr11_baseline_fit.txt', 4,
fidTheory, 2.34)
def __init__(self):
obh2 = 0.022
Om = 0.31
h = 0.676
mnu = 0.06
fidTheory = LCDMCosmology(obh2, Om, h, mnu)
TabulatedBAOLikelihood.__init__(self,
"CombBAOzb3",
'data/combxiFBBAOlikgrid_zb3.dat',
-2,
fidTheory,
0.61,
order_aper_apar=False)
def __init__(self, varyw=True, varyOde=True, userd_DE=True):
# two parameters: Om and h
self.userd_DE = userd_DE
print('userd', userd_DE)
self.varyw = varyw
self.varyOde = varyOde
self.w0 = w_par.value
self.Ode = Ode_par.value
self.oC = LCDMCosmology()
LCDMCosmology.__init__(self)
def __init__(self):
# fiducial cosmology for Lya data.
# Taken from https://github.com/igmhub/picca/blob/master/data/deSainteAgatheetal2019/auto_alone_stdFit/auto_alone_stdFit..ap.at.scan.dat
# fiducial model -- see Table 2 of Victoria's paper
obh2 = 0.02222
h = 0.6731
Om = 0.1426 / h**2
mnu = 0.06 # rd=147.33
fidTheory = LCDMCosmology(obh2, Om, h, mnu)
TabulatedBAOLikelihood.__init__(
self,
"DR14LyaAuto",
'data/deSainteAgatheetal2019_ap_at_scan.dat',
2,
fidTheory,
2.34,
aperp_col=1,
apar_col=0,
skiprows=1)
def __init__(self):
# fiducial cosmology for Lya data.
# Taken from https://github.com/igmhub/picca/tree/master/data/Blomqvistetal2019/cross_alone_stdFit
# fiducial model -- double check
obh2 = 0.02222
h = 0.6731
Om = 0.1426 / h**2
mnu = 0.06 # rd=147.33
fidTheory = LCDMCosmology(obh2, Om, h, mnu)
# File from 5/16 from Nicolas.
TabulatedBAOLikelihood.__init__(
self,
"DR14LyaCross",
'data/Blomqvistetal2019_ap_at_scan.dat',
2,
fidTheory,
2.34,
aperp_col=1,
apar_col=0,
skiprows=1)
def __init__(self, varymu=True, varyAmp=True, varysig=True, varyCos=True):
## two parameters: Om and h
## we start with false here...
varyOk = False
## this is my "original cosmology" -- outside gaussian not much will change.
self.varyOk = varyOk
self.varymu = varymu
self.varyAmp = varyAmp
self.varysig = varysig
self.varyCos = varyCos
self.Ok = Ok_par.value
self.mu = mu_par.value
self.Amp = Amp_par.value
self.sig = sig_par.value
## auxiliary self cosmology
self.oC = LCDMCosmology()
LCDMCosmology.__init__(self)
def ParseModel(model):
if model == "LCDM":
T = LCDMCosmology()
elif model == "LCDMmasslessnu":
T = LCDMCosmology(mnu=0)
elif model == "nuLCDM":
T = LCDMCosmology()
T.setVaryMnu()
elif model == "NeffLCDM":
LCDMCosmology.rd_approx = "CuestaNeff"
T = LCDMCosmology()
T.setVaryNnu()
elif model == "noradLCDM":
T = LCDMCosmology(disable_radiation=True)
elif model == "oLCDM":
T = oLCDMCosmology()
elif model == "nuoLCDM":
T = oLCDMCosmology()
T.setVaryMnu()
elif model == "wCDM":
T = wLCDMCosmology()
elif model == "nuwCDM":
T = wLCDMCosmology()
T.setVaryMnu()
elif model == "waCDM":
T = owa0CDMCosmology(varyOk=False)
elif model == "owCDM":
T = owa0CDMCosmology(varywa=False)
elif model == "owaCDM":
T = owa0CDMCosmology()
elif model == "JordiCDM":
T = JordiCDMCosmology()
elif model == "WeirdCDM":
T = WeirdCDMCosmology()
elif model == "TLight":
T = TiredLightDecorator(PolyCDMCosmology())
elif model == "StepCDM":
T = StepCDMCosmology()
elif model == "Spline":
T = SplineLCDMCosmology()
elif model == "DecayFrac":
T = DecayLCDMCosmology()
elif model == "Decay":
T = DecayLCDMCosmology(varyxfrac=False, xfrac=1.0)
elif model == "Decay01":
T = DecayLCDMCosmology(varyxfrac=False, xfrac=0.1)
elif model == "Decay05":
T = DecayLCDMCosmology(varyxfrac=False, xfrac=0.5)
elif model == "PolyCDM":
T = PolyCDMCosmology()
elif model == "fPolyCDM":
T = PolyCDMCosmology(varyOk=False)
elif model == "EarlyDE":
T = EarlyDECosmology(varyw=False, userd_DE=False)
elif model == "EarlyDE_rd_DE":
T = EarlyDECosmology(varyw=False)
elif model == "SlowRDE":
T = SlowRDECosmology(varyOk=False)
elif model == "Quint_last":
T = QuintCosmology()
elif model == 'wDM':
T = wDMCosmology()
else:
print(("Cannot recognize model", model))
sys.exit(1)
return T
def ParseModel(model):
if model == "LCDM":
T = LCDMCosmology()
elif model == "LCDMmasslessnu":
T = LCDMCosmology(mnu=0)
elif model == "nuLCDM":
T = LCDMCosmology()
T.setVaryMnu()
elif model == "NeffLCDM":
LCDMCosmology.rd_approx = "CuestaNeff"
T = LCDMCosmology()
T.setVaryNnu()
elif model == "noradLCDM":
T = LCDMCosmology(disable_radiation=True)
elif model == "oLCDM":
T = oLCDMCosmology()
elif model == "nuoLCDM":
T = oLCDMCosmology()
T.setVaryMnu()
elif model == "wCDM":
T = wLCDMCosmology()
elif model == "nuwCDM":
T = wLCDMCosmology()
T.setVaryMnu()
elif model == "waCDM":
T = owa0CDMCosmology(varyOk=False)
elif model == "owCDM":
T = owa0CDMCosmology(varywa=False)
elif model == "owaCDM":
T = Lowa0CDMCosmology()
elif model == 'FCDM':
T = OzcosCosmology()
elif model == 'HornFco_c_cubic':
T = HornFcoCosmology(varyc=True)
elif model == 'HornFco_f_cubic':
T = HornFcoCosmology(varyf=True)
elif model == "JordiCDM":
T = JordiCDMCosmology()
elif model == "WeirdCDM":
T = WeirdCDMCosmology()
elif model == "TLight":
T = TiredLightDecorator(PolyCDMCosmology())
elif model == "StepCDM":
T = StepCDMCosmology()
elif model == "Spline":
T = SplineLCDMCosmology()
elif model == "DecayFrac":
T = DecayLCDMCosmology()
elif model == "Decay":
T = DecayLCDMCosmology(varyxfrac=False)
elif model == "Decay01":
T = DecayLCDMCosmology(varyxfrac=False, xfrac=0.1)
elif model == "Decay05":
T = DecayLCDMCosmology(varyxfrac=False, xfrac=0.5)
elif model == "PolyCDM":
T = PolyCDMCosmology()
elif model == "fPolyCDM":
T = PolyCDMCosmology(varyOk=False)
elif model == "EarlyDE":
T = EarlyDECosmology(varyw=False, userd_DE=False)
elif model == "EarlyDE_rd_DE":
T = EarlyDECosmology(varyw=False)
elif model == "SlowRDE":
T = SlowRDECosmology(varyOk=False)
elif model == "Quint":
T = QuintCosmology_try()
elif model == "Quintess":
T = QuintomCosmology(varymquin=True)
elif model == "Phantom":
T = QuintomCosmology(varymphan=True)
elif model == "Quintom":
T = QuintomCosmology(varymquin=True, varymphan=True, varyiniphi=True)
elif model == "Anisotropic":
T = AnisotropicCosmology()
elif model == "Fourier":
T = FourierCosmology()
elif model == 'wDM':
T = wDMCosmology()
else:
print "Cannot recognize model", model
sys.exit(1)
return T
