Syntax:
-SoftenedTransit <bls | blsfixper | P0 T00 eta0 delta0 mconst0 cval0> fitephem fiteta fitcval fitdelta fitmconst correctlc omodel [model_outdir] fit_harm [<"aov" | "ls" | "bls" | "spec" Pharm> nharm nsubharm] Example 1.$ ./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.
Fig 1. Softened Transit Model fit to the light curve in Example 1.