- 1 Experiment Definitions
- 2 Input Sky Components
- 3 Sky Models
- 4 Sim Data Sets
- 5 NERSC Details
- 6 Data Challenge 1
- 6.1 Instrument specification
- 6.2 Masks: Clem
- 6.3 Scalar, tensor & non-Gaussian CMB: Borrill
- 6.4 Extragalactic foregrounds: Alvarez / Battaglia / Bond / Stein
- 6.5 Galactic foregrounds
- 6.6 Noise: homogeneous isotropic matched to N_ell: Victor / Clem / Colin / John K
- 6.7 Combined noise + CMB + dust + sync maps
- 6.8 Phase 1 Analysis
- 6.9 Schedule
Input Sky Components
Sky model 00 is simply LCDM + Gaussian dust + Gaussian Sync. Since the dust and sync levels are set to the levels found in the BICEP/Keck field it wouldn't make any sense to use this model for a large sky fraction.
Sky Model 01 is PySM run in a1d1f1s1 mode - i.e. with the default settings for AME, dust, free-free and synchrotron.
Sim Data Sets
01.00 is the 01 experiment definition plus the 00 sky model.
Shared space is available on NERSC, where everyone is welcome to sign-up for an account under the "Data Analysis for Post-Planck CMB Experiments" allocation (PI Borrill). The files are located in the shared CMB-S4 file space (/project/projectdirs/cmbs4). In order to manage this space, all files stored there should be owned by the cmbs4 project account but accessible to the cmbs4 group. To do this:
- log in to the cmbs4 project account following the instructions [here]
- sync the files to an appropriate subdirectory in the project space
- ensure that the permissions are set appropriately (g+rX,o-rwx)
Remember to include a README and to post the details on this wiki page.
For any NERSC issues, including access to the filegroup and/or project account, please contact Julian Borrill.
Data Challenge 1
Intended as a simple test to develop our practice, compare spectral- and map-based forecasting, contrast different map-based codes and methods, and coordinate r and non-r science forecasting.
Maps pager for data challenge 1 (version 3, posted 2017-02-10)
- All maps are in Healpix format with nside = 512 and ℓmax = 1024.
- Maps have been filtered to remove all signal and noise below ℓ = 30.
|Beam FWHM (arcmin)||76.6||57.5||27.0||24.2||15.9||14.8||10.7||8.5|
- Tophat bandpasses described in this posting, added to above table --Cbischoff (talk) 15:51, 8 November 2016 (UTC)
- Beam widths copied from Victor's posting here --Cbischoff (talk) 19:20, 14 December 2016 (UTC)
Beam window functions for each frequency can be found on NERSC in:
In addition to beam smoothing at high ell, these window functions drop sharply to zero below ell=30 (maps for the data challenge contain no information below ell=30).
Noise inverse variance map is defined as a flat circular region with radius = 12°, surrounded by a tapered region that falls to zero with cos profile over an additional 15 degrees. This map is defined to range from 0 (in unobserved pixels) to 1 (in the flat 12° circular region); noise simulations are generated over the full sky from a C_l spectrum, then multiplied by 1 / sqrt(inverse variance) so that the noise amplitude blows up near the edge of the field.
- (Sum of inverse variance map) / (# of pixels in full sky) = 0.0293, so total noise power should be equivalent to fsky = 3%.
- Fraction of pixels that are observed is 5.4%.
File located on NERSC at:
Scalar, tensor & non-Gaussian CMB: Borrill
Lensed scalar CMB, consistent with Planck 2015 and low τ parameters
- Lensed CMB maps are posted on NERSC in directory:
- Unlensed CMB maps are posted on NERSC in directory:
Corresponding Kappa Maps?
Extragalactic foregrounds: Alvarez / Battaglia / Bond / Stein
Not included in Data Challenge 1
- Amplitude (in BB spectrum) is 4.25 μK2 at ν = 353 GHz, ℓ = 80. This is the value that was used for Science Book forecasts and is listed in Victor's 2016-05-31 posting (Section 2). It also corresponds to the best fit dust amplitude from the BK14 result.
- Dust amplitude in EE is 2x larger than BB. Dust amplitude in TT is 10x larger than EE (20x larger than BB). There is no TE, TB, or EB correlation.
- Dust scaling in frequency follows a greybody spectrum with βdust = 1.6 and Tdust = 19.6 K. This choice of parameters comes from PIPXXX.
- Dust Dℓ scaling follows a power law in ell with exponent αdust = -0.4. This parameter also comes from PIPXXX.
- Dust maps are posted on NERSC in directory:
- Amplitude (in BB spectrum) is 3.8 μK2 at ν = 23 GHz, ℓ = 80. This is the value that was used for Science Book forecasts and is listed in Victor's 2016-05-31 posting (Section 2). It also corresponds to the 95% upper limit from the BK14 result.
- Sync amplitude in EE is 2x larger than BB. Sync amplitude in TT is 10x larger than EE (20x larger than BB). There is no TE, TB, or EB correlation.
- Sync scaling in frequency follows a power-law spectrum with βsync = -3.1.
- Sync Dℓ scaling follows a power law in ell with exponent αsync = -0.6.
- Sync maps are posted on NERSC in directory:
Noise: homogeneous isotropic matched to N_ell: Victor / Clem / Colin / John K
White + 1/ℓ noise
- Nℓ functional form and parameters listed in Victor's 2016-12-20 posting (Table 2).
- Noise aℓm are calculated out to ℓ = 4096, but maps are rendered at nside = 512 with ℓmax = 1024. No beam smoothing or pixel window function is applied.
- Full sky noise maps are generated from Nℓ, then scaled by the square root of variance map to produce noise that is lowest in the center of the observed field and blows up towards the edge.
- Noise maps are posted on NERSC in director:
- Maps with filename like r0000_f030_n512.fits are full sky noise maps. Filenames like r0000_f030_n512_b00p0_mfsky03.fits have been multiplied by the square root of variance map.
Combined noise + CMB + dust + sync maps
- Four versions of these maps, with different levels of residual lensing signal. Example file names are:
- AL = 1 : r0000_f030_AL1p0_n512_mfsky03.fits
- AL = 0.3 : r0000_f030_AL0p3_n512_mfsky03.fits
- AL = 0.1 : r0000_f030_AL0p1_n512_mfsky03.fits
- AL = 0.03 : r0000_f030_AL0p03_n512_mfsky03.fits
- Maps with partial lensing are constructed by the following combination:
partially_lensed_map = sqrt(A_L) * lensed_map + (1 - sqrt(A_L)) * unlensed_map
- Maps are posted on NERSC in directory:
Phase 1 Analysis
Please add your main results here and link to your more detailed posting in the usual logbook. In the logbook, please specify which maps you made use of in your analysis, i.e. whether you used the combined maps, CMB simulations, noise simulations, etc.
- Dec 7: Beta version of Phase 1 maps
- Dec 14: Telecon
- Dec 17: Noise power (N_ell) prescriptions to Clem from Victor and Colin
- Dec 21: A few realizations of Phase 1 maps delivered by Clem
- Dec 21: Telecon (note phase shift)
- Jan 4: Telecon
- Jan 12: 100 realizations of Phase 1 maps posted to NERSC
- Jan 18: Telecon
- Feb 1: Telecon
- Feb 15: Telecon
- Feb 27, 28, Mar 1: S4 general meeting at SLAC
- Mar 2, 3: S4 CDT meeting at SLAC