Spectralindex

Return to Analysis Utilities

This page documents the spectralindex function of Python module analysisUtils.

This function is designed to fit a power-law spectral index to the results output from casa's fluxscale task. Currently, it also requires the output from the listobs task to determine the center frequencies of each spectral window. It runs a brief Monte-Carlo series of fits to determine the uncertainty on the fitted slope on the basis of the error bars on each flux density. Finally, it produces a plot. Note: spectralIndex is a synonym for spectralindex.

Usage:

spectralIndex(filename, yfilename, source=, verbose=False, maxpoints=0, trials=2000, spw=, help=False)

   filename: contains a listobs output file
   yfilename: contains a fluxscale output file
   source: sourcename to choose from the (possibly) multi-source fluxscale file
   maxpoints: the maximum number of spws to select for the fit (0=no max)
   trials: number of Monte-Carlo fits to run to estimate the fit uncertainties
   spw: the spws to use, e.g. =all, '1~3,5,6~8'=[1,2,3,5,6,7,8]