Question

我一直在尝试使用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模式？
- 类似于以下链接： http://cosmologist.info/polar/EBsupport.html
在除了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，有没有办法在healpy中这样做？由于POLARIZATION参数似乎未包含在此版本中。

Answer 1

请参阅synalm http://healpy.readthedocs.org/en/1.x/generated/healpy.sphtfunc.synalm.html#healpy.sphtfunc.synalm的帮助，如果您hp.synfast(cl)正在创建温度为cl的温度图。您应该为每个组件创建一个Cl列表，即TT，EE，BB，TE，然后将所有Cl归零并且只有一个非零。
您可以使用高斯光束来平滑输出地图，例如1度，以模拟真实乐器如何看待该信号。
不支持在healpy

使用CAMB和Healpy生成E和B模式映射

1 个答案: