modify paramfiles

Author: sgiardie

project/SO/pISO/paramfiles/global_dr6v4xlegacy.dict

@@ -1,21 +1,47 @@
+run_name = "dr6xlegacy"
+data_dir = "/scratch/c.spxsg6/dr6xlegacy/"
+
 surveys = ["dr6", "Planck"]
 arrays_dr6 = ["pa4_f220", "pa5_f090", "pa5_f150", "pa6_f090", "pa6_f150"]
 arrays_Planck = ["f100", "f143", "f217"]
 planck_version = "legacy"
 
+# path to chain of dr6 lcdm (from which we extract the covmat of FG params)
 p_act_chain_filename = "/scratch/c.spxsg6/chains/p-actbase_lcdm_camb/p-actbase_lcdm_camb"
 
-release_dir = "/scratch/c.spxsg6/dr6/"
+release_dir = "/global/cfs/cdirs/cmb/data/act_dr6/dr6.02/"
 
 binning_file = release_dir + "binning/binning_50"
 apply_kspace_filter  = True
 k_filter_dr6 = {"type":'binary_cross', "vk_mask":[-90, 90], "hk_mask":[-50, 50], "weighted":False}
 k_filter_Planck = {"type":'binary_cross', "vk_mask":[-90, 90], "hk_mask":[-50, 50], "weighted":False}
-kspace_tf_path  = release_dir + "kspace_dr6xplanck_binning_50/"
-sim_spec_dir = release_dir + f"dr6xlegacy_newbestfit/sim_spectra/dr6xplanck/{planck_version}/"
-planck_mc_correction = release_dir + f"dr6xlegacy_newbestfit/for_planck/mc_correction/{planck_version}"
+
+kspace_tf_path  = release_dir  + "pspipe/kspace_tf/dr6xplanck_binning_50/"
 source_catalog = release_dir + "pspipe/source_catalog/cat_skn_multifreq_20220526_nightonly.txt"
 
+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/{run_name}"
+plots_base_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}"
+catalog_dir = data_dir + f"catalogs/{run_name}"
+sim_spec_dir = data_dir + f"sim_spectra/dust_binning_50/{run_name}"
+planck_projected_dir = data_dir + f"planck_projected/{run_name}"
+sim_spec_dir = data_dir + f"sim_spectra/{run_name}"
+
+# if you want to compute the leakage contribution and end-to-end sim corrections
+
+sim_spectra_planck_noise_and_syst_dir = data_dir + f"sim_spectra_planck_noise_and_syst/{run_name}"
+noise_alms_dir = data_dir + f"noise_alms/{run_name}"
+planck_mc_correction_dir = data_dir + f"montecarlo_planck/{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}"
+
 deconvolve_pixwin = True
 pixwin_dr6 = {"pix": "CAR", "order": 0}
 pixwin_Planck = {"pix": 'HEALPIX', "nside": 2048}
@@ -46,9 +72,9 @@ maps_dr6_pa6_f090 = [release_dir + "maps/act_dr6.02_std_AA_night_pa6_f090_4way_s
 maps_dr6_pa6_f150 = [release_dir + "maps/act_dr6.02_std_AA_night_pa6_f150_4way_set%d_map_srcfree.fits" % (i) for i in range(4)]
 
 # These are produced by the pipeline
-maps_Planck_f100 = [release_dir + "maps/HFI_SkyMap_2048_R3.01_halfmission-%d_f100_map_srcfree.fits" % i for i in [1,2]]
-maps_Planck_f143 = [release_dir + "maps/HFI_SkyMap_2048_R3.01_halfmission-%d_f143_map_srcfree.fits" % i for i in [1,2]]
-maps_Planck_f217 = [release_dir + "maps/HFI_SkyMap_2048_R3.01_halfmission-%d_f217_map_srcfree.fits" % i for i in [1,2]]
+maps_Planck_f100 = [ data_dir + f"{planck_version}/maps/HFI_SkyMap_2048_R3.01_halfmission-%d_f100_map_srcfree.fits" % i for i in [1,2]]
+maps_Planck_f143 = [ data_dir + f"{planck_version}/maps/HFI_SkyMap_2048_R3.01_halfmission-%d_f143_map_srcfree.fits" % i for i in [1,2]]
+maps_Planck_f217 = [ data_dir + f"{planck_version}/maps/HFI_SkyMap_2048_R3.01_halfmission-%d_f217_map_srcfree.fits" % i for i in [1,2]]
 
 
 cal_dr6_pa4_f150 = 1.0
@@ -77,19 +103,19 @@ pol_eff_Planck_f217 = 0.9817
 
 do_bandpass_integration = True
 
-passband_dir_dr6 = release_dir + "passbands/"
+passband_dir_dr6 = release_dir + "passbands/processed/"
 freq_info_dr6_pa4_f220 = {"freq_tag": 220, "passband": passband_dir_dr6 + "passband_dr6_pa4_f220.dat"}
 freq_info_dr6_pa5_f090 = {"freq_tag": 90, "passband": passband_dir_dr6 + "passband_dr6_pa5_f090.dat"}
 freq_info_dr6_pa5_f150 = {"freq_tag": 150, "passband": passband_dir_dr6 + "passband_dr6_pa5_f150.dat"}
 freq_info_dr6_pa6_f090 = {"freq_tag": 90, "passband": passband_dir_dr6 + "passband_dr6_pa6_f090.dat"}
 freq_info_dr6_pa6_f150 = {"freq_tag": 150, "passband": passband_dir_dr6 + "passband_dr6_pa6_f150.dat"}
 
-passband_dir_planck = release_dir + "passbands/"
+passband_dir_planck = release_dir + "pspipe/for_planck/passbands/"
 freq_info_Planck_f100 = {"freq_tag": 100, "passband": passband_dir_planck + "passband_npipe_f100.dat"}
 freq_info_Planck_f143 = {"freq_tag": 143, "passband": passband_dir_planck + "passband_npipe_f143.dat"}
 freq_info_Planck_f217 = {"freq_tag": 217, "passband": passband_dir_planck + "passband_npipe_f217.dat"}
 
-beam_file_dir_dr6 = release_dir + "dr6xlegacy_newbestfit/beams/"
+beam_file_dir_dr6 = release_dir + "beams/main_beams/nominal/"
 beam_T_dr6_pa4_f220 = beam_file_dir_dr6 + "coadd_pa4_f220_night_beam_tform_jitter_cmb.txt"
 beam_T_dr6_pa5_f090 = beam_file_dir_dr6 + "coadd_pa5_f090_night_beam_tform_jitter_cmb.txt"
 beam_T_dr6_pa5_f150 = beam_file_dir_dr6 + "coadd_pa5_f150_night_beam_tform_jitter_cmb.txt"
@@ -109,20 +135,20 @@ beam_mono_dr6_pa5_f150 = beam_file_dir_dr6 + 'coadd_pa5_f150_night_beam_tform_ji
 beam_mono_dr6_pa6_f090 = beam_file_dir_dr6 + 'coadd_pa6_f090_night_beam_tform_jitter_mono.txt'
 beam_mono_dr6_pa6_f150 = beam_file_dir_dr6 + 'coadd_pa6_f150_night_beam_tform_jitter_mono.txt'
 
-planck_fits_beam_path = release_dir + "dr6xlegacy_newbestfit/for_planck/beams/"
+planck_fits_beam_path = data_dir + f"{run_name}"
 
 # These are produced by the pipeline
-beam_T_Planck_f100 = "beams/legacy/bl_T_legacy_100hm1x100hm2.dat"
-beam_T_Planck_f143 = "beams/legacy/bl_T_legacy_143hm1x143hm2.dat"
-beam_T_Planck_f217 = "beams/legacy/bl_T_legacy_217hm1x217hm2.dat"
+beam_T_Planck_f100 = planck_fits_beam_path + "/beams/legacy/bl_T_legacy_100hm1x100hm2.dat"
+beam_T_Planck_f143 = planck_fits_beam_path + "/beams/legacy/bl_T_legacy_143hm1x143hm2.dat"
+beam_T_Planck_f217 = planck_fits_beam_path + "/beams/legacy/bl_T_legacy_217hm1x217hm2.dat"
 
-beam_pol_Planck_f100 = "beams/legacy/bl_pol_legacy_100hm1x100hm2.dat"
-beam_pol_Planck_f143 = "beams/legacy/bl_pol_legacy_143hm1x143hm2.dat"
-beam_pol_Planck_f217 = "beams/legacy/bl_pol_legacy_217hm1x217hm2.dat"
+beam_pol_Planck_f100 = planck_fits_beam_path + "/beams/legacy/bl_pol_legacy_100hm1x100hm2.dat"
+beam_pol_Planck_f143 = planck_fits_beam_path + "/beams/legacy/bl_pol_legacy_143hm1x143hm2.dat"
+beam_pol_Planck_f217 = planck_fits_beam_path + "/beams/legacy/bl_pol_legacy_217hm1x217hm2.dat"
 
 include_beam_chromaticity_effect_in_best_fit = True
 
-leakage_file_dir_dr6 = release_dir + "dr6xlegacy_newbestfit/beams/"
+leakage_file_dir_dr6 = release_dir + "beams/leakage_beams/nominal/"
 
 leakage_beam_dr6_pa4_f220_TE = [leakage_file_dir_dr6 + "pa4_f220_gamma_t2e.txt" for i in range(4)]
 leakage_beam_dr6_pa5_f090_TE = [leakage_file_dir_dr6 + "pa5_f090_gamma_t2e.txt" for i in range(4)]
@@ -146,7 +172,7 @@ leakage_beam_Planck_f143_TB = ["beams/legacy/gamma_legacy_143%s_t2b.dat" % split
 leakage_beam_Planck_f217_TB = ["beams/legacy/gamma_legacy_217%s_t2b.dat" % split for split in ["hm1", "hm2"]]
 
 
-mask_dir = release_dir + "masks/"
+mask_dir = release_dir + "pspipe/masks/"
 
 ps_mask_dr6_pa4_f220 = mask_dir + "source_mask_15mJy_and_dust_rad12.fits"
 ps_mask_dr6_pa5_f090 = mask_dir + "source_mask_15mJy_and_dust_rad12.fits"
@@ -187,7 +213,7 @@ cross_link_threshold = 0.97
 n_med_ivar = 3
 
 
-window_dir = release_dir + "dr6xlegacy_newbestfit/windows/"
+window_dir = release_dir + "pspipe/windows/dr6xplanck/"
 
 window_kspace_dr6_pa4_f220 = window_dir + "window_dr6_pa4_f220_kspace.fits"
 window_T_dr6_pa4_f220 =  window_dir +"window_dr6_pa4_f220_baseline.fits"