current working status

cosmicfishpie/configs/config.py

@@ -209,6 +209,8 @@ def init(
         ),
     )
     settings.setdefault("specs_dir", settings["specs_dir_default"])
+    settings.setdefault("external_data_dir", os.path.join(os.path.dirname(os.path.realpath(__file__)), 
+                                                          "external_data"))
     settings.setdefault("survey_name", surveyName)
     settings.setdefault("survey_specs", "ISTF-Optimistic")
     settings.setdefault("survey_name_photo", "Euclid-Photometric-ISTF-Pessimistic")

cosmicfishpie/configs/other_survey_specifications/Euclid-Photometric-DeboleR1.yaml

@@ -0,0 +1,45 @@
+specifications:
+  'z_bins_GCph': [0.2, 0.39, 0.5, 0.65, 0.83, 1.03, 2.50]
+  'z_bins_WL': [0.2, 0.44, 0.69, 0.91, 1.18, 1.51, 2.50]
+  'zm': 0.9
+  'ngamma': 1.5
+  'ngal_sqarmin_GCph': 10.13
+  'ngal_sqarmin_WL': 35.88
+  photo_z_params:
+    'fout': 0.1
+    'co': 1
+    'cb': 1
+    'sigma_o': 0.05
+    'sigma_b': 0.05
+    'zo': 0.1
+    'zb': 0.0
+  'area_survey_GCph': 2500
+  'lmax_GCph': 750
+  'lmin_GCph': 10
+  'ellipt_error': 0.3
+  'area_survey_WL': 2500
+  'lmax_WL': 3000
+  'lmin_WL': 10
+  'kmax': 30.
+  'vary_ph_bias': 0.06
+  'vary_IA_pars': 0.035
+  ph_bias_model: 'binned'
+  ph_bias_parametrization: 
+    binned:
+      keystr: 'b'
+      generate_bias_keys: true
+    binned_constant:
+      keystr: 'b'
+      generate_bias_keys: true
+    sqrt:
+      b0: 1.0
+    flagship:
+      A: 1.0
+      B: 2.5
+      C: 2.8
+      D: 1.6
+  IA_params:
+    IA_model: 'eNLA'
+    AIA: 1.72
+    betaIA: 2.17
+    etaIA: -0.41

cosmicfishpie/cosmology/nuisance.py

@@ -48,6 +48,7 @@ def __init__(
         self.specs = self.config.specs
         self.settings = self.config.settings
         self.specsdir = self.settings["specs_dir"]
+        self.extdir = self.settings["external_data_dir"]
         self.surveyname = self.specs["survey_name"]
         if "GCsp" in self.observables or "IM" in self.observables:
             self.sp_zbins = self.gcsp_zbins()
@@ -220,9 +221,8 @@ def luminosity_ratio(self):
         """
         try:
             # Lumratio file for IA
-            lum = np.loadtxt(os.path.join(self.specsdir, "lumratio_file.dat"))
+            lum = np.loadtxt(os.path.join(self.extdir, "lumratio_file.dat"))
             lumratio = interp1d(lum[:, 0], lum[:, 1], kind="linear")
-            logger.info("Successfully loaded luminosity ratio file")
         except (FileNotFoundError, OSError) as e:
             logger.warning(f"Could not load luminosity ratio file: {e}. Using default value of 1.0")