Syntax:

$ ./vartools -i EXAMPLES/3.transit -oneline \
    -BLS q 0.01 0.1 0.5 5.0 20000 200 7 1 0 0 0 \
    -SoftenedTransit bls 1 1 1 1 1 0 1 EXAMPLES/OUTDIR1 0

Name                         = EXAMPLES/3.transit
BLS_Period_1_0               =     2.12312625
BLS_Tc_1_0                   = 53727.297046247397
BLS_SN_1_0                   =  38.39425
BLS_SR_1_0                   =   0.00237
BLS_SDE_1_0                  =   4.77204
BLS_Depth_1_0                =   0.01136
BLS_Qtran_1_0                =   0.03000
BLS_i1_1_0                   =   0.98500
BLS_i2_1_0                   =   1.01000
BLS_deltaChi2_1_0            = -24130.93833
BLS_fraconenight_1_0         =   0.42662
BLS_Npointsintransit_1_0     =   146
BLS_Ntransits_1_0            =     4
BLS_Npointsbeforetransit_1_0 =   106
BLS_Npointsaftertransit_1_0  =   120
BLS_Rednoise_1_0             =   0.00156
BLS_Whitenoise_1_0           =   0.00490
BLS_SignaltoPinknoise_1_0    =  12.89679
BLS_Period_invtransit_0      =     1.14599569
BLS_deltaChi2_invtransit_0   = -3289.67397
BLS_MeanMag_0                =  10.16740
SoftenedTransit_Period_1     =     2.12322112
SoftenedTransit_T0_1         = 53727.29783160
SoftenedTransit_eta_1        =     0.06171206
SoftenedTransit_cval_1       =   -10.87159958
SoftenedTransit_delta_1      =    -0.01206461
SoftenedTransit_mconst_1     =    10.16686817
SoftenedTransit_chi2perdof_1 =    27.04335183

Fit a Protopapas et al. 2005 softened transit model to the light curve EXAMPLES/3.transit. We first identify the transit signal using -BLS and then use the results from BLS to initialize the parameters for -SoftenedTransit. We vary the ephemeris, eta, cval, delta and mconst. We do not subtract the best-fit model from the light curve. The model is output to the directory EXAMPLES/OUTDIR1 (the filename will be 3.transit.softenedtransit.model). We do not fit a harmonic series to the light curve simultaneously with the transit.