imagename = '18A-229_Q_singlefield_selfcal_iter1_wterms'
myclean(vis=cont_vis,
fields="Sgr B2 N Q,Sgr B2 NM Q,Sgr B2 MS Q".split(","),
spws='',
imsize=1000,
phasecenters={"Sgr B2 N Q":'J2000 17h47m19.897 -28d22m17.340',
"Sgr B2 NM Q":'J2000 17h47m20.166 -28d23m04.968',
"Sgr B2 MS Q":'J2000 17h47m20.166 -28d23m04.968',
"Sgr B2 S Q":'J2000 17h47m20.461 -28d23m45.059',
},
cell='0.01arcsec',
name=imagename,
gridder='wproject',
wprojplanes=64,
niter=10000,
threshold='3mJy',
scales=[0,3,9],
robust=0.5,
savemodel='modelcolumn',
mask=mask,
)
caltable = '18A-229_Q_concatenated_cal_iter1_wterms.cal'
if not os.path.exists(caltable):
gaincal(vis=cont_vis,
caltable=caltable,
field='Sgr B2 N Q,Sgr B2 NM Q,Sgr B2 MS Q',
assert os.getenv('SCRIPT_DIR') is not None
sys.path.append(os.getenv('SCRIPT_DIR'))
from continuum_imaging_general import myclean, tclean, makefits
from continuum_windows import all_Qmses, Kamses, Kmses
Qmses = all_Qmses
calibrators = ['J1744-3116', 'J1733-1304', "1331+305=3C286",]
for ms in Qmses:
myclean('../'+ms,
name=ms[:22],
threshold='5mJy',
spws=Qmses[ms],
niter=1000,
imsize=200,
cell='0.01arcsec',
fields=calibrators,
noneg=False,
)
for ms in Kamses:
myclean('../'+ms,
name=ms[:22],
threshold='5mJy',
cell='0.01arcsec',
niter=1000,
imsize=200,
spws=Kamses[ms],
fields=calibrators,
noneg=False,
for gridder in ('standard', 'wproject', 'widefield', 'mosaic', 'awproject'):
try:
imagename = '18A-229_Q_singlefield_imaging_smallfield_test_{0}'.format(
gridder)
myclean(vis=cont_vis,
fields=selfcal_fields,
spws='',
imsize=1000,
phasecenters={
"Sgr B2 N Q": 'J2000 17h47m19.897 -28d22m17.340',
"Sgr B2 NM Q": 'J2000 17h47m20.166 -28d23m04.968',
"Sgr B2 MS Q": 'J2000 17h47m20.166 -28d23m04.968',
"Sgr B2 S Q": 'J2000 17h47m20.461 -28d23m45.059',
},
cell='0.01arcsec',
name=imagename,
gridder=gridder,
niter=10000,
threshold='0.5mJy',
scales=[0, 3, 9],
robust=0.5,
savemodel='none',
mask=mask,
noneg=False,
**extrapars[gridder])
except Exception as ex:
myprint(ex)
for gridder in ('standard', 'wproject', 'widefield', 'mosaic', 'awproject'):
interp=[],
gaintable=[caltable],
calwt=[False],
vis=cont_ms,
applymode='calonly',
antenna='*&*',
spwmap=[],
parang=True)
if False: # slow version
myclean(
vis=cont_ms,
name=name,
imsize=8000,
cell='0.01arcsec',
fields=['Sgr B2 N Q', 'Sgr B2 NM Q', 'Sgr B2 MS Q'],
threshold='2mJy',
savemodel='modelcolumn',
spws='', # all windows are continuum now
scales=[0, 3, 9, 27],
)
else:
# Diagnostic: make sure the before & after are different
myclean(
vis=cont_ms,
name="before_cutout_" + name,
fields=[
'Sgr B2 MS Q',
],
spws='',
imsize=1000,
Overall imaging script to make combined continuum images for Ka, Q, and K bands
"""
import sys
sys.path.append('.')
print(sys.path)
from continuum_imaging_general import myclean, tclean, makefits
from continuum_windows import Qmses, Kamses, Kmses
for ms in Qmses:
myclean('../'+ms,
name=ms[:22],
threshold='5mJy',
spws=Qmses[ms],
niter=1000,
imsize=2000,
cell='0.04arcsec',
)
for ms in Kamses:
myclean('../'+ms,
name=ms[:22],
threshold='5mJy',
cell='0.05arcsec',
niter=1000,
imsize=2000,
fields=['Sgr B2 MN Ka', 'Sgr B2 MS Ka', 'Sgr B2 S Ka', 'Sgr B2 DS1 Ka', 'Sgr B2 DS2 Ka'],
spws=Kamses[ms],
)
"""
Overall imaging script to make combined continuum images for Ka, Q, and K bands
"""
from continuum_imaging_general import myclean, tclean, makefits
from continuum_windows import Qmses, Kamses, Kmses
myclean(
['../' + x for x in Qmses],
name='18A-229_combined',
threshold='5mJy',
spws=[x for x in Qmses.values()],
)
myclean(
['../' + x for x in Kamses],
name='18A-229_combined',
threshold='5mJy',
cell='0.015arcsec',
fields=[
'Sgr B2 MN Ka', 'Sgr B2 MS Ka', 'Sgr B2 S Ka', 'Sgr B2 DS1 Ka',
'Sgr B2 DS2 Ka'
],
spws=[x for x in Kamses.values()],
)
myclean(
['../' + x for x in Kmses],
name='18A-229_combined',
threshold='5mJy',
cell='0.02arcsec',
myprint("Beginning iteration 1")
imagename = '18A-229_Q_singlefield_selfcal_iter1'
myclean(vis=cont_vis,
fields="Sgr B2 NM Q,Sgr B2 MS Q".split(","),
spws='',
imsize=1000,
phasecenters={"Sgr B2 N Q":'J2000 17h47m19.897 -28d22m17.340',
"Sgr B2 NM Q":'J2000 17h47m20.166 -28d23m04.968',
"Sgr B2 MS Q":'J2000 17h47m20.166 -28d23m04.968',
"Sgr B2 S Q":'J2000 17h47m20.461 -28d23m45.059',
},
cell='0.01arcsec',
name=imagename,
niter=10000,
threshold='3mJy',
scales=[0,3,9],
robust=0.5,
savemodel='modelcolumn',
mask=mask,
)
caltable = '18A-229_Q_concatenated_cal_iter1.cal'
if not os.path.exists(caltable):
mygaincal(vis=cont_vis,
caltable=caltable,
field='Sgr B2 NM Q,Sgr B2 MS Q',
calmode='p',
gaintables = [x for x in pipeline_tables if os.path.exists(x)
] + [caltable_nocombine]
gainfield = ['' for x in pipeline_tables if os.path.exists(x)] + []
calwt = [False for x in pipeline_tables if os.path.exists(x)] + [False]
spwmap = [[] for x in pipeline_tables if os.path.exists(x)] + []
interp = [
'linear,nearestflag' if 'BPcal' in x else ''
for x in pipeline_tables if os.path.exists(x)
] + ['']
applycal(flagbackup=False,
gainfield=gainfield,
interp=interp,
gaintable=gaintables,
calwt=calwt,
vis=ms,
applymode='calonly',
antenna='*&*',
spwmap=spwmap,
parang=True)
myclean(
vis=ms,
name=name + '_taper',
imsize=2000,
cell='0.04arcsec',
fields=['Sgr B2 N Q', 'Sgr B2 NM Q', 'Sgr B2 MS Q'],
uvrange='0~2000klambda',
threshold='5mJy',
)