- Small changes to CAMB sideband model to match method described in t…

…he paper draft. - Corresponding small change to CAMB BBand1 models. - Allow multiple model-config options that are applied in turn (INI file first). - Add EH98 models (see comments in config/flux_Aln2.ini for an example of using these).
igmhub · Aug 27, 2012 · 3b8af16 · 3b8af16
2 parents 1edbe31 + b7d63fc
commit 3b8af16
Show file tree

Hide file tree

Showing 32 changed files with 8,790 additions and 1,427 deletions.
diff --git a/config/flux_Aln2.ini b/config/flux_Aln2.ini
@@ -1,33 +1,42 @@
-# Homogeneous cosmology for converting (ll,sep,z) into co-moving coords
+## Homogeneous cosmology for converting (ll,sep,z) into co-moving coords
 omega-matter = 0.27
 hubble-constant = 0.7
 
-# Tabulated models to use
+## Tabulated models to use: CAMB / EH98
 modelroot = ../models/
 fiducial =  DR9LyaMocks
 nowiggles = DR9LyaMocksSB
 broadband = DR9LyaMocksBBand
+#fiducial = EH98
+#nowiggles = EH98NoWiggles
+#broadband = EH98BBand
 
-# Model parameter configuration
-# fiducial model with only BAO scale floating
-##model-config = fix[alpha]=3.8; fix[beta]=1.4; fix[(1+beta)*bias]=-0.336; fix[BAO amplitude]; fix[BBand*]=0;
-# perform linear bias fit with no broadband power
-##model-config = fix[BAO*]=1; fix[BBand*]=0;
-# perform linear bias fit with floating broadband power
-##model-config = fix[BAO*]=1;
-# performm BAO fit with fixed beta, bb
-
-#model-config = fix[alpha-bias]=3.8; fix[alpha-beta]=0; fix[beta]=1.4; fix[(1+beta)*bias]=-0.336; fix[BBand1*]; fix[BAO scale *];
-model-config = fix[alpha-bias]=3.8; fix[alpha-beta]=0; fix[beta]=1.4; fix[(1+beta)*bias]=-0.336; fix[BBand1*]; fix[BAO scale];
-anisotropic=yes
-
-# Data to analyze
-plateroot = Method3/data_Aln2.SN/
-platelist = plates.list
+## redshift-space distortion: fix / float
+model-config = fix[beta]=1.4; fix[alpha-beta]=0; fix[(1+beta)*bias]=-0.336; fix[alpha-bias]=3.8;
+#model-config = value[beta]=1.4; value[alpha-beta]=0; value[(1+beta)*bias]=-0.336; value[alpha-bias]=3.8;
+
+## broadband distortion: none (fixed) / BBand1 / BBand2 / BBand12
+#model-config = fix[BBand*]=0;
+model-config = fix[BBand2*]=0;
+#model-config = fix[BBand1*]=0;
+#model-config = fix[BBand2*]=0; error[BBand2 quad*]=0.1;
+
+## BAO amplitude: fix / float
+#model-config = fix[BAO amplitude]=1;
+model-config = value[BAO amplitude]=1;
+
+## BAO scale: isotropic / anisotropic
+model-config = value[BAO scale]=1; fix[BAO scale *];
+#model-config = fix[BAO scale]; value[BAO scale *]=1;
+#anisotropic = yes
+
+## Data to analyze
+plateroot = /Volumes/Data/boss/lya/Method3/delta_pk_Aln2_d3d/
+platelist = mocks_by_realization.list
 weighted = yes
-#reuse-cov = yes
+#reuse-cov = 1
 
-# Cosmolib binning of input correlation data in...
+## Cosmolib binning of input correlation data in...
 # ...redshift
 nz = 3
 minz = 1.75
@@ -42,7 +51,7 @@ maxll = 0.27 # cosmolib's xi3d_llend
 dll = 0.02   # cosmolib's xi3d_llstep
 dll2 = 0.002 # cosmolib's Delta_pix_3D (uses two_step if this is non-zero)
 
-# Final pruning cuts to apply before fitting
+## Final pruning cuts to apply before fitting
 rmin = 20
 rmax = 200
 llmin = -1

diff --git a/config/flux_Aln2_xi.ini b/config/flux_Aln2_xi.ini
@@ -9,15 +9,15 @@ nowiggles = DR9LyaMocksSB
 broadband = DR9LyaMocksBBand
 
 # Model parameter configuration
-xi-model = yes
-num-xi = 14 # number points spanning [rmin,rmax] set below
-model-config = fix[alpha]=3.8; fix[beta]=1.4; fix[(1+beta)*bias]=-0.336;
+xi-points = 20,50,80,95,110,125,140,170,200
+xi-method = linear
+model-config = fix[beta]=1.4; fix[(1+beta)*bias]=-0.336; fix[alpha-*];
 
 # Data to analyze
-plateroot = /Users/david/Cosmo/LyAlpha/cosmolib/mocks/
-platelist = flux_Aln2_1.list
+plateroot = /Volumes/Data/boss/lya/Method3/delta_pk_Aln2_d3d/
+platelist = mocks_all.list
 weighted = yes
-reuse-cov = yes
+reuse-cov = 1
 
 # Cosmolib binning of input correlation data in...
 # ...redshift
@@ -35,9 +35,11 @@ dll = 0.02   # cosmolib's xi3d_llstep
 dll2 = 0.002 # cosmolib's Delta_pix_3D (uses two_step if this is non-zero)
 
 # Final pruning cuts to apply before fitting
-rmin = 50
-rmax = 180
+rmin = 20
+rmax = 200
 llmin = -1
 
 # Bootstrap analysis options
 ##bootstrap-trials = 100
+mcmc-save = 100
+mcmc-interval = 25
diff --git a/models/DR9LyaMocksBBand1.0.dat b/models/DR9LyaMocksBBand1.0.dat
@@ -76,73 +76,73 @@
 80	-0.000142521
 81	-0.0001286498
 82	-0.00011586736
-83	-0.00010552853594877503
-84	-0.00009473457193030779
-85	-0.00008480373822367375
-86	-0.00007567890288170751
-87	-0.00006730680057992811
-88	-0.00005963774300114567
-89	-0.00005262535342359644
-90	-0.00004622632326886853
-91	-0.000040400188594650204
-92	-0.00003510912472067511
-93	-0.000030317757357191912
-94	-0.000025992988766535273
-95	-0.00002210383763222043
-96	-0.00001862129143848814
-97	-0.000015518170278132522
-98	-0.000012769001109328298
-99	-0.000010349901574391051
-100	-8.23847257614952e-6
-101	-6.413698881940922e-6
-102	-4.855857092080015e-6
-103	-3.54643036982189e-6
-104	-2.4680293840539067e-6
-105	-1.604318964854132e-6
-106	-9.399500162070828e-7
-107	-4.604962700657583e-7
-108	-1.52395502044482e-7
-109	-2.894861703284929e-9
+83	-0.00010410822
+84	-0.00009339252
+85	-0.0000837175
+86	-0.00007509024
+87	-0.00006694115318396415
+88	-0.000059281749173460116
+89	-0.000052280380137649414
+90	-0.00004589359684391429
+91	-0.000040080799339152305
+92	-0.00003480403196013978
+93	-0.00003002779474599853
+94	-0.00002571886979894099
+95	-0.00002184616128111578
+96	-0.00001838054786198679
+97	-0.000015294746543954513
+98	-0.000012563186895412654
+99	-0.000010161894811438068
+100	-8.068385004020557e-6
+101	-6.261561497169243e-6
+102	-4.721625468301627e-6
+103	-3.4299898368191996e-6
+104	-2.369200054417415e-6
+105	-1.522860600079658e-6
+106	-8.755667264313989e-7
+107	-4.128410436598194e-7
+108	-1.2107456298108588e-7
+109	1.2528145960791563e-8
 110	0.
-111	-1.3242770334501598e-7
-112	-3.895617680777473e-7
-113	-7.614118713374475e-7
-114	-1.2385752144466675e-6
-115	-1.8122007332590237e-6
-116	-2.4739556866215734e-6
-117	-3.215994449322856e-6
-118	-4.03092934978224e-6
-119	-4.91180340541332e-6
-120	-5.852064820172651e-6
-121	-6.845543119395687e-6
-122	-7.886426806706648e-6
-123	-8.969242436660504e-6
-124	-0.000010088835004898935
-125	-0.000011240349565057978
-126	-0.000012419213988492315
-127	-0.000013621122789167252
-128	-0.00001484202194182448
-129	-0.000016078094626835814
-130	-0.000017325747840027427
-131	-0.000018581599810248504
-132	-0.000019842468171588747
-133	-0.000021105358840955108
-134	-0.000022367455555233616
-135	-0.000023626110025496533
-136	-0.000024878832668709385
-137	-0.000026123283880151605
-138	-0.00002735726581231832
-139	-0.000028578714628429545
-140	-0.000029785693200859267
-141	-0.000030976384226814003
-142	-0.00003214908373546474
-143	-0.000033302194962471223
-144	-0.00003443422256944126
-145	-0.000035543767187365206
-146	-0.000036629520264442494
-147	-0.00003769025920000893
-148	-0.000038568404
-149	-0.000039732186
+111	-1.473410607594119e-7
+112	-4.188477788944194e-7
+113	-8.045022849888098e-7
+114	-1.2948775789830072e-6
+115	-1.881101582420163e-6
+116	-2.5548235730656327e-6
+117	-3.3081828285170093e-6
+118	-4.1337793192388135e-6
+119	-5.024646304073641e-6
+120	-5.9742246928140915e-6
+121	-6.976339050965012e-6
+122	-8.025175131480005e-6
+123	-9.115258827095529e-6
+124	-0.000010241436444989777
+125	-0.000011398856212923115
+126	-0.000012582950932838108
+127	-0.000013789421704159907
+128	-0.000015014222644792038
+129	-0.00001625354654309701
+130	-0.000017503811379016974
+131	-0.000018761647656975727
+132	-0.000020023886497332357
+133	-0.00002128754843696154
+134	-0.00002254983289304832
+135	-0.00002380810824742448
+136	-0.000025059902511764892
+137	-0.000026302894536726454
+138	-0.000027534905730665867
+139	-0.000028753892255937128
+140	-0.00002995793767295459
+141	-0.0000311452460042327
+142	-0.000032314135192486745
+143	-0.000033463030928620106
+144	-0.00003459046082703263
+145	-0.00003569504892718028
+146	-0.000036775510501702725
+147	-0.00003783064715272754
+148	-0.0000388593421791518
+149	-0.00003986055619881953
 150	-0.00004082368
 151	-0.000041851816
 152	-0.000042823955999999995