from msaexp import spectrum
import msaexp.utils

sp = spectrum.SpectrumSampler("https://s3.amazonaws.com/msaexp-nirspec/extractions/smacs0723-ero-v4/smacs0723-ero-v4_g395m-f290lp_2736_6355.spec.fits")

z = 7.665

fig, axes = sp.multiplot(
    ny=4,
    sx=8, sy=2,
    z=z,
    color='k', alpha=0.5
)

# Overplot line list
li = msaexp.utils.lines.LineList()
for ax in axes:
    li.add_to_axis(ax, alpha=0.5)

axes[-1].set_xlabel(f'rest wavelength, z={z:.3f}')
fig.tight_layout(pad=1)