Merge pull request #1 from SBU-COSMOLIKE/main

v3.1
CosmoLike · Aug 24, 2023 · b8d661f · b8d661f
2 parents 1a4c158 + 52ae290
commit b8d661f
Show file tree

Hide file tree

Showing 67 changed files with 16,434 additions and 6,872 deletions.
diff --git a/EXAMPLE_EVALUATE1.yaml b/EXAMPLE_EVALUATE1.yaml
@@ -5,12 +5,20 @@ stop_at_error: True
 likelihood:
   lsst_y1.lsst_cosmic_shear:
     path: ./external_modules/data/lsst_y1
+    data_file: LSST_Y1_M1_GGL0.05.dataset  # 705 non-masked elements  (EE2 delta chi^2 ~ 11.8)
+    #data_file: LSST_Y1_M2_GGL0.05.dataset  # 600 non-masked elements  (EE2 delta chi^2 ~ 6.7)
+    #data_file: LSST_Y1_M3_GGL0.05.dataset  # 495 non-masked elements  (EE2 delta chi^2 ~ 3.2)
+    #data_file: LSST_Y1_M4_GGL0.05.dataset  # 405 non-masked elements  (EE2 delta chi^2 ~ 1.4)
+    #data_file: LSST_Y1_M5_GGL0.05.dataset  # 315 non-masked elements  (EE2 delta chi^2 ~ 0.5)
+    #data_file: LSST_Y1_M6_GGL0.05.dataset  # 780 non-masked elements  (EE2 delta chi^2 ~ 12.96)
     print_datavector: False
     print_datavector_file: "./projects/lsst_y1/chains/example1_lsst_y1_theory.modelvector"
     accuracyboost: 1.0
     integration_accuracy: 0
     kmax_boltzmann: 5.0
-
+    # 1 EE2, 2 Halofit
+    non_linear_emul: 2
+
 params:
   As_1e9:
     prior:
@@ -153,7 +161,7 @@ theory:
     use_renames: True
     extra_args:
       halofit_version: takahashi
-      AccuracyBoost: 1.05
+      AccuracyBoost: 1.1
       lens_potential_accuracy: 1.0
       num_massive_neutrinos: 1
       nnu: 3.046

diff --git a/EXAMPLE_EVALUATE2.yaml b/EXAMPLE_EVALUATE2.yaml
@@ -5,19 +5,29 @@ stop_at_error: True
 likelihood:
   lsst_y1.lsst_3x2pt:
     path: ./external_modules/data/lsst_y1
+    data_file: LSST_Y1_M1_GGL0.05.dataset 
+    #data_file: LSST_Y1_M2_GGL0.05.dataset  
+    #data_file: LSST_Y1_M3_GGL0.05.dataset  
+    #data_file: LSST_Y1_M4_GGL0.05.dataset  
+    #data_file: LSST_Y1_M5_GGL0.05.dataset  
+    #data_file: LSST_Y1_M6_GGL0.05.dataset 
     print_datavector: False
-    #print_datavector_file: "./projects/lsst_y1/chains/example2_lsst_y1_theory.modelvector"
     accuracyboost: 1.0
     integration_accuracy: 0 
     kmax_boltzmann: 5.0
+    # 1 EE2, 2 Halofit
+    non_linear_emul: 2
     # ------------------------------------------------------------------------
     # Settings adopted when printing the lsst_y1_theory.modelvector
     # ------------------------------------------------------------------------
+    #data_file: LSST_Y1_M6_GGL0.05.dataset
     #print_datavector: True
     #print_datavector_file: "./projects/lsst_y1/chains/lsst_y1_theory.modelvector"
-    #accuracyboost: 1.0
-    #integration_accuracy: 1 
+    #accuracyboost: 3.0
+    #integration_accuracy: 10 
     #kmax_boltzmann: 20.0
+    # 1 EE2, 2 Halofit
+    #non_linear_emul: 2
 
 params:
   As_1e9:
@@ -152,6 +162,23 @@ params:
   omegan:
     derived: 'lambda omegan2, H0: omegan2/((H0/100)**2)'
     latex: \Omega_\mathrm{\\nu}
+# FIXED LSST PARAMETERS
+  LSST_BMAG_1:
+    value: 0.0
+    latex: b_\mathrm{BMAG-LSST}^1
+  LSST_BMAG_2:
+    value: 0.0
+    latex: b_\mathrm{BMAG-LSST}^2
+  LSST_BMAG_3:
+    value: 0.0
+    latex: b_\mathrm{BMAG-LSST}^3
+  LSST_BMAG_4:
+    value: 0.0
+    latex: b_\mathrm{BMAG-LSST}^4
+  LSST_BMAG_5:
+    value: 0.0
+    latex: b_\mathrm{BMAG-LSST}^5
+
 
 theory:
   camb:
@@ -160,17 +187,14 @@ theory:
     use_renames: True
     extra_args:
       halofit_version: takahashi
+      AccuracyBoost: 1.1
       lens_potential_accuracy: 1.0
       num_massive_neutrinos: 1
       nnu: 3.046
       dark_energy_model: ppf
       accurate_massive_neutrino_transfers: false
       k_per_logint: 20
-      AccuracyBoost: 1.05
-      # ------------------------------------------------------------------------
-      # Settings adopted when printing the lsst_y1_theory.modelvector
-      # ------------------------------------------------------------------------
-      #AccuracyBoost: 1.15
+      #AccuracyBoost: 2.0
 
 sampler:
   evaluate:
@@ -198,7 +222,7 @@ sampler:
       LSST_DZ_L2: 0.000309875
       LSST_DZ_L3: 0.00855907
       LSST_DZ_L4: -0.00316269
-      LSST_DZ_L5: -0.0146753
+      LSST_DZ_L5: -0.0146753 
       LSST_B1_1: 1.72716
       LSST_B1_2: 1.65168
       LSST_B1_3: 1.61423

diff --git a/EXAMPLE_EVALUATE3.yaml b/EXAMPLE_EVALUATE3.yaml
@@ -9,7 +9,9 @@ likelihood:
     accuracyboost: 1.0
     integration_accuracy: 0
     kmax_boltzmann: 5.0
-
+    # 1 EE2, 2 Halofit
+    non_linear_emul: 2
+
 params:
   As_1e9:
     prior:

diff --git a/EXAMPLE_EVALUATE4.yaml b/EXAMPLE_EVALUATE4.yaml
@@ -9,7 +9,9 @@ likelihood:
     accuracyboost: 1.0
     integration_accuracy: 0
     kmax_boltzmann: 5.0
-
+    # 1 EE2, 2 Halofit
+    non_linear_emul: 2
+
 params:
   As_1e9:
     prior:

diff --git a/EXAMPLE_EVALUATE5.yaml b/EXAMPLE_EVALUATE5.yaml
@@ -9,7 +9,9 @@ likelihood:
     accuracyboost: 1.0
     integration_accuracy: 0
     kmax_boltzmann: 5.0
-
+    # 1 EE2, 2 Halofit
+    non_linear_emul: 2
+
 params:
   As_1e9:
     prior:

diff --git a/EXAMPLE_MCMC1.covmat b/EXAMPLE_MCMC1.covmat
@@ -0,0 +1,19 @@
+# As_1e9 ns H0 omegab omegam w LSST_DZ_S1 LSST_DZ_S2 LSST_DZ_S3 LSST_DZ_S4 LSST_DZ_S5 LSST_A1_1 LSST_A1_2 LSST_M1 LSST_M2 LSST_M3 LSST_M4 LSST_M5
+8.147056434834848426e-02 4.250444435033904657e-03 -1.946004033369889141e+00 -1.621784588170465053e-03 6.569860462970322085e-04 1.257382254757036591e-02 -6.514781257094201877e-06 -3.173247923400945088e-05 -8.963020365848351420e-06 2.723323747596670863e-05 1.286809047296887000e-05 -6.182013764548815841e-03 -4.161936073088944726e-02 -3.037141872562700462e-05 -8.441536121298864859e-05 -4.532632193962631744e-05 3.738689725912665491e-05 8.162022253457763301e-05
+4.250444435033904657e-03 1.395610298376724189e-03 -1.935747802432359888e-01 -1.208874930062493736e-04 -1.656492325073450819e-04 -3.518466284966141003e-03 1.063589463192090629e-06 1.376674000082042774e-06 -1.304833964955558142e-06 -1.419360813318523520e-06 6.134637304927333259e-07 8.861829888951925195e-04 1.244595171303405521e-03 2.480374504287038556e-06 1.585004447932609113e-06 -6.324942628166939441e-06 -1.117980793124111911e-05 6.756030743808678041e-06
+-1.946004033369889141e+00 -1.935747802432359888e-01 7.593365706488788192e+01 6.660845623087062461e-02 -4.398076177042117840e-02 -7.172148972691005930e-02 3.762927466548532558e-04 8.469106372599704679e-04 6.280585332831347973e-04 -1.244249542540311945e-04 -4.260436939618453180e-04 2.421646621564214896e-02 2.209873938491425938e-01 1.258523112864857324e-03 2.574008940964310835e-03 2.817004038729869694e-03 7.637679288374588180e-04 -5.538935403315359447e-03
+-1.621784588170465053e-03 -1.208874930062493736e-04 6.660845623087062461e-02 1.099930910443445373e-04 3.986676156670956433e-05 5.290860905449195533e-05 -1.357033377808488756e-07 -5.217070829260963597e-07 -6.317030974096768386e-07 -4.066099155546820885e-07 4.903819784784667332e-07 1.953868488418605124e-05 9.421725313187949414e-05 -5.609264239707860776e-07 -1.479691643120089858e-06 -1.402818145767477712e-06 -2.573176357883316328e-07 4.639809768378610909e-06
+6.569860462970322085e-04 -1.656492325073450819e-04 -4.398076177042117840e-02 3.986676156670956433e-05 3.666916581700847746e-04 1.875813482051728659e-03 -2.681560719919659457e-06 -5.062639327422614876e-06 -3.702270732728634503e-06 -1.133579047104119779e-06 2.655507920246068324e-06 -1.840575917286123742e-04 5.474513929148928015e-04 -7.932368041742284689e-06 -1.501007636135896810e-05 -1.279564303242187595e-05 -7.603303715069855376e-07 3.525945247520988353e-05
+1.257382254757036591e-02 -3.518466284966141003e-03 -7.172148972691005930e-02 5.290860905449195533e-05 1.875813482051728659e-03 2.262729709868453037e-02 -1.480868957396114730e-05 -3.388503225471228641e-05 -8.696836200988197261e-06 9.787700349773876198e-06 1.548438486195624468e-05 -5.034895673604316592e-03 -1.579702037317391727e-02 -4.455498499923921715e-05 -8.843223806756485771e-05 -2.463397255434189652e-05 5.926323694008239356e-05 1.628103585465690028e-04
+-6.514781257094201877e-06 1.063589463192090629e-06 3.762927466548532558e-04 -1.357033377808488756e-07 -2.681560719919659457e-06 -1.480868957396114730e-05 3.745896751730360265e-06 3.392143010809439666e-07 -2.614123498590056113e-09 -9.397337279913465347e-08 -2.697920122504008479e-08 -1.602302074233713755e-06 -3.294090302001190067e-04 -4.457452854321212434e-07 5.314889440720499303e-07 2.019267890999811298e-07 8.042057305909256749e-08 -3.875011038586829060e-07
+-3.173247923400945088e-05 1.376674000082042774e-06 8.469106372599704679e-04 -5.217070829260963597e-07 -5.062639327422614876e-06 -3.388503225471228641e-05 3.392143010809439666e-07 3.317838019184899857e-06 3.149186175161560198e-07 1.743397516059241798e-07 -9.709345073427335288e-08 2.328035583147670113e-05 1.557671925905976728e-04 3.692300935113772405e-07 -1.144357356875239888e-06 6.444042628104154237e-07 3.535510988316960909e-07 -3.481052607250834130e-07
+-8.963020365848351420e-06 -1.304833964955558142e-06 6.280585332831347973e-04 -6.317030974096768386e-07 -3.702270732728634503e-06 -8.696836200988197261e-06 -2.614123498590056113e-09 3.149186175161560198e-07 3.491052490727321512e-06 2.662315332280460696e-07 3.123110797441771725e-08 8.820052920498775299e-06 1.295076492768112758e-04 7.181289791655983966e-08 5.192720063828188049e-07 -1.369457362174064736e-06 5.794234387741672692e-07 1.298695340157753149e-07
+2.723323747596670863e-05 -1.419360813318523520e-06 -1.244249542540311945e-04 -4.066099155546820885e-07 -1.133579047104119779e-06 9.787700349773876198e-06 -9.397337279913465347e-08 1.743397516059241798e-07 2.662315332280460696e-07 3.630630936216653247e-06 9.685998380315394033e-08 1.542954609794885866e-06 8.094082404904401568e-05 -6.295560736855802647e-09 2.532532214356249740e-07 5.456114338635830018e-07 -1.308915476765712806e-06 5.892544166965478468e-07
+1.286809047296887000e-05 6.134637304927333259e-07 -4.260436939618453180e-04 4.903819784784667332e-07 2.655507920246068324e-06 1.548438486195624468e-05 -2.697920122504008479e-08 -9.709345073427335288e-08 3.123110797441771725e-08 9.685998380315394033e-08 3.931127619020217686e-06 -1.155349565554091907e-05 -9.258458850131971772e-05 9.780403635346436906e-09 -5.749969029409392803e-08 9.760597382413269880e-08 3.281460782297497434e-07 -3.779833131714117822e-07
+-6.182013764548815841e-03 8.861829888951925195e-04 2.421646621564214896e-02 1.953868488418605124e-05 -1.840575917286123742e-04 -5.034895673604316592e-03 -1.602302074233713755e-06 2.328035583147670113e-05 8.820052920498775299e-06 1.542954609794885866e-06 -1.155349565554091907e-05 5.633427949569455064e-03 3.680120827936887024e-02 -6.856867302296777180e-06 9.228786642592550458e-06 -9.185716445689369870e-06 -3.124616192443146958e-05 5.225148099304805070e-06
+-4.161936073088944726e-02 1.244595171303405521e-03 2.209873938491425938e-01 9.421725313187949414e-05 5.474513929148928015e-04 -1.579702037317391727e-02 -3.294090302001190067e-04 1.557671925905976728e-04 1.295076492768112758e-04 8.094082404904401568e-05 -9.258458850131971772e-05 3.680120827936887024e-02 4.108061244573737447e-01 -2.658533885315810712e-04 8.278955173626539529e-05 3.568001267738926658e-05 -1.068328260022480280e-04 5.979694946226428453e-05
+-3.037141872562700462e-05 2.480374504287038556e-06 1.258523112864857324e-03 -5.609264239707860776e-07 -7.932368041742284689e-06 -4.455498499923921715e-05 -4.457452854321212434e-07 3.692300935113772405e-07 7.181289791655983966e-08 -6.295560736855802647e-09 9.780403635346435252e-09 -6.856867302296777180e-06 -2.658533885315810712e-04 2.420181950125698599e-05 8.609914435324823678e-07 6.335513729061734817e-07 2.590337076277333225e-07 -8.198036456101131296e-07
+-8.441536121298864859e-05 1.585004447932609113e-06 2.574008940964310835e-03 -1.479691643120089858e-06 -1.501007636135896810e-05 -8.843223806756485771e-05 5.314889440720499303e-07 -1.144357356875239888e-06 5.192720063828188049e-07 2.532532214356249740e-07 -5.749969029409392803e-08 9.228786642592550458e-06 8.278955173626539529e-05 8.609914435324823678e-07 2.166578982468595320e-05 1.916179168234973924e-06 1.267205659521347302e-06 -1.034049099167444953e-06
+-4.532632193962631744e-05 -6.324942628166939441e-06 2.817004038729869694e-03 -1.402818145767477712e-06 -1.279564303242187595e-05 -2.463397255434189652e-05 2.019267890999811298e-07 6.444042628104154237e-07 -1.369457362174064736e-06 5.456114338635830018e-07 9.760597382413271204e-08 -9.185716445689369870e-06 3.568001267738926658e-05 6.335513729061734817e-07 1.916179168234973924e-06 1.785276143788537746e-05 3.047284087323038894e-06 1.803289887217451155e-06
+3.738689725912666169e-05 -1.117980793124111911e-05 7.637679288374588180e-04 -2.573176357883315799e-07 -7.603303715069854317e-07 5.926323694008239356e-05 8.042057305909256749e-08 3.535510988316960909e-07 5.794234387741672692e-07 -1.308915476765712806e-06 3.281460782297497434e-07 -3.124616192443146958e-05 -1.068328260022480280e-04 2.590337076277333225e-07 1.267205659521347302e-06 3.047284087323038894e-06 1.625734873801078332e-05 4.768116354103480280e-06
+8.162022253457763301e-05 6.756030743808678041e-06 -5.538935403315359447e-03 4.639809768378610909e-06 3.525945247520988353e-05 1.628103585465690028e-04 -3.875011038586829060e-07 -3.481052607250834130e-07 1.298695340157753149e-07 5.892544166965478468e-07 -3.779833131714117822e-07 5.225148099304805070e-06 5.979694946226427098e-05 -8.198036456101131296e-07 -1.034049099167444953e-06 1.803289887217451155e-06 4.768116354103480280e-06 1.976127575503067640e-05