Commits
Ville Suoranta authored 8baa7e0646d Merge
2312 2312 | def test_widefield_wbaproj_mfs(self): |
2313 2313 | """ [widefield] Test_Widefield_wbaproj_mfs : MFS with wideband AWProjection (wbawp=T, allspw) and nt=1 stokes I """ |
2314 2314 | |
2315 2315 | # casalog.post("EMPTY TEST") |
2316 2316 | # return |
2317 2317 | |
2318 2318 | self.prepData("refim_mawproject.ms") |
2319 2319 | ret = tclean(vis=self.msfile,field='*',imagename=self.img,imsize=512,cell='10.0arcsec',phasecenter="J2000 19:59:28.500 +40.44.01.50", |
2320 2320 | niter=30,gridder='awproject',cfcache=self.img+'.cfcache',wbawp=True,conjbeams=True,psterm=False,computepastep=360.0, |
2321 2321 | rotatepastep=360.0,deconvolver='hogbom',pblimit=0.3,parallel=self.parallel) |
2322 - | # report=self.th.checkall(imexist=[self.img+'.image', self.img+'.psf', self.img+'.weight'],imval=[(self.img+'.image',1.0,[256,256,0,0]),(self.img+'.weight',0.493,[256,256,0,0]) ] ) |
2322 + | report=self.th.checkall(imexist=[self.img+'.image', self.img+'.psf', self.img+'.weight'],imval=[(self.img+'.image',1.0,[256,256,0,0]),(self.img+'.weight',0.493,[256,256,0,0]) ] ) |
2323 2323 | # |
2324 2324 | # Changed to the following for 5.5.0 release of AWP. Will revisit and replace the test MS later. |
2325 2325 | # |
2326 - | report=self.th.checkall(imexist=[self.img+'.image', self.img+'.psf', self.img+'.weight'],imval=[(self.img+'.image',0.698,[256,256,0,0]),(self.img+'.weight',0.493,[256,256,0,0]) ] ) |
2326 + | #report=self.th.checkall(imexist=[self.img+'.image', self.img+'.psf', self.img+'.weight'],imval=[(self.img+'.image',0.698,[256,256,0,0]),(self.img+'.weight',0.493,[256,256,0,0]) ] ) |
2327 2327 | self.checkfinal(report) |
2328 2328 | |
2329 2329 | |
2330 2330 | #do stokes V too.. |
2331 2331 | |
2332 2332 | def test_widefield_aproj_mtmfs(self): |
2333 2333 | """ [widefield] Test_Widefield_aproj_mtmfs : MFS with AWProjection (wbawp=T,conjbeams=F, allspw) and nt=2 stokes I """ |
2334 2334 | |
2335 2335 | # casalog.post("EMPTY TEST") |
2336 2336 | # return |
2337 2337 | |
2338 2338 | self.prepData("refim_mawproject.ms") |
2339 2339 | ret = tclean(vis=self.msfile,spw='*',field='*',imagename=self.img,imsize=512,cell='10.0arcsec',phasecenter="J2000 19:59:28.500 +40.44.01.50", |
2340 2340 | niter=0,gridder='awproject',cfcache=self.img+'.cfcache',wbawp=True,conjbeams=False,psterm=False,computepastep=360.0, |
2341 2341 | rotatepastep=360.0,deconvolver='mtmfs',parallel=self.parallel) |
2342 - | # report=self.th.checkall(imexist=[self.img+'.image.tt0', self.img+'.psf.tt0', self.img+'.weight.tt0'],imval=[(self.img+'.image.tt0',0.96,[256,256,0,0]),(self.img+'.weight.tt0',0.48,[256,256,0,0]),(self.img+'.alpha',0.04,[256,256,0,0]) ] ) |
2342 + | report=self.th.checkall(imexist=[self.img+'.image.tt0', self.img+'.psf.tt0', self.img+'.weight.tt0'],imval=[(self.img+'.image.tt0',0.96,[256,256,0,0]),(self.img+'.weight.tt0',0.48,[256,256,0,0]),(self.img+'.alpha',0.04,[256,256,0,0]) ] ) |
2343 2343 | # |
2344 2344 | # Changed to the following for 5.5.0 release of AWP. Will revisit and replace the test MS later. |
2345 2345 | # |
2346 - | report=self.th.checkall(imexist=[self.img+'.image.tt0', self.img+'.psf.tt0', self.img+'.weight.tt0'],imval=[(self.img+'.image.tt0',0.705,[256,256,0,0]),(self.img+'.weight.tt0',0.48,[256,256,0,0]),(self.img+'.alpha',0.04,[256,256,0,0]) ] ) |
2346 + | #report=self.th.checkall(imexist=[self.img+'.image.tt0', self.img+'.psf.tt0', self.img+'.weight.tt0'],imval=[(self.img+'.image.tt0',0.705,[256,256,0,0]),(self.img+'.weight.tt0',0.48,[256,256,0,0]),(self.img+'.alpha',0.04,[256,256,0,0]) ] ) |
2347 2347 | self.checkfinal(report) |
2348 2348 | ## alpha should represent that of the mosaic PB (twice).. -0.1 doesn't look right. Sigh.... well.. it should converge to zero. |
2349 2349 | ## alpha keeps increasing in magnitude with niter.... not right. |
2350 2350 | ## restricting this check to niter=0 only. Alpha should represent that of the PB. |
2351 2351 | |
2352 2352 | def test_widefield_wbaproj_mtmfs(self): |
2353 2353 | """ [widefield] Test_Widefield_wbaproj_mtmfs : MFS with wideband AWProjection (wbawp=T,conjbeams=T, allspw) and nt=2 stokes I """ |
2354 2354 | |
2355 2355 | # casalog.post("EMPTY TEST") |
2356 2356 | # return |
2357 2357 | |
2358 2358 | self.prepData("refim_mawproject.ms") |
2359 2359 | ret = tclean(vis=self.msfile,field='*',imagename=self.img,imsize=512,cell='10.0arcsec',phasecenter="J2000 19:59:28.500 +40.44.01.50", |
2360 2360 | niter=30,gridder='awproject',cfcache=self.img+'.cfcache',wbawp=True,conjbeams=True,psterm=False, |
2361 2361 | computepastep=360.0,rotatepastep=360.0,deconvolver='mtmfs',pblimit=0.1,parallel=self.parallel) |
2362 - | # report=self.th.checkall(imexist=[self.img+'.image.tt0', self.img+'.psf.tt0', self.img+'.weight.tt0'],imval=[(self.img+'.image.tt0',0.96,[256,256,0,0]),(self.img+'.weight.tt0',0.486,[256,256,0,0]),(self.img+'.alpha',0.0,[256,256,0,0]) ] ) |
2362 + | report=self.th.checkall(imexist=[self.img+'.image.tt0', self.img+'.psf.tt0', self.img+'.weight.tt0'],imval=[(self.img+'.image.tt0',0.96,[256,256,0,0]),(self.img+'.weight.tt0',0.486,[256,256,0,0]),(self.img+'.alpha',0.0,[256,256,0,0]) ] ) |
2363 2363 | # |
2364 2364 | # Changed to the following for 5.5.0 release of AWP. Will revisit and replace the test MS later. |
2365 2365 | # |
2366 - | report=self.th.checkall(imexist=[self.img+'.image.tt0', self.img+'.psf.tt0', self.img+'.weight.tt0'],imval=[(self.img+'.image.tt0',0.696,[256,256,0,0]),(self.img+'.weight.tt0',0.486,[256,256,0,0]),(self.img+'.alpha',0.0,[256,256,0,0]) ] ) |
2366 + | #report=self.th.checkall(imexist=[self.img+'.image.tt0', self.img+'.psf.tt0', self.img+'.weight.tt0'],imval=[(self.img+'.image.tt0',0.696,[256,256,0,0]),(self.img+'.weight.tt0',0.486,[256,256,0,0]),(self.img+'.alpha',0.0,[256,256,0,0]) ] ) |
2367 2367 | ## alpha should be ZERO as the pb spectrum has been taken out. |
2368 2368 | self.checkfinal(report) |
2369 2369 | |
2370 2370 | |
2371 2371 | # def test_widefield_wbaproj_subsets(self): |
2372 2372 | # """ [widefield] Test_Widefield_wbaproj_subsets : MFS with the AWProjection gridder and A,W turned off """ |
2373 2373 | # self.prepData("refim_mawproject.ms") |
2374 2374 | # ## PS only |
2375 2375 | # ret = tclean(vis=self.msfile,spw='*',field='*',imagename=self.img,imsize=512,cell='10.0arcsec',phasecenter="J2000 19:59:28.500 +40.44.01.50",niter=30,gridder='awproject',psterm=True,aterm=False,wprojplanes=1,computepastep=360.0,rotatepastep=360.0,deconvolver='hogbom',pblimit=0.3) |
2376 2376 | # #report=self.th.checkall(imexist=[self.img+'.image', self.img+'.psf', self.img+'.weight'],imval=[(self.img+'.image',1.0,[256,256,0,0]),(self.img+'.weight',0.493,[256,256,0,0]) ] ) |