Syntax:

-restricttimes
    ["exclude"]
    < "JDrange" minJD maxJD |
      "JDrangebylc"
        < "fix" minJD | "list" ["column" col] | "fixcolumn" <colname | colnum> |
          "expr" expression >
        < "fix" maxJD | "list" ["column" col] | "fixcolumn" <colname | colnum> |
          "expr" expression >
      "JDlist" JDfilename |
      "imagelist" imagefilename |
      "expr" expression >

$ ./vartools -i EXAMPLES/3 -stats t min,max \
             -restricttimes JDrange 53740 53750 \
             -stats t min,max -oline

Name                  = EXAMPLES/3
STATS_t_MIN_0         = 53725.173920000001
STATS_t_MAX_0         = 53756.281021000003
RestrictTimes_MinJD_1 = 53740
RestrictTimes_MaxJD_1 = 53750
STATS_t_MIN_2         = 53740.336210000001
STATS_t_MAX_2         = 53745.478681000001

Use the vartools -restricttimes command to remove points outside the range 53740 < JD < 5370 from the light curve. The calls to -stats before and after the call to -restricttimes show how the command affects the timespan of the light curve.

$ ./vartools -i EXAMPLES/3 -stats mag min,max \
             -restricttimes expr '(mag>10.16311)&&(mag<10.17027)' \
             -stats mag min,max -oneline

Name            = EXAMPLES/3
STATS_mag_MIN_0 = 10.141400000000001
STATS_mag_MAX_0 = 10.1921
STATS_mag_MIN_2 = 10.163119999999999
STATS_mag_MAX_2 = 10.170260000000001

Filter the light curve to remove any points outside the range 10.16311 < mag < 10.17027. The calls to -stats before and after the call to -restricttimes show how the command affects the light curve range.

$ ./vartools -i EXAMPLES/3 -stats mag pct20.0,pct80.0 \
             -restricttimes expr \
               '(mag>STATS_mag_PCT20_00_0)&&(mag<STATS_mag_PCT80_00_00)' \
             -stats mag min,max -oneline

Name                 = EXAMPLES/3
STATS_mag_PCT20.00_0 = 10.16311
STATS_mag_PCT80.00_0 = 10.17027
STATS_mag_MIN_2      = 10.163119999999999
STATS_mag_MAX_2      = 10.170260000000001

Same as Example 2, but here we make the cut in terms of the magnitude percentiles, which are calculated with the -stats command, rather than giving the magnitude range explicitly.

$ ./vartools -i EXAMPLES/3.transit -oneline \
             -BLS q 0.01 0.1 0.1 20.0 100000 200 0 1 \
                  0 0 0 fittrap nobinnedrms \
             -expr 'minph=0.5-BLS_Qtran_1_0/2.0' \
             -expr 'maxph=0.5+BLS_Qtran_1_0/2.0' \
             -expr 'ph=t' \
             -changevariable t ph \
             -Phase bls T0 bls 0.5 \
             -restricttimes exclude JDrangebylc expr minph expr maxph \
             -o EXAMPLES/3.cliptransit columnformat t,mag,err,ph

Name                         = EXAMPLES/3.transit
BLS_Period_1_0               =     2.12334706
BLS_Tc_1_0                   = 53727.297293937358
BLS_SN_1_0                   =   7.26127
BLS_SR_1_0                   =   0.00238
BLS_SDE_1_0                  =   6.34195
BLS_Depth_1_0                =   0.01220
BLS_Qtran_1_0                =   0.03576
BLS_Qingress_1_0             =   0.19618
BLS_OOTmag_1_0               =  10.16686
BLS_i1_1_0                   =   0.98213
BLS_i2_1_0                   =   1.01790
BLS_deltaChi2_1_0            = -24217.21939
BLS_fraconenight_1_0         =   0.43155
BLS_Npointsintransit_1_0     =   165
BLS_Ntransits_1_0            =     4
BLS_Npointsbeforetransit_1_0 =   127
BLS_Npointsaftertransit_1_0  =   143
BLS_Rednoise_1_0             =   0.00151
BLS_Whitenoise_1_0           =   0.00489
BLS_SignaltoPinknoise_1_0    =  14.38935
BLS_Period_invtransit_0      =     1.14594782
BLS_deltaChi2_invtransit_0   = -3301.69183
BLS_MeanMag_0                =  10.16740
RestrictTimes_MinJD_6        = 0.48211760476606086
RestrictTimes_MaxJD_6        = 0.5178823952339392

Use the vartools -restricttimes command together with the -BLS, -expr, -changevariable, and -Phase commands to remove points in transit from the light curve of 3.transit. First the -BLS command is used to detect the transit signal and determine its ephemeris. The following two -expr commands calculate the minimum and maximum transit phase assuming the transit is at phase 0.5. Here the variable "BLS_Qtran_1_0" stores the fractional transit duration = T14/P determined from BLS. The third -expr command initializes the ph variable to be equal to the times. The -changevariable command switches the variable to use as storing times in the subsequent commands to ph. This is done so that the -Phase and -restricttimes commands will operate on ph, leaving the original time measurements in t. The -Phase command phases the light curve using the previous BLS command for the period and for the transit epoch, and setting the transit center to occur at phase 0.5. We then issue the -restricttimes command using the "exclude" keyword to specify the points to remove from the light curve rather than the points to keep in the light curve. The "JDrangebylc" keyword is needed to let the range be taken from analytic expressions. In this case we use the contents of the variables minph and maxph, computed above, define the range of points that will be excluded from the light curves. The resulting light curve is output to EXAMPLES/3.cliptransit using the columnformat keyword to output the times, magnitudes, errors and phases of the points. Note that the light curve will be sorted by phase rather than time. To re-sort by time, one could use a "-changevariable t t" command before the "-o" command.