HXI能谱分析教程 — Zhentong Li

本文面向已有 FITS、IDL/SolarSoft 或太阳高能数据分析基础的学生和研究者。目标是完成一次可检查、可复现的 HXI 能谱分析：从选择观测区间开始，得到背景方案明确、响应文件匹配、参数含义清楚的拟合结果。

操作流程依据 2026 年 5 月 12 日培训报告；2023 年用户手册用于补充数据产品与界面说明。软件界面可能随版本变化，但数据筛选、背景检验与前向拟合的逻辑相同。

开始之前：理解拟合的对象

HXI 记录的是探测器的计数谱，而我们想讨论的热等离子体与非热电子通常要通过光子谱模型来描述。探测器响应矩阵（DRM）把入射光子谱映射为预期计数谱：

C_i^{\rm model}=\sum_j R_{ij}F_j + B_i.

这里 $F_j$ 是模型光子谱， $R_{ij}$ 描述能量响应， $B_i$ 是背景。由模型算计数相对直接；反过来从计数唯一地“还原”光子谱并不成立。因此实际分析是选定物理模型、经 DRM 前向折叠，再与观测计数比较。下面按 2026 年 5 月 12 日 HXI 培训报告的流程整理；早期软件手册用来补充界面与数据产品的背景信息。

实际文件可能使用计数、计数率或单位能量计数率，响应也可能包含相应的归一化。上式表达的是前向关系，不是可以忽略文件单位直接相乘的通用计算式。

准备一个独立的事件目录，将原始数据、生成的能谱、响应矩阵、背景文件、拟合配置及输出图表分开保存。至少记录事件日期、目标时间、探测器、数据级别和软件版本；之后改变时间 bin 或背景时，为结果使用不同名称。

1. 先判断这段数据能不能拟合

先看 HXI quicklook 的全开流量探测器光变和动态能谱，再对照背景探测器、仪器状态标记及 GOES 等独立信息。需要回答四个问题：目标时段有没有有效数据？增强是否与耀斑一致？目标能区的计数是否足以做能谱？背景能否可靠估计？

背景探测器并非全开探测器的完美替身；简单相减后仍可能留有残差。尤其要检查粒子环境、轨道变化和高能端是否存在真实耀斑辐射。高能端不能未经核对就当成纯背景。

首次处理一个事件，建议先画出覆盖耀斑前、上升、峰值与衰减阶段的多能段光变，并标记准备拟合的时段。不要只截取峰值附近几秒，否则很难判断背景变化和数据完整性。

检查	具体操作	决策
数据覆盖	查看目标区间及前后有效数据、状态标记	避开缺测和异常区间
事件归属	对照全开探测器、背景探测器、动态能谱及 GOES	排查粒子增强或轨道背景变化
信噪比	比较目标能段的总谱、背景谱与误差	确定可以拟合的能量上限和时间长度
快速演化	比较不同能段光变及硬度变化	过快变化时缩短区间，计数不足时谨慎合并

2. 生成能谱文件和 DRM

下载目标时段的 HXI 数据，确认 HXI/SSW 环境与数据目录。培训示例中的用户数据路径设置在 hxi/setup/setup.hxi_env：

setenv HXI_DATA_USER "your path"

在 HXI GUI 中设置时间范围、探测器、时间 bin 与能量 bin，查看光变和能谱后保存能谱 FITS及与之匹配的响应矩阵 FITS。2023 年用户手册以 D94 全开探测器和邻近的 D99 背景探测器为示例；实际应按研究目标、可用计数和背景行为选择探测器，不应把这个组合当成所有事件的固定答案。

保存前记录数据级别、所选探测器、时间/能量 bin 和背景方案。过细的时间或能量 bin 会让统计不确定度增大；过粗又可能抹掉演化。图像分析的 bin 建议不能直接套用于能谱拟合。

一次生成后，逐项检查：

能谱文件覆盖计划拟合的时间段，背景所需的时段或独立背景文件也已准备好。
能量通道边界与预期一致；不要只根据文件名猜测 bin 设置。
DRM 对应同一探测器和所生成能谱的能量设置。重新生成能谱后，一并检查响应，不沿用未经确认的旧文件。
如选择了 Use AI BKG，在文件记录中明确“已扣背景”，后续不得再次扣除。

时间长度没有适合所有耀斑的固定值。可以先选取形态相对稳定、计数较充分的一段完成试拟合，再逐步缩短到所需的时间分辨率。高能端以可用统计量和背景不确定度为依据，不应为了画出更长的谱线而延伸拟合区间。

3. 背景方案要和数据一起审查

2026 年培训报告给出三类思路：

用耀斑前后相对平静的区间，在 OSPEX 内建模背景；前提是背景随轨道变化的趋势可被合理描述。
用约 48 小时前或后的相近轨道数据估计背景。需要保留相应的 Level 1.0 文件，将其时间轴对齐后，通过脚本替换 OSPEX 背景；这不是只在 GUI 中选一个背景区间。
在生成能谱时选择 Use AI BKG。这一方案生成的能谱已扣除预测背景，进入 OSPEX 后不要再重复扣除。2026 年培训材料说明 AI 背景数据当时尚未常规公开发布，测试数据通过培训资料提供；实际使用前应向数据团队确认对应时段的数据和软件版本是否可用。

比较候选背景时，应同时查看低能和高能端的计数、拟合残差，以及不同背景选择对关键参数的影响。背景误差可能比形式上的拟合误差更重要。

使用约 48 小时前后轨道背景时，应先保存当前事件的能谱，再读取参照轨道数据。核对探测器、能道、曝光和有效区间，将参照数据按相应轨道阶段映射并插值到目标时间轴，之后才在 OSPEX 中替换背景。培训报告给出的接口是 o->replacedata；具体调用参数应以所用 HXI 示例脚本和当前 OSPEX 方法定义为准。不要仅把日期平移 48 小时就认定背景一定相同。

若不同背景方案使关键参数明显变化，可以把这种变化作为背景敏感性单独报告。不能只保留统计误差最小的一次拟合，也不应把扣背景后的负值全部截成零来获得更“好看”的谱线。

4. 在 OSPEX 中前向拟合

更新 SSW 的 spex，启动 OSPEX 并选用 HXI 文件读取器：

o = ospex()
o->set, spex_file_reader='hxi'

导入匹配的能谱和 DRM，设置拟合时间段、能量范围与背景，再选择模型函数和参数范围。常见的起点是热成分 vth 与非热成分；bpow 描述光子幂律，thick2 则是电子厚靶模型，两者参数不能混读。培训报告建议用 thick2 时先测试较简单的单幂律电子分布，确有需要再增加自由度。

GUI 适合先检查数据和单个时间段；确认流程后，再用脚本保存和复用配置。下面是培训报告中涉及的主要设置入口，右侧变量需先由实际文件和研究目标定义，不能原样直接运行：

; IDL / OSPEX：配置骨架，先为右侧变量赋实际值。
o = ospex()
o->set, spex_file_reader='hxi'
o->set, spex_specfile=specfile
o->set, spex_drmfile=drmfile
o->set, spex_erange=energy_range
o->set, spex_fit_time_interval=fit_interval
o->set, fit_function=model_name
o->set, fit_comp_params=initial_parameters
o->set, fit_comp_minima=lower_bounds
o->set, fit_comp_maxima=upper_bounds
o->set, fit_comp_free_mask=free_parameters

先在界面中确认能谱读取正确，再设置背景和执行拟合。spex_uncert、迭代次数以及 albedo 等选项，应依据当前软件说明和物理目标设置并记录；不要复制一套未说明含义的默认数值。

5. 选择模型与自由参数

模型或选项	主要用途	使用时注意
`vth`	光学薄热辐射成分	温度与发射度相关；参数单位以当前 OSPEX 定义为准
`bpow` / `bpow_ep`	描述光子谱的折断幂律形状	光子谱指数不等于电子谱指数
`thick2`	由注入电子分布计算厚靶轫致辐射	低能截止与热成分、电子通量常有强相关
`thick2_vnorm`	用参考电子能量处的分布值归一化	参考能量应处于截止能量之间，拟合时保持固定
albedo 修正	考虑光球反照对模型的影响	核对源位置和各向异性假设，按科学目的使用

建议先用尽可能少的成分解释谱形。如果只在高能端出现显著且稳定的剩余计数，再评估非热成分；若残差在连续能段呈系统性结构，先检查背景、响应和能区选择，再考虑增加模型自由度。

thick2 的参数描述的是电子通量分布，组件图展示的是该电子分布产生的 X 射线光子谱。其积分归一化与截止能量相互关联，不能把归一化参数当作与其他参数完全独立的常数。设置成单幂律形式的方法应核对当前模型的折断能量和两段指数定义，不凭参数序号猜测。

6. 判断结果是否可信

拟合完成后不要只看一个 $\chi^2$ 数字。应核对计数谱与前向模型、各能道残差、能量范围和时间范围是否合理；再改变背景方案与初值，检查温度、谱指数和低能截止等结论是否稳定。仅凭能谱拟合往往不足以排除不同物理模型。

逐项完成下面的检查，再开始讨论物理结论：

收敛与边界： 是否达到迭代上限？参数是否顶住上下界？换一个合理初值是否落在相近解上？
残差形状： 残差是否在某段连续偏正或偏负？高能端的模型成分是否依赖少数低计数通道？
区间敏感性： 略微调整时间或能量范围后，关键参数是否发生超出误差的大幅变化？
背景敏感性： 采用另一合理背景后，非热成分与低能截止是否仍受约束？
跨仪器约束： 热成分是否与 GOES、AIA DEM 或其他 X 射线观测相容？比较时注意时间、视场及柱/体积发射度定义。

若热成分覆盖了低能截止附近的辐射，能谱可能只给出截止能量的范围或上限，而不是精确测量。此时应报告约束的局限，不把一次拟合输出的小误差当作充分证据。

7. 常见问题与处理顺序

现象	优先检查
无法读取能谱或响应	HXI 读取器、文件路径、SSW/spex 版本、文件是否成对生成
扣背景后大量负值	背景是否过高、重复扣除、轨道或时间轴匹配是否错误
高能端异常平坦或隆起	背景统计、粒子环境、是否存在真实高能成分
拟合不收敛或参数撞边界	数据单位、初值、拟合区间、模型复杂度与计数质量
减小误差后拟合更差	误差模型是否合理；不能通过随意放大不确定度掩盖系统偏差
相邻时段参数剧烈跳动	真实演化、计数不足、背景变化以及参数退化分别排查

8. 保存一份可复现的分析结果

至少保留原始/背景数据来源、所用探测器、能谱和 DRM 文件名、时间与能量 bin、拟合区间、模型组件、初值与参数边界、自由参数掩膜、统计量和不确定度设定。输出总计数谱、背景、前向模型、残差，以及主要参数随时间变化的图表。

对于批量时间段分析，先验证一个代表性时段，再把它的设置推广到其他时段。可以用前一时段结果作为初值，但仍要检查每一段的拟合质量，避免错误解连续传递。

最终报告应能让另一个人仅凭你的配置与文件清单复现结果，并理解哪些结论来自数据，哪些依赖背景、模型或几何假设。

资料与版本

李振同，《HXI 能谱数据处理实例》，2026 年 5 月 12 日培训报告（本文操作流程与背景处理的主要依据）。
ASO-S/HXI software User Guide（2023-06-01）：数据产品、GUI 与 OSPEX 的早期说明；如与 2026 年培训材料不一致，以后者为准。
OSPEX 拟合模型组件说明：拟合函数及参数定义。

本文不分发事件数据和培训测试版软件。需要对应示例数据、AI 背景或背景替换脚本时，可通过联系方式联系作者，并说明事件时间、软件版本和拟合设置。