UMICH-2015: Simulations & Data Analysis Break-Out Session 2
What S&DA do we need to do to support the science sections of the CMB-S4 Science Book?
What are the right inputs to the simulations?
- Sky model
- scalar, tensor, non-Gaussian
- Mission model
- beams, band-passes, noise PSDs
- detector electronics, data acquisition
- scanning strategy, flags
- atmosphere, ground pickup
What are the required output(s) of the simulations?
What are the required analysis outputs/metrics/figures of merit?
- Power spectra
What is the right balance between the veracity of the models used in the simulations and the tractability of their production and analysis?
- Spectral domain
- Map domain
- Time domain
What resources are available over the coming year?
- Cycles & storage
Forecasting using simulations with varying balance between complexity and feasibility?
Right balance issue: Large scale simulations with full complexity require long investment and significant resources. How can we go beyond very basic forecasting that neglects foregrounds and systematics and defines an instrument with only a beam size, a sky coverage, and a sensitivity, without implementing the full sims? How many levels of complexity should we consider, e.g.
- basic: sensitivity with Gaussian white stationary noise; gaussian beams; sky fraction
- moderate: inhomogeneous noise with non-stationarity and low-frequency excess
- advanced: model of atmospheric noise based on map differences from existing data on patches or on model; scanning and filtering + map-making for detector sets, scaling-up as a function of number of detectors; simple foreground subtraction with, e.g. decorrelation, NILC, or simple model fitting per pixel
- full: as realistic as possible
A few additional ideas:
- 1/ell noise: usually experiments do not perform as well on large scales as extrapolated from noise rms. Can we work-out an empirical law, possibly based on existing experiments?
- foregrounds: separate polarization issues (for primordial and lensing B-modes mostly) from temperature issues (on small scales, for SZ and extragalactic astrophysics). Should we worry about EE and TE foregrounds for parameter estimations?
- to what extend can we consider that we have a trustable model of the level of foreground residuals (to treat them as a noise in forecasting)?
- for cluster science we probably should consider the issues of
- relativistic corrections
- point source contamination in clusters