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.

Example Injected Transit

Fig 1. Light Curve with Transit Injected from in Example 1.