コード例 #1
0
    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,
コード例 #2
0
    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)
コード例 #3
0
            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",