4. Transmission spectra: full¶
4.1. Setup¶
%load_ext autoreload
%autoreload 2
import glob as glob
import matplotlib as mpl
import matplotlib.patheffects as PathEffects
import matplotlib.pyplot as plt
import matplotlib.transforms as transforms
import numpy as np
import pandas as pd
import seaborn as sns
import corner
import json
import pathlib
import pickle
import utils
import warnings
from astropy import constants as const
from astropy import units as uni
from astropy.io import ascii, fits
from astropy.time import Time
from mpl_toolkits.axes_grid1 import ImageGrid
# Default figure dimensions
FIG_WIDE = (11, 5)
FIG_LARGE = (8, 11)
# Figure style
sns.set(style="ticks", palette="colorblind", color_codes=True, context="talk")
params = utils.plot_params()
plt.rcParams.update(params)
4.2. Dowload data¶
Unzip this into a folder named data in the same level as this notebook
4.3. Plot¶
base_dir = "data/tspec_full"
fig, ax = plt.subplots(figsize=FIG_WIDE)
utils.plot_tspec_IMACS(ax, base_dir)
# Save
fig.tight_layout()
fig.set_size_inches(FIG_WIDE)
#utils.savefig(f"../paper/figures/tspec_full/tspec_full.pdf")
offsets: [-381.66382639 -250.19806149 604.72304087 238.023367 -134.88932493]
offsets (% mean wlc depth): [-2.94346912 -1.92957838 4.66374719 1.83568466 -1.04029393]
Saving tspec to: data/tspec_full/tspec_c.csv
mean WLC depth: 12966.462729191497 228.463440286341
Rp (Rj): 1.2765205003564386 jupiterRad
Rs (Rsun): 1.152 solRad
gp (m/s^2): 2920.5267508681527 cm / s2
base_dir = "data/tspec_full"
fig, ax = plt.subplots(figsize=FIG_WIDE)
utils.plot_tspec_IMACS(ax, base_dir)
# Save
fig.tight_layout()
fig.set_size_inches(FIG_WIDE)
utils.savefig(f"../paper/figures/tspec_full/tspec_full.pdf")
offsets: [-381.66382639 -250.19806149 604.72304087 238.023367 -134.88932493]
offsets (% mean wlc depth): [-2.94346912 -1.92957838 4.66374719 1.83568466 -1.04029393]
Saving tspec to: data/tspec_full/tspec_c.csv
mean WLC depth: 12966.462729191497 228.463440286341
Rp (Rj): 1.2765205003564386 jupiterRad
Rs (Rsun): 1.152 solRad
gp (m/s^2): 2920.5267508681527 cm / s2
4.4. Table¶
pd.read_csv("data/tspec_full/tspec_c.csv")#.to_clipboard(index=False)
| Wavelength (Å) | Transit 1 | Transit 2 | Transit 3 | Transit 4 | Transit 5 | Combined | |
|---|---|---|---|---|---|---|---|
| 0 | 5400.0 - 5600.0 | 13100^{+1300}_{-1400} | 13900^{+300}_{-300} | 13600^{+600}_{-600} | 13300^{+600}_{-600} | 14000^{+600}_{-600} | 13800 \pm 200 |
| 1 | 5600.0 - 5800.0 | 14200^{+1400}_{-1400} | 12000^{+1200}_{-1100} | 13300^{+700}_{-700} | 13300^{+600}_{-600} | 13400^{+500}_{-500} | 13300 \pm 300 |
| 2 | 5800.0 - 5985.8 | 13600^{+1100}_{-1200} | 11900^{+1100}_{-1000} | 14000^{+700}_{-800} | 13700^{+400}_{-400} | 13000^{+800}_{-700} | 13500 \pm 300 |
| 3 | 5985.8 - 6185.8 | 12800^{+1000}_{-1000} | 12300^{+1000}_{-500} | 12800^{+1200}_{-1300} | 13200^{+500}_{-500} | 14000^{+700}_{-700} | 13200 \pm 300 |
| 4 | 6655.0 - 6855.0 | 14000^{+900}_{-800} | 13000^{+600}_{-700} | 12600^{+500}_{-600} | 13200^{+600}_{-600} | 13800^{+400}_{-400} | 13300 \pm 200 |
| 5 | 6855.0 - 7055.0 | 13100^{+900}_{-900} | 13400^{+400}_{-500} | 13000^{+600}_{-600} | 14000^{+600}_{-600} | 14100^{+400}_{-400} | 13700 \pm 200 |
| 6 | 7055.0 - 7255.0 | 13000^{+700}_{-700} | 12600^{+900}_{-900} | 12900^{+700}_{-700} | 13800^{+600}_{-600} | 13800^{+300}_{-300} | 13600 \pm 200 |
| 7 | 7255.0 - 7455.0 | 13400^{+800}_{-700} | 12500^{+600}_{-600} | 13500^{+600}_{-600} | 13300^{+500}_{-400} | 14300^{+200}_{-300} | 13800 \pm 200 |
| 8 | 7455.0 - 7655.0 | 12900^{+1000}_{-1000} | 12200^{+600}_{-600} | 12800^{+600}_{-600} | 13400^{+400}_{-400} | 14000^{+500}_{-500} | 13200 \pm 200 |
| 9 | 7655.0 - 7709.0 | 12200^{+900}_{-800} | 12300^{+700}_{-700} | 13800^{+1000}_{-1000} | 13500^{+700}_{-700} | 13200^{+500}_{-500} | 13000 \pm 300 |
| 10 | 7709.0 - 7909.0 | 13100^{+800}_{-800} | 12200^{+600}_{-600} | 13000^{+500}_{-400} | 14000^{+700}_{-700} | 14100^{+400}_{-400} | 13400 \pm 200 |
| 11 | 7909.0 - 8109.0 | 12400^{+600}_{-600} | 12200^{+700}_{-700} | 12600^{+500}_{-500} | 13600^{+500}_{-600} | 13600^{+400}_{-400} | 13000 \pm 200 |
| 12 | 8109.0 - 8269.0 | 12500^{+400}_{-300} | 12600^{+700}_{-800} | 12300^{+500}_{-500} | 13200^{+500}_{-500} | 13700^{+500}_{-500} | 12900 \pm 200 |
| 13 | 8269.0 - 8469.0 | 12600^{+400}_{-400} | 12600^{+800}_{-900} | 12400^{+400}_{-400} | 13400^{+500}_{-500} | 13800^{+400}_{-400} | 13100 \pm 200 |
| 14 | 8469.0 - 8669.0 | 12300^{+600}_{-600} | 12400^{+600}_{-500} | 11800^{+500}_{-500} | 13600^{+500}_{-500} | 13400^{+400}_{-400} | 12700 \pm 200 |
| 15 | 8669.0 - 8869.0 | 12600^{+700}_{-600} | 12300^{+400}_{-500} | 12800^{+400}_{-400} | 12900^{+400}_{-400} | 12800^{+700}_{-700} | 12700 \pm 200 |
| 16 | 8869.0 - 9069.0 | 13100^{+800}_{-900} | 12800^{+400}_{-500} | 12300^{+600}_{-500} | 12700^{+400}_{-300} | 13000^{+600}_{-600} | 12700 \pm 200 |