Tutorial 3: Transits with GPs

See on GitHub.

Imagine the star from the previous example would have been a bit more active, or your telescope would have suffered from some systematics. Then your data might have looked more like this:

You can download the data file here: Leonardo.csv

As an experienced allesfitter, the team asks you to model the signal. The discovery report gives you a first guess for the transit signal:

  • Epoch: 1.09 +- 0.01 days after start of observations
  • Period: 3.41 +- 0.01 days
  • R_planet / R_star: 0.10 +- 0.01
  • (R_star + R_planet) / semi-major axis: somewhere between 0.1 and 0.3
  • R_star: 1.00+-0.01 R_sun
  • M_star: 1.00+-0.01 M_sun
  • T_eff: 5700+-100 K

Now, time to use a GP baseline! Start up the GUI -- or simply update the baseline settings and params in your template settings.csv and params.csv files from the last tutorial to make the following changes:




baseline_gp_offset_flux_Leonardo,0,1,uniform -0.01 0.01,$\mathrm{gp \ln \sigma (Leonardo)}$,
baseline_gp_matern32_lnsigma_flux_Leonardo,-5,1,uniform -15 0,$\mathrm{gp \ln \sigma (Leonardo)}$,
baseline_gp_matern32_lnrho_flux_Leonardo,0,1,uniform -1 15,$\mathrm{gp \ln \rho (Leonardo)}$,

Now run the fit. Success? Brilliant! You got great transit parameters despite all the red noise, and that allowed TMNT to schedule follow-up observations. You are now the hero of your team! They even baked you a cake! And even better: they nicely prepared all the old and new data for you in tutorial 04_transits_and_rvs:

  • Leonardo.csv (discovery photometry)
  • Michelangelo.csv (follow-up photometry)
  • Donatello.csv (decent RV data)
  • Raphael.csv (good RV data)

Now let’s allesfit all that data! Move on to the next tutorial and learn how to tackle all these data sets from scratch.