modify fit_dust code

Author: Serena Giardiello

project/SO/pISO/python/dust/fit_dust_amplitude.py

@@ -1,6 +1,6 @@
 """
-This script uses 143/353 GHz spectra from Planck to fit dust amplitude within ACT survey
-if --use-220 is used , we will also fit the high ell 220 GHz ACT channel, the idea being to try to break the degeneracy between dust and CIB
+This script uses 143/353 GHz spectra from Planck to fit dust amplitude within the deep56 patch
+if --use-220 is used , we will also fit the high ell 220 GHz ACT channel, the idea being to try to break the degeneracy between dust and CIB. We will not use this in the first round of dust fits for ISO, will modify the relevant code later on if needed.
 example use:
 python fit_dust_amplitude.py global_dust.dict --mode BB
 python fit_dust_amplitude.py global_dust.dict --mode TT --use-220  --dr6-result-path ./dr6
@@ -160,9 +160,9 @@ def get_fg_model(a_p, a_c, a_gtt, a_gte, a_gee, a_gbb, a_gtb, bandint_shift_dr6_
 Rminus1_cl_stop = 0.1
 mc_cov = True
 
-result_dir = f"chains/dust_from_planck353_{mode}" # YOU MAY WANT TO CHANGE TO FULL PATH, AND/OR FROM PARAMFILE
-chain_name = f"{result_dir}/dust" # YOU MAY WANT TO CHANGE TO FULL PATH, AND/OR FROM PARAMFILE
-plot_dir = f"plots/dust_from_planck353_{mode}" # YOU MAY WANT TO CHANGE TO FULL PATH, AND/OR FROM PARAMFILE
+result_dir = d["chain_dir"] + f"/dust_from_planck353_{mode}" 
+chain_name = f"{result_dir}/dust" 
+plot_dir = d["plots_base_dir"] + f"/dust_from_planck353_{mode}" 
 
 if use_220:
     chain_name += "_with_220"
@@ -186,8 +186,8 @@ def get_fg_model(a_p, a_c, a_gtt, a_gte, a_gee, a_gbb, a_gtb, bandint_shift_dr6_
 lmin, lmax = 300, 2000
 idx_planck = np.where((bin_low >= lmin) & (bin_high <= lmax))[0]
 lb, res_planck, cov_res_planck = dust_utils.get_residual_and_cov(["Planck_f143", "Planck_f353"],
-                                                                 "spectra", # YOU MAY WANT TO CHANGE TO FULL PATH, AND/OR FROM PARAMFILE
-                                                                 "covariances", # YOU MAY WANT TO CHANGE TO FULL PATH, AND/OR FROM PARAMFILE
+                                                                 d["spec_dir"], 
+                                                                 d["cov_dir"], 
                                                                  mode,
                                                                  spectra,
                                                                  op="aa+bb-2ab",

project/SO/pISO/python/dust/global_dust.dict

@@ -7,7 +7,7 @@ binning_file = release_dir + "pspipe/binning/binning_50"
 apply_kspace_filter  = True
 k_filter_Planck = {"type":'binary_cross', "vk_mask":[-90, 90], "hk_mask":[-50, 50], "weighted":False}
 kspace_tf_path  = release_dir + "pspipe/kspace_tf/dust_binning_50/"
-sim_spec_dir = release_dir + "pspipe/sim_spectra/dust_binning_50/"
+sim_spec_dir = "/global/cfs/cdirs/sobs/users/sgiardie/iso/dust/sim_spectra/dust_binning_50/"
 source_catalog = release_dir + "pspipe/source_catalog/cat_skn_multifreq_20220526_nightonly.txt"
 
 mcm_dir = "/scratch/c.spxsg6/iso/dust/mcms"
@@ -21,7 +21,7 @@ cov_dir = "/scratch/c.spxsg6/iso/dust/covariances"
 montecarlo_dir = "/scratch/c.spxsg6/iso/dust/montecarlo"
 noise_alms_dir = "/scratch/c.spxsg6/iso/dust/noise_alms"
 catalog_dir = "/scratch/c.spxsg6/iso/dust/catalogs"
-
+chain_dir = "/scratch/c.spxsg6/iso/dust/chains"
 
 deconvolve_pixwin = True
 pixwin_Planck = {"pix": 'HEALPIX', "nside": 2048}