fix dictionary paths + add option to save diagonal kspace filter for cov computation

Author: Serena Giardiello

project/SO/pISO/paramfiles/global_legacy_deep56.dict

@@ -18,23 +18,23 @@ legacy_map_dir = "/global/cfs/cdirs/cmb/data/planck2018/pr3/frequencymaps/"
 npipe_sim_dir = "/global/cfs/cdirs/cmb/data/planck2020/npipe"
 legacy_sim_dir = "/global/cfs/cdirs/sobs/users/sgiardie/iso/legacy/legacy_sims"
 
-mcm_dir =  data_dir + f"mcms/{run_name}"
+mcm_dir =  data_dir + f"mcms/{run_name}/"
 alms_dir =  data_dir + f"alms/{run_name}/"
 spec_dir =  data_dir + f"spectra/{run_name}/"
 best_fits_dir =  data_dir + f"best_fits"
-plots_dir =  data_dir + f"plots/{run_name}"
-noise_model_dir =  data_dir + f"noise_model/{run_name}"
-sq_win_alms_dir =  data_dir + f"sq_win_alms/{run_name}"
-cov_dir =  data_dir + f"covariances/{run_name}"
-montecarlo_dir =  data_dir + f"montecarlo/{run_name}"
-noise_alms_dir =  data_dir + f"noise_alms/{run_name}"
-chain_dir =  data_dir + f"chains/{run_name}"
-sim_spectra_planck_noise_and_syst_dir =  data_dir + f"sim_spectra_planck_noise_and_syst/{run_name}"
-planck_mc_correction_dir =  data_dir + f"montecarlo/{run_name}"
-spectra_leak_corr_dir =  data_dir + f"spectra_leak_corr/{run_name}"
-spectra_leak_corr_planck_bias_corr_dir =  data_dir + f"spectra_leak_corr_planck_bias_corr/{run_name}"
-montecarlo_beam_leakage_dir =  data_dir + f"montecarlo_beam_leakage/{run_name}"
-sim_spec_dir =  data_dir + f"sim_spectra/dust_binning_50/{run_name}"
+plots_dir =  data_dir + f"plots/{run_name}/"
+noise_model_dir =  data_dir + f"noise_model/{run_name}/"
+sq_win_alms_dir =  data_dir + f"sq_win_alms/{run_name}/"
+cov_dir =  data_dir + f"covariances/{run_name}/"
+montecarlo_dir =  data_dir + f"montecarlo/{run_name}/"
+noise_alms_dir =  data_dir + f"noise_alms/{run_name}/"
+chain_dir =  data_dir + f"chains/{run_name}/"
+sim_spectra_planck_noise_and_syst_dir =  data_dir + f"sim_spectra_planck_noise_and_syst/{run_name}/"
+planck_mc_correction_dir =  data_dir + f"montecarlo/{run_name}/"
+spectra_leak_corr_dir =  data_dir + f"spectra_leak_corr/{run_name}/"
+spectra_leak_corr_planck_bias_corr_dir =  data_dir + f"spectra_leak_corr_planck_bias_corr/{run_name}/"
+montecarlo_beam_leakage_dir =  data_dir + f"montecarlo_beam_leakage/{run_name}/"
+sim_spec_dir =  data_dir + f"sim_spectra/dust_binning_50/{run_name}/"
 
 k_filter_Planck = {'type': 'binary_cross', 'vk_mask': [-90, 90], 'hk_mask': [-50, 50], 'weighted': False}
 kspace_tf_path = "analytical" 

project/SO/pISO/paramfiles/global_npipe_deep56.dict

@@ -19,23 +19,23 @@ npipe_sim_dir = "/global/cfs/cdirs/cmb/data/planck2020/npipe"
 legacy_sim_dir = "/global/cfs/cdirs/sobs/users/sgiardie/iso/legacy/legacy_sims"
 
 
-mcm_dir =  data_dir + f"mcms/{run_name}"
+mcm_dir =  data_dir + f"mcms/{run_name}/"
 alms_dir =  data_dir + f"alms/{run_name}/"
 spec_dir =  data_dir + f"spectra/{run_name}/"
 best_fits_dir =  data_dir + f"best_fits"
-plots_dir =  data_dir + f"plots/{run_name}"
-noise_model_dir =  data_dir + f"noise_model/{run_name}"
-sq_win_alms_dir =  data_dir + f"sq_win_alms/{run_name}"
-cov_dir =  data_dir + f"covariances/{run_name}"
-montecarlo_dir =  data_dir + f"montecarlo/{run_name}"
-noise_alms_dir =  data_dir + f"noise_alms/{run_name}"
-chain_dir =  data_dir + f"chains/{run_name}"
-sim_spectra_planck_noise_and_syst_dir =  data_dir + f"sim_spectra_planck_noise_and_syst/{run_name}"
-planck_mc_correction_dir =  data_dir + f"montecarlo/{run_name}"
-spectra_leak_corr_dir =  data_dir + f"spectra_leak_corr/{run_name}"
-spectra_leak_corr_planck_bias_corr_dir =  data_dir + f"spectra_leak_corr_planck_bias_corr/{run_name}"
-montecarlo_beam_leakage_dir =  data_dir + f"montecarlo_beam_leakage/{run_name}"
-sim_spec_dir =  data_dir + f"sim_spectra/dust_binning_50/{run_name}"
+plots_dir =  data_dir + f"plots/{run_name}/"
+noise_model_dir =  data_dir + f"noise_model/{run_name}/"
+sq_win_alms_dir =  data_dir + f"sq_win_alms/{run_name}/"
+cov_dir =  data_dir + f"covariances/{run_name}/"
+montecarlo_dir =  data_dir + f"montecarlo/{run_name}/"
+noise_alms_dir =  data_dir + f"noise_alms/{run_name}/"
+chain_dir =  data_dir + f"chains/{run_name}/"
+sim_spectra_planck_noise_and_syst_dir =  data_dir + f"sim_spectra_planck_noise_and_syst/{run_name}/"
+planck_mc_correction_dir =  data_dir + f"montecarlo/{run_name}/"
+spectra_leak_corr_dir =  data_dir + f"spectra_leak_corr/{run_name}/"
+spectra_leak_corr_planck_bias_corr_dir =  data_dir + f"spectra_leak_corr_planck_bias_corr/{run_name}/"
+montecarlo_beam_leakage_dir =  data_dir + f"montecarlo_beam_leakage/{run_name}/"
+sim_spec_dir =  data_dir + f"sim_spectra/dust_binning_50/{run_name}/"
 
 k_filter_Planck = {'type': 'binary_cross', 'vk_mask': [-90, 90], 'hk_mask': [-50, 50], 'weighted': False}
 kspace_tf_path = "analytical" 

project/SO/pISO/python/get_pseudo2datavec.py

@@ -22,7 +22,9 @@
 
 mcm_dir = d['mcm_dir']
 plot_dir = opj(d['plots_dir'], 'mcms')
+spec_dir = d["spec_dir"]
 pspy_utils.create_directory(plot_dir)
+pspy_utils.create_directory(spec_dir)
 
 surveys = d["surveys"]
 lmax = d['lmax']
@@ -70,6 +72,13 @@
         if np.count_nonzero(v.diagonal() == 0):
             log.info(f'WARNING: 0 in kspace_transfer_matrix {k}')
 
+    # this will be used in the old covariance computation
+    for spec_name in spec_name_list:
+        one_d_tf = kspace_transfer_matrix[spec_name].diagonal()
+        cov_T_E_only = d["cov_T_E_only"]
+        if cov_T_E_only == True: one_d_tf = one_d_tf[:4 * n_bins]
+        np.savetxt(f"{spec_dir}/one_dimension_kspace_tf_{spec_name}.dat", one_d_tf)
+
 if d[f"pixwin_{sv}"]["pix"] == "HEALPIX" and deconvolve_pixwin:
     pixwins = {}
     for sv in surveys: