wls = np.linspace(1.5, 16, 500)
plt.plot(wls, qcl.full_gain_spectrum(wls))
plt.axhline(0, ls='--', lw=0.5)
plt.xlabel('Wavelength (μm)')
plt.ylabel('Gain (cm$^{-1}$)')
plt.show()