Analysis Utilities

Analysis Utilities (or analysisUtils for short) is a small set of Python source code files that provide a number of analysis and plotting utilities. The utilities were developed for ALMA data reduction and are, in many cases, also useful for EVLA data reduction. This CASA Guide documents some of the most useful functions contained in the analysisUtils Python module from a user's perspective.

Getting Started

If you are working on a machine at NRAO Charlottesville, Socorro, or Santiago, then you can skip Step 1, because the latest modules are already available at: /users/thunter/AIV/science/analysis_scripts.

Step 1: Download Analysis Utilities from here and extract the tar ball. From a Unix command line this can be done with

$ tar xvf analysis_scripts.tar

Seven files will be extracted -- a README file, and six python scripts. The README file contains the time and date that the tar ball was generated, which can be useful for reporting bugs. The history file available at the ftp site contains a list of major changes in each version.

Step 2: Edit your existing casapy initialization file in $HOME/.casa/init.py or create a new empty file if it does not already exist, and add the following

import sys
sys.path.append("/PATH_TO_ANALYSIS_SCRIPTS/")
import analysisUtils as au

where /PATH_TO_ANALYSIS_SCRIPTS/ is the path to the directory you just extracted from the tar ball. When this is done, start casapy and you will have access to all the functions contained in the analysisUtils module.

Support

For help with analysisUtils installation or usage, you may submit a helpdesk ticket to the Data Reduction category at https://help.almascience.org/. Be advised that analysisUtils is a commissioning tool, so the level of support available is not as extensive as for CASA, but we do welcome feedback. If your question is with regard to plotbandpass, then please provide the version number of this program (it is printed to the screen at execution time and on the plots). In all cases, please provide the release date of your analysisUtils, which can be viewed from inside casapy as follows:

au.version()

Key Function

plotbandpass

a faster version of plotcal for bandpass tables, with useful overlay capabilities. This function has been incorporated into CASA as a task.

Some Other Functions

Up-to-date, online help is available for all of these functions (and more) using the command: help(au.functionName). To search for functions, you can type au.help('phrase'), where phrase is case-insensitive and minimum-match.

angularSeparationOfFields

reads an ms and computes the angular separation of all fields, or a subset

buildConfigurationFile

reads the list of antenna stations in an .ms and creates a .cfg file suitable for simobserve

editIntents

change the intents for a specified field in an ms (based on John Lightfoot's pipeline script)

getBaselineLengths

returns a dictionary of the baseline lengths in your ms, by default sorted by length

getBaselineStats

returns a bunch of statistics on the baseline lengths in your ms

gjincBeam

compute the effective restoring beam obtained from the casa command sdimaging when using the GJINC gridding kernel

lstrange

lists the range of LST, UT, MJD for the whole ms, and for scans with OBSERVE_TARGET intent (including the elevation range)

obslist

prints the antenna station coordinates in local offsets from the Center of Array, and computes longest/shortest baselines

planet

contacts the JPL Horizons telnet service and returns the J2000 position, velocity, angular diameter, range and range rate of a solar system object

planetFlux

calls Bryan Butler's solar_system_setjy module in CASA to compute flux density vs. time or flux density vs. frequency

planetPlots

generates a plots of uv amplitude vs. uv distance for a grid of ALMA configurations and observing frequencies

plotSunDuringTrack

generates a plot of the Sun's az vs. el during an ms or ASDM

plotconfig

plot any standard telescope configuration of observatories known to casa, and return an array of the sorted baseline lengths

plotmosaic

shows relative location of pointings in an .ms

plotPWV

reads the PWV from each antenna's WVR from the ASDM_CALWVR table in an ms, and creates a plot vs. time

plotweather

plot weather conditions vs. time for your ms

plotWVRSolutions

plot the WVR solutions in terms of baseline-based phase corrections

spectralindex

performs a least-squares fit to the multi-spw output from fluxscale

timeOnSource

compute the total time spent integrating on-source for each specified field

Documentation of additional functions can be found at Todd's wikipage