Tutorial 3: Automated Self-Calibration#

Link to Script

This script was written for CASA 6.6.4.34
Prepared for ALMA Data Processing Workshop @ University of Victoria
J. Speedie (jspeedie@uvic.ca) October 2024

Additional scripts needed: *.py from
https://github.com/jjtobin/auto_selfcal

Dataset: TW Hya calibrated continuum measurement set (435 MB)

twhya_calibrated.ms.tar (Download)

Download auto_selfcal scripts#

Additional scripts needed!

Download *.py from jjtobin/auto_selfcal

Your working directory should look like:

TUTORIAL3 >> ls
auto_selfcal-main/
script_auto_selfcal_setup.py
twhya_calibrated.ms.tar

Untar with:

tar -xvzf twhya_calibrated.ms.tar

Prepare MS for auto_selfcal#

Here, we generate what the auto_selfcal.py script expects.

os.system('rm -rf twhya_target.ms') # We're about to create this

Split out the science target:

split(vis='twhya_calibrated.ms', # Old measurement set (435 MB)
      field='5', # TW Hya
      width='8', # For tutorial purposes; in practice, better not to
      outputvis='twhya_target.ms',  # New measurement set (43 MB)
      datacolumn='data')

Inspect the new MS:

listobs(vis='twhya_target.ms',
        listfile='twhya_target_listobs.txt')

Exit CASA.

Run automated self-calibration#

Copy the auto_selfcal/ scripts into the same directory as twhya_target.ms.

Run the auto_selfcal.py script:

casa -c auto_selfcal.py

Open and inspect the resulting Weblog! weblog/index.html

Open & inspect resulting Target_TW_Hya_Band7* files in CARTA!

Slides#

SRDP_autoselfcal

Download Slides