[daip] Aips task SAD fails to work when In2name is used for rms

[Eric Greisen]

Subhashis Roy wrote:
> Hi Eric,
> 
> 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
> ---------------------------------------

> 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

You have asked for all components < 5 Jy to be discarded and the task 
has done that for you.  DPARM(9) >- 2 makes some things (those with 
units of Jy/beam) change to units of S/N ratio - CPARM, ICUT, DPARM(1), 
DPARM(3).  The Units of DPARM(2) are Jy and so cannot be switched.

Eric Greisen