MNGuenther / allesfitter

allesfitter is a convenient wrapper around the packages ellc (light curve and RV models), dynesty (static and dynamic nested sampling) emcee (Markov Chain Monte Carlo sampling) and celerite (Gaussian Process models).
MIT License
60 stars 36 forks source link

Double Lined (SB2) Binaries in Allesfitter #37

Closed JackActon96 closed 3 years ago

JackActon96 commented 3 years ago

Is there a way to fit double lined binaries within allesfitter, and if so what extra parameters/settings are needed to do this - apologies if this is in the documentation, but I couldn't find anything.

MNGuenther commented 3 years ago

Hi @JackActon96! Yes, this is readily possible and super easy.

All you need to do are the following touches:

Data files: save the secondary's RV series as a separate data file. So you would have, for example, "HARPS_Primary.csv" and "HARPS_Secondary.csv".

Settings file: Pass these infos into the settings file as "inst_rv,HARPS_Primary" and "inst_rv2,HARPS_Secondary".

Params file: Make sure you add two lines to fit for, namely "B_K" and "B_q", where K is the usual semi-amplitude of the primary, and "q" refers to the mass ratio between the two bodies.

I have made this more visible now in the website's tutorials tabs and with an explicit GitHub example.. Hope this helps and happy fitting!

(Sorry for the delay, this fell through the cracks while I was moving from MIT to ESA.)

MNGuenther commented 3 years ago

(Edited for clarity)