ldetect -- determine source positions from an input qpoe
ldetect qpoe tab cell_size eband
LDETECT runs the detect tasks (CELLMAP, SNRMAP, LPEAKS, BKDEN, and BEPOS) for a single detect cell size and threshold. Intermediate files are saved. An overview of the detect algorithms and implementation is available in the EXPLAIN_DETECT on-line documentation.
Name of the input qpoe file.
Required output table root. The is specified by supplying the rootname and the intermediate files will use the root in their names. See helps for lpeaks and bepos for for a list of what's included in their tables.
6 | 9 | 12 | 24 | 36 | 48Cellsizes for ROSAT HRI qpoe.
30 [h/b]| 45 [s]| 60 [s/h/b]| 120 [s/h/b]Cellsizes for ROSAT PSPC qpoe.
6 | 9 | 12 | 24 | 48Cellsizes for Einstein HRI qpoe.
144 [h/b]| 240 [s]Cellsizes for Einstein IPC qpoe.
Energy band. Prompted when instrument is PSPC or IPC. Choice of Soft, hard, or broad. The PI filter appropriate to the band choice is appended to the qpoe filename. PI bands must be respecified when running BEPOS.
Required Region from which the background density is computed. If "default" is specified, a default region is defined. The default region is (circle img_cenx img_ceny radius), where radius is the specified bkden parameter below specified in arc-minutes.
circle radius for default background density region
Acceptable Fluctuation, F above.
Maximum iterations for computing bkden.
Output parameters. It is the computed background density, from the last run of the task.
Minimum SNR threshold for local detect (lpeaks). Sources with computed SNR's below this threshold are thrown out.
Minimum SNR threshold for Maximum likelihood determination (bepos). Sources with computed SNR's below this threshold are thrown out.
PRF Sigma of the field in arcseconds. When 0.0, the PRF sigma is computed by keying on Telescope, instrument (from input header) & ldetect.energy in a PRF lookup table (ldetect.prf_table) and using the coefficients that are returned to compute the PRF.
Energy (KeV), key look up in PRF coefficients table, coeffs are used to compute the Sigma for a Gaussian Point Response Function.
PRF coefficients table. Entries for ROSAT PSPC/HRI & Einstein IPC/HRI. Lookup based on value of bepos.energy. Returns coefficients assiciated with a default lookup when energy = 0.0, returns coefficients closest to the energy when exact match is not found.
Maximum iterations for convergence to best positions values.
Maximum iterations for confidence error calculations.
Number of confidence levels that confidence errors are determined for below.
The 68 percent confidence level for minimum c-statistic calculation.
The 90 percent confidence level for minimum c-statistic calculation.
Pixel limit between LPEAKS and BEPOS position. When the difference between the 2 positions is greater that this limit the source is thrown out.
Factor in calculating the radius of the bin over which we will search for the position. (bin_radius = optimum_radius * src_sigma)
Convergence minimum for c-statistic
Convergence minimum for source
X center header keyword for off-axis angle calculation
Y center header keyword for off-axis angle calculation
Boolean flag specifying whether or not the output Table file can be overwritten, if it already exists.
Display info to screen. Integer of range 0-5. The higher the number more verbose the display.
xp> ldetect qpoe file : xdata$rp110590.qp Output table root: parlac45 ROSAT PSPC Cellsize in arcsecs - 30 [h/b]| 45 [s]| 60 [s/h/b]| 120 [s/h/b]: 45 Energy Band (soft|hard|broad) (broad): soft ... Running Cellmap ... xdata$rp110590.qp[bl=30,pi=11:41] -> parlac45.imh add parlac45.imh,cellsize = 45 parlac45.imh updated xdata$rp110590.qp[pi=11:41] ... Running Snrmap ... Convolving 3x3 cells Convolving 5x5 cells Calculating Snrmap with equation: snrmap=(1.7777*cell3-frame)/((1.7777**2*cell3+oframe)**0.5) Writing Snrmap Image: parlac45_snr.imh ... Running Lpeaks ... src# logx logy physx physy cellx/y snr --------------------------------------------------------- 1 233.00 260.00 6990.00 7800.00 45 45 3.527 2 232.00 312.00 6960.00 9360.00 45 45 17.798 Writing to Output file: parlac45_ruf.tab ... Running Bkden ... Bkden region is circle 256.00 256.00 7.50' Bkden = 1.1024632 counts/arcmin**2 ... Running Bepos ... *** PI bands must be Respecified with input QPOE *** Input Qpoe: xdata$rp110590.qp[pi=11:41] bkden = 1.1024632 (cts/sq arc-min), 0.0000766 (cts/sq pixel) 2 Rough Position(s) Input src# ra dec physx physy cellcnts framecnts snr --------------------------------------------------------------------------- 1 22:08:38.3 45:31:50.5 6996.9 7802.5 32 18 3.4 2 22:08:41.4 45:44:49.5 6934.5 9360.6 2046 2086 16.8 Writing 2 Sources to Output file: parlac45_pos.tab
Documentation giving an overview of the detect package and run order of tasks. (help using_detect )
Documentation on the detect algorithms and philosophy along with limitations in (help explain_detect ).
Documentation Match Source task. (help lmatchsrc )