ciao的简单使用，基于heasoft中xspec的能谱拟合

#ciao的简单使用，基于heasoft中xspec的能谱拟合| 来源: 网络整理| 查看: 265

前言

为了记录正在学习使用钱德拉望远镜的天文软件包ciao，并学习使用heasoft中xspec对数据进行能谱拟合，写下这篇博客. 如有错误之处及建议，欢迎到评论区讨论. PS：heasoft的安装请参看博客：https://blog.csdn.net/qq_19905867/article/details/109039709

一、下载样本数据

1.下载有专属的命令，本文以下载观测号19464为例.当然，你也可以官网先下载样本数据，之后解压就行.

下载样本数据的详细步骤请参看我的博客：https://blog.csdn.net/qq_46753404/article/details/109283789.这里不在赘述. 二、建立源区域和背景区域

可以先查看repro目录下都有哪些文件. 1.切换目录到repro：

cd 19464/repro

2.我们需要定义两个区域，一个用于源，另一个用于背景.要做到这一点，首先用ds9显示图像:

ds9 acisf19464_repro_evt2.fits -cmap b -scale log -zoom 4

3.选源区域和背景区域（PS：源区域要略小于背景区域，注意记得保存） 3.1源区域在这里插入图片描述 3.2背景区域

三、使用xspec进行能谱拟合

1.仍然在ciao初始化的环境下.依次输入以下命令.

punlearn specextract pset specextract infile="acisf19464_repro_evt2.fits[sky=region(src.reg)]" pset specextract bkgfile="acisf19464_repro_evt2.fits[sky=region(bkg.reg)]" pset specextract outroot=19464 pset specextract clobber=yes pset specextract grouptype=NUM_CTS pset specextract binspec=10 pset specextract weight=no specextract

注意第7行的代码数字默认是15，你可以根据需要改变，以下：

pset specextract binspec=10

2.进入xspec前，需要初始化，我的初始化快捷是在终端敲入h，打开xspec：

usrname$:xspec 效果如下： XSPEC version: 12.11.1 Build Date/Time: Wed Oct 14 18:55:19 2020 XSPEC12>

2.1 导入数据

data 19464_grp.pi

2.2 打开画图工具

cpd /xw plot data

效果如下：在这里插入图片描述

2.3将横坐标的能道（Channel）换为能量（Energy）

setplot energy plot data

效果如下：在这里插入图片描述 2.4 忽略不好的数据点（ps：小数点不能省略）

ignore **-1.0,8.0-** renorm

2.5 加入模型

model phabs*powerlaw /*

2.6 拟合

fit 1000

2.7 画出有残差的图

pl eu de

效果如下：在这里插入图片描述 2.8 保存

#先进入PLT iplot hardcopy 19464_fit.eps/cps quit save all 19464.xcm quit 总结

提示：这里对文章进行总结：

一、下载样本数据二、建立源区域和背景区域三、使用xspec进行能谱拟合

【本文地址】

公司简介

联系我们