Difference between revisions of "ForecastingSims"

From CMB-S4 wiki
Jump to navigationJump to search
(Created page with "===Simulations for forecasting r performance=== Here is a selection of codes and sim maps: ---- ====PySM==== Ben Thorne, David Alonso, Jo Dunkley https://github.com/bthorne93...")
 
Line 10: Line 10:
 
https://github.com/bthorne93/PySM
 
https://github.com/bthorne93/PySM
  
Spit out Nside=256 fits maps by running 'python main.py main_config.ini'. See the readme for any required dependencies.
+
Spits out Nside=256 fits maps by running 'python main.py main_config.ini'. See the readme for any required dependencies.
 +
 
 +
Code already includes a few options for each component, and can be easily extended to add new inputs.
 +
 
  
 
Suggested starting points and ini files for each:
 
Suggested starting points and ini files for each:
Line 23: Line 26:
  
 
----
 
----
 
  
 
====(Other available codes - TBD)====
 
====(Other available codes - TBD)====

Revision as of 04:59, 13 May 2016

Simulations for forecasting r performance

Here is a selection of codes and sim maps:


PySM

Ben Thorne, David Alonso, Jo Dunkley

https://github.com/bthorne93/PySM

Spits out Nside=256 fits maps by running 'python main.py main_config.ini'. See the readme for any required dependencies.

Code already includes a few options for each component, and can be easily extended to add new inputs.


Suggested starting points and ini files for each:

  • over-simple = single component dust, power-law synchrotron, minimal spatial variation of index, no polarised AME
  • towards more realism = 2 component dust, spatial steepening of the synchrotron index away from plane, 2% polarised AME.
  • towards more realism = 1 component dust with spatially varying dust index, spatial steepening of the synchrotron index away from plane, 2% polarised AME.

Current issues: The version of code on github (v0.3) doesn't yet have small-scale realizations added, so the polarization in low S/N regions is still noisy and all maps are smoothed to 1 degree resolution. This means that power spectra of e.g. the B2 region are noisy and so Victor's method using power spectra won't instantly work.

Solution: could try Victor's method on higher S/N regions as a test of how well things work. Or wait for the small-scales (~2 weeks). Pixel-based methods will work a little better on these sims.


(Other available codes - TBD)