fig12, ax12 = pl.subplots(num=12)
ax12.errorbar(cat['v_rms'][hot]*np.sqrt(8*np.log(2))*gcorfactor[hot], [149]*hot.sum(),
lolims=True, linestyle='none', capsize=0, alpha=0.3,
marker='^', color='r')
for mask,color,alpha,markersize in masks_colors:
ax12.errorbar(cat['v_rms'][mask]*np.sqrt(8*np.log(2))*gcorfactor[mask], cat['temperature_chi2'][mask],
#yerr=[cat['elo_t'][mask], cat['ehi_t'][mask]],
markersize=10 if any(mask & is_leaf) else 5,
markeredgecolor='none',
linestyle='none', capsize=0, alpha=alpha, marker='.', color=color)
ax12.set_xlabel(r"Line FWHM (km s$^{-1}$)")
ax12.set_ylabel("H$_2$CO Temperature (K)")
linewidths = np.linspace(0,2.35*cat['v_rms'].max())*u.km/u.s
ax12.plot(linewidths*2.35, [heating.tkin_all(10**4*u.cm**-3, sigma, 10*u.pc,
5*u.km/u.s/u.pc, 30*u.K)
for sigma in linewidths],
linestyle='--', color='k', label='$n=10^4$ cm$^{-3}$', zorder=-5)
ax12.plot(linewidths*2.35, [heating.tkin_all(10**4*u.cm**-3, sigma, 10*u.pc,
1*u.km/u.s/u.pc, 30*u.K)
for sigma in linewidths],
linestyle='--', color='r', label='$n=10^4$ cm$^{-3}$, $dv/dr=1$', zorder=-5, linewidth=2, alpha=0.5)
ax12.plot(linewidths*2.35, [heating.tkin_all(10**4*u.cm**-3, sigma, 20*u.pc,
5*u.km/u.s/u.pc, 30*u.K)
for sigma in linewidths],
linestyle='--', color='b', label='$n=10^4$ cm$^{-3}$, $L=20$ pc', zorder=-5, alpha=0.5, linewidth=2)
ax12.plot(linewidths*2.35, [heating.tkin_all(10**5*u.cm**-3, sigma, 10*u.pc,
5*u.km/u.s/u.pc, 30*u.K)
for sigma in linewidths],
linestyle=':', color='k', label='$n=10^5$ cm$^{-3}$', zorder=-5)
ax12.plot(linewidths*2.35, [heating.tkin_all(10**6*u.cm**-3, sigma, 10*u.pc,
linewidth=0.5,
alpha=0.6,
color='r')
mask = maps & lolim_conservative
ax7.plot(pcfittable['width'][mask] * (8 * np.log(2))**0.5,
pcfittable['tmin1sig_chi2'][mask],
marker='^',
markersize=10,
markeredgecolor='none',
color='r',
alpha=0.4,
linestyle='none')
linewidths = np.linspace(0, pcfittable['width'].max()) * u.km / u.s
ax7.plot(linewidths * 2.35, [
heating.tkin_all(10**4 * u.cm**-3, sigma, 10 * u.pc, 5 * u.km / u.s / u.pc,
30 * u.K) for sigma in linewidths
],
linestyle='--',
color='k',
label='$n=10^4$ cm$^{-3}$',
zorder=-5)
ax7.plot(linewidths * 2.35, [
heating.tkin_all(10**4 * u.cm**-3, sigma, 10 * u.pc, 1 * u.km / u.s / u.pc,
30 * u.K) for sigma in linewidths
],
linestyle='--',
color='r',
label='$n=10^4$ cm$^{-3}$, $dv/dr=1$',
zorder=-5,
linewidth=2,
alpha=0.5)
markersize=10 if any(mask & is_leaf) else 5,
markeredgecolor="none",
linestyle="none",
capsize=0,
alpha=alpha,
marker=".",
color=color,
)
ax12.set_xlabel(r"Line FWHM (km s$^{-1}$)")
ax12.set_ylabel("H$_2$CO Temperature (K)")
linewidths = np.linspace(0, 2.35 * cat["v_rms"].max()) * u.km / u.s
ax12.plot(
linewidths * 2.35,
[
heating.tkin_all(10 ** 4 * u.cm ** -3, sigma, 10 * u.pc, 5 * u.km / u.s / u.pc, 30 * u.K)
for sigma in linewidths
],
linestyle="--",
color="k",
label="$n=10^4$ cm$^{-3}$",
zorder=-5,
)
ax12.plot(
linewidths * 2.35,
[
heating.tkin_all(10 ** 4 * u.cm ** -3, sigma, 10 * u.pc, 1 * u.km / u.s / u.pc, 30 * u.K)
for sigma in linewidths
],
linestyle="--",
color="r",
