@@ -36,6 +36,20 @@ use_toeplitz_cov = True
3636# maps
3737src_free_maps_so = False
3838
39+ maps_dir = '/scratch/gpfs/SIMONSOBS/users/zatkins/projects/lat-iso/piso/tiger_dr6xdeep56_20251019/maps/'
40+ maps_so_i1_f090 = [maps_dir + 'so/sky_map_i1_f090_set%d_20251019.fits' % (i) for i in range(4)]
41+ maps_so_i1_f150 = [maps_dir + 'so/sky_map_i1_f150_set%d_20251019.fits' % (i) for i in range(4)]
42+ maps_so_i3_f090 = [maps_dir + 'so/sky_map_i3_f090_set%d_20251019.fits' % (i) for i in range(4)]
43+ maps_so_i3_f150 = [maps_dir + 'so/sky_map_i3_f150_set%d_20251019.fits' % (i) for i in range(4)]
44+ maps_so_i4_f090 = [maps_dir + 'so/sky_map_i4_f090_set%d_20251019.fits' % (i) for i in range(4)]
45+ maps_so_i4_f150 = [maps_dir + 'so/sky_map_i4_f150_set%d_20251019.fits' % (i) for i in range(4)]
46+ maps_so_i6_f090 = [maps_dir + 'so/sky_map_i6_f090_set%d_20251019.fits' % (i) for i in range(4)]
47+ maps_so_i6_f150 = [maps_dir + 'so/sky_map_i6_f150_set%d_20251019.fits' % (i) for i in range(4)]
48+ maps_so_c1_f220 = [maps_dir + 'so/sky_map_c1_f220_set%d_20251019.fits' % (i) for i in range(4)]
49+ maps_so_c1_f280 = [maps_dir + 'so/sky_map_c1_f280_set%d_20251019.fits' % (i) for i in range(4)]
50+ maps_so_i5_f220 = [maps_dir + 'so/sky_map_i5_f220_set%d_20251019.fits' % (i) for i in range(4)]
51+ maps_so_i5_f280 = [maps_dir + 'so/sky_map_i5_f280_set%d_20251019.fits' % (i) for i in range(4)]
52+
3953cal_so_i1_f090 = 1.02
4054cal_so_i1_f150 = 1.02
4155cal_so_i3_f090 = 1.02
@@ -132,6 +146,49 @@ leakage_beam_so_c1_f280_TB = leakage_dir + 'gamma_t2b_c1.txt'
132146leakage_beam_so_i5_f220_TB = leakage_dir + 'gamma_t2b_i5.txt'
133147leakage_beam_so_i5_f280_TB = leakage_dir + 'gamma_t2b_i5.txt'
134148
149+ mask_dir = data_dir + 'mask/'
150+
151+ xtra_mask_ivar_quantile = 0.1
152+
153+ ps_mask_so_i1_f090 = mask_dir + 'source_mask_15mJy_and_dust_rad12.fits'
154+ ps_mask_so_i1_f150 = mask_dir + 'source_mask_15mJy_and_dust_rad12.fits'
155+ ps_mask_so_i3_f090 = mask_dir + 'source_mask_15mJy_and_dust_rad12.fits'
156+ ps_mask_so_i3_f150 = mask_dir + 'source_mask_15mJy_and_dust_rad12.fits'
157+ ps_mask_so_i4_f090 = mask_dir + 'source_mask_15mJy_and_dust_rad12.fits'
158+ ps_mask_so_i4_f150 = mask_dir + 'source_mask_15mJy_and_dust_rad12.fits'
159+ ps_mask_so_i6_f090 = mask_dir + 'source_mask_15mJy_and_dust_rad12.fits'
160+ ps_mask_so_i6_f150 = mask_dir + 'source_mask_15mJy_and_dust_rad12.fits'
161+ ps_mask_so_c1_f220 = mask_dir + 'source_mask_15mJy_and_dust_rad12.fits'
162+ ps_mask_so_c1_f280 = mask_dir + 'source_mask_15mJy_and_dust_rad12.fits'
163+ ps_mask_so_i5_f220 = mask_dir + 'source_mask_15mJy_and_dust_rad12.fits'
164+ ps_mask_so_i5_f280 = mask_dir + 'source_mask_15mJy_and_dust_rad12.fits'
165+
166+ gal_mask_so_i1_f090 = mask_dir + 'HFI_Mask_GalPlane-apo0_2048_R2.00_GAL070_fejer1.fits'
167+ gal_mask_so_i1_f150 = mask_dir + 'HFI_Mask_GalPlane-apo0_2048_R2.00_GAL070_fejer1.fits'
168+ gal_mask_so_i3_f090 = mask_dir + 'HFI_Mask_GalPlane-apo0_2048_R2.00_GAL070_fejer1.fits'
169+ gal_mask_so_i3_f150 = mask_dir + 'HFI_Mask_GalPlane-apo0_2048_R2.00_GAL070_fejer1.fits'
170+ gal_mask_so_i4_f090 = mask_dir + 'HFI_Mask_GalPlane-apo0_2048_R2.00_GAL070_fejer1.fits'
171+ gal_mask_so_i4_f150 = mask_dir + 'HFI_Mask_GalPlane-apo0_2048_R2.00_GAL070_fejer1.fits'
172+ gal_mask_so_i6_f090 = mask_dir + 'HFI_Mask_GalPlane-apo0_2048_R2.00_GAL070_fejer1.fits'
173+ gal_mask_so_i6_f150 = mask_dir + 'HFI_Mask_GalPlane-apo0_2048_R2.00_GAL070_fejer1.fits'
174+ gal_mask_so_c1_f220 = mask_dir + 'HFI_Mask_GalPlane-apo0_2048_R2.00_GAL070_fejer1.fits'
175+ gal_mask_so_c1_f280 = mask_dir + 'HFI_Mask_GalPlane-apo0_2048_R2.00_GAL070_fejer1.fits'
176+ gal_mask_so_i5_f220 = mask_dir + 'HFI_Mask_GalPlane-apo0_2048_R2.00_GAL070_fejer1.fits'
177+ gal_mask_so_i5_f280 = mask_dir + 'HFI_Mask_GalPlane-apo0_2048_R2.00_GAL070_fejer1.fits'
178+
179+ extra_mask_so_i1_f090 = mask_dir + 'xtra_mask_so_i1_f090.fits'
180+ extra_mask_so_i1_f150 = mask_dir + 'xtra_mask_so_i1_f150.fits'
181+ extra_mask_so_i3_f090 = mask_dir + 'xtra_mask_so_i3_f090.fits'
182+ extra_mask_so_i3_f150 = mask_dir + 'xtra_mask_so_i3_f150.fits'
183+ extra_mask_so_i4_f090 = mask_dir + 'xtra_mask_so_i4_f090.fits'
184+ extra_mask_so_i4_f150 = mask_dir + 'xtra_mask_so_i4_f150.fits'
185+ extra_mask_so_i6_f090 = mask_dir + 'xtra_mask_so_i6_f090.fits'
186+ extra_mask_so_i6_f150 = mask_dir + 'xtra_mask_so_i6_f150.fits'
187+ extra_mask_so_c1_f220 = mask_dir + 'xtra_mask_so_c1_f220.fits'
188+ extra_mask_so_c1_f280 = mask_dir + 'xtra_mask_so_c1_f280.fits'
189+ extra_mask_so_i5_f220 = mask_dir + 'xtra_mask_so_i5_f220.fits'
190+ extra_mask_so_i5_f280 = mask_dir + 'xtra_mask_so_i5_f280.fits'
191+
135192apod_pts_source_degree = 0.3
136193apod_survey_degree = 2
137194edge_skip_rescale = 1
0 commit comments