STScI Logo

using_timcor xray.xtiming.timcor



timcor -- package of ROSAT timing correction tasks


This document describes the tasks available in the timcor package.


This is not a conventional help file, but a brief writeup to introduce users to the set of prototype timing correction tasks recently installed in PROS.

The timing correction package is currently composed of three tasks:

	1) apply_bary - apply the barycenter time corrections to the 
	                time-sorted QPOE file

	2) calc_bary - calculate the barycenter timing corrections
	3) scc_to_utc - convert from ROSAT space-craft clock to UTC

PLEASE NOTE: The original code for this package was written at MPE. It was ported and modified from a FORTRAN/EXSAS/MIDAS environment to an SPP/PROS/IRAF environment.

These tasks are written SPECIFICALLY for ROSAT data. They will not work correctly on data from other missions, i.e. Einstein. There are important differences in the data:

	1) Einstein data doesn't have the split orbit ephemeris file 
           needed by "calc_bary".

	2) The Einstein data format differs from the ROSAT data format

	3) There are spacecraft dependencies in "calc_bary" in that 
	   the times are adjusted to account for the spacecraft clock
We look forward to any comments/feedback that you might have and also welcome suggestions for both short and long term improvements. (email

The following is an example of how to run barycenter timing corrections
from start to finish:

STEP 1: Run the task "timsort"

cl> xray
xr> xtiming
ti> timsort
input qpoe file name: xdata$rp110590n00.qp[time=(1606720.00:1607479.00)]
Output time-ordered file name: parlac
source region descriptor: c 6942.0 9375.0 1.5
Output time-ordered bkgd file name (NONE): 
starting source creation - parlac_sti.qp

STEP 2: Run the task "calc_bary"

xr> xtiming
xt> timcor
ti> calc_bary
input ephemeris table: xdata$
root name for output file(s) [root_cor(ephem).tab]: parlac
input RA (hh:mm:ss.s): 22:08:41.60
input DEC (dd:mm:ss.s): 45:44:47.30

RDF input ... No conversion!!
Output ephemeris table:

Computing Correction table ...
Opening file:
Output Correction table:

	The output table "" will be used as input to the task
"apply_bary".  The ra/dec input to calc_bary is written to the table header, 
and is read by apply_bary. The correction table contains four columns:

            1) integer part of uncorrected orbit time UT time
            2) real part of uncorrected orbit time in UT time
            3) integer part of corrected orbit time UT time
            4) real part of corrected orbit time in UT time

STEP 3: Run the task "apply_bary"

cl> xray
xr> xtiming
xt> timcor
ti> apply_bary
input qpoe file name: parlac_sti.qp
input correction table filename:
output qpoe file name: parlac_tcor.qp
event definition (large): 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. 2448062.98482851,  SCC 1607479.00000000
first corr table time = J.D. 2448062.97500000
last corr table time  = J.D. 2448063.04722222

Barycenter Corrections applied to event times in parlac_tcor.qp
using correction table 
with reference position at RA =  22:08:41.60, DEC =  45:44:47.30
Creating QPOE output file: parlac_tcor.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:

        MJDREFI + MJDREFF + (photon_time / 86400)

The values MJDREFI and MJDREFF can be found in the header of the
output QPOE file. (MJDREFF should be 0.5)

The file parlac_tcor.qp is barycenter corrected and can be used
as input to the xtiming analysis tasks in the xray package.  See
the help files for individual tasks for further details.

Search Form · STSDAS