Self-calibration usage modes for initial capability:
- Phase-only, continuum-only
- Amplitude attempts as future development possibility
- Single pointing observations (no mosaics)
- Aggregate continuum self-calibration (gaincal combine=spw, applycal spwmap)
- Long and short baseline observations
Later Usage modes to consider:
- Support for self-calibration within mosaics (interpolation over entire scans).
- Per-spw calibration
- allowing applymode=calonly for some data to attempt to solve for those data by cleaning deeper holding solint constant
Selfcal Prototype:
Obtain via 'git clone https://github.com/jjtobin/auto_selfcal'
- Copy scripts to working directory where _target.ms files reside
- Start CASA in MPI mode (i.e., mpicasa -n 9 /path/to/casa/bin/casa); otherwise edit script to set parallel=False
- run via execfile('auto_selfcal.py')
- will make use of cont.dat to flag line channels if cont.dat is present, otherwise will use all channels
- Can make use of spectral averaging by setting spectral_average=True in auto_selfcal.py
- averaging depends on frequency for VLA, will simply average to 15.625 MHz channels for ALMA
Test Suite:
| Telescope | Project Code | MOUS/FSID | Band | Conf. | EBs | Tgts | Version | Notes | Original Peak SNR | Selfcal Peak SNR | min solint | Original RMS (mJy/bm or uJy/beam) | Final RMS (mJy/bm or uJy/beam) |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| ALMA | 2019.1.01709.S | uid://A001/X1465/X435 | 6 | C-3 | 2 | 1 | Jan 5 | All cont spw Strong point-like continuum source and faint radio jet separated by 3x beam size Final selfcal image captures the faint radio jet better than the original PL image does | PKS1030-357: 109 | PKS1030-357: 254 | PKS1030-357: 10.08s | PKS1030-357: 0.053 | PKS1030-357: 0.030 |
| ALMA | 2019.1.01709.S | uid://A001/X1465/X431 | 3 | C-5 | 1 | 1 | Jan 5 | All cont spw strong point source continuum Selfcal failed due to increased final RMS. But actual images (original PL image and selfcal image) are not very different in terms of RMS and peak SNR. This dataset needs some investigation. | 4C-03.79: 5 | 4C-03.79: 715 Seems to be increased by erroneous estimate of the original RMS | 4C-03.79: None | 4C-03.79: 0.009 This reported value (9.8e-6 Jy) from pickle file seems wrong. Original PL image has RMS of 1 mJy/beam | 4C-03.79: 1 |
| ALMA | 2019.1.01197.S | uid://A001/X1465/X15aa | 6 | C-3 | 1 | 1 | Jan 5 | All cont spw Continuum source with lensed morphology | PJ011646.8: 34 | PJ011646.8: 38 | PJ011646.8: 52.42s | PJ011646.8: 0.1 | PJ011646.8: 0.09 |
| ALMA | 2019.1.00463.S | uid://A001/X146a/Xf | 6 | C-6 | 1 | 1 | Jan 5 | mixture of narrow/wide band width continuum source off phase center | BHR7: 322 | BHR7: 742 | BHR7: 42.34s | BHR7: 0.072 | BHR7: 0.043 |
| ALMA | 2019.1.01425.S | uid://A001/X1465/Xd60 | 6 | C-10 | 1 | 5 | Nov | uses full MS low-SB extended emission surrounding point-like sources | L1448IRS3B: 85 L1448IRS1: 97 L1448IRS3C: 90 Per-emb-17: 100 Per-emb-22: 107 | L1448IRS3B: 93 L1448IRS1: 101 L1448IRS3C: 97 Per-emb-17: 108 Per-emb-22: 150 | L1448IRS3B: inf_EB L1448IRS1: inf_EB L1448IRS3C: inf_EB Per-emb-17: inf_EB Per-emb-22: inf | L1448IRS3B: 0.033 L1448IRS1: 0.031 L1448IRS3C: 0.031 Per-emb-17: 0.031 Per-emb-22: 0.031 | L1448IRS3B: 0.032 L1448IRS1: 0.031 L1448IRS3C: 0.031 Per-emb-17: 0.031 Per-emb-22: 0.031 |
| ALMA | 2017.1.01413.S | uid://A001/X1273/X7bc | 6 | C-4 | 2 | 5 | Dec 3 | uses full MSes, cont.dat, test case for spectral averaging compact, disk sources | L1527-IRS Band_6: 471 TMC1A Band_6: 411 IRAS_04302+2247 Band_6: 158 L1489IRS Band_6: 115 TMC1 Band_6: 211 | L1527-IRS Band_6: 1729 TMC1A Band_6: 2002 IRAS_04302+2247 Band_6: 157 L1489IRS Band_6: 135 TMC1 Band_6: 295 | L1527-IRS Band_6: 30.24s TMC1A Band_6: 30.24s IRAS_04302+2247 Band_6: None L1489IRS Band_6: inf TMC1 Band_6: 30.24s | L1527-IRS Band_6: 0.214 TMC1A Band_6: 0.299 IRAS_04302+2247 Band_6: 0.131 L1489IRS Band_6: 0.058 TMC1 Band_6: 0.073 | L1527-IRS Band_6: 0.061 TMC1A Band_6: 0.065 IRAS_04302+2247 Band_6: 0.131 L1489IRS Band_6: 0.054 TMC1 Band_6: 0.057 |
| ALMA | '' | '' | '' | '' | '' | '' | Nov | does not use spectral averaging compact, disk sources | L1527-IRS: 451 TMC1A: 398 IRAS_04302+2247: 151 L1489IRS: 110 TMC1: 227 | L1527-IRS: 1681 TMC1A: 1909 IRAS_04302+2247: 231 L1489IRS: 139 TMC1: 324 | L1527-IRS: inf TMC1A: 30.24s IRAS_04302+2247: None L1489IRS: inf TMC1: 12.10s | L1527-IRS: 0.222 TMC1A: 0.307 IRAS_04302+2247: 0.137 L1489IRS: 0.061 TMC1: 0.068 | L1527-IRS: 0.062 TMC1A: 0.068 IRAS_04302+2247: 0.099 L1489IRS: 0.052 TMC1: 0.051 |
| ALMA | 2019.1.00261.L | uid://A001/X146c/X26f | 6 | C-5 | 2 | 2 | Dec 3 | spectrally averaged, 2 spws, used for fast runs | IRAS32: 507 IRS5N: 588 | IRAS32: 1244 IRS5N: 1039 | IRAS32: 18.14s IRS5N: int | IRAS32: 0.092 IRS5N: 0.061 | IRAS32: 0.039 IRS5N: 0.036 |
| ALMA | 2018.1.00427.S | uid://A001/X133d/X1b87 | 6 | C-9 | 2 | 1 | Dec 3 | uses only continuum spw compact sources | IRAS2A Band_6: 224 | IRAS2A Band_6: 995 | IRAS2A Band_6: inf | IRAS2A Band_6: 0.109 | IRAS2A Band_6: 0.041 |
| ALMA | 2015.1.01503.S | uid://A001/X2d8/X287 | 7 | C-4 | 6 | 11 | Dec 3 | polarization project, all data in a single combined MS, uses TDM sources range from very extended to somewhat extended | Per26 Band_7: 224 Per18 Band_7: 174 Per21 Band_7: 200 Per11 Band_7: 270 Per14 Band_7: 266 Per50 Band_7: 223 Per2 Band_7: 161 Per5 Band_7: 292 Per41 Band_7: 164 Per29 Band_7: 251 | Per26 Band_7: 1162 Per18 Band_7: 583 Per21 Band_7: 704 Per11 Band_7: 1029 Per14 Band_7: 794 Per50 Band_7: 1157 Per2 Band_7: 390 Per5 Band_7: 1087 Per41 Band_7: 283 Per29 Band_7: 1149 | Per26 Band_7: int Per18 Band_7: int Per21 Band_7: int Per11 Band_7: int Per14 Band_7: int Per50 Band_7: int Per2 Band_7: int Per5 Band_7: int Per41 Band_7: 14.11s Per29 Band_7: int | Per26 Band_7: 0.976 Per18 Band_7: 0.493 Per21 Band_7: 0.368 Per11 Band_7: 0.601 Per14 Band_7: 0.329 Per50 Band_7: 0.553 Per2 Band_7: 0.485 Per5 Band_7: 0.759 Per41 Band_7: 0.148 Per29 Band_7: 0.643 | Per26 Band_7: 0.231 Per18 Band_7: 0.170 Per21 Band_7: 0.128 Per11 Band_7: 0.180 Per14 Band_7: 0.130 Per50 Band_7: 0.135 Per2 Band_7: 0.220 Per5 Band_7: 0.233 Per41 Band_7: 0.098 Per29 Band_7: 0.163 |
| ALMA | 2018.1.00962.S | uid://A001/X133d/X2727 | 9 | C-3 | 1 | 1 | Nov | uses full MSes, cont.dat, extended emission | NGC6334I: 46 | NGC6334I: 106 | NGC6334I: int | NGC6334I: 62.618 | NGC6334I: 28.644 |
| ALMA | 2017.1.00586.S | uid://A001/X1284/X20cf | 8 | C-2? | 3 | 2 | Dec 3 | uses full MSes, cont.dat, some extended emission, collections of point sources | NGC_1068 Band_8: 61 NGC_1068_SBring Band_8: 49 NGC_1068_CH3OH Band_8: 11 | NGC_1068 Band_8: 77 NGC_1068_SBring Band_8: 61 NGC_1068_CH3OH Band_8: 11 | NGC_1068 Band_8: inf_EB NGC_1068_SBring Band_8: inf_EB NGC_1068_CH3OH Band_8: None | NGC_1068 Band_8: 0.322 NGC_1068_SBring Band_8: 0.301 NGC_1068_CH3OH Band_8: 0.237 | NGC_1068 Band_8: 0.266 NGC_1068_SBring Band_8: 0.257 NGC_1068_CH3OH Band_8: 0.237 |
| ALMA | 2018.1.00510.S | 6 | C-9 | 2 | 1 | Dec 3 | uses 3 TDM windows fails due to beam, calonly would help here well-resolved in C-10 | NGC1333_IRAS_4A Band_6: 129 | NGC1333_IRAS_4A Band_6: 129 | NGC1333_IRAS_4A Band_6: None | NGC1333_IRAS_4A Band_6: 0.143 | NGC1333_IRAS_4A Band_6: 0.143 | |
| '' | '' | '' | '' | '' | '' | '' | Dec 3 | '', calonly | NGC1333_IRAS_4A Band_6: 129 | NGC1333_IRAS_4A Band_6: 460 | NGC1333_IRAS_4A Band_6: int | NGC1333_IRAS_4A Band_6: 0.143 | NGC1333_IRAS_4A Band_6: 0.054 |
| ALMA-ACA | 2017.1.00857.S | uid://A001/X12c4/X217 | 8 | ACA | 8 | 4 | Dec 3 | uses full MSes, cont.dat, point sources | DR_Tau Band_8: 126 DO_Tau Band_8: 38 DL_Tau Band_8: 38 FZ_Tau Band_8: 21 | DR_Tau Band_8: 272 DO_Tau Band_8: 227 DL_Tau Band_8: 279 FZ_Tau Band_8: 21 | DR_Tau Band_8: inf DO_Tau Band_8: inf DL_Tau Band_8: int FZ_Tau Band_8: None | DR_Tau Band_8: 5.514 DO_Tau Band_8: 12.132 DL_Tau Band_8: 17.709 FZ_Tau Band_8: 2.403 | DR_Tau Band_8: 2.578 DO_Tau Band_8: 2.227 DL_Tau Band_8: 2.648 FZ_Tau Band_8: 2.403 |
| ALMA-ACA | 2017.1.00857.S | uid://A001/X12c4/X21f | 8 | ACA | 1 | 1 | Dec 3 | full MSes, cont.dat point source | AS_205A: 63 | AS_205A: 385 | AS_205A: inf | AS_205A: 24.770 | AS_205A: 4.102 |
| VLA | 13B-454 | 13B-454.sb29509255.eb29947371.56956.79185517361 | C | C | 1 | 1 | Dec 6 | M82, high DR, extended | SN2014J EVLA_C: 1506 | SN2014J EVLA_C: 6393 | SN2014J EVLA_C: 17.00s | SN2014J EVLA_C: 0.105 | SN2014J EVLA_C: 0.030 |
VLA | 13B-454 | 13B-454.sb29105214.eb29105303.56758.17242523148 | C | A | 1 | 1 | Nov | M82, high DR, extended Very long run time | SN2014J EVLA_C: 5933 | SN2014J EVLA_C: 6078 | SN2014J EVLA_C: 35.00s | SN2014J EVLA_C: 0.005 | SN2014J EVLA_C: 0.005 |
| VLA | 15A-370 | 15A-370.sb30976913.eb30980293.57204.36680753472 15A-370.sb30976913.eb30978455.57203.32521244213 | Ka, Q | A | 2 | 1 | Dec 3 | compact, low DR | L1157 EVLA_KA: 28 L1157 EVLA_Q: 19 | L1157 EVLA_KA: 27 L1157 EVLA_Q: 19 | L1157 EVLA_KA: None L1157 EVLA_Q: None | L1157 EVLA_KA: 0.011 L1157 EVLA_Q: 0.020 | L1157 EVLA_KA: 0.011 L1157 EVLA_Q: 0.020 |
| VLA | 17B-188 | 17B-188.sb35564398.eb35590549.58363.10481791667 | S,C,Ku,K,Ka | C | 1 | 2 | Dec 3 | compact at HF, compact+extended at LF, multiband | 14C_KuQ_C EVLA_K: 110 14C_KuQ_C EVLA_KA: 73 14C_KuQ_C EVLA_KU: 185 14C_SC_C EVLA_C: 1196 14C_SC_C EVLA_S: 1558 | 14C_KuQ_C EVLA_K: 212 14C_KuQ_C EVLA_KA: 116 14C_KuQ_C EVLA_KU: 257 14C_SC_C EVLA_C: 1384 14C_SC_C EVLA_S: 1635 | 14C_KuQ_C EVLA_K: int 14C_KuQ_C EVLA_KA: inf 14C_KuQ_C EVLA_KU: 33.00s 14C_SC_C EVLA_C: inf 14C_SC_C EVLA_S: int | 14C_KuQ_C EVLA_K: 0.079 14C_KuQ_C EVLA_KA: 0.080 14C_KuQ_C EVLA_KU: 0.068 14C_SC_C EVLA_C: 0.021 14C_SC_C EVLA_S: 0.025 | 14C_KuQ_C EVLA_K: 0.050 14C_KuQ_C EVLA_KA: 0.058 14C_KuQ_C EVLA_KU: 0.052 14C_SC_C EVLA_C: 0.018 14C_SC_C EVLA_S: 0.024 |
| VLA | TDRW0001 | TDRW0001.sb35624494.eb35628826.58395.23719237269 | S | D | 1 | 1 | Dec 3 | extended, high DR | 3C75 EVLA_S: 3815 | 3C75 EVLA_S: 6733 | 3C75 EVLA_S: 25.00s | 3C75 EVLA_S: 0.038 | 3C75 EVLA_S: 0.022 |
| VLA | 16B-427 | 16B-427.sb33651951.eb34041883.57963.57173581018 | L,S,C,X,K, Q | C | 1 | 1 | Dec 3 | multiband | S255 NIRS3 EVLA_C: 1107 S255 NIRS3 EVLA_K: 473 S255 NIRS3 EVLA_L: 570 S255 NIRS3 EVLA_Q: 250 S255 NIRS3 EVLA_S: 596 S255 NIRS3 EVLA_X: 1849 | S255 NIRS3 EVLA_C: 1173 S255 NIRS3 EVLA_K: 681 S255 NIRS3 EVLA_L: 550 S255 NIRS3 EVLA_Q: 564 S255 NIRS3 EVLA_S: 622 S255 NIRS3 EVLA_X: 2267 | S255 NIRS3 EVLA_C: 80.00s S255 NIRS3 EVLA_K: 21.00s S255 NIRS3 EVLA_L: inf S255 NIRS3 EVLA_Q: 36.00s S255 NIRS3 EVLA_S: inf_EB S255 NIRS3 EVLA_X: 9.00s | S255 NIRS3 EVLA_C: 0.014 S255 NIRS3 EVLA_K: 0.053 S255 NIRS3 EVLA_L: 0.059 S255 NIRS3 EVLA_Q: 0.109 S255 NIRS3 EVLA_S: 0.025 S255 NIRS3 EVLA_X: 0.011 | S255 NIRS3 EVLA_C: 0.013 S255 NIRS3 EVLA_K: 0.041 S255 NIRS3 EVLA_L: 0.061 S255 NIRS3 EVLA_Q: 0.069 S255 NIRS3 EVLA_S: 0.024 S255 NIRS3 EVLA_X: 0.010 |
| VLA | 20A-390 | 20A-390.sb37850360.eb37919292.58894.32864972222 20A-390.sb37850360.eb37920910.58895.323193819444 20A-390.sb37850360.eb37926293.58898.356132071756 20A-390.sb37850360.eb37933285.58900.24178115741 20A-390.sb37850360.eb37933287.58900.32599280092 | Ka | C | 5 | 1 | Dec 3 | point source, low dec | TWHya EVLA_KA: 164 | TWHya EVLA_KA: 278 | TWHya EVLA_KA: 129.00s | TWHya EVLA_KA: 0.005 | TWHya EVLA_KA: 0.005 |
| VLA | 21A-178 | 21A-178.sb39472833.eb39589695.59326.113018460645 | X | D | 1 | 21 | Nov | variety of S/N | SN 2014L EVLA_X: 113 SN 2004gk EVLA_X: 37 SN 1990B EVLA_X: 173 SN 2004cc EVLA_X: 176 SN 1987K EVLA_X: 661 PTF10bfz EVLA_X: 421 SN 1997dq EVLA_X: 17 SN 2011bm EVLA_X: 209 SN 2015bn EVLA_X: 217 SN 2003dg EVLA_X: 15 SN 2007I EVLA_X: 32 SN 2004gn EVLA_X: 456 SN 2007rw EVLA_X: 18 SN 2000ew EVLA_X: 21 SN 2014ad EVLA_X: 33 SN 2002bm EVLA_X: 60 SN 2008aq EVLA_X: 14 SN 2005bh EVLA_X: 6 SN 2004gj EVLA_X: 34 SN 2002ji EVLA_X: 76 LSQ13abf EVLA_X: 9 | SN 2014L EVLA_X: 114 SN 2004gk EVLA_X: 37 SN 1990B EVLA_X: 174 SN 2004cc EVLA_X: 175 SN 1987K EVLA_X: 732 PTF10bfz EVLA_X: 429 SN 1997dq EVLA_X: 17 SN 2011bm EVLA_X: 216 SN 2015bn EVLA_X: 240 SN 2003dg EVLA_X: 15 SN 2007I EVLA_X: 32 SN 2004gn EVLA_X: 464 SN 2007rw EVLA_X: 18 SN 2000ew EVLA_X: 21 SN 2014ad EVLA_X: 34 SN 2002bm EVLA_X: 61 SN 2008aq EVLA_X: 14 SN 2005bh EVLA_X: 6 SN 2004gj EVLA_X: 35 SN 2002ji EVLA_X: 76 LSQ13abf EVLA_X: 9 | SN 2014L EVLA_X: 9.00s SN 2004gk EVLA_X: 81.00s SN 1990B EVLA_X: inf_EB SN 2004cc EVLA_X: inf_EB SN 1987K EVLA_X: 18.00s PTF10bfz EVLA_X: 39.00s SN 1997dq EVLA_X: inf_EB SN 2011bm EVLA_X: int SN 2015bn EVLA_X: 39.00s SN 2003dg EVLA_X: None SN 2007I EVLA_X: inf_EB SN 2004gn EVLA_X: inf_EB SN 2007rw EVLA_X: inf_EB SN 2000ew EVLA_X: inf_EB SN 2014ad EVLA_X: 81.00s SN 2002bm EVLA_X: 39.00s SN 2008aq EVLA_X: None SN 2005bh EVLA_X: None SN 2004gj EVLA_X: 81.00s SN 2002ji EVLA_X: 162.00s LSQ13abf EVLA_X: None | SN 2014L EVLA_X: 0.014 SN 2004gk EVLA_X: 0.011 SN 1990B EVLA_X: 0.013 SN 2004cc EVLA_X: 0.013 SN 1987K EVLA_X: 0.046 PTF10bfz EVLA_X: 0.012 SN 1997dq EVLA_X: 0.013 SN 2011bm EVLA_X: 0.012 SN 2015bn EVLA_X: 0.015 SN 2003dg EVLA_X: 0.012 SN 2007I EVLA_X: 0.012 SN 2004gn EVLA_X: 0.014 SN 2007rw EVLA_X: 0.012 SN 2000ew EVLA_X: 0.012 SN 2014ad EVLA_X: 0.013 SN 2002bm EVLA_X: 0.012 SN 2008aq EVLA_X: 0.012 SN 2005bh EVLA_X: 0.012 SN 2004gj EVLA_X: 0.012 SN 2002ji EVLA_X: 0.012 LSQ13abf EVLA_X: 0.011 | SN 2014L EVLA_X: 0.014 SN 2004gk EVLA_X: 0.011 SN 1990B EVLA_X: 0.013 SN 2004cc EVLA_X: 0.013 SN 1987K EVLA_X: 0.042 PTF10bfz EVLA_X: 0.012 SN 1997dq EVLA_X: 0.013 SN 2011bm EVLA_X: 0.012 SN 2015bn EVLA_X: 0.014 SN 2003dg EVLA_X: 0.012 SN 2007I EVLA_X: 0.012 SN 2004gn EVLA_X: 0.014 SN 2007rw EVLA_X: 0.012 SN 2000ew EVLA_X: 0.013 SN 2014ad EVLA_X: 0.013 SN 2002bm EVLA_X: 0.012 SN 2008aq EVLA_X: 0.012 SN 2005bh EVLA_X: 0.012 SN 2004gj EVLA_X: 0.012 SN 2002ji EVLA_X: 0.012 LSQ13abf EVLA_X: 0.011 |