|
|
|
|
|
SynopsisSearch for variability using Gregory-Loredo algorithm. Syntaxglvary infile outfile lcfile effile [frac] [seed] [mmax] [mmin] [nbin] [mintime] [clobber] [verbose] Description`glvary' implements the Gregory-Loredo variability test algorithm. It splits the events into multiple time bins and looks for significant deviations. Input data includes an event file with good time data and a normalized effective area. Two output files are created: a table of odds ratios and a light curve file with includes +/-3 sigma curves. The odds ratio file includes information on the total odds ratio, the corresponding probability of a variable signal, the mth bin with the maximum odds ratio and the odds-weighted first moment of the mth bin, as well as the characteristic time scales for these two values. The light curve file consists of the binnings weighted by the odd ratios and shows the optimal binning for the curve. The standard deviation is provided for each point on the light curve. The range of probabilities: 0.5 < P < 0.9 (above 0.9 all is variable and below 0.5 all is non-variable) is found to be ambiguous and therefore additional criterion is required based on the light curve, its average standard deviation, and the average count rate. The fractions f3 and f5 of the light curve that are within 3 sigma and 5 sigma respectively, of the average count rate. If f3 > 0.997 AND f5 = 1.0 for cases in the ambiguous interval, the source is considered to be non-variable. The variability index is calculated: Table 1.1
Exampleglvary infile=acis_evt.fits outfile=gl_out.fits lc=gl_lc.fits eff="dither.fits[cols time,dtf=psffrac]" Search the events in 'acis_evt.fits' (TIME column) and output both a probability-weighted light curve and a table of odds. The PSFFRAC in the 'dither.fits' file is used to correct the time bins for intervals when the aperture moved on and off the chip. Parameters
Detailed Parameter DescriptionsParameter=infile (file required filetype=input)Input event file The event times and GTI are used to search for time varability. If there is more than one GTI for ACIS data, the first GTI listed is used. Parameter=outfile (file required filetype=output)Output probabilities A probability is output for each division of the time interval. Parameter=lcfile (file required filetype=output)Probability-weighted light curve Parameter=effile (file required filetype=input)Efficiency file The file must have TIME and DTF columns, e.g. an HRC dead time factor (DTF) file. Parameter=frac (real not required default=1.0)Fraction of events to use (1.0 == 100 percent) Parameter=seed (integer not required default=1)Random seed used to discard events if "frac" value is less than 1.0 Parameter=mmax (integer not required default=INDEF)Maximum number of bins to split time range into. INDEF == tool will determine when to stop. Parameter=mmin (integer not required default=INDEF)Minimum number of bins to split time range into. INDEF == tool will determine when to stop. Parameter=nbin (integer not required default=0)Number of bins in the output lightcurve; 0 == tool will determine the optimal value. Parameter=mintime (float not required default=50)Smallest time bin to allow. Parameter=clobber (boolean not required filetype=input default=no)Overwrite existing output dataset with same name? BugsSee the bugs page for this tool on the CIAO website for an up-to-date listing of known bugs. See Also
|
 |
The Chandra X-Ray
Center (CXC) is operated for NASA by the Smithsonian Astrophysical Observatory. 60 Garden Street, Cambridge, MA 02138 USA. Email: cxcweb@head.cfa.harvard.edu Smithsonian Institution, Copyright © 1998-2004. All rights reserved. |
![[CIAO Logo]](../imgs/ciao_logo_navbar.gif){kind=logo}