The barycentric corrections computed in 8.7.1 are applied to the event times in a time-sorted QPOE file. In addition, the times are converted from spacecraft-clock units to UTC.
input files: rp90_sti.qp, rp110590n00_cor.tab
output files: rp90_out.qp
ti> apply_bary input qpoe file name: rp90_sti.qp[time=(1606720.0:1612809.8)] # see Note 1&2 input correction table filename: rp110590_cor.tab output qpoe file name: rp90_out event definition: medium SCC->UTC conversion table is calibrated to SCC 110883482.0000 first qpoe gti time = J.D. 2448062.97604379, SCC 1606720.00000000 last qpoe gti time = J.D. 2448063.04652758, SCC 1612809.80000000 first corr table time = J.D. 2448062.97569444 last corr table time = J.D. 2448063.04652778 Barycenter Corrections applied to event times in rp90_out.qp using correction table rp110590_cor.tab with reference position at RA = 22:08:41.60, DEC = 45:44:47.30 Creating QPOE output file: rp90_out.qp
NOTE: The output QPOE times are in seconds from noon of the MJD of the start of the first OBI.
To convert the times to MJD use the following formula:
The values MJDREFI and MJDREFF can be found in the header of the output QPOE file. (MJDREFF should be 0.5)
NOTE 1: The orbit data records in the root_eph.tab file are generated every 60 seconds during the observation. The apply_bary task calculates the barycenter timing corrections using a linear interpolation between the two orbit records for which an event time falls. Occasionally there is not a final orbit record, so that there may be events with arrival times occurring after the last orbit record. When calculating the barycenter timing corrections for these events, the apply_bary task uses a linear extrapolation from the previous two orbit records and gives the user a warning that an extrapolation is being used. This extrapolation should effect at most 60 seconds of data. The apply_bary task prints the times of the first and last event and the first and last orbit record, so that the user may assess the appropriateness of the extrapolation being performed.
For this dataset, there are 62 events that arrive up to 7 seconds after the last orbit record. The example above uses a time filter on the input qpoe file to filter out these events. The user may also use the qpcopy task to perform the filtering.
NOTE 2: The apply_bary task takes as input a time-sorted qpoe file, and allows the user to specify a region descriptor. There is a known bug when using large region descriptors with time-sorted qpoe files. The error usually manifests itself in the form of ``ERROR: PLIO: reference out of bounds on mask'' or ``ERROR: PLIO: stack overflow (plr_open)''. Since the standard qpoe file is spatially sorted, the user is encouraged to apply region descriptors when running the timsort task (as in example 8.3.1).