Difference between revisions of "ForecastingSims"

From CMB-S4 wiki
Jump to navigationJump to search
 
(11 intermediate revisions by the same user not shown)
Line 10: Line 10:
 
https://github.com/bthorne93/PySM
 
https://github.com/bthorne93/PySM
  
Spits out Nside=512 fits maps by running 'python main.py main_config.ini'. See the readme for any required dependencies.
+
Spits out Healpix fits maps up to Nside=512 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.
 
Code already includes a few options for each component, and can be easily extended to add new inputs.
  
 +
'''(Friday 13: Inputs not quite ready for use as need to check code before others use it - will send note out when ready)'''
  
 
Suggested starting points and ini files for each (noise level to be adjusted):
 
Suggested starting points and ini files for each (noise level to be adjusted):
*over-simple 's1d1a1' = power-law synchrotron, single component dust, minimal spatial variation of indices, no polarised AME  
+
*over-simple 's1d1a1' = power-law synchrotron, single component dust, minimal spatial variation of indices, no polarised AME. [[File:config_s1d1a1.txt]]
*towards more realism 's4d3a2'  = spatial steepening of the synchrotron index away from plane and curvature, 1 component dust with spatially varying dust index sigma=0.3, 2% polarised AME.   
+
*towards more realism 's4d3a2'  = spatial steepening of the synchrotron index away from plane and curvature, 1 component dust with spatially varying dust index sigma=0.3, 2% polarised AME.  [[File:config_s4d3a2.txt]]
*towards more realism  's4d4a2' = spatial steepening of the synchrotron index away from plane and curvature, 2 component dust,  2% polarised AME
+
*towards more realism  's4d4a2' = spatial steepening of the synchrotron index away from plane and curvature, 2 component dust,  2% polarised AME [[File:config_s4d4a2.txt]]
  
  
 
Current issues:
 
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.
+
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 in temp and 2 degrees in pol. This means that power spectra of e.g. the B2 region are noisy and so Victor's method using power spectra won't work right away.
  
 
Solution: could try Victor's method on higher S/N regions as a test of how well things work. Or wait for the small-scale Gaussian sims (~2 weeks). Pixel-based methods will work a little better on these current v0.3 sims.  
 
Solution: could try Victor's method on higher S/N regions as a test of how well things work. Or wait for the small-scale Gaussian sims (~2 weeks). Pixel-based methods will work a little better on these current v0.3 sims.  

Latest revision as of 07:17, 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 Healpix fits maps up to Nside=512 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.

(Friday 13: Inputs not quite ready for use as need to check code before others use it - will send note out when ready)

Suggested starting points and ini files for each (noise level to be adjusted):

  • over-simple 's1d1a1' = power-law synchrotron, single component dust, minimal spatial variation of indices, no polarised AME. File:Config s1d1a1.txt
  • towards more realism 's4d3a2' = spatial steepening of the synchrotron index away from plane and curvature, 1 component dust with spatially varying dust index sigma=0.3, 2% polarised AME. File:Config s4d3a2.txt
  • towards more realism 's4d4a2' = spatial steepening of the synchrotron index away from plane and curvature, 2 component dust, 2% polarised AME File:Config s4d4a2.txt


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 in temp and 2 degrees in pol. This means that power spectra of e.g. the B2 region are noisy and so Victor's method using power spectra won't work right away.

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


N.B. - this code is still under development. Please let us know about any issues/bugs you spot! And we welcome contributors!


(Other available codes - TBD)