sdfit – Fit a spectral line – single dish task

Description

Task sdfit is a basic line-fitter for single-dish spectra. It assumes that the spectra have been calibrated in tsdcal or sdreduce.

Parameters

Title
Parameter	Default	Description
infile	`''`	name of input SD dataset
datacolumn	`'data'`	name of data column to be used [“data”, “float_data”, or “corrected_data”]
antenna	`''`	select data by antenna name or ID, e.g. “PM03”
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)
timerange	`''`	select data by time range, e.g. “09:14:0~09:54:0” (“”=all) (see examples in help)
scan	`''`	select data by scan numbers, e.g. “21~23” (“”=all)
pol	`''`	select data by polarization IDs, e.g. “XX,YY” (“”=all)
intent	`''`	select data by observational intent, e.g. “ON_SOURCE” (“”=all)
timebin	`''`	bin width for time averaging
timespan	`''`	span the timebin across “scan”, “state”, “field”, or a combination of them (e.g., “scan,state”)
polaverage	`''`	polarization averaging mode (“”, “stokes” or “geometric”).
fitfunc	`'gaussian'`	function for fitting [“gaussian”, “lorentzian”]
fitmode	`'list'`	mode for setting additional channel masks. “list” and “auto” are available now.
nfit	`numpy.array( [ int(0) ] )`	list of number of lines to fit in maskline region.
thresh	`float(5.0)`	S/N threshold for linefinder
avg_limit	`int(4)`	channel averaging for broad lines
minwidth	`int(4)`	the minimum channel width to detect as a line
edge	`numpy.array( [ int(0),int(0) ] )`	channels to drop at beginning and end of spectrum
outfile	`''`	name of output file
overwrite	`False`	overwrite the output file if already exists [True, False]

Parameter Explanations

infile

''

name of input SD dataset

datacolumn

'data'

name of data column to be used [“data”, “float_data”, or “corrected_data”]

antenna

''

select data by antenna name or ID, e.g. “PM03”

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)

timerange

''

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

scan

''

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

pol

''

select data by polarization IDs, e.g. “XX,YY” (“”=all)

intent

''

select data by observational intent, e.g. “ON_SOURCE” (“”=all)

timebin

''

bin width for time averaging

timespan

''

span the timebin across “scan”, “state”, “field”, or a combination of them (e.g., “scan,state”)

polaverage

''

polarization averaging mode (“”, “stokes” or “geometric”).

fitfunc

'gaussian'

function for fitting

fitmode

'list'

mode for setting additional channel masks.

nfit

numpy.array( [ int(0) ] )

list of number of lines to fit in maskline region.

thresh

float(5.0)

S/N threshold for linefinder

avg_limit

int(4)

channel averaging for broad lines

minwidth

int(4)

the minimum channel width to detect as a line

edge

numpy.array( [ int(0),int(0) ] )

channels to drop at beginning and end of spectrum

outfile

''

name of output file

overwrite

False

overwrite the output file if already exists