label='$%s~\mathrm{Gyr}$~$\mathrm{M\&H~(2008)}$' %ages[i], zorder=0) xs, ys = iso_calc(pars3, ages[i]) xs -= .45 pl.plot(xs, ys, color='k', linestyle='-', linewidth=lw, label = '$%s~\mathrm{Gyr}$~$\mathrm{Angus~\emph{et~al.}~(2015)}$' % ages[i], zorder=0) xs, ys1 = iso_calc(pars3-pars3_err, ages[i]) xs -= .45 xs, ys2 = iso_calc(pars3+pars3_err, ages[i]) xs -= .45 # pl.fill_between(xs, ys1, ys2, facecolor='0.5', alpha=0.3, edgecolor='None', # zorder=0) from prob_uncerts import sigmas samples = sigmas("ACHF") for j in range(np.shape(samples)[1]): par_samp = np.zeros(4) par_samp[:3] = samples[:, j] par_samp[3] = 0.45 xs, ys = iso_calc(par_samp, ages[i]) xs -= 0.45 pl.plot(xs, ys, color="k", zorder=0, alpha=.04) pl.xlabel("$\mathrm{B-V-}~c$") pl.ylabel("$\mathrm{P_{rot} (days)}$") pl.xlim(10**-3, 1.4) pl.ylim(1, 10**2) pl.legend(loc='upper left') pl.loglog() pl.savefig("/Users/angusr/Python/Gyro/gyro_paper/p_vs_bv%s.pdf" % i)
pl.clf() sun = a==4.568 bv -= 0.45 print bv_err[l], "\n" print t_err[l] print logg_err[l] print p_err[l], "\n" pl.errorbar(bv[l], p[l], xerr=bv_err[l], yerr=p_err[l], color='k', \ fmt='o', mec='k', capsize=0, markersize=5, ecolor='.6') pl.errorbar(bv[l2], p[l2], xerr=bv_err[l2], yerr=p_err[l2], color='k', \ fmt='o', mec='.7', capsize=0, markersize=5, ecolor='.6') pl.errorbar(bv[sun], p[sun], xerr=bv_err[sun], yerr=p_err[sun], color='r', \ fmt='o', mec='r', capsize=0, markersize=10, ecolor='.6') from prob_uncerts import sigmas samples = sigmas(fname) for i in range(np.shape(samples)[1]): par_samp = np.zeros(4) par_samp[:3] = samples[:, i] par_samp[3] = 0.45 xs, ys = iso_calc(par_samp, age) xs -= 0.45 pl.plot(xs, ys, color=c, zorder=0, alpha=.1) # Add Isochrones xs, ys = iso_calc(pars, age) xs -= 0.4 c = "k" pl.plot(xs, ys, color=c, linestyle='-.', linewidth=lw, \ label='$%s~\mathrm{Gyr}$~$\mathrm{(Barnes~2007)}$' %age, zorder=0) xs, ys = iso_calc(pars2, age)