Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Execution Block ID uid://A001/X15a0/Xb2 Sgr_A_st_d_03_TM1 #247

Open
6 of 15 tasks
keflavich opened this issue Sep 12, 2022 · 29 comments
Open
6 of 15 tasks

Execution Block ID uid://A001/X15a0/Xb2 Sgr_A_st_d_03_TM1 #247

keflavich opened this issue Sep 12, 2022 · 29 comments
Assignees
Labels
Delivered EB Execution Block Needs Reimaging: Spectral Resolution size mitig Needs to be reimaged with full spectral resolution TM1

Comments

@keflavich
Copy link
Contributor

keflavich commented Sep 12, 2022

Sgr_A_st_d_03_TM1
uid://A001/X15a0/Xb2

Product Links:

Reprocessed Product Links:

@keflavich
Copy link
Contributor Author

Had a very hard crash in the reimaging step of the pipeline run:

2022-10-19 09:05:33 INFO: Executing tclean(vis=['uid___A002_Xfe3986_X9083_target.ms', 'uid___A002_Xfe62c1_X1871_target.ms'], field='Sgr_A_star', spw=['25:85.9547018882~86.0120749351GHz;86.0919089194~86.1590475913GHz;86.3697409507~86.4158835288GHz', '25:85.9547088863~86.0120819332GHz;86.0919159176~86.1590545894GHz;86.3697479488~86.4158905269GHz'], antenna=['0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43&', '0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43&'], scan=['6,9,12,15,18', '6,9,12,15,18'], intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', imagename='uid___A001_X15a0_Xb2.s8_0.Sgr_A_star_sci.spw25.mfs.I.iter1', imsize=[1728, 2160], cell=['0.23arcsec'], phasecenter='ICRS 17:46:45.9499 -028.42.46.695', stokes='I', specmode='mfs', nchan=-1, outframe='LSRK', perchanweightdensity=False, gridder='mosaic', mosweight=True, usepointing=False, pblimit=0.2, deconvolver='hogbom', restoration=True, restoringbeam='common', pbcor=True, weighting='briggs', robust=0.5, npixels=0, niter=20000000, threshold='0.000697Jy', nsigma=0.0, interactive=0, usemask='auto-multithresh', sidelobethreshold=2.0, noisethreshold=4.25, lownoisethreshold=1.5, negativethreshold=0.0, minbeamfrac=0.3, growiterations=75, dogrowprune=True, minpercentchange=1.0, fastnoise=False, restart=True, savemodel='none', calcres=False, calcpsf=False, parallel=True)
*** Error in `/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3': double free or corruption (fasttop): 0x00002b4ad8000b20 ***
======= Backtrace: =========
/lib64/libc.so.6(+0x816b9)[0x2b4a6bdb26b9]
/lib64/libstdc++.so.6(_ZNSs6assignERKSs+0x9e)[0x2b4a81b6bdbe]
/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/casatools/__casac__/lib/libcasatools.cpython-36m-x86_64-linux-gnu.so(_ZN9casa6core9LogOrigin9classNameERKNS_6StringE+0xd)[0x2b4a8783ef1d]
/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/casatools/__casac__/_logsink.cpython-36m-x86_64-linux-gnu.so(_ZN5casac7logsink11postLocallyERKSsS2_S2_+0x6b)[0x2b4a95c0cceb]
/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/casatools/__casac__/_logsink.cpython-36m-x86_64-linux-gnu.so(+0x1344d)[0x2b4a95c1544d]
/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3(PyCFunction_Call+0xc6)[0x4f35e6]
/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3(_PyEval_EvalFrameDefault+0x59df)[0x54a17f]

...
7ffcc40d9000-7ffcc412a000 rw-p 00000000 00:00 0                          [stack]
7ffcc4180000-7ffcc4182000 r-xp 00000000 00:00 0                          [vdso]
ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0                  [vsyscall]
[c0703a-s11:42337] *** Process received signal ***
[c0703a-s11:42337] Signal: Aborted (6)
[c0703a-s11:42337] Signal code:  (-6)
[c0703a-s11:42337] [ 0] /lib64/libpthread.so.0(+0xf630)[0x2b4a6b213630]
[c0703a-s11:42337] [ 1] /lib64/libc.so.6(gsignal+0x37)[0x2b4a6bd67377]
[c0703a-s11:42337] [ 2] /lib64/libc.so.6(abort+0x148)[0x2b4a6bd68a68]
[c0703a-s11:42337] [ 3] /lib64/libc.so.6(+0x78ec7)[0x2b4a6bda9ec7]
[c0703a-s11:42337] [ 4] /lib64/libc.so.6(+0x816b9)[0x2b4a6bdb26b9]
[c0703a-s11:42337] [ 5] /lib64/libstdc++.so.6(_ZNSs6assignERKSs+0x9e)[0x2b4a81b6bdbe]
[c0703a-s11:42337] [ 6] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/casatools/__casac__/lib/libcasatools.cpython-36m-x86_64-linux-gnu.so(_ZN9casa6core9LogOrigin9classNameERKNS_6StringE+0xd)[0x2b4a8783ef1d]
[c0703a-s11:42337] [ 7] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/casatools/__casac__/_logsink.cpython-36m-x86_64-linux-gnu.so(_ZN5casac7logsink11postLocallyERKSsS2_S2_+0x6b)[0x2b4a95c0cceb]
[c0703a-s11:42337] [ 8] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/casatools/__casac__/_logsink.cpython-36m-x86_64-linux-gnu.so(+0x1344d)[0x2b4a95c1544d]
[c0703a-s11:42337] [ 9] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3(PyCFunction_Call+0xc6)[0x4f35e6]
[c0703a-s11:42337] [10] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3(_PyEval_EvalFrameDefault+0x59df)[0x54a17f]
[c0703a-s11:42337] [11] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3[0x543cb7]
[c0703a-s11:42337] [12] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3[0x54d6eb]
[c0703a-s11:42337] [13] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3[0x54c61c]
[c0703a-s11:42337] [14] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3(_PyEval_EvalFrameDefault+0x2ca)[0x544a6a]
[c0703a-s11:42337] [15] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3[0x543cb7]
[c0703a-s11:42337] [16] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3[0x54d6eb]
[c0703a-s11:42337] [17] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3[0x54c61c]
[c0703a-s11:42337] [18] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3(_PyEval_EvalFrameDefault+0x2ca)[0x544a6a]
[c0703a-s11:42337] [19] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3[0x543cb7]
[c0703a-s11:42337] [20] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3[0x54d6eb]
[c0703a-s11:42337] [21] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3[0x54c61c]
[c0703a-s11:42337] [22] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3(_PyEval_EvalFrameDefault+0x1244)[0x5459e4]
[c0703a-s11:42337] [23] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3[0x54d62d]
[c0703a-s11:42337] [24] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3[0x54c61c]
[c0703a-s11:42337] [25] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3(_PyEval_EvalFrameDefault+0x2ca)[0x544a6a]
[c0703a-s11:42337] [26] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3[0x543cb7]
[c0703a-s11:42337] [27] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3[0x54d6eb]
[c0703a-s11:42337] [28] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3[0x54c61c]
[c0703a-s11:42337] [29] /blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/bin/python3(_PyEval_EvalFrameDefault+0x2ca)[0x544a6a]
[c0703a-s11:42337] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 42337 on node c0703a-s11 exited on signal 6 (Aborted).
--------------------------------------------------------------------------

@keflavich
Copy link
Contributor Author

These files, from nov 9, are uncalibrated:

 ls -dlhrt member.uid___A001_X15a0_Xb2/calibrated/working/*ms
drwxr-sr-x 28 adamginsburg adamginsburg 4.0K Nov  9 00:35 member.uid___A001_X15a0_Xb2/calibrated/working/uid___A002_Xfe3986_X9083.ms
drwxr-sr-x 28 adamginsburg adamginsburg 4.0K Nov  9 00:58 member.uid___A001_X15a0_Xb2/calibrated/working/uid___A002_Xfe3986_X9083_target.ms
drwxr-sr-x 28 adamginsburg adamginsburg 4.0K Nov  9 01:45 member.uid___A001_X15a0_Xb2/calibrated/working/uid___A002_Xfe62c1_X1871.ms

There were no errors in the pipeline log.

The imaging pipeline log has:

2022-11-09 05:59:01     INFO    mstransform::::casa     ##########################################
2022-11-09 05:59:01     INFO    mstransform::::casa     ##### Begin Task: mstransform        #####
2022-11-09 05:59:01     INFO    mstransform::::casa     mstransform( vis='uid___A002_Xfe62c1_X1871.ms', outputvis='uid___A002_Xfe62c1_X1871_target.ms', createmms=False, separationaxis='auto', numsubms='auto', tileshape=[0], field='3,4,5,6,
7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87
,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121', spw='25,27,29,31,33,35', scan='', antenna='', correlation='', timerange='', intent='OBSERVE_TARGET#ON_SOURCE', a
rray='', uvrange='', observation='', feed='', datacolumn='corrected', realmodelcol=False, keepflags=True, usewtspectrum=False, combinespws=False, chanaverage=False, chanbin=1, hanning=False, regridms=False, mode='channel', nchan=-1, start=
0, width=1, nspw=1, interpolation='linear', phasecenter='', restfreq='', outframe='', veltype='radio', preaverage=False, timeaverage=False, timebin='0s', timespan='', maxuvwdistance=0.0, docallib=False, callib='', douvcontsub=False, fitspw
='', fitorder=0, want_cont=False, denoising_lib=True, nthreads=1, niter=1, disableparallel=False, ddistart=-1, taql='', monolithic_processing=False, reindex=False )
2022-11-09 05:59:01     INFO    MSTransformManager::parseMsSpecParams   Input file name is uid___A002_Xfe62c1_X1871.ms
2022-11-09 05:59:01     INFO    MSTransformManager::parseMsSpecParams   Data column is CORRECTED
2022-11-09 05:59:01     INFO    MSTransformManager::parseMsSpecParams   Output file name is uid___A002_Xfe62c1_X1871_target.ms
2022-11-09 05:59:01     INFO    MSTransformManager::parseMsSpecParams   Re-index is disabled
2022-11-09 05:59:01     INFO    MSTransformManager::parseMsSpecParams   Tile shape is [0]
2022-11-09 05:59:01     INFO    MSTransformManager::parseDataSelParams  field selection is 3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54
,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121
2022-11-09 05:59:01     INFO    MSTransformManager::parseDataSelParams  spw selection is 25,27,29,31,33,35
2022-11-09 05:59:01     INFO    MSTransformManager::parseDataSelParams  scan intent selection is OBSERVE_TARGET#ON_SOURCE
2022-11-09 05:59:01     WARN    MSTransformManager::checkDataColumnsToFill      CORRECTED_DATA column requested but not available in input MS
2022-11-09 05:59:01     INFO    MSTransformManager::initDataSelectionParams     Selected SPWs Ids are Axis Lengths: [6, 4]  (NB: Matrix in Row/Column order)
2022-11-09 05:59:01     INFO    MSTransformManager::initDataSelectionParams+    [25, 0, 1919, 1
2022-11-09 05:59:01     INFO    MSTransformManager::initDataSelectionParams+     27, 0, 1919, 1
2022-11-09 05:59:01     INFO    MSTransformManager::initDataSelectionParams+     29, 0, 1919, 1
2022-11-09 05:59:01     INFO    MSTransformManager::initDataSelectionParams+     31, 0, 1919, 1
2022-11-09 05:59:01     INFO    MSTransformManager::initDataSelectionParams+     33, 0, 3839, 1
2022-11-09 05:59:01     INFO    MSTransformManager::initDataSelectionParams+     35, 0, 3839, 1]
2022-11-09 05:59:01     INFO    MSTransformManager::open        Select data
2022-11-09 05:59:01     INFO    MSTransformManager::createOutputMSStructure     Create output MS structure
2022-11-09 05:59:02     SEVERE  mstransform::::casa     Task mstransform raised an exception of class RuntimeError with the following message: Desired column (CORRECTED_DATA) not found in the input MS (/orange/adamginsburg/ACES/rawdata/202
1.1.00172.L/science_goal.uid___A001_X1590_X30a8/group.uid___A001_X1590_X30a9/member.uid___A001_X15a0_Xb2/calibrated/working/uid___A002_Xfe62c1_X1871.ms).
2022-11-09 05:59:02     INFO    mstransform::::casa     Task mstransform complete. Start time: 2022-11-09 00:59:01.119776 End time: 2022-11-09 00:59:01.634003
2022-11-09 05:59:02     INFO    mstransform::::casa     ##### End Task: mstransform          #####
2022-11-09 05:59:02     INFO    mstransform::::casa     ##########################################

@keflavich
Copy link
Contributor Author

@keflavich
Copy link
Contributor Author

The MSes have both data & corrected data columns that are identical, with values in the range 0.00211 to 0.00275

@keflavich
Copy link
Contributor Author

I've deleted the calibrated/ folder and will try again.

@keflavich
Copy link
Contributor Author

#134 (comment), #134 (comment), #134 (comment) all belong here.

Re: @pyhsiehATalma you are correct, that log file refers to Xb2 and all of the related files belong here. I must have had two tabs open.

@pyhsiehATalma
Copy link
Contributor

Original product shows divergence in spw 33.

uid___A001_X15a0_Xb2 s38_0 Sgr_A_star_sci spw33 cube I iter1 image I mom8 sky

The divergence seems to be resolved with newer pipeline run. I check the cube in /calibrated/working/Sgr_A_star_sci.spw33.cube I.iter1.image.pbcor.
The mom 8 map of new cube is shown below in the right panel.

So I think there is no need to reclean.
uid___A001_X15a0_Xb2 s38_0 Sgr_A_star_sci spw33 cube I iter1 image pbcor fits-uid___A001_X15a0_Xb2 s38_0 Sgr_A_star_sci spw33 cube I iter1 image pbcor fits moment maximum-image-2022-12-29-18-36-49

@keflavich
Copy link
Contributor Author

terminal failure of pipeline was:

2022-11-21 05:34:30 INFO: Executing mstransform(vis='uid___A002_Xfe3986_X9083.ms', outputvis='uid___A002_Xfe3986_X9083_target.ms', field='3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63
,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121', spw='25,27,29,31,33,35', intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='corrected', chanaverage=False, chanbin=1,
timeaverage=False, timebin='0s', reindex=False)
2022-11-21 05:34:31     WARN    MSTransformManager::checkDataColumnsToFill      CORRECTED_DATA column requested but not available in input MS
2022-11-21 05:34:31     SEVERE  mstransform::::casa     Task mstransform raised an exception of class RuntimeError with the following message: Desired column (CORRECTED_DATA) not found in the input MS (/orange/adamginsburg/ACES/rawdata/2021.1.00172.L/science_goal.uid___A001_X1590_X30a8/group.uid___A001_X1590_X30
a9/member.uid___A001_X15a0_Xb2/calibrated/working/uid___A002_Xfe3986_X9083.ms).
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa     Error executing pipeline task hif_mstransform.
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa     Traceback (most recent call last):
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+      File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 640, in execute
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+        return self._handle_multiple_vis(dry_run, **parameters)
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+      File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 762, in _handle_multiple_vis
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+        single_result = self.execute(dry_run=dry_run, **parameters)
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+      File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 40, in attach_timestamp_to_results
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+        result = method(self, *args, **kw)
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+      File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 120, in handle_matplotlibrc
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+        result = method(self, *args, **kwargs)
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+      File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 88, in capture
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+        result = method(self, *args, **kw)
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+      File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 61, in finalise_pipeline_result
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+        result = method(self, *args, **kw)
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+      File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 660, in execute
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+        result = self.prepare(**prepare_parameters)
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+      File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/hif/tasks/mstransform/mstransform.py", line 167, in prepare
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+        self._executor.execute(mstransform_job)
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+      File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 88, in capture
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+        result = method(self, *args, **kw)
2022-11-21 05:34:31     SEVERE  mstransform::pipeline.infrastructure.basetask::casa+      File "/blue/adaESC[31m2022-11-21 05:34:30 ERROR: Error executing pipeline task hif_mstransform.ESC[0m
ESC[31m2022-11-21 05:34:30 ERROR: Traceback (most recent call last):ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 640, in executeESC[0m
ESC[31m    return self._handle_multiple_vis(dry_run, **parameters)ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 762, in _handle_multiple_visESC[0m
ESC[31m    single_result = self.execute(dry_run=dry_run, **parameters)ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 40, in attach_timestamp_to_resultsESC[0m
ESC[31m    result = method(self, *args, **kw)ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 120, in handle_matplotlibrcESC[0m
ESC[31m    result = method(self, *args, **kwargs)ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 88, in captureESC[0m
ESC[31m    result = method(self, *args, **kw)ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 61, in finalise_pipeline_resultESC[0m
ESC[31m    result = method(self, *args, **kw)ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 660, in executeESC[0m
ESC[31m    result = self.prepare(**prepare_parameters)ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/hif/tasks/mstransform/mstransform.py", line 167, in prepareESC[0m
ESC[31m    self._executor.execute(mstransform_job)ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 88, in captureESC[0m
ESC[31m    result = method(self, *args, **kw)ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 818, in executeESC[0m
ESC[31m    result = job.execute(dry_run=self._dry_run, **kwargs)ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/jobrequest.py", line 243, in executeESC[0m
ESC[31m    return self.fn(*self.args, **self.kw)ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/casatasks/mstransform.py", line 620, in __call__ESC[0m
ESC[31m    task_result = _mstransform_t( _pc.document['vis'], _pc.document['outputvis'], _pc.document['createmms'], _pc.document['separationaxis'], _pc.document['numsubms'], _pc.document['tileshape'], _pc.document['field'], _pc.document['spw'], _pc.document['scan'], _pc.document['antenna'], _pc.document['correlation'], _pc.document['timerange'], _pc.document['intent'], _pc.document['array'], _pc.document['uvrange'], _pc.document['observation'], _pc.document['feed'], _pc.document['datacolumn'], _pc.document['realmodelcol'], _pc.document['keepflags'], _pc.document['usewtspectrum'], _pc.document['combinespws'], _pc.document['chanaverage'], _pc.document['chanbin'], _pc.document['hanning'], _pc.document['regridms'], _pc.document['mode'], _pc.document['nchan'], _pc.document['start'], _pc.document['width'], _pc.document['nspw'], _pc.document['interpolation'], _pc.document['phasecenter'], _pc.document['restfreq'], _pc.document['outframe'], _pc.document['veltype'], _pc.document['preaverage'], _pc.document['timeaverage'], _pc.document['timebin'], _pc.document['timespan'], _pc.document['maxuvwdistance'], _pc.document['docallib'], _pc.document['callib'], _pc.document['douvcontsub'], _pc.document['fitspw'], _pc.document['fitorder'], _pc.document['want_cont'], _pc.document['denoising_lib'], _pc.document['nthreads'], _pc.document['niter'], _pc.document['disableparallel'], _pc.document['ddistart'], _pc.document['taql'], _pc.document['monolithic_processing'], _pc.document['reindex'] )ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/casatasks/private/task_mstransform.py", line 307, in mstransformESC[0m
ESC[31m    mtlocal.open()ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/casatools/mstransformer.py", line 44, in openESC[0m
ESC[31m    return self._swigobj.open()ESC[0m
ESC[31m  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/casatools/__casac__/mstransformer.py", line 176, in openESC[0m
ESC[31m    return _mstransformer.mstransformer_open(self)ESC[0m
ESC[31mRuntimeError: Desired column (CORRECTED_DATA) not found in the input MS (/orange/adamginsburg/ACES/rawdata/2021.1.00172.L/science_goal.uid___A001_X1590_X30a8/group.uid___A001_X1590_X30a9/member.uid___A001_X15a0_Xb2/calibrated/working/uid___A002_Xfe3986_X9083.ms).ESC[0m
ESC[31mESC[0m
mginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 818, in execute

so, basically same as above.

@keflavich
Copy link
Contributor Author

In the pipeline rerun, this has failed again in step 2:
image

@keflavich
Copy link
Contributor Author

In the latest run, CORRECTED_DATA exists:

CASA <3>: ms.open('uid___A002_Xfe3986_X9083.ms')
CASA <9>: ms.selectinit(0)
CASA <10>: ms.getdata('CORRECTED_DATA')['corrected_data'].shape
Out[10]: (2, 1, 110880)
CASA <11>: ms.open('uid___A002_Xfe62c1_X1871.ms')
Out[11]: True

CASA <12>: ms.selectinit(0)
Out[12]: True

CASA <13>: ms.getdata('CORRECTED_DATA')['corrected_data'].shape
Out[13]: (2, 1, 110880)

so... hm.

@keflavich
Copy link
Contributor Author

@ashleythomasbarnes and I are trying to image the data to see if they make sense. The data seem to have reasonable amplitudes, maybe?

CASA <37>: ms.open('/orange/adamginsburg/ACES/rawdata/2021.1.00172.L/science_goal.uid___A001_X1590_X30a8/group.uid___A001_X1590_X30a9/member.uid___A001_X15a0_Xb2/calibrated/uid___A002_Xfe62c1_X1871.ms')
Out[37]: True

CASA <38>: ms.selectinit(0)
Out[38]: True

CASA <39>: data = ms.getdata(['amplitude', 'corrected_amplitude'])

CASA <40>: data['amplitude'].mean()
Out[40]: 0.0023897540722077433

@keflavich
Copy link
Contributor Author

These cube images are still junk. I'm trying to run this pseudo-manually to see if there's any more useful information in an 'interactive session' run (as compared to a batch job run in the background)

@keflavich
Copy link
Contributor Author

Manual pipeline run failed:

2023-02-08 23:59:17 ERROR: Error executing pipeline task hifa_restoredata.2023-02-08 23:59:17     SEVERE  flagmanager::pipeline.infrastructure.basetask::casa     Error executing pipeline task hifa_restoredata.
2023-02-08 23:59:17 ERROR: Traceback (most recent call last):
  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 660, in execute    result = self.prepare(**prepare_parameters)
  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/h/tasks/restoredata/restoredata.py", line 262, in prepare
    self._do_restore_calstate(pipemanifest)  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/h/tasks/restoredata/restoredata.py", line 488, in _do_restore_calstate
    append=append)
  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/callibrary.py", line 2180, in import_state
    calapp = CalApplication.from_export(line)
  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/callibrary.py", line 102, in from_export
    ''.format(os.path.gaintable))
AttributeError: module 'posixpath' has no attribute 'gaintable'

?!?!?!

@keflavich
Copy link
Contributor Author

I think that's a straight up bug in CASA; the code is a clear typo. How did I reach that corner case?

@keflavich
Copy link
Contributor Author

I hack-fixed that bug, but it crashed a few lines lower with the real error:

2023-02-09 04:19:58 INFO: Importing calibration state from /scratch/local/56787088/tmpusxjcgls
2023-02-09 04:19:58 WARNING: Could not access uid___A002_Xfe3986_X9083.ms.hif_uvcontfit.s5_1.Sgr_A_star.uvcont.tbl. Using heuristics to determine caltable type
2023-02-09 04:19:59     WARN    flagmanager::pipeline.infrastructure.callibrary::casa   Could not access uid___A002_Xfe3986_X9083.ms.hif_uvcontfit.s5_1.Sgr_A_star.uvcont.tbl. Using heuristics to determine caltable type
2023-02-09 04:19:59 ERROR: Error executing pipeline task hifa_restoredata.
2023-02-09 04:19:59     SEVERE  flagmanager::pipeline.infrastructure.basetask::casa     Error executing pipeline task hifa_restoredata.
2023-02-09 04:19:59 ERROR: Traceback (most recent call last):
  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/basetask.py", line 660, in execute
    result = self.prepare(**prepare_parameters)
  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/h/tasks/restoredata/restoredata.py", line 262, in prepare
    self._do_restore_calstate(pipemanifest)
  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/h/tasks/restoredata/restoredata.py", line 488, in _do_restore_calstate
    append=append)
  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/callibrary.py", line 2180, in import_state
    calapp = CalApplication.from_export(line)
  File "/blue/adamginsburg/adamginsburg/casa/casa-6.4.3-2-pipeline-2021.3.0.17/lib/py/lib/python3.6/site-packages/Pipeline-2021.3.0.17-py3.6.egg/pipeline/infrastructure/callibrary.py", line 114, in from_export
    raise ValueError(gaintable)
ValueError: uid___A002_Xfe3986_X9083.ms.hif_uvcontfit.s5_1.Sgr_A_star.uvcont.tbl

i.e., this is not a recognized calibration table type.

I'm going to delete the entire directory structure and try again.

@claire-cook claire-cook self-assigned this Feb 10, 2023
@claire-cook claire-cook added Needs Reimaging: Spectral Resolution size mitig Needs to be reimaged with full spectral resolution Needs Reimaging: Divergence Needs reimaging b/c of divergence labels Feb 10, 2023
@claire-cook claire-cook removed their assignment Feb 10, 2023
@claire-cook claire-cook self-assigned this Feb 10, 2023
@keflavich
Copy link
Contributor Author

Still junk. Needs the full delete-everything treatment.

It looks like the ASDM was only 2GB! That was never right.

$ ls -lh *_Xb2_* *X484c*
-rw-r--r-- 1 adamginsburg adamginsburg 167G Feb  9 12:37 2021.1.00172.L_uid___A001_X15a0_Xb2_001_of_001.tar
-rw-r--r-- 1 adamginsburg adamginsburg 654M Feb  9 12:38 2021.1.00172.L_uid___A001_X15a0_Xb2_auxiliary.tar
-rw-r--r-- 1 adamginsburg adamginsburg 2.0G Jul  5  2022 2021.1.00172.L_uid___A002_Xf4df6f_X484c.asdm.sdm.tar
$ md5sum *_Xb2_* *X484c*
9865e5801a0c44d7b6ebb78aa77ab501  2021.1.00172.L_uid___A001_X15a0_Xb2_001_of_001.tar
974ff1bbb73ae5e8b17cfc41ab6d28fc  2021.1.00172.L_uid___A001_X15a0_Xb2_auxiliary.tar
c3abf65fb6b5b6224cba38d0cf9eefae  2021.1.00172.L_uid___A002_Xf4df6f_X484c.asdm.sdm.tar

@keflavich
Copy link
Contributor Author

Cubes are good now! Continuum is junk. Needs the same treatment as x

@ashleythomasbarnes
Copy link
Collaborator

Any update @claire-cook?

@claire-cook
Copy link
Contributor

Any update @claire-cook?

I'm actively running tclean on the entire divergent SPW 33 cube on a KU machine! I raised the cyclefactor from 2 to 3.0 to attempt to fix the divergence, so it's taking a while.

@claire-cook
Copy link
Contributor

At a cyclefactor of 2.0, my tclean run of SPW 33 remained divergent solely in channel 615. After increasing the cyclefactor to 3.0, this divergence is resolved. The non-divergent channel 615 in the recleaned SPW 33 cube:
member uid_A001_x15a0_xb2 SPW 33 channel 615 RECLEAN NO DIVERGENCE

There should be no need for further recleaning.

@pyhsiehATalma
Copy link
Contributor

@claire-cook thanks for the re-imaging!

could you please upload the new cubes to globus under /upload? Then I could copy them for mosaic.
let me know if you have any questions.

https://github.com/ACES-CMZ/reduction_ACES/wiki/Using-Globus-to-transfer-ACES-data-(Linux)

@keflavich
Copy link
Contributor Author

I moved the uploaded files into /orange/adamginsburg/ACES/rawdata/2021.1.00172.L/science_goal.uid___A001_X1590_X30a8/group.uid___A001_X1590_X30a9/member.uid___A001_X15a0_Xb2/reclean/

@claire-cook claire-cook removed Needs Reimaging: Divergence Needs reimaging b/c of divergence Needs Reimaging: Spectral Resolution size mitig Needs to be reimaged with full spectral resolution labels Jun 22, 2023
@d-l-walker
Copy link
Contributor

To-do: check whether the products are still size mitigated. @claire-cook can you look into this? I see that you removed the label for size mitigation, but I don't see any updated parameters in override_teclean_commands.json.

(Tagging #179 for tracking)

@claire-cook
Copy link
Contributor

To-do: check whether the products are still size mitigated. @claire-cook can you look into this? I see that you removed the label for size mitigation, but I don't see any updated parameters in override_teclean_commands.json.

(Tagging #179 for tracking)

The recleaned product for the previously-divergent SPW 33 should no longer contain size mitigation. It appears that raising the cyclefactor to 3.0 not only fixed the divergence, but also seemed to fix the size mitigation (see number of channels in the following screenshots). I can check if there's size mitigation in the reprocessed products of the other SPWs.

Original (size mitigated) product for SPW 33:
Original Xb2 SPW 33 size mitigation

Recleaned SPW 33 product without size mitigation:
Reprocessed Xb2 SPW 33 no size mitigation

The only updated tclean parameter used was a raised cyclefactor of 3.0
Updated tclean parameters Xb2 SPW 33

@claire-cook
Copy link
Contributor

Checking other SPWs for size mitigation: spectral resolution size mitigation is still present for everything but SPW 33.
Recleaning SPW 33 resolved divergence and size mitigation, but (though I haven't downloaded everything) the other SPWs appear to also have size mitigation in both the original and reprocessed products:
Here's the spectral resolution/number of channels for the original and the reprocessed SPW 25 cubes, both with the same spectral resolution size mitigation:
Original Xb2 SPW 25 size mitgated (1x resolution + channel number)
Reprocessed Xb2 SPW 25 size mitigated (1x resolution + channel number)

While I haven't physically downloaded the other cubes, a quick check shows that the file sizes of the original vs. reprocessed products on Globus are identical, indicating that size mitigation is still present for those as well. When size mitigation is solved, we'd expect the unmitigated cube to be ~twice the physical file size of the mitigated one.

@claire-cook claire-cook added the Needs Reimaging: Spectral Resolution size mitig Needs to be reimaged with full spectral resolution label Dec 6, 2023
d-l-walker added a commit to d-l-walker/reduction_ACES that referenced this issue Dec 13, 2023
keflavich pushed a commit that referenced this issue Dec 13, 2023
* Undo size mitigation for aa (#251)

* Undo size mitigation for l (#137)

* Undo size mitigation for r (#263)

* Undo size mitigation for d (#247)

* Undo size mitigation for g (#134)
@keflavich
Copy link
Contributor Author

Moved files to get spw35 unmitigated.
/orange/adamginsburg/ACES/data/2021.1.00172.L/science_goal.uid___A001_X1590_X30a8/group.uid___A001_X1590_X30a9/member.uid___A001_X15a0_Xb2/calibrated/working$ mv *spw35.cube.* spw35_mitigated/

@d-l-walker
Copy link
Contributor

QA for latest continuum selection/imaging:

  • Note: new SPW 33 & 35 image missing
  • Note: apparently the individual SPW MFS images do not exist for this region
  • Possible signs of contamination in the lower frequency windows. I'll check the cubes.

region_d_continuum

@d-l-walker
Copy link
Contributor

Reminder that SPWs 25,27,29,31 all need to be un-size-mitigated and re-stat-cont-ed @keflavich

@keflavich
Copy link
Contributor Author

Moved files to sizemitigated/: mv *spw{25,27,29,31}* sizemitigated/. Restart forthcoming

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Delivered EB Execution Block Needs Reimaging: Spectral Resolution size mitig Needs to be reimaged with full spectral resolution TM1
Projects
None yet
Development

No branches or pull requests

5 participants