juliet icon indicating copy to clipboard operation
juliet copied to clipboard

possibility to know evry single transit depth (each r_1 and r_2)

Open MohammadEftekhar opened this issue 2 years ago • 7 comments

Dear Nestor

Is it possible to reach every single transit depth, in order to every single r_1 and r_2?

thank you Mohammad

MohammadEftekhar avatar Oct 26 '21 16:10 MohammadEftekhar

Hi Mohammad,

Not yet, but I've thought of implementing this though. Should be very similar (and easier than) the TTV implementation.

N.

nespinoza avatar Mar 16 '22 21:03 nespinoza

This would be very useful, indeed. I also think that there should be a way of outputting the duration of the fitted transit model as well! (at some point, even the duration of single transits for TTV, TDVs..)

LucaNap avatar Jun 01 '24 09:06 LucaNap

Hi folks! We've already implemented in the latest version of juliet the possibility to fit different transit depths simultaneously (thanks so @Jayshil's efforts!). You would do so by doing something like:

p_p1_TESS     uniform 0., 1.
p_p1_JWST    uniform 0., 1.

And so on. I haven't had time to document this, but likely will do so in the next few months! Will be closing this once docs are updated, but will leave the issue open for the time being.

N.

nespinoza avatar Jun 01 '24 19:06 nespinoza

Great, @nespinoza ! Did you mean:

p_p1_TESS_0
p_p1_TESS_1
...

like the transit times, or did you mean to fit different depths only for different instruments? And if it's exactly like the transit times (T_p1_TESS_0... etc), what about the r1, r2 parametrization?

By the way, having at least the transit duration as an optional output of any transit fitting, would be extremely valuable.

LucaNap avatar Jun 02 '24 17:06 LucaNap

Hi @LucaNap,

No, the implementation I described above fits depths on different instruments, which might have common transit parameters. As in:

P_p1                 normal 4.1,0.01
t0_p1                normal 2455242.,0.01
.
.
.
p_p1_TESS1     uniform 0., 1.
p_p1_TESS2    uniform 0., 1.

In theory, with this framework you can still do what you want to do --- the t0's can come from different instruments too, so you can divide first manually a lightcurve into different chunks, and fit a different depth and time-of-transit for each lightcurve (assuming they are "different instruments"). You can link, e.g., mflux, q1, q2, etc., to be the same for different instruments (e.g., mflux_TESS1_TESS2) and that would be equivalent to what you want to do.

N.

nespinoza avatar Jun 03 '24 21:06 nespinoza

As for the r1 and r2 parameters --- no, only p can be done like those.

nespinoza avatar Jun 03 '24 21:06 nespinoza

And as for the question about outputting the transit duration, it is useful, but I prefer to leave that to the user for now :).

nespinoza avatar Jun 03 '24 21:06 nespinoza