Howdy,
Problems in paradise: I am trying to generate a new level=1.5 event file
for my HETG/ACIS data by following the CIAO science thread http://asc.har-
vard.edu/ciao/threads/level1.5.hetg.thread.html. When I try to get the
position of the zero order image with tgdetect, it fails with this cryptic
comment:
/data/patrickw/CXCDS_ROOT_SOLARIS/bin/tgdetect: test: argument expected
Thinking that this might be a problem with the new level 1 event file I
had produced, I tried and failed to run tgdetect on the original level 1
event file. The "diagnostic" information returned to the screen follows.
If it matters, I am running CIAO 1.1.5. I didn't run into this problem
last time I tried this (when I was using CIAO 1.1.3). I don't see a bug
report for this tool. What's the deal?
ich bein das chrome dinnette,
Chris Mauche
--------------------------------------------------------------------------
| Christopher W. Mauche Phone: (925) 422-7017 |
| Lawrence Livermore Nat'l Lab. Fax: (925) 423-5998 |
| L-43, 7000 East Avenue e-mail: mauche@cygnus.llnl.gov |
| Livermore, CA 94550 mauche@llnl.gov |
|------------------------------------------------------------------------|
| www: http://www-phys.llnl.gov/Staff/MaucheChris/ |
--------------------------------------------------------------------------
> punlearn tgdetect
> pset tgdetect infile=acisf01706_001N001_evt1.fits [original evt1 file]
> pset tgdetect outfile=GC_evt1_aspflt_src1a.fits
> pset tgdetect eenergy=0.8 [as in example, fails either way]
> pset tgdetect verbose=2
> tgdetect
Input L1 event file or @list (root_evt1.fits or .lis) (acisf01706_001N001_evt1.fits):
Input source position(s) file from previous OBI or @list or NONE (NONE):
Output source position(s) file name (root_src1a.fits) (GC_evt1_aspflt_src1a.fits):
RUNNING tool: tgdetect
#################### PARAMETERS SECTION ####################
Input parameters are ...
infile=acisf01706_001N001_evt1.fits
OBI_srclist_file=NONE
outfile=GC_evt1_aspflt_src1a.fits
zo_pos_x=default
zo_pos_y=default
zo_sz_filt_x=default
zo_sz_filt_y=default
snr_thresh=40
thresh=3
findpeaks=yes
celldetect_log=no
psftable=/data/patrick1/CXCDS_ROOT_SOLARIS/data/psfsize.fits
fixedcell=0
fixedcell_cc_mode=99
expfile=none
bkgfile=none
bkgvalue=0
bkgerrvalue=0
eband=1.4967
eenergy=0.8
snrfile=none
convolve=no
xoffset=INDEF
yoffset=INDEF
cellfile=none
centroid=yes
kernel=fits
snr_ratio_limit=1
setsrcid=yes
max_separation=3
clobber=yes
verbose=2
mode=ql
Root name is GC_evt1_aspflt
#################### TG_FILT SECTION ####################
RUNNING dmmakepar...
dmmakepar input=acisf01706_001N001_evt1.fits[EVENTS] output=/tmp/tgdetect_4413.par
Wrote hdr file from acisf01706_001N001_evt1.fits -> /tmp/tgdetect_4413.par
Reading hdr parameters: ra_nom dec_nom instrume detnam grating
Values of hdr parameters ...
ra_nom = 193.099257
dec_nom = -29.253873
instrume = ACIS
detnam = ACIS-456789
grating = HETG
acsys = SKY:ASC-FP-1.1
readmode = TIMED
CENTER of the sky pixel plane is: 4096.5, 4096.5
SETTING the default gzo filter x = 400
SETTING the default gzo filter y = 400
Creating image filter acisf01706_001N001_evt1.fits[EVENTS][bin x=3896:4296,y=3896:4296] for celldetect input
#################### CELLDETECT SECTION ####################
RUNNING celldetect...
celldetect infile=acisf01706_001N001_evt1.fits[EVENTS][bin x=3896:4296,y=3896:4296] outfile=GC_evt1_aspflt_cell_out.fits log=no verbose=2 clobber=yes psftable=/data/patrick1/CXCDS_ROOT_SOLARIS/data/psfsize.fits fixedcell=0 expfile=none bkgfile=none bkgvalue=0 bkgerrvalue=0 eband=1.4967 eenergy=0.8 thresh=3 findpeaks=yes snrfile=none convolve=no xoffset=INDEF yoffset=INDEF cellfile=none centroid=yes kernel=fits
input file: acisf01706_001N001_evt1.fits[EVENTS][bin x=3896:4296,y=3896:4296]
output file: GC_evt1_aspflt_cell_out.fits
Warning: Keyword TLMVER is not in the header.
Warning: Keyword HDUSPEC is not in the header.
Warning: Keyword HDUDOC is not in the header.
Warning: Keyword HDUVERS is not in the header.
creating cell size table
reading data image
creating candidate source list
cell size 6
Warning: Keyword HISTNUM is not in the header.
OUTPUT file from celldetect is: GC_evt1_aspflt_cell_out.fits
#################### SRC SELECT SNR SECTION ####################
RUNNING dmcopy w/ snr threshold filter...
Using snr_thresh = 40
dmcopy infile=GC_evt1_aspflt_cell_out.fits[SRCLIST][snr=40:] outfile=GC_evt1_aspflt_snr_select.fits kernel=fits verbose=2 clobber=yes
The debug level is 2
Open GC_evt1_aspflt_cell_out.fits[SRCLIST][snr=40:] as table or image
CXC Data Model: DM V1.98
DM parse clause snr=40:
DM: Opened Table block: GC_evt1_aspflt_cell_out.fits[SRCLIST][snr=40:]
Output file is GC_evt1_aspflt_snr_select.fits
Clobbering old output file: GC_evt1_aspflt_snr_select.fits
DM: Opened Image block: PRIMARY
DM: Opened Image block: PRIMARY
DM: Opened Table block: SRCLIST
HIST = dmcopy GC_evt1_aspflt_cell_out.fits[SRCLIST][snr=40:] GC_evt1_aspflt_snr_select.fits kernel=FITS
OUTPUT file from dmcopy snr select is: GC_evt1_aspflt_snr_select.fits
/data/patrick1/CXCDS_ROOT_SOLARIS/bin/tgdetect: test: argument expected
Check tgdetect parameters:
> plist tgdetect
Parameters for /home/mauche/cxcds_param/tgdetect.par
##
## TGDETECT -- Create filter; run celldetect; narrow down detected
## 'zero order' source list; set source id's; match
## sources to previous OBI source list.
##
## Note: if either "infile" or "OBI_srclist_file" are @lists, only
## the first item on the list is read in; this tool only works on
## one set of input files; if more than one file is listed,
## everything but the first are ignored.
##
infile = acisf01706_001N001_evt1.fits Input L1 event file or @list (root_evt1.fits or .lis)
OBI_srclist_file = NONE Input source position(s) file from previous OBI or @list or NONE
outfile = GC_evt1_aspflt_src1a.fits Output source position(s) file name (root_src1a.fits)
(zo_pos_x = default) Center GZO filter sky X position (default=pixel(ra_nom))
(zo_pos_y = default) Center GZO filter sky Y position (default=pixel(dec_nom))
(zo_sz_filt_x = default) Size of GZO filter in X pixels (ACIS=200; HRC=1800)
(zo_sz_filt_y = default) Size of GZO filter in Y pixels (ACIS=200; HRC=1800)
(snr_thresh = 40) SNR threshold to select the detected sources
#
# celldetect parameters
#
(thresh = )celldetect.thresh -> 3) celldetect source threshold
(findpeaks = yes) find local peaks for celldetect
(celldetect_log = )celldetect.log -> no) make a celldetect log file?
(psftable = )celldetect.psftable -> /data/patrick1/CXCDS_ROOT_SOLARIS/data/psfsize.fits) table of PSF size data, for celldetect
(fixedcell = )celldetect.fixedcell -> 0) celldetect fixed cell size to use
(fixedcell_cc_mode = 99) celldetect fixed cell size to use for CC mode ACIS data
(expfile = none) exposure map file name, for celldetect
(bkgfile = none) background file, for celldetect
(bkgvalue = )celldetect.bkgvalue -> 0) background count/pixel, for celldetect
(bkgerrvalue = )celldetect.bkgerrvalue -> 0) background error, for celldetect
(eband = )celldetect.eband -> 1.4967) energy band, for celldetect
(eenergy = 0.8) encircled energy of PSF, for celldetect
(snrfile = none) celldetect snr output file (for convolution only)
(convolve = )celldetect.convolve -> no) use convolutions for celldetect
(xoffset = INDEF) celldetect offset of x axis from optical axis
(yoffset = INDEF) celldetect offset of y axis from optical axis
(cellfile = none) output cell size image file
(centroid = yes) compute source centroids in celldetection?
(kernel = fits) output file format for celldetect
#
# tgidselectsrc parameters
#
(snr_ratio_limit = )tgidselectsrc.snr_ratio_limit -> 1) Value of SNR ratio to use as lower limit
(setsrcid = )tgidselectsrc.setsrcid -> yes) Set src ids in output file?
#
# tgmatchsrc parameters
#
(max_separation = )tgmatchsrc.max_separation -> 3) Maximum allowed separation (arcsec) for sources to match
#
#
(clobber = yes) OK to overwrite existing output file(s)?
(verbose = 2) Verbosity level (0 = no display)
(mode = ql)
This archive was generated by hypermail 2b29 : Thu Dec 13 2012 - 01:00:03 EST