new version of pspipe

project/data_analysis/paramfiles/global_dr6_v3_4pass.dict

@@ -2,16 +2,17 @@ surveys = ["dr6"]
 
 arrays_dr6 = ["pa4_f150", "pa4_f220", "pa5_f090", "pa5_f150", "pa6_f090", "pa6_f150"]
 
+
 data_dir = '/project/projectdirs/act/data/tlouis/s17s18s19/'
 map_dir = '/project/projectdirs/act/data/tlouis/s17s18s19/maps_dr6_v3/4pass/'
 npass_dr6 = 4
 
 deconvolve_pixwin = True
-binning_file = data_dir + "binning/binning_dr6_small_bin.dat"
+binning_file = data_dir + "binning/BIN_ACTPOL_50_4_SC_large_bin_at_low_ell"
 niter = 0
 remove_mean = False
-binned_mcm = True
-lmax = 7000
+binned_mcm = False
+lmax = 5400
 type = 'Dl'
 write_splits_spectra = True
 multistep_path = data_dir
@@ -41,8 +42,9 @@ cross_link_threshold = 0.97
 n_med_ivar = 3
 
 # kspace filter parameters
-
-k_filter_dr6 = {"apply":True, "type":"binary_cross","vk_mask":[-90, 90], "hk_mask":[-50, 50], "weighted":False, "tf": "analytic"}
+apply_kspace_filter  = True
+kspace_tf_path  = "/project/projectdirs/act/data/tlouis/s17s18s19/transfer_fcns/transfer_functions_choi_binning"
+k_filter_dr6 = {"type":"binary_cross","vk_mask":[-90, 90], "hk_mask":[-50, 50], "weighted":False}
 
 deconvolve_map_maker_tf_dr6 = False
 
@@ -80,19 +82,20 @@ pol_eff_dr6_pa5_f150 = 1.024
 pol_eff_dr6_pa6_f090 = 1.0
 pol_eff_dr6_pa6_f150 = 1.0
 
+
 nu_eff_dr6_pa4_f150 = 150
 nu_eff_dr6_pa4_f220 = 220
 nu_eff_dr6_pa5_f090 = 90
 nu_eff_dr6_pa5_f150 = 150
 nu_eff_dr6_pa6_f090 = 90
 nu_eff_dr6_pa6_f150 = 150
 
-beam_dr6_pa4_f150 = data_dir + 'beams/coadd_pa4_f150_night_beam_tform_jitter_cmb.txt'
-beam_dr6_pa4_f220 = data_dir + 'beams/coadd_pa4_f220_night_beam_tform_jitter_cmb.txt'
-beam_dr6_pa5_f090 = data_dir + 'beams/coadd_pa5_f090_night_beam_tform_jitter_cmb.txt'
-beam_dr6_pa5_f150 = data_dir + 'beams/coadd_pa5_f150_night_beam_tform_jitter_cmb.txt'
-beam_dr6_pa6_f090 = data_dir + 'beams/coadd_pa6_f090_night_beam_tform_jitter_cmb.txt'
-beam_dr6_pa6_f150 = data_dir + 'beams/coadd_pa6_f150_night_beam_tform_jitter_cmb.txt'
+beam_dr6_pa4_f150 = data_dir + 'beams/20220517_beams/coadd_pa4_f150_night_beam_tform_jitter_cmb.txt'
+beam_dr6_pa4_f220 = data_dir + 'beams/20220517_beams/coadd_pa4_f220_night_beam_tform_jitter_cmb.txt'
+beam_dr6_pa5_f090 = data_dir + 'beams/20220517_beams/coadd_pa5_f090_night_beam_tform_jitter_cmb.txt'
+beam_dr6_pa5_f150 = data_dir + 'beams/20220517_beams/coadd_pa5_f150_night_beam_tform_jitter_cmb.txt'
+beam_dr6_pa6_f090 = data_dir + 'beams/20220517_beams/coadd_pa6_f090_night_beam_tform_jitter_cmb.txt'
+beam_dr6_pa6_f150 = data_dir + 'beams/20220517_beams/coadd_pa6_f150_night_beam_tform_jitter_cmb.txt'
 
 window_T_dr6_pa4_f150 =  "windows/window_dr6_pa4_f150.fits"
 window_pol_dr6_pa4_f150 =  "windows/window_dr6_pa4_f150.fits"
@@ -116,13 +119,16 @@ window_pol_dr6_pa6_f150 =  "windows/window_dr6_pa6_f150.fits"
 # best fit params (only used for sim generation and covariances computation)
 cosmo_params = {"cosmomc_theta":0.0104085, "logA": 3.044, "ombh2": 0.02237, "omch2": 0.1200, "ns": 0.9649, "Alens": 1.0, "tau": 0.0544}
 fg_norm = {"nu_0": 150.0, "ell_0": 3000, "T_CMB": 2.725}
-fg_components = {"tt": ["tSZ_and_CIB", "cibp", "kSZ", "radio", "dust"], "te": ["radio", "dust"], "ee": ["radio", "dust"]}
-fg_params = {"a_tSZ": 3.30, "a_kSZ": 1.60, "a_p": 6.90, "beta_p": 2.08, "a_c": 4.90, "beta_c": 2.20, "a_s": 3.10, "a_gtt": 8.7, "a_gte": 0.355, "a_gee": 0.13, "a_psee": 0.05, "a_pste": 0, "xi": 0.1, "T_d": 9.60}
+fg_components = {'tt': ['tSZ_and_CIB', 'cibp', 'kSZ', 'radio', 'dust'], 'te': ['radio', 'dust'], 'ee': ['radio', 'dust'], 'bb': ['radio', 'dust'], 'tb': ['radio', 'dust'], 'eb': []} 
+fg_params = {"a_tSZ": 3.30, "a_kSZ": 1.60, "a_p": 6.90, "beta_p": 2.08, "a_c": 4.90, "beta_c": 2.20, "a_s": 3.10, "a_gtt": 8.7,  "xi": 0.1, "T_d": 9.60, "a_gte": 0, "a_gtb": 0, "a_gee": 0,  "a_gbb": 0, "a_pste": 0, "a_pstb": 0, "a_psee": 0, "a_psbb": 0}
 
 #sim
 iStart = 0
-iStop = 79
+iStop = 39
 sim_alm_dtype = "complex64"
+read_noise_alms_from_disk = False
+noise_sim_type = "fdw"
+noise_model_parameters = {"downgrade": 4, "mask_est_name": "dr6v3_20220316_baseline_union_mask", "mask_obs_name": "dr6v3_xlink_union_mask_0.001", "union_sources": "regular_20220316", "notes": "20220619"}
 
 
 #plot