IGM transmission example

Plot IGM transmission curves (Inoue+14 & Madau+96

/home/runner/work/synthesizer/synthesizer/src/synthesizer/emission_models/transformers/igm.py:317: RuntimeWarning: divide by zero encountered in power
  * (lam_obs[cond0 & ~cond1] / lam_l) ** (-0.3)
/home/runner/work/synthesizer/synthesizer/src/synthesizer/emission_models/transformers/igm.py:319: RuntimeWarning: divide by zero encountered in power
  - 0.2905 * (lam_obs[cond0 & ~cond1] / lam_l) ** (-0.3)
/home/runner/work/synthesizer/synthesizer/src/synthesizer/emission_models/transformers/igm.py:312: RuntimeWarning: invalid value encountered in cast
  tau_lc_dla_value[cond0 & ~cond1] = (
/home/runner/work/synthesizer/synthesizer/src/synthesizer/emission_models/transformers/igm.py:434: RuntimeWarning: overflow encountered in exp
  transmission = np.exp(-tau)

import cmasher as cmr
import matplotlib.pyplot as plt
import numpy as np
from unyt import angstrom

from synthesizer.emission_models.attenuation import Inoue14, Madau96

lam = np.arange(0, 20000) * angstrom


redshifts = [3.0, 5.0, 7.0]
colors = cmr.take_cmap_colors("cmr.guppy", len(redshifts))

for IGM, ls in zip([Inoue14, Madau96], ["-", ":"]):
    igm = IGM()
    for z, color in zip(redshifts, colors):
        plt.plot(
            lam,
            igm.get_transmission(z, lam),
            ls=ls,
            c=color,
            label=f"{igm.name} z={z}",
        )

plt.legend()
plt.xlabel(r"$\lambda_{obs}/\AA$")
plt.ylabel(r"$T$")
plt.ylim([0, 1.1])
plt.show()
# plt.savefig('../docs/source/images/img.png', bbox_inches='tight', dpi=200)