我一直在尝试使用 CAMB 和 healpy 制作 E 和 B 模式地图,但是到目前为止我制作的东西看起来并不多。
我正在使用 CAMB 在以下链接生成模拟地图:http: //lambda.gsfc.nasa.gov/toolbox/tb_camb_form.cfm,除了选中张量 cl 的框和不同的 r 之外,我保留了所有默认参数.
有了这个,我一直在使用以下代码来生成地图:
# read in the power spectra produced by camb, ordering is: [TT, EE, BB, TE]
cl = hp.read_cl('scalar_and_tensor_power_spectra_produced_by_CAMB.fits')
# set the random seed for phase information
numpy.random.seed(some number) #Im using 1 for now
#Produce E mode and B mode maps
EEmap = hp.synfast(cl[1], NSIDE)
BBmap = hp.synfast(cl[2], NSIDE)
快速的谷歌搜索显示与我得到的结果截然不同的结果,您可以在以下链接中看到: http://tinypic.com/view.php ?pic=ifvrpt &s=9# .VmY7Jx8So8o http://tinypic.com /view.php?pic=a5ngx&s=9#.VmY7KR8So8o
我期待在 E 和 B 中看到某些图案,因为它们是渐变和卷曲成分,这在我在网上找到的图像中很明显,但在我的产品中却没有。
所以我有三个问题:
我用来生成地图的逻辑是否有错误?
是否需要额外的过滤才能看到预期的 E 和 B 模式?
在 healpy 以外的 healpix 版本中,mollview 中有一个名为 POLARIZATION 的参数(http://healpix.jpl.nasa.gov/html/idlnode33.htm#idl:mollview:polarization),其中设置 POLARIZATION = 3 允许您过度绘制极化信息作为无头矢量。这就是我认为他们在这里所做的事情:http: //bicepkeck.org/media/b_over_b_rect_BICEP2.png 和这里的第一张图片: https ://inspirehep.net/record/787860/plots ,有没有办法做到这一点在愈合?由于 POLARIZATION 参数似乎不包含在此版本中。