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

前言

为了记录正在学习使用钱德拉望远镜的天文软件包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进行能谱拟合