GJ 1214b MIRI phase curve analysis

Outputs of the fiducial Gj 1214b analysis. Files: result.txt: eclipse depth, transit depth, and nightside Fp/Fs per wavelength bin. This file also reports amplitude and phase offset, but these are calculated ONLY from the first-order component of the phase curve, ignoring the second-order component. See amplitudes_and_phases.txt for the (relative) amplitude and phase offset while taking the second-order component into account. lightcurves.txt: best fit spectroscopic light curves Generate plots.ipynb: a Jupyter notebook that generates the *.png plots in the folder. If you're confused about the format of the data files, consult this notebook. *.png files: generated by the aforementioned Jupyter notebook. If you're confused about what a plot is showing, consult the notebook. white_lightcurve.txt: best fit white light curve white_envelope.txt, white_envelope_3sigma.txt: 1 sigma and 3 sigma envelopes of the white light curve white_light_result.txt: system parameters inferred from white light curve white_chain.npy: the emcee chain for the white light curve. Parameters: delta_t_transit, delta_t_eclipse, Fp, C1, D1, C2, D2, rp/rs, a/Rs, b, error_factor, Fstar, Aramp, tau, cy, slope. white_lnprobs.npy: ln(posterior) for every sample of the white chain chain_6000_6500.npy: the chain for the spectroscopic light curve spanning 6000--6500 Angstroms. Parameters: Fp, C1, D1, C2, D2, error_factor, rp/rs, Fstar, Aramp, tau, cy, slope lnprobs_6000_6500.npy: ln(posterior) for the spectroscopic chain mentioned above amplitudes_and_phases.txt: relative amplitude, (max - min) / max, of the phase curve, along with the phase offset, defined as the position of the phase curve maximum relative to mid-eclipse. Negative means east. calc_amplitude_offset.py: script that calculates the amplitudes and phase offsets in the above file, starting from the spectroscopic chain. Usage: python calc_amplitude_offset.py 6000, for the bin 6000-6500 Angstrom. fiducial_stellar_spectrum.fits: the stellar spectrum we adopted data.pkl: the reduced data that we fit for with emcee positions.txt: x and y positions inferred for the trace

本数据集包含GJ 1214b基准分析（fiducial analysis）的全部输出结果，各文件详情如下：
1. `result.txt`：存储各波长区间的食深（eclipse depth）、凌星深度（transit depth）以及夜面行星光度与恒星光度之比（Fp/Fs）。该文件同时报告了振幅与相位偏移，但二者仅通过相位曲线的一阶分量计算得到，未考虑二阶分量。如需获取考虑二阶分量时的（相对）振幅与相位偏移，请参考`amplitudes_and_phases.txt`文件。
2. `lightcurves.txt`：经最优拟合的光谱测光曲线（spectroscopic light curves）。
3. `Generate plots.ipynb`：用于生成本文件夹内所有*.png格式图表的Jupyter Notebook。若对数据文件的格式存在疑问，可查阅该Notebook。
4. *.png格式文件：由上述Jupyter Notebook生成。若对图表展示的内容存在疑问，可查阅该Notebook。
5. `white_lightcurve.txt`：经最优拟合的白光测光曲线（white light curve）。
6. `white_envelope.txt`与`white_envelope_3sigma.txt`：白光测光曲线的1σ与3σ置信包络线。
7. `white_light_result.txt`：通过白光测光曲线反演得到的系统参数。
8. `white_chain.npy`：用于白光测光曲线拟合的emcee采样链，包含参数包括：凌星时间偏移量（delta_t_transit）、食时间偏移量（delta_t_eclipse）、Fp、C1、D1、C2、D2、行星半径与恒星半径之比（rp/rs）、半长轴与恒星半径之比（a/Rs）、冲击参数（b）、误差因子（error_factor）、恒星光度（Fstar）、振幅斜坡（Aramp）、弛豫时间（tau）、基线偏移（cy）、斜率（slope）。
9. `white_lnprobs.npy`：白光采样链中每个样本的对数后验概率（ln(posterior)）。
10. `chain_6000_6500.npy`：对应6000~6500埃（Angstroms）波段的光谱测光曲线采样链，包含参数包括：Fp、C1、D1、C2、D2、误差因子（error_factor）、行星半径与恒星半径之比（rp/rs）、恒星光度（Fstar）、振幅斜坡（Aramp）、弛豫时间（tau）、基线偏移（cy）、斜率（slope）。
11. `lnprobs_6000_6500.npy`：上述光谱测光采样链对应每个样本的对数后验概率。
12. `amplitudes_and_phases.txt`：相位曲线的相对振幅（计算方式为(max - min)/max），以及相位偏移——即相位曲线峰值相对于食甚时刻的位置，负值代表东侧。
13. `calc_amplitude_offset.py`：用于从光谱测光采样链出发，计算上述文件中振幅与相位偏移的Python脚本。使用方法：以6000为例，执行`python calc_amplitude_offset.py 6000`即可对应6000~6500埃的波长区间。
14. `fiducial_stellar_spectrum.fits`：本次分析采用的恒星光谱（stellar spectrum）。
15. `data.pkl`：经过预处理的观测数据，用于通过emcee进行拟合。
16. `positions.txt`：通过迹线反演得到的x、y坐标位置。