Difference between revisions of "Data Challenges"
(→Galactic foregrounds: Added PIPXXX link for dust spectrum parameters)
|Line 104:||Line 104:|
=== Analysis ===
=== Analysis ===
=== Schedule ===
=== Schedule ===
Revision as of 13:13, 22 February 2017
At the Chicago meeting it was agreed that it would be useful to provide common inputs for all science forecasting, to include extra-galactic foreground, galactic foreground and noise maps at a range of frequencies. It was further agreed that, in order to support the widest range of uses, these would be provided as all-sky HEALPix maps at Nside-8192.
Possible sky-model inputs are:
- Extra-galactice: Alvarez/Battaglia/Bond, PSM
- Galactic: pySM, PSM
Please add to these lists as appropriate.
Inputs will be hosted at 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 will then be 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.
- 1 Data Challenge 1
- 1.1 Instrument specification
- 1.2 Masks: Clem
- 1.3 Scalar, tensor & non-Gaussian CMB: Borrill
- 1.4 Extragalactic foregrounds: Alvarez / Battaglia / Bond / Stein
- 1.5 Galactic foregrounds
- 1.6 Noise: homogeneous isotropic matched to N_ell: Victor / Clem / Colin / John K
- 1.7 Combined noise + CMB + dust + sync maps
- 1.8 Phase 1 Analysis
- 1.9 Schedule
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