GJ 1214b MIRI phase curve analysis

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