夏尔巴射出包(星系团、星系团和星系的x射线分析)
deproject的Python项目详细描述
Deproject是一个Sherpa扩展包 便于二维环形x射线能谱的回注 三维源特性。对于典型的热模型,这将 包括径向温度和密度分布。这个基本方法 广泛用于x射线聚类分析,是xspec的基础。 型号projct。deproject模块将此功能 sherpa是一个易于使用和理解的python模块。
deproject的基本物理假设是 发射率在半径为 对应于用于提取specta的环隙。假设一 为每个环形光谱构造一个线性体积加权模型 外壳模型的组合。
deproject的0.2版仅限于环形。
更多文档可在https://deproject.readthedocs.io/
许可证
deproject模块在 BSD 2-Clause license, 作为源发行版中的文件LICENSE提供。
要求
安装假定您正在将deproject安装到 CIAO environment(CIAO 4.11或 因为这是获得xspec模型的最简单方法 和夏尔巴一起。第standalone Sherpa 可以使用版本,但在这种情况下,您需要build Sherpa with XSPEC support。
需要以下python包:
安装
deproject模块应使用以下命令安装 (假设已安装CIAO 4.11):
echo “numpy==1.12.1” > constraints.txt pip install -c constraints.txt ‘astropy<3.1’ deproject
示例
如果您有一组称为src<;n>;.pi的x射线pha光谱,其中<;n>;是 一个整数,表示环空数,文件包含 XFLT0001到XFLT0005头关键字由 XSPEC projct model, 然后a Deproject object 可以使用 deproject_from_xflt 带有命令的助手例程:
>>> from deproject import deproject_from_xflt >>> from astropy import units as u >>> dep = deproject_from_xflt('src*.pi', 0.492 * u.arcsec)
在本例中,XFLT0001和XFLT0002关键字, 规定环空内外半径的 ACIS像素,因此需要乘以0.492 ArcSeconds 转换为角度(第二个参数)。
这将自动将光谱加载到单独的夏尔巴数据集中, 哪个可以单独安装,但通常更易于使用 由deproject_from_xflt返回的对象。例如, 下面将为每个光谱设置适合的数据范围,并确保 在拟合之前减去背景:
>>> dep.ignore(None, 0.5) >>> dep.ignore(7.0, None) >>> dep.subtract()
夏尔巴函数用于更改统计和优化程序:
>>> from sherpa.astro import ui >>> ui.set_stat('chi2xspecvar') >>> ui.set_method('levmar')
数据可以拟合,所有参数的误差估计,使用 洋葱皮脱毛方法,使用以下命令:
>>> onion = dep.fit() >>> errs = dep.conf()
返回值包括密度(和误差,如果合适的话),如 一个Astropy Quantity。
>>> print(onion['density']) print(onion['density']) density 1 / cm3 -------------------- 0.1100953546292787 0.07736622021374819 0.04164827967805805 0.03630168106524076 0.025221797991301052 0.021845331641349316 ... 0.012396857131392835 0.01336640115325031 0.012303975980575187 0.013631563529090736 0.013996131292837352 0.010843683594144967 0.023067220584935984 Length = 20 rows
那on-line documentation 包含更多信息,包括创建Deproject对象 直接(不需要XFLTxxxx关键字)。