PERDGRM (Jul96) ftools.xte PERDGRM (Jul96)
NAME
perdgrm -- compute the periodogram (a discrete fourier analysis) on
light curves which are comprised of binned intensities at unevenly
spaced time intervals. This routine uses the definition and
normalization of the periodogram given in Scargle, J. D. 1982,
ApJ, 263, 835.
USAGE
perdgrm infile outfile times rms
DESCRIPTION
This task performs fourier analysis using the algorithm and
normalization given by Scargle, J. D. 1982, ApJ, 163, 835. The
algorithm works with any binned time series, but is particularly
useful for series in which the bins are unevenly spaced in time.
Rather than fill data gaps with the mean intensity, the algorithm
uses a periodogram technique for determining the strength of
sinusoidal variations present in the time series. The algorithm
utilizes a normalization which allows for the computation of "false
alarm probabilities" of observed peaks in the power spectrum.
This task computes the power spectrum and false alarm probabilities
for an input light curve. The light curve may be broken up into
intervals, and the intervals separately analyzed. Because in
general the frequency sampling of the intervals may differ,
averaging the power spectra from different intervals is not
straight-forward and is not done in this routine. Instead, power
spectra from diffent intervals are recorded in separate FITS
extensions. The user may also modify the frequency resolution,
determine the criteria for subtracting a polynomial from the input
light curve, and modify the upper and lower period limits in which
to compute the power spectrum. The user may also specify whether
the window function is to be computed.
PARAMETERS
infile [file name]
The name of the input FITS file. A file following the OGIP
Timing FITS File format is expected. If no TIME column exists,
the routine constructs the times from the TIMEZERO and TIMEDEL
keywords.
outfile [file name]
The name of the output FITS file. Results from different time
intervals are placed in separate extensions.
times [string]
A list of paired start and stop time ranges for intervals to be
analyzed. A "-" will result in all the data being analyzed in a
single interval. Omitting a start or stop time in the first or
last interval will default to the earliest or latest time in
the file. Up to 15 time intervals may be specified.
(inres = 1) [integer]
The frequency resolution for the fourier transform. Higher
values of inres proportionally decreases the size of the
frequency steps, resulting in an increase in the number of
frequencies. The number of frequency steps is given by
inres * (max time - min time) * (mean freq. step)
where
(mean freq. step) = (1/(N-1)) [sum over i] 1 / T(i) - T(i-1))
rms [real*4]
The rms error to be achieved by subtracting a polynomial from
the light curve. By specifiying a negative value, the highest
degree of polynomial to be tried is abs(rms).
(inminP = 0) [real*4]
The minimum period to be searched, in the same time unit as the
time series. The INDEF value is 0.
(inmaxP = INDEF) [real*4]
The maximum period to be searched, in the same time unit as the
time series. The INDEF value is determined by the sampling in
the light curve.
(window = no) [boolean]
Whether a window function is to be computed. If so, it is
added as a column in the output extension.
(chatter = 9) [integer]
Amount of information to tell the user.
(timename = TIME) [string]
Name for the Time column in the input FITS file. Note that the
input string is case sensitive.
(ratename = RATE) [string]
Name for the Rate column in the input FITS file. Note that the
input string is case sensitive.
(errname = ERROR) [string]
Name for the Error column in the input FITS file. Note that the
input string is case sensitive.
(copyprime = yes) [boolean]
Whether to copy the primary header and array to the output.
(Not supported at this time.)
(copyall = no) [boolean]
If true, all other extensions in the input file will be copied
to the output (only in effect when 1 input file is specified).
However, this has not yet been implemented for SCARGLE.
EXAMPLES
1. Analyze an input light curve, subtracting a
polynomial of order 5 or less. Also compute the window function.
perdgrm sinetest2.lc sinetest2.ft - -5 window=yes
2. Analyze the input light curve in 3 intervals, subtracting a
polynomial of order 2 or less in each interval. The results are output to 3
separate extensions in the output file.
perdgrm sinetest.lc sinetest.ft 49000.001-49000.003,49000.003-49000.006,
49000.006- -2
NOTES:
The implementation of Scargle's algorithm embodied here is based on
routines developed by Koji Mukai and Alan Smale.
This routine can also be used with event lists, if the intensity at
each time stamp is given a value of 1.
BUGS
Version PERDGRM_V3.5.2.
Please report problems to xtehelp@athena.gsfc.nasa.gov.
SEE ALSO