plotbandpass – Makes detailed plots of Tsys and bandpass solutions. – visualization, calibration task

Description

Developed at the NAASC, this is a generic task to display CASA: Tsys and bandpass solution tables with options to overlay them in various combinations, and/or with an atmospheric transmission or sky temperature model. It works with both the ‘new’ (casa 3.4) and ‘old’ calibration table formats, and allows for mixed mode spws (e.g. TDM and FDM for ALMA). It uses the new msmd tool to access the information about an ms. This task is still being developed as new ALMA observing modes are commissioned. So if you encounter problems, please report them.

Parameters

Title
Parameter	Default	Description
caltable	`''`
antenna	`''`
field	`''`
spw	`''`
yaxis	`'amp'`
xaxis	`'chan'`
figfile	`''`
plotrange	`numpy.array( [ ] )`
caltable2	`''`
overlay	`''`
showflagged	`False`
timeranges	`''`
buildpdf	`False`
caltable3	`''`
markersize	`int(3)`
density	`int(108)`
interactive	`True`
showpoints	`'auto'`
showlines	`'auto'`
subplot	`'22'`
zoom	`''`
poln	`''`
showatm	`False`
pwv	`'auto'`
gs	`'gs'`
convert	`'convert'`
chanrange	`''`
solutionTimeThresholdSeconds	`float(30.0)`
debug	`False`
phase	`''`
vis	`''`
showtsky	`False`
showfdm	`False`
showatmfield	`''`
lo1	`''`
showimage	`False`
showatmpoints	`False`
parentms	`''`
pdftk	`'pdftk'`
channeldiff	`False`
edge	`int(8)`
resample	`int(1)`
platformingThreshold	`float(10.0)`
platformingSigma	`float(10.0)`
basebands	`''`
showBasebandNumber	`False`
scans	`''`
figfileSequential	`False`
chanrangeSetXrange	`False`

Parameter Explanations

caltable

''

Input table name, either a bandpass solution or a Tsys solution

antenna

''

A comma-delimited string list of antennas (either names or integer indices) for which to display solutions. Default = all antennas.

field

''

A comma-delimited string list of fields (either names or integer indices) for which to display solutions. Default = all fields.

spw

''

A comma-delimited string list of spws for which to display solutions. Default = all spws.

yaxis

'amp'

The quantity to plot on the y-axis (“amp”, “phase”, “both”, “tsys”, append “db” for dB).

xaxis

'chan'

The quantity to plot on the x-axis (“chan” or “freq”).

figfile

''

The name of the plot file to produce.

plotrange

numpy.array( [ ] )

The axes limits to use [x0,x1,y0,y1].

caltable2

''

A second cal table, of type BPOLY or B, to overlay on a B table

overlay

''

Show multiple solutions in same frame in different colors (time, antenna, spw, baseband, or time,antenna)

showflagged

False

Show the values of the solution, even if flagged

timeranges

''

Show only these timeranges, the first timerange being 0

buildpdf

False

If True, assemble all the pngs into a pdf

caltable3

''

A third cal table, of type BPOLY, to overlay on the first two tables

markersize

int(3)

Size of points

density

int(108)

dpi to use in creating PNGs and PDFs (default=108)

interactive

True

if False, then run to completion automatically without pause

showpoints

'auto'

Draw points for the data (default=F for amp, T for phase)

showlines

'auto'

Draw lines connecting the data (default=T for amp, F for phase)

subplot

'22'

11..81,22,32 or 42 for RowsxColumns (default=22), any 3rd digit is ignored

zoom

''

“intersect” will zoom to overlap region of caltable with caltable2

poln

''

Polarizations to plot: “” = all, or “RR”,”RL”,”LR”,”LL”,”XX”,”XY”,”YX”,”YY”,”RR,LL”,”XX,YY”

showatm

False

Compute and overlay the atmospheric transmission curve

pwv

'auto'

Define the pwv to use for the showatm option: “auto” or value in mm

gs

'gs'

For buildpdf=T, full path for ghostscript command (in case it is not found)

convert

'convert'

For buildpdf=T, full path for the ImageMagick convert command (in case it is not found)

chanrange

''

Set xrange (“5~100”) over which to autoscale y-axis for xaxis=”freq”

solutionTimeThresholdSeconds

float(30.0)

Consider 2 solutions simultaneous if within this interval in seconds

debug

False

Print verbose messages for debugging purposes

phase

''

The y-axis limits to use for phase plots when yaxis=”both”

vis

''

name of the ms for this table, in case it does not match the string in the caltable

showtsky

False

Compute and overlay the sky temperature curve instead of transmission

showfdm

False

when showing TDM spws, draw the locations of the corresponding FDM spws

showatmfield

''

for overlay=”time”, use first observation of this fieldID or name

lo1

''

specify the LO1 setting (in GHz) for the observation (’’ = automatic)

showimage

False

also show the atmospheric curve for the image sideband (in black)

showatmpoints

False

Draw atmospheric curve with points instead of a line

parentms

''

if showimage=T, name of the parent ms (only needed if the ms has been previously split)

pdftk

'pdftk'

For buildpdf=T, full path for pdftk command (in case it is not found)

channeldiff

False

Set to a value > 0 (sigma) to plot derivatives of the solutions

edge

int(8)

The number of edge channels to ignore in finding outliers (for channeldiff>0)

resample

int(1)

The channel expansion factor to use when computing MAD of derivative (for channeldiff>0)

platformingThreshold

float(10.0)

if platformingSigma=0, then declare platforming if the amplitude derivative exceeds this percentage of the median

platformingSigma

float(10.0)

declare platforming if the amplitude derivative exceeds this many times the MAD

basebands

''

A baseband number or list of baseband numbers for which to display solutions. Default = all.

showBasebandNumber

False

Put the baseband converter number (BBC_NO) in the title of each plot

scans

''

A scan or list of scans for which to display solutions. Default = all. Does not work with overlay=”time”.

figfileSequential

False

naming scheme for pngs: False: name by spw/antenna (default), True: figfile.000.png, figfile.001.png, etc.

chanrangeSetXrange

False

If True, then chanrange also sets the xrange to display