#################################################################### ##################### General information ########################### ########## This script contains the steps that have already been performed to create the measurement sets SDP.81_Band4_continuum.ms and SDP.81_Band4.ms, which can be found in the Imaging directory. This script does not need to run by the students, it is purely informative. #################################################################### ##################### Startup ########################### ########## importing packages and aliases import re import os import analysisUtils as au import analysisUtils as aU es = aU.stuffForScienceDataReduction() ################################################################### # NOTE spectral averaging was already done (every 16 channels for # continuum and 20 for the line data) in the individual calibration # scripts. # Concatenate datasets: band4_calibrated_data = ['uid___A002_X91068d_Xa1e.ms.split.cal.science', 'uid___A002_X91068d_Xc50.ms.split.cal.science', 'uid___A002_X91068d_Xead.ms.split.cal.science', 'uid___A002_X915f1c_X5d0.ms.split.cal.science', 'uid___A002_X924c91_X12c.ms.split.cal.science', 'uid___A002_X924c91_X34f.ms.split.cal.science', 'uid___A002_X924c91_X572.ms.split.cal.science', 'uid___A002_X924c91_X847.ms.split.cal.science', 'uid___A002_X924c91_Xa47.ms.split.cal.science', 'uid___A002_X93514d_X1634.ms.split.cal.science', 'uid___A002_X93514d_X1857.ms.split.cal.science', 'uid___A002_X93514d_Xfc6.ms.split.cal.science'] os.system('rm -rf SDP.81_Band4.ms') concat(vis = band4_calibrated_data, concatvis = 'SDP.81_Band4.ms') # Define spectral windows for the continuum. spw_cont = '0~2,4~6,8~10,12~14,16~18,20~22,24~26,28~30,32~34,36~38,40~42,44~46' os.system('rm -rf SDP.81_Band4_continuum.ms') split(vis='SDP.81_Band4.ms',outputvis='SDP.81_Band4_continuum.ms', spw=spw_cont,datacolumn='data') ### prepare a dataset with only the bandpass calibrator, from only one execution: split(vis='SDP81_B4_uncalibrated.ms.split',outputvis='bandpass.ms',field='0',datacolumn='corrected',spw='0,1,2')