夏尔巴射出包(星系团、星系团和星系的x射线分析)

deproject的Python项目详细描述


Deproject是一个Sherpa扩展包 便于二维环形x射线能谱的回注 三维源特性。对于典型的热模型,这将 包括径向温度和密度分布。这个基本方法 广泛用于x射线聚类分析,是xspec的基础。 型号projctdeproject模块将此功能 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包:

  • sherpa
  • Astropy (restricted to version 3.0 when using CIAO 4.11)
  • SciPy.

安装

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>;是 一个整数,表示环空数,文件包含 XFLT0001XFLT0005头关键字由 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)

在本例中,XFLT0001XFLT0002关键字, 规定环空内外半径的 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关键字)。

欢迎加入QQ群-->: 979659372 Python中文网_新手群

推荐PyPI第三方库


热门话题
java计时器结束得太早?   java通过行选择编辑jtable内容   java如何将数组的值传递给另一个类?   java对象。hashCode()算法   代号为1且lambda不工作的java Android本机代码   泛型Java方法签名中是否真的需要“?扩展”?   异常处理声纳问题“main”不应“抛出”任何JAVA 7   java无法使用MVEL导入静态方法   java什么是类路径容器?   java从InputStream解码字节   Eclipse重构仅限Java   java使用Qpid通过SSL与AMQP1.0代理通信   java有没有一种方法可以在Oracle中使用带有时区的时间戳作为函数输入   保存文件时运行安卓 lint的java出错   尽管添加了JScrollPane,java JList视图区域仍在增长   java计算第n个浮点   java如何在mysql中插入日期?   GWT中异步回调的java问题