From 171674f72bdaae19f0341ea9a0e64b65c6711e4e Mon Sep 17 00:00:00 2001 From: Victor de Souza magalhaes Date: Wed, 23 Aug 2023 14:19:34 -0600 Subject: [PATCH] Locit: CASA pointing data reduction script --- etc/casa/point-redux.py | 158 ++++++++++++++++++++++++++++++++++++++++ 1 file changed, 158 insertions(+) create mode 100644 etc/casa/point-redux.py diff --git a/etc/casa/point-redux.py b/etc/casa/point-redux.py new file mode 100644 index 00000000..e8dfe064 --- /dev/null +++ b/etc/casa/point-redux.py @@ -0,0 +1,158 @@ +################################################################################ +### ### +### Suggested path forward ### +### ### +################################################################################ + +# Fringe fitting a single source i.e. BlLac seems to be the best +# compromise between flagging and having reasonable solutions. +# However a few crap scans remain, specially 57, 117 & 120 have 360 phase +# spread because of their low elevation, 120 even being below the horizon + +# The pointing ms name +msname = 'pointing-2015-10-04.ms' +# Root for the products of this script +nameroot = '2015-10-04' +# Scans to flag, for example due to low elevation +scanstoflag = '' # here no scans below 5 degrees +# Reference antenna +ref_ant = 'ea04' # ea02 would also be a good choice(both are very close) +# Fringe fit source +frg_src = '0319+415' # Seems to be the brightest, old friend from the 30m pnt +# Number of channels +n_chan = 64 + +# Moved antennas are: 1, 6, 7, 12, 17, 20, 26, 28 + +# Extensions +ext_pnt = '-pnt.ms' +ext_frg = '-frg.cal' +ext_avg = '-avg.ms' +ext_cal = '-pha.cal' + +# The desired intent +pnt_intent = 'CALIBRATE_POINTING#ON_SOURCE' + + +################################################################################ +# Here starts the script, DO NOT CHANGE BELOW +################################################################################ + +# file names +point_only = nameroot+ext_pnt +fring_cal = nameroot+ext_frg +avg_data = nameroot+ext_avg +gaincal_tab = nameroot+ext_cal + +split(vis = msname, # Name of input visibility file + outputvis = point_only, # Name of output visibility file + keepmms = True, # keep Multi MS + field = '', # Select field using field id(s) or field name(s) + spw = '', # Select spectral window/channels + scan = '', # Scan number range + antenna = '', # Select data based on antenna/baseline + correlation = '', # Select data based on correlation + timerange = '', # Select data based on time range + intent = pnt_intent, # Select observing intent + array = '', # Select (sub)array(s) by array ID number. + uvrange = '', # Select data by baseline length. + observation = '', # Select by observation ID(s) + feed = '', # Multi-feed numbers: Not yet implemented. + datacolumn = 'data', # Which data column(s) to process. + keepflags = True, # + width = 1, # Number of channels to average + timebin = '0s', # Bin width for time averaging +) + +if len(scanstoflag) > 0: + flagdata(vis = point_only, # Name of input visibility file + mode = 'manual', # Flagging mode + autocorr = False, # Flag only the auto-correlations? + spw = '', # Select spectral window/channels + field = '', # Select field id(s) or field name(s) + antenna = '', # Select data based on antenna/baseline + uvrange = '', # Select data by baseline length. + timerange = '', # Select data based on time range + correlation = '', # Select data based on correlation + scan = scanstoflag, # Scan number range + intent = '', # Select observing intent + array = '', # (Sub)array numbers + observation = '', # Select by observation ID(s) + action = 'apply', # Action to perform + display = 'report', # Display (data/report/both). + flagbackup = False, # Save in flagversions + savepars = False, # Save parameters to the FLAG_CMD + writeflags = True, # Do not modify. + ) + # Create our new flag state under the name 'baseflags' + flagmanager(vis=point_only, mode='save', versionname='baseflags') + + + flagmanager(vis=point_only, mode='restore', versionname='baseflags') + + +fringefit(vis = point_only, # Name of input visibility file + caltable = fring_cal, # Name of output gain calibration table + field = frg_src, # field names + solint = 'inf', # Solution interval: egs. 'inf', '60s' + refant = ref_ant, # Reference antenna name(s) + minsnr = 3.0, # Reject solutions below this snr + zerorates = True, # Zero delay-rates in solution table + globalsolve = True, # Refine estimates with global lst-sq solver + niter = 100, # Maximum number of iterations + corrdepflags= False, # Respect correlation-dependent flags + paramactive = [], # Control which parameters are solved for + parang = False, # Apply para angle correction on the fly +) + +applycal(vis=point_only, + gaintable=[fring_cal], + interp=['nearest'], + parang=False) + + +# Now we create a new dataset that is colapsed on the channel axis +# within each spw, also create a flagversion to store current flag +# state on the averaged MS +split(vis = point_only, # Name of input visibility file + outputvis = avg_data, # Name of output visibility file + keepmms = True, # keep Multi-ms + field = '', # Select field field name(s) + spw = '', # Select spectral window/channels + scan = '', # Scan number range + antenna = '', # Select data based on antenna/baseline + correlation = '', # Select data based on correlation + timerange = '', # Select data based on time range + intent = '', # Select observing intent + array = '', # Select (sub)array(s) by array ID number. + uvrange = '', # Select data by baseline length. + observation = '', # Select by observation ID(s) + feed = '', # Multi-feed numbers: Not yet implemented. + datacolumn = 'corrected', # Which data column(s) to process. + keepflags = False, # + width = n_chan, # Number of channels to average + timebin = '0s', # Bin width for time averaging +) +flagmanager(vis=avg_data, mode='save', versionname='original') + + +flagmanager(vis=avg_data, mode='restore', versionname='original') +gaincal(vis = avg_data, # Name of input visibility file + caltable = gaincal_tab, # Name of output gain cal table + field = '', # Select field using field name(s) + solint = '10min', # Solution interval + refant = ref_ant, # Reference antenna name(s) + refantmode = 'flex', # Reference antenna mode + minblperant= 3, # Minimum baselines _per antenna_ + minsnr = 3.0, # Reject solutions below this SNR + gaintype = 'G', # G is gain + calmode = 'p', # Type of solution" ('ap', 'p', 'a') + parang = False, # Apply parallactic angle correction + solmode = 'L1' # Solution mode L1 => least squares +) + +applycal(vis=avg_data, + gaintable=[gaincal_tab], + interp=['nearest'], + parang=False) +