2

我一直在尝试使用 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 中看到某些图案，因为它们是渐变和卷曲成分，这在我在网上找到的图像中很明显，但在我的产品中却没有。

所以我有三个问题：

  1. 我用来生成地图的逻辑是否有错误？

  2. 是否需要额外的过滤才能看到预期的 E 和 B 模式？

  3. 在 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 参数似乎不包含在此版本中。

4

1 回答 1

1

  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 设为零，只有一个非零。

  2. 您可以使用假设为 1 度的高斯光束来平滑输出图，以模拟真实仪器如何看待该信号。

  3. 不支持在healpy

于 2015-12-08T07:06:29.797 回答