Plotbandpass: Difference between revisions

From CASA Guides
Jump to navigationJump to search
No edit summary
No edit summary
Line 3: Line 3:
This page documents the plotbandpass function of Python module [[Analysis Utilities|analysisUtils]].
This page documents the plotbandpass function of Python module [[Analysis Utilities|analysisUtils]].


The '''plotbandpass''' function provides extended capabilities beyond those of CASA's plotcal. It was written to facilitate rapid visualization of ALMA's bandpass and spectral Tsys measurements in Science Verification and Early Science data, but it can be used on EVLA data as well. It is distributed with the analysisUtils Python module which provides numerous other Python extensions callable from within CASA. The function parameters and some relevant examples are given below. The current list of parameters can always be obtained by giving parameter: '''help=True'''. This function has been tested on dual-polarization ALMA data, and single-, dual-, and full-polarization EVLA data.  
The '''plotbandpass''' function provides extended capabilities beyond those of CASA's plotcal. It was written to facilitate rapid visualization of ALMA's bandpass and spectral Tsys measurements in Science Verification and Early Science data, but it can be used on EVLA data as well. It is distributed with the analysisUtils Python module which provides numerous other Python extensions callable from within CASA. Analysis Utilities has been tested on dual-polarization ALMA data, and single-, dual-, and full-polarization EVLA data. The function parameters can always be obtained by giving parameter: '''help=True'''.  For example,


Regarding the '''showatm=True''' parameter, since most people reducing ALMA data do not load the ASDM_CALWVR table into the MS, the atmospheric transmission curve will be computed by default for 1.0 mm, but you can set the PWV to any value using the '''pwv''' option. If you want the fully correct curve for ALMA, then copy the <tt>CalWVR.xml</tt> file from your ASDM to the present working directory and specify '''pwv='auto' '''.
<source lang="python">
# In CASA
import analysisUtils as au
au.plotbandpass(help=True)
</source>
<code><pre>
$Id: plotbandpass.py,v 1.118 2011/12/13 20:12:47 thunter Exp $
Usage: plotbandpass(caltable='', antenna='', field='', spw='', yaxis='amp',
  xaxis='chan', figfile='', plotrange=[0,0,0,0], help=False, caltable2='',
  overlay='', showflagged=False, timeranges='', buildpdf=False, caltable3='',
  markersize=3, density=108, interactive=True, showpoints='auto',
  showlines='auto', subplot='22', zoom='', poln='', showatm=False, pwv='1.0',
  gs='gs', convert='convert', chanrange='', debug=False, vm='',
  solutionTimeThresholdSeconds=60.0, phase='')
antenna: must be ID (int or string or list), or a single antenna name or list
atm: overlay the atmospheric transmission curve for the weather and elevation
buildpdf: True/False, if True and figfile is set, assemble pngs into a pdf
caltable: a bandpass table, of type B or BPOLY
caltable2: a second cal table, of type BPOLY or B, to overlay on a B table
caltable3: a third cal table, of type BPOLY, to overlay on the first two
chanrange: set xrange (e.g. '5~100') over which to autoscale y-axis for xaxis='freq'
convert: full path for convert command (in case it's not found)
density: dpi to use in creating PNGs and PDFs (default=108)
field: must be an ID, source name, or list thereof; can use trailing *: 'J*'
figfile: the base_name of the png files to save: base_name.antX.spwY.png
gs: full path for ghostscript command (in case it's not found)
help: print this message
interactive: if False, then figfile will run to completion automatically
overlay: 'antenna' or 'time', make 1 plot with different items in colors
showflagged:  show the values of data, even if flagged
markersize: size of points (default=3)
phase: the y-axis limits to use for phase plots when yaxis='both'
plotrange: define axis limits: [x0,x1,y0,y1] where 0,0 means auto
poln: polarizations to plot (e.g. 'XX','YY','RR','LL' or '' for both)
pwv: define the pwv to use for the showatm option: 'auto' or value in mm
  showatm: compute and overlay the atmospheric transmission curve
showlines: draw lines connecting the data (default=T for amp, F for phase)
showpoints: draw points for the data (default=F for amp, T for phase)
solutionTimeThresholdSeconds: consider 2 solutions simultaneous if within this interval (default=60)
spw: must be single ID or list or range (e.g. 0~4, not the original ID)
subplot: 11,22,32 or 42 for rows x columns; options: 11,22,32 or 42 (default=22)
timeranges: show only these timeranges, the first timerange being 1
vm: a result from ValueMapping, returned from a previous call
xaxis: 'chan' or 'freq'
yaxis: 'amp', 'phase', or 'both' == 'ap'; append 'db' for dB
zoom: 'intersect' will zoom to overlap region of caltable with caltable2
</pre></code>


Most of the regular CASA list syntax is accepted (commas and tildes) but negation via exclamation point is not supported, nor is channel selection within an spw. You are welcome to submit a ticket to the ALMA helpdesk if you want either of these capabilities. If you find any problems or unexpected behavior, please report the following to the [http://help.almascience.org helpdesk]:  
''Regarding the '''showatm=True''' parameter, since most people reducing ALMA data do not load the ASDM_CALWVR table into the MS, the atmospheric transmission curve will be computed by default for 1.0 mm, but you can set the PWV to any value using the '''pwv''' option.  If you want the fully correct curve for ALMA, then copy the <tt>CalWVR.xml</tt> file from your ASDM to the present working directory and specify '''pwv='auto' '''. ''
 
Most of the regular CASA list syntax is accepted (commas and tildes) but negation via exclamation point is not supported, nor is channel selection within an spw. You are welcome to submit a ticket to the ALMA helpdesk if you want either of these capabilities. If you find any problems or unexpected behavior, please include the following information in your bug or enhancement request and send to [http://help.almascience.org helpdesk]:  
# the version number you are using, which is printed as the first line of output;  
# the version number you are using, which is printed as the first line of output;  
# the exact command; and  
# the exact command; and  
# a pointer to the dataset.
# a pointer to the dataset.

Revision as of 12:16, 16 December 2011

Return to Analysis Utilities

This page documents the plotbandpass function of Python module analysisUtils.

The plotbandpass function provides extended capabilities beyond those of CASA's plotcal. It was written to facilitate rapid visualization of ALMA's bandpass and spectral Tsys measurements in Science Verification and Early Science data, but it can be used on EVLA data as well. It is distributed with the analysisUtils Python module which provides numerous other Python extensions callable from within CASA. Analysis Utilities has been tested on dual-polarization ALMA data, and single-, dual-, and full-polarization EVLA data. The function parameters can always be obtained by giving parameter: help=True. For example,

# In CASA
import analysisUtils as au
au.plotbandpass(help=True)
$Id: plotbandpass.py,v 1.118 2011/12/13 20:12:47 thunter Exp $
Usage: plotbandpass(caltable='', antenna='', field='', spw='', yaxis='amp',
   xaxis='chan', figfile='', plotrange=[0,0,0,0], help=False, caltable2='',
   overlay='', showflagged=False, timeranges='', buildpdf=False, caltable3='',
   markersize=3, density=108, interactive=True, showpoints='auto',
   showlines='auto', subplot='22', zoom='', poln='', showatm=False, pwv='1.0',
   gs='gs', convert='convert', chanrange='', debug=False, vm='',
   solutionTimeThresholdSeconds=60.0, phase='')
 antenna: must be ID (int or string or list), or a single antenna name or list
 atm: overlay the atmospheric transmission curve for the weather and elevation
 buildpdf: True/False, if True and figfile is set, assemble pngs into a pdf
 caltable: a bandpass table, of type B or BPOLY
 caltable2: a second cal table, of type BPOLY or B, to overlay on a B table
 caltable3: a third cal table, of type BPOLY, to overlay on the first two
 chanrange: set xrange (e.g. '5~100') over which to autoscale y-axis for xaxis='freq'
 convert: full path for convert command (in case it's not found)
 density: dpi to use in creating PNGs and PDFs (default=108)
 field: must be an ID, source name, or list thereof; can use trailing *: 'J*'
 figfile: the base_name of the png files to save: base_name.antX.spwY.png
 gs: full path for ghostscript command (in case it's not found)
 help: print this message
 interactive: if False, then figfile will run to completion automatically
 overlay: 'antenna' or 'time', make 1 plot with different items in colors
 showflagged:  show the values of data, even if flagged
 markersize: size of points (default=3)
 phase: the y-axis limits to use for phase plots when yaxis='both'
 plotrange: define axis limits: [x0,x1,y0,y1] where 0,0 means auto
 poln: polarizations to plot (e.g. 'XX','YY','RR','LL' or '' for both)
 pwv: define the pwv to use for the showatm option: 'auto' or value in mm
 showatm: compute and overlay the atmospheric transmission curve
 showlines: draw lines connecting the data (default=T for amp, F for phase)
 showpoints: draw points for the data (default=F for amp, T for phase)
 solutionTimeThresholdSeconds: consider 2 solutions simultaneous if within this interval (default=60)
 spw: must be single ID or list or range (e.g. 0~4, not the original ID)
 subplot: 11,22,32 or 42 for rows x columns; options: 11,22,32 or 42 (default=22)
 timeranges: show only these timeranges, the first timerange being 1
 vm: a result from ValueMapping, returned from a previous call
 xaxis: 'chan' or 'freq'
 yaxis: 'amp', 'phase', or 'both' == 'ap'; append 'db' for dB
 zoom: 'intersect' will zoom to overlap region of caltable with caltable2

Regarding the showatm=True parameter, since most people reducing ALMA data do not load the ASDM_CALWVR table into the MS, the atmospheric transmission curve will be computed by default for 1.0 mm, but you can set the PWV to any value using the pwv option. If you want the fully correct curve for ALMA, then copy the CalWVR.xml file from your ASDM to the present working directory and specify pwv='auto' .

Most of the regular CASA list syntax is accepted (commas and tildes) but negation via exclamation point is not supported, nor is channel selection within an spw. You are welcome to submit a ticket to the ALMA helpdesk if you want either of these capabilities. If you find any problems or unexpected behavior, please include the following information in your bug or enhancement request and send to helpdesk:

  1. the version number you are using, which is printed as the first line of output;
  2. the exact command; and
  3. a pointer to the dataset.