[daip] Aips task SAD fails to work when In2name is used for rms
Subhashis Roy
roy at ncra.tifr.res.in
Tue Feb 7 15:00:34 EST 2012
Hi,
As shown below, I find that SAD (both in 31DEC11 and 31DEC10 version)
finds all the fit to be rejected due to minimum flux while using SNR in
CPARM and DPARM with DPARM(9) >2 and valid In2name for rms noise.
-----------------------------------------
roy.nc> SAD 1: Task SAD (release of 31DEC11) begins
roy.nc> SAD 1: Create G358.7+0.DBC.MODRES. 2 (MA) on disk 3 cno 390
roy.nc> SAD 1: Fit to histogram gives RMS = 2.026E-03
roy.nc> SAD 1: Looking for up to 500 sources in subfield:
roy.nc> SAD 1: BLC: 1. 1. 1. 1. 1. 1. 1.
roy.nc> SAD 1: TRC: 512. 512. 1. 1. 1. 1. 1.
roy.nc> SAD 1: Using Clean beam 14.35781 10.72465 16.51 from CG table
roy.nc> SAD 1: Searching down to 50.000000 (Sig/Nois) pass 1
roy.nc> SAD 1: Found 1 Islands in pass 1
roy.nc> SAD 1: Debug display of rejection for, in order,
roy.nc> SAD 1: peak, flux, rms, width, Xpos, Ypos, Xout, Yout, residual flux
roy.nc> SAD 1: Options selected are '1 '
roy.nc> SAD 1: Flux RMS Width BLC TRC 'Reasons '
roy.nc> SAD 1: Searching down to 15.000000 (Sig/Nois) pass 2
roy.nc> SAD 1: Found 1 Islands in pass 2
roy.nc> SAD 1: Searching down to 4.500000 (Sig/Nois) pass 3
roy.nc> SAD 1: Found 1 Islands in pass 3
roy.nc> SAD 1: Successfully solved 3 islands, all converged
roy.nc> SAD 1: 3 total components had their fits rejected
roy.nc> SAD 1: 3 components rejected for minimum flux
roy.nc> SAD 1: Resulting number of sources fit is 0
roy.nc> SAD 1: Destroyed residual image as requested
roy.nc> SAD 1: Appears to have ended successfully
roy.nc> SAD 1: roy 31DEC11 TST: Cpu= 0.0 Real= 0 IO= 10
---------------------------------------
As an input file, I generated a simple Gaussian model with width of 4
pixels and a peak flux density of 1 Jy at the centre with a noise (flux)
of 2 mJy/beam using IMMOD. This model was used as inname in SAD.
(However, SAD does work when dparm(9)=0 and actual flux density cutoffs
are used in cparm and dparm.)
For your verification I am attaching the input Fits file (gzipped), and
the noise map (generated from the above using RMSD).
The SAD input is also shown below.
----------------------------------
AIPS 1: SAD: Task to fit Gaussian models to an image by least-squares
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME 'G358.7+0.DBC' Input image name (name)
AIPS 1: INCLASS 'IMMOD1' Input image name (class)
AIPS 1: INSEQ 1 Input image name (seq. #)
AIPS 1: INDISK 3 Input image disk drive #
AIPS 1: INVERS -1 Input image MF ext. version
AIPS 1: <= 0 => new
AIPS 1: BLC *all 0 Bottom Left corner of fit
AIPS 1: TRC *all 0 Top Right corner of fit
AIPS 1: DORESID 0 >0 => 1 Catalog residual map
AIPS 1: OUTNAME 'G358.7+0.DBC' Output image name (name)
AIPS 1: OUTCLASS 'MODRES' Output image name (class)
AIPS 1: OUTSEQ 0 Output image name (seq. #)
AIPS 1: OUTDISK 3 Output image disk drive #
AIPS 1: NGAUSS 500 Max. Number of components
AIPS 1: CPARM 50 15 Flux cutoff levels in
AIPS 1: 4.5 *rest 0 descending order
AIPS 1: ICUT 0 Retry level
AIPS 1: BWSMEAR 0.003 Bandwidth smearing corr.
AIPS 1: SORT ' ' Sort Order of output ' '=RA
AIPS 1: DOCRT 1 >0 => print fits on CRT
AIPS 1: > 72 => terminal width
AIPS 1: =0 => print suppressed
AIPS 1: <0 => line printer
AIPS 1: FITOUT ' '
AIPS 1: Disk file to save fit info
AIPS 1: OUTVERS -1 CC file version # for
AIPS 1: solutions: -1 => none, 0 new
AIPS 1: STVERS -1 STars output file version
AIPS 1: -1 => none, 0 => new
AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
#
AIPS 1: DOALL 1 >0 => fit multiple peaks
AIPS 1: DOWIDTH 1 1 >0 => fit widths
AIPS 1: 1 1 1 1
AIPS 1: 1 1 1 1
AIPS 1: 1 1
AIPS 1: GAIN 0.1 Amp-dependent part of retry
AIPS 1: and warning levels
AIPS 1: DPARM 4 5 Reject components with
AIPS 1: 0 8 (1) peak < D(1)
AIPS 1: 0 0 (2) flux < D(2)
AIPS 1: 5 0 (3) rms > D(3) + GAIN * flux
AIPS 1: 5 0 (in quadrature)
AIPS 1: (4) widths > D(4) cells
AIPS 1: (5) peaks > D(5) cells
AIPS 1: outside interior area
AIPS 1: (6) peaks > D(6) cells
AIPS 1: outside image
AIPS 1: (7) total residual flux in
AIPS 1: the fitting box
AIPS 1: [no default]
AIPS 1: (9) Use IN2NAME etc to set
AIPS 1: image of rms
AIPS 1: >= 2 => CPARM is in S/N
AIPS 1: (see HELP)
AIPS 1: (10) Expand fit box pixels
AIPS 1: IN2NAME 'G358.7+0.DBC' RMS image name (name)
AIPS 1: IN2CLASS 'TSTRMS' RMS image name (class)
AIPS 1: IN2SEQ 1 RMS image name (seq. #)
AIPS 1: IN2DISK 3 RMS image disk drive #
AIPS 1: PRTLEV 1 Debug print level
AIPS 1: PBPARM *all 0 Primary beam parameters:
AIPS 1: (1) level to believe - <= 0
AIPS 1: means do not apply a primary
AIPS 1: beam (2) > 0 use (3)-(7)
----------------
With regards,
Subhashis
-------------- next part --------------
A non-text attachment was scrubbed...
Name: G358.7+0.DBC.IMMOD1.FIT.gz
Type: application/x-gzip
Size: 296560 bytes
Desc:
URL: <http://listmgr.nrao.edu/pipermail/daip/attachments/20120208/34f06650/attachment.gz>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: G358.7+0.DBC.TSTRMS.FIT.gz
Type: application/x-gzip
Size: 212651 bytes
Desc:
URL: <http://listmgr.nrao.edu/pipermail/daip/attachments/20120208/34f06650/attachment-0001.gz>
More information about the Daip
mailing list