sdreduceold – ASAP SD task [DEPRECATED]: do sdcal, sdaverage, and sdbaseline in one task – single dish task

Description

### DEPRECATION WARNING ################################################# This task will be removed in CASA 5.1. The functionality of this task with MeasurementSet format is replicated with sdcal, sdbaseline, sdsmooth, and mstransform. #########################################################################

Task sdreduceold performs data selection, calibration, spectral averaging and/or baseline fitting for single-dish spectra. This task internally calls the tasks, sdcal, sdaverage, and sdbaseline and it can be used to run all the three steps in one task execution. This task has better performance than invoking the three tasks separately because it runs all three steps without writing intermediate data to disk.

It is possible to skip arbitrary operations by setting calmode = ‘none’ (for calibration), average=False (for time and polarization averaging), kernel = ‘none’ (for smoothing), and/or blfunc=’none’ (for baseline fitting).

Please take a look at descriptions of tasks, sdcal, sdaverage, and sdbalseline, for more information.

Parameters

Title
Parameter	Default	Description
infile	`''`
antenna	`int(0)`
fluxunit	`''`	units of the flux ['K' or 'Jy'] (''=current)
telescopeparam	`''`
field	`''`
spw	`''`
restfreq	`''`
frame	`''`	frequency reference frame ['LSRK', 'TOPO', 'LSRD', 'BARY', 'GALACTO', 'LGROUP', or 'CMB'] (''=current)
doppler	`''`	doppler convention ['RADIO', 'OPTICAL', 'Z', 'BETA', or 'GAMMA'] (''=current)
timerange	`''`
scan	`''`
pol	`''`
calmode	`'none'`	SD calibration mode ['ps', 'nod', 'otf', 'otfraster', 'fs', 'fsotf' or 'none'] ('none' = skip calibration)
fraction	`'10%'`
noff	`int(-1)`	number of the OFF data to mark, e.g., 10 (-1 = use fraction instead of number)
width	`float(0.5)`
elongated	`False`	the observed area is elongated in one direction [True, False]
markonly	`False`
plotpointings	`False`	plot pointing directions for ON and OFF [True, False]
tau	`float(0.0)`
average	`False`
timeaverage	`False`
tweight	`'tintsys'`	weighting for time averaging ['tintsys', 'tsys', 'tint', 'var', or 'median']
scanaverage	`False`
averageall	`False`
polaverage	`False`
pweight	`'tsys'`	weighting for polarization averaging ['tsys' or 'var']
kernel	`'none'`	type of spectral smoothing kernel ['hanning', 'gaussian', 'boxcar', 'regrid'] (''=no smoothing)
kwidth	`int(5)`
chanwidth	`'5'`
maskmode	`'auto'`	mode of setting additional channel masks ['auto', 'list', or 'interact']
thresh	`float(5.0)`
avg_limit	`int(4)`
edge	`numpy.array( [ int(0) ] )`
blfunc	`'none'`	baseline model function ['poly', 'chebyshev', 'cspline', 'sinusoid' or 'none'] ('none' = skip baseline fit)
order	`int(5)`
npiece	`int(2)`
applyfft	`True`
fftmethod	`'fft'`	method for automatically set wave numbers of sinusoids (only 'fft' is available now)
fftthresh	`float(3.0)`
addwn	`UNKNOWN`
rejwn	`UNKNOWN`
clipthresh	`float(3.0)`
clipniter	`int(0)`
verifycal	`False`
verifysm	`False`
verifybl	`False`	interactively verify the results of operation for each spectrum [True, False] (see description in help)
verbosebl	`True`
bloutput	`True`
blformat	`''`	format of the text file specified with bloutput ['' or 'csv']
showprogress	`True`
minnrow	`int(1000)`
outfile	`''`
outform	`'ASAP'`	output file format ['ASAP', 'MS2', 'ASCII', or 'SDFITS']
overwrite	`False`
plotlevel	`int(0)`

Parameter Explanations

infile

''

name of input SD dataset

antenna

int(0)

select an antenna name or ID, e.g. 'PM03' (only effective for MS input)

fluxunit

''

units of the flux (''=current)

telescopeparam

''

parameters of telescope for flux conversion (see description in help of sdcal)

field

''

select data by field IDs and names, e.g. '3C2*' (''=all)

spw

''

select data by IF IDs (spectral windows), e.g. '3,5,7' (''=all)

restfreq

''

the rest frequency, e.g. '1.41GHz' (default unit: Hz) (see examples in help)

frame

''

frequency reference frame (''=current)

doppler

''

doppler convention (''=current). Effective only when spw selection is in velocity unit.

timerange

''

select data by time range, e.g. '09:14:0~09:54:0' (''=all) (see examples in help of sdcal)

scan

''

select data by scan numbers, e.g. '21~23' (''=all)

pol

''

select data by polarization IDs, e.g. '0,1' (''=all)

calmode

'none'

SD calibration mode ('none' = skip calibration)

fraction

'10%'

fraction of the OFF data to mark as OFF spectra, e.g., '10%'

noff

int(-1)

number of the OFF data to mark (-1 = use fraction instead of number)

width

float(0.5)

width of the pixel for edge detection

elongated

False

the observed area is elongated in one direction

markonly

False

do calibration (False) or just mark OFF (True)

plotpointings

False

plot pointing direction for ON and OFF

tau

float(0.0)

the zenith atmospheric optical depth for correction (0. = no correction)

average

False

data averaging [True, False]

timeaverage

False

average spectra over time [True, False] (see examples in help of sdaverage)

tweight

'tintsys'

weighting for time averaging

scanaverage

False

average spectra within a scan number [True, False] (see examples in help of sdaverage)

averageall

False

set True only when averaging spectra with different spectral resolutions

polaverage

False

average spectra over polarizations [True, False]

pweight

'tsys'

weighting for polarization averaging

kernel

'none'

type of spectral smoothing kernel ('none'=no smoothing)

kwidth

int(5)

width of smoothing kernel in channels

chanwidth

'5'

width of regridded channels

maskmode

'auto'

mode of setting additional channel masks

thresh

float(5.0)

S/N threshold for linefinder

avg_limit

int(4)

channel averaging for broad lines

edge

numpy.array( [ int(0) ] )

channels to drop at beginning and end of spectrum

blfunc

'none'

baseline model function ('none' = skip baseline fit)

order

int(5)

order of baseline model function

npiece

int(2)

number of element polynomials for cubic spline curve

applyfft

True

automatically set wave numbers of sinusoids [True, False]

fftmethod

'fft'

method for automatically set wave numbers of sinusoids

fftthresh

float(3.0)

threshold to select wave numbers of sinusoids

addwn

*UNKNOWN*

additional wave numbers to use

rejwn

*UNKNOWN*

wave numbers NOT to use

clipthresh

float(3.0)

clipping threshold for iterative fitting

clipniter

int(0)

maximum iteration number for iterative fitting

verifycal

False

interactively verify the results of calibration [True, False] (see description in sdcal)

verifysm

False

interactively verify the results of smoothing for each spectrum. [not available for kernel=”regrid”]

verifybl

False

interactively verify the results of baseline fitting for each spectrum (only for blfunc=”poly”. see description in help)

verbosebl

True

output baseline fitting results to logger [True, False]

bloutput

True

output baseline fitting results to a text file [True, False]

blformat

''

format of the text file specified with bloutput

showprogress

True

show progress status for large data [True, False]

minnrow

int(1000)

minimum number of input spectra to show progress status in baseline fitting

outfile

''

name of output file (See a WARNING in help)

outform

'ASAP'

output file format (See a WARNING in help)

overwrite

False

overwrite the output file if already exists [True, False]

plotlevel

int(0)

plot and summarize results (0=none). See description in each task