Syntax:
-Injecttransit <"Pspec" | "Pfix" per | "Prand" minp maxp | "Plogrand" minp maxp | "randfreq" minf maxf | "lograndfreq" minf maxf> <"Rpspec" | "Rpfix" Rp | "Rprand" minRp maxRp | "Rplogrand" minRp maxRp> <"Mpspec" | "Mpfix" Mp | "Mprand" minMp maxMp | "Mplogrand" minMp maxMp> <"phasespec" | "phasefix" phase | "phaserand> <"sinispec" | "sinifix" sin_i | "sinirand"> <"eomega" <"espec" | "efix" e | "erand"> <"ospec" | "ofix" omega | "orand"> | "hk" <"hspec" | "hfix" h | "hrand"> <"kspec" | "kfix" k | "krand">> <"Mstarspec" | "Mstarfix" Mstar> <"Rstarspec" | "Rstarfix" Rstar> <"quad" | "nonlin"> <"ldspec" | "ldfix" ld1 ... ldn> ["dilute" <"spec" | "fix" dilute>] omodel [modeloutdir] Example 1.$ ./vartools -i EXAMPLES/3 -oneline -randseed 1 \ -Injecttransit lograndfreq 0.2 2.0 \ Rpfix 1.0 Mpfix 1.0 \ phaserand sinirand \ eomega efix 0. ofix 0. \ Mstarfix 1.0 Rstarfix 1.0 \ quad ldfix 0.3471 0.3180 \ 1 EXAMPLES/OUTDIR1 \ -BLS q 0.01 0.1 0.5 5.0 20000 200 7 1 0 1 EXAMPLES/OUTDIR1 \ 1 fittrap Name = EXAMPLES/3 Injecttransit_Period_0 = 0.72240757 Injecttransit_Rp_0 = 1.00000 Injecttransit_Mp_0 = 1.00000 Injecttransit_phase_0 = 0.39438 Injecttransit_sin_i_0 = 0.96564 Injecttransit_h_0 = 0.00000 Injecttransit_k_0 = 0.00000 Injecttransit_Mstar_0 = 1.00000 Injecttransit_Rstar_0 = 1.00000 Injecttransit_ld_1_0 = 0.34710 Injecttransit_ld_2_0 = 0.31800 BLS_Period_1_1 = 0.72239761 BLS_Tc_1_1 = 53725.889333375417 BLS_SN_1_1 = 37.02975 BLS_SR_1_1 = 0.00159 BLS_SDE_1_1 = 5.93220 BLS_Depth_1_1 = 0.00935 BLS_Qtran_1_1 = 0.05315 BLS_Qingress_1_1 = 0.19369 BLS_OOTmag_1_1 = 10.16683 BLS_i1_1_1 = 0.96375 BLS_i2_1_1 = 1.01691 BLS_deltaChi2_1_1 = -10851.70035 BLS_fraconenight_1_1 = 0.30928 BLS_Npointsintransit_1_1 = 179 BLS_Ntransits_1_1 = 10 BLS_Npointsbeforetransit_1_1 = 227 BLS_Npointsaftertransit_1_1 = 186 BLS_Rednoise_1_1 = 0.00164 BLS_Whitenoise_1_1 = 0.00490 BLS_SignaltoPinknoise_1_1 = 14.75474 BLS_Period_invtransit_1 = 1.25145489 BLS_deltaChi2_invtransit_1 = -2806.72319 BLS_MeanMag_1 = 10.16718
Inject a transit into the light curve "EXAMPLES/3" and recover it using -BLS. We draw the period from a uniform-log random distribution in frequency with limits of 0.2 c/d and 2.0 c/d. We fix the planet radius to 1.0 R_J and the planet mass to 1.0 M_J. We adopt a random phase, and draw sin(i) from a random orientation distribution (subject to the constraint that a transit must occur). We fix the eccentricity and argument of periastron to 0, we fix the star mass and radius to 1.0 M_sun and 1.0 R_sun, respectively, and we adopt a quadratic limb darkening law, with coefficients 0.3471 and 0.3180. We output the model to the directory EXAMPLES/OUTDIR1 (the filename will be EXAMPLES/OUTDIR1/3.injecttransit.model). After injecting the transit we use the -BLS command to recover the transit.
Fig 1. Light Curve with Transit Injected from in Example 1.