NAME

fakelc -- creates a fake light curve given a mean count rate. The output may be either an event list or a binned light curve.


USAGE

        fakelc outfile type binsize binunit span spanunit cts ctsunit noise

DESCRIPTION

This task creates a single light curve with a given mean count rate. The output light curve may take the form of either an event list or a binned time series. The ouput may be either an uninterrupted light curve, or one containing data gaps. The length of the continuous segments is random, given the average duration of the data gaps. Continuous segments and data gaps resulting from data sampling typcially obtained from satellites in low earth orbit may also be chosen.

The user may also specify sampling appropriate for the sampling from one HEXTE cluster on XTE. The output FITS file conforms to the OGIP standards for a timing file. No energy information is included in the output.

The mean intensity may be specified in either counts or counts/s. Noise may be added in the form of Poisson statistics or Gaussian statistics. If the latter, a sigma for the distribution is input.


PARAMETERS

outfile [file name]
The name of the output FITS file. This file complies with the OGIP Timing FITS File Format.

type = E or B [string]
The type of light curve to be produced: E for an event list; B for a binned light curve.

binsize [double precision real]
The sampling to be used for the fake light curve, in units of binunit. If a binned light curve, then this is the bin size. For an event list, it is the accuracy of the time stamp. For an event list, multiple events occuring within the time stamp accuracy are retained. Hence, there may be multiple time stamps in the event list.

binunit = s or d [string]
Time unit to be used in the light curve - seconds (s) or days (d). Note that for an event list, TUNIT also determines the time unit in the count rate: count/s or count/d.

span [double precision real]
The total live time for the fake observation.

spanunit = s or d [string]
Time unit to be used for the total live time of the fake observation, in seconds (s) or days (d). This time unit is also used for gaplength.

cts [real*4]
The mean intensity for the light curve, in units determined by the ctsunit parameter.

ctsunit = count or count/s [string]
The units on the input intensity. For an event list, if count is used, a rate is computed using the BINSIZE. Also, for an event list the time unit for a rate is determined by BINUNIT, so that the rate may be either count/s or count/d.

noise [integer]
The noise option to be used in producing a binned light curve. The integer values take on the following meanings: 0 = no noise; 1 = Poisson noise; 2 = Gaussian noise. The Poisson statistics are governed by the mean intensity as expressed in counts. (If CTSUNIT = count/s, the conversion to counts is done internally using the sample size.) If Gaussian noise is specified, the sigma for the Gaussian distribution is specified using the sigma parameter. For event lists, this parameter is ignored.

(sigma) [real*4]
The sigma to be used for Gaussian noise.

(tstart) [double precision real]
The absolute start time of the fake observation, nominally in MJD, but may be input as zero. The default value is MJD 49352.000696574074 (Jan 1, 1994, 0 hr, 1 m, 0.184 s), the reference time for the XTE mission.

(hexte = no) [logical]
A switch to utilize the XTE HEXTE standard modulation of 32-s on source and 32-s off source for a single cluster.

(nsegs = 1) [integer]
The number of data segments in the light curve. Each data segment is a continuous time series. If nsegs > 1, then gaps will be generated in the time axis. The output FITS file has at least a TIME and RATE column (and an ERROR column if noise > 0).

(gaplength) [double precision real]
The average length of each data gap for nsegs > 1. The units are those given by the SPANUNIT parameter. This average length is used to generate segments of random duration and gaps of random size. For nsegs = 1, this parameter has no effect.

(orbit = no) [boolean]
An option to assign times according to the sampling for a typical XTE orbit.

(seed = -3267) [integer]
The initial seed for the random number generator. Note that this must be a negative integer.


EXAMPLES

1. Create a binned fake light curve using a constant value intensity of 50 count/s. The bin size is set to 0.25 s, and the total time span is 864 s in one segment of data. No noise is added, and output is directed to test.lc

     fakelc outfile=test.lc type=B binsize=0.25 binunit=s cts=50
     ctsunit=count/s noise=0 span=864 spanunit=s 

2. Create a binned fake light curve with a bin size of 86.4 s. The total live time is 1 day, but the data is broken up into 10 segments having an average data gap of 0.05 d. The mean intensity is 50, with units in counts. Poisson noise is added. Output is directed to the file test2.lc.

     fakelc outfile=test2.lc type=B binsize=86.4 binunit=s cts=50
     ctsunit=count noise=1 span=1 spanunit=d  nsegs=10 gaplength=0.05

3. Create an event list in which the average count rate is 0.5 count/s, and the timing accuracy is 1 s. A total time span of 1 day is divided into 5 intervals with an average data gap length of 0.1 days. Output is directed to test3.lc

     fakelc outfile=test3.evt type=E binsize= 1.0 binunit=s cts=0.5
     ctsunit=count/s noise=0 span=0.5 spanunit=d nsegs=5 gaplength=0.1


NOTES:


BUGS

Version FAKELC_V3.3

Please report problems to xtehelp@athena.gsfc.nasa.gov.


SEE ALSO

addsine addshots

CATEGORY

Jan95 ftools.xte