Build: #123 was successful
Job: Build and Package RHEL6 was successful
Build log
The build generated 2,109 lines of output.The output is too long and has been truncated to the last 1,000 lines. Download or view full build log.
| 28-Jul-2019 12:47:03 | - pointingspacing = "%fPB" % (gridratio_tp/pbcoeff) |
| 28-Jul-2019 12:47:03 | + pointingspacing = "%fPB" % (gridratio_tp/pbcoeff) |
| 28-Jul-2019 12:47:03 | if str.upper(pointingspacing)=="NYQUIST": |
| 28-Jul-2019 12:47:03 | pointingspacing="%fPB" % nyquist |
| 28-Jul-2019 12:47:03 | q = re.compile('(\d+.?\d+)\s*PB') |
| 28-Jul-2019 12:47:03 | @@ -617,7 +622,7 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | # find imcenter - phase center |
| 28-Jul-2019 12:47:03 | - imcenter , offsets = util.median_direction(pointings) |
| 28-Jul-2019 12:47:03 | + imcenter , offsets = util.median_direction(pointings) |
| 28-Jul-2019 12:47:03 | epoch, ra, dec = util.direction_splitter(imcenter) |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | # model is centered at model_refdir, and has model_size; |
| 28-Jul-2019 12:47:03 | @@ -650,12 +655,12 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | shift = [ (qa.convert(ra,'deg')['value'] - |
| 28-Jul-2019 12:47:03 | qa.convert(mra,'deg')['value'])*pl.cos(qa.convert(dec,'rad')['value'] ), |
| 28-Jul-2019 12:47:03 | (qa.convert(dec,'deg')['value'] - qa.convert(mdec,'deg')['value']) ] |
| 28-Jul-2019 12:47:03 | - if verbose: |
| 28-Jul-2019 12:47:03 | + if verbose: |
| 28-Jul-2019 12:47:03 | msg("pointings are shifted relative to the model by %g,%g arcsec" % (shift[0]*3600,shift[1]*3600),origin='simobserve') |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | xmax = qa.convert(model_size[0],'deg')['value']*0.5 |
| 28-Jul-2019 12:47:03 | ymax = qa.convert(model_size[1],'deg')['value']*0.5 |
| 28-Jul-2019 12:47:03 | - # add PB halfwidth (relmargin=0.5) |
| 28-Jul-2019 12:47:03 | + # add PB halfwidth (relmargin=0.5) |
| 28-Jul-2019 12:47:03 | # for mosaics of small model images |
| 28-Jul-2019 12:47:03 | xmax=xmax+pb*relmargin/3600 |
| 28-Jul-2019 12:47:03 | ymax=ymax+pb*relmargin/3600 |
| 28-Jul-2019 12:47:03 | @@ -716,13 +721,13 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | ################################################################## |
| 28-Jul-2019 12:47:03 | - # create one figure for model and pointings - need antenna diam |
| 28-Jul-2019 12:47:03 | + # create one figure for model and pointings - need antenna diam |
| 28-Jul-2019 12:47:03 | # to determine primary beam |
| 28-Jul-2019 12:47:03 | if grfile: |
| 28-Jul-2019 12:47:03 | file = fileroot + "/" + project + ".skymodel.png" |
| 28-Jul-2019 12:47:03 | else: |
| 28-Jul-2019 12:47:03 | file = "" |
| 28-Jul-2019 12:47:03 | - |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | if grscreen or grfile: |
| 28-Jul-2019 12:47:03 | util.newfig(show=grscreen) |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | @@ -734,7 +739,7 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | pl.axis("equal") |
| 28-Jul-2019 12:47:03 | else: |
| 28-Jul-2019 12:47:03 | discard = util.statim(modelflat,plot=True,incell=model_cell) |
| 28-Jul-2019 12:47:03 | - |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | lims = pl.xlim(),pl.ylim() |
| 28-Jul-2019 12:47:03 | if pb <= 0 and verbose: |
| 28-Jul-2019 12:47:03 | msg("unknown primary beam size for plot",priority="warn") |
| 28-Jul-2019 12:47:03 | @@ -764,7 +769,7 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | ylim = max(abs(pl.array(lims[1]))) |
| 28-Jul-2019 12:47:03 | # show entire pb: (statim doesn't by default) |
| 28-Jul-2019 12:47:03 | pl.xlim([max([xlim,pb/2]),min([-xlim,-pb/2])]) |
| 28-Jul-2019 12:47:03 | - pl.ylim([min([-ylim,-pb/2]),max([ylim,pb/2])]) |
| 28-Jul-2019 12:47:03 | + pl.ylim([min([-ylim,-pb/2]),max([ylim,pb/2])]) |
| 28-Jul-2019 12:47:03 | pl.xlabel("resized model sky",fontsize="x-small") |
| 28-Jul-2019 12:47:03 | util.endfig(show=grscreen,filename=file) |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | @@ -817,10 +822,10 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | if hourangle=="transit": |
| 28-Jul-2019 12:47:03 | haoffset=0.0 |
| 28-Jul-2019 12:47:03 | - else: |
| 28-Jul-2019 12:47:03 | + else: |
| 28-Jul-2019 12:47:03 | haoffset="no" |
| 28-Jul-2019 12:47:03 | # is this a time quantity? |
| 28-Jul-2019 12:47:03 | - if qa.isquantity(str(hourangle)+"h"): |
| 28-Jul-2019 12:47:03 | + if qa.isquantity(str(hourangle)+"h"): |
| 28-Jul-2019 12:47:03 | if qa.compare(str(hourangle)+"h","s"): |
| 28-Jul-2019 12:47:03 | haoffset=qa.convert(qa.quantity(str(hourangle)+ |
| 28-Jul-2019 12:47:03 | "h"),'s')['value'] |
| 28-Jul-2019 12:47:03 | @@ -882,7 +887,7 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | sm.setconfig(telescopename=telescopename, x=stnx, y=stny, z=stnz, |
| 28-Jul-2019 12:47:03 | dishdiameter=diam.tolist(), |
| 28-Jul-2019 12:47:03 | - mount=[mounttype], antname=antnames, padname=padnames, |
| 28-Jul-2019 12:47:03 | + mount=[mounttype], antname=antnames, padname=padnames, |
| 28-Jul-2019 12:47:03 | coordsystem='global', referencelocation=posobs) |
| 28-Jul-2019 12:47:03 | if str.upper(telescopename).find('VLA') >= 0: |
| 28-Jul-2019 12:47:03 | sm.setspwindow(spwname=fband, freq=qa.tos(model_start), |
| 28-Jul-2019 12:47:03 | @@ -894,12 +899,12 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | else: |
| 28-Jul-2019 12:47:03 | sm.setspwindow(spwname=fband, freq=qa.tos(model_start), |
| 28-Jul-2019 12:47:03 | deltafreq=qa.tos(model_width), |
| 28-Jul-2019 12:47:03 | - freqresolution=qa.tos(model_width), |
| 28-Jul-2019 12:47:03 | + freqresolution=qa.tos(model_width), |
| 28-Jul-2019 12:47:03 | nchannels=model_nchan, refcode=outframe, |
| 28-Jul-2019 12:47:03 | stokes='XX YY') |
| 28-Jul-2019 12:47:03 | sm.setfeed(mode='perfect X Y',pol=['']) |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | - if verbose: |
| 28-Jul-2019 12:47:03 | + if verbose: |
| 28-Jul-2019 12:47:03 | msg(" spectral window set at %s" % qa.tos(model_specrefval),origin='simobserve') |
| 28-Jul-2019 12:47:03 | sm.setlimits(shadowlimit=0.01, elevationlimit='10deg') |
| 28-Jul-2019 12:47:03 | if uvmode: |
| 28-Jul-2019 12:47:03 | @@ -910,7 +915,7 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | mereftime = me.epoch('UTC', refdate) |
| 28-Jul-2019 12:47:03 | # integration is a scalar quantity, etime is a vector of seconds |
| 28-Jul-2019 12:47:03 | - sm.settimes(integrationtime=integration, usehourangle=usehourangle, |
| 28-Jul-2019 12:47:03 | + sm.settimes(integrationtime=integration, usehourangle=usehourangle, |
| 28-Jul-2019 12:47:03 | referencetime=mereftime) |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | for k in range(0,nfld): |
| 28-Jul-2019 12:47:03 | @@ -922,7 +927,7 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | else: |
| 28-Jul-2019 12:47:03 | sourcefieldlist = sourcefieldlist + ',' + src |
| 28-Jul-2019 12:47:03 | if docalibrator: |
| 28-Jul-2019 12:47:03 | - sm.setfield(sourcename="phase calibrator", |
| 28-Jul-2019 12:47:03 | + sm.setfield(sourcename="phase calibrator", |
| 28-Jul-2019 12:47:03 | sourcedirection=caldirection,calcode='C', |
| 28-Jul-2019 12:47:03 | distance='0m') |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | @@ -972,9 +977,9 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | if kfld == nfld: |
| 28-Jul-2019 12:47:03 | if docalibrator: |
| 28-Jul-2019 12:47:03 | - endtime = sttime + qa.convert(integration,'s')['value'] |
| 28-Jul-2019 12:47:03 | + endtime = sttime + qa.convert(integration,'s')['value'] |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | - # need to observe cal singly to get new row in obs table, so |
| 28-Jul-2019 12:47:03 | + # need to observe cal singly to get new row in obs table, so |
| 28-Jul-2019 12:47:03 | # first observemany the on-source pointing(s) |
| 28-Jul-2019 12:47:03 | sm.observemany(sourcenames=srces,spwname=fband,starttimes=starttimes,stoptimes=stoptimes,project=project) |
| 28-Jul-2019 12:47:03 | # and clear the list |
| 28-Jul-2019 12:47:03 | @@ -1000,7 +1005,7 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | sm.setdata(fieldid=range(0,nfld)) |
| 28-Jul-2019 12:47:03 | if uvmode or components_only: #Interferometer only |
| 28-Jul-2019 12:47:03 | sm.setvp(dovp=True,usedefaultvp=False) |
| 28-Jul-2019 12:47:03 | - # only use mosaic gridding for Het arrays for now - |
| 28-Jul-2019 12:47:03 | + # only use mosaic gridding for Het arrays for now - |
| 28-Jul-2019 12:47:03 | # the standard gridding with a VPSkyJones is less susceptible |
| 28-Jul-2019 12:47:03 | # to issues if the image is too small which can happen a lot |
| 28-Jul-2019 12:47:03 | # in Simulation. |
| 28-Jul-2019 12:47:03 | @@ -1010,7 +1015,7 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | else: |
| 28-Jul-2019 12:47:03 | msg("Heterogeneous array only supported for mosaics (nfld>1), and make sure that your image is larger than the primary beam or results may be unstable",priority="error") |
| 28-Jul-2019 12:47:03 | else: |
| 28-Jul-2019 12:47:03 | - # checks have to be manual since there's no way to |
| 28-Jul-2019 12:47:03 | + # checks have to be manual since there's no way to |
| 28-Jul-2019 12:47:03 | # get the "diam" out of PBMath AFAIK |
| 28-Jul-2019 12:47:03 | if telescopename=="ALMA": |
| 28-Jul-2019 12:47:03 | if (diam[0]<10)|(diam[0]>13): |
| 28-Jul-2019 12:47:03 | @@ -1083,13 +1088,13 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | restfreq[i]=qa.convert(qa.quantity(model_specrefval),'Hz')['value'] |
| 28-Jul-2019 12:47:03 | tb.putcol("REST_FREQUENCY",restfreq) |
| 28-Jul-2019 12:47:03 | tb.done() |
| 28-Jul-2019 12:47:03 | - |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | ############################################ |
| 28-Jul-2019 12:47:03 | - # create figure |
| 28-Jul-2019 12:47:03 | + # create figure |
| 28-Jul-2019 12:47:03 | if grfile: |
| 28-Jul-2019 12:47:03 | file = fileroot + "/" + project + ".observe.png" |
| 28-Jul-2019 12:47:03 | else: |
| 28-Jul-2019 12:47:03 | @@ -1152,7 +1157,7 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | # if obs is unknown, casalog will send a warning to screen - temporarily suppress that |
| 28-Jul-2019 12:47:03 | if not telescopename in me.obslist(): |
| 28-Jul-2019 12:47:03 | - casalog.filter("ERROR") |
| 28-Jul-2019 12:47:03 | + casalog.filter("ERROR") |
| 28-Jul-2019 12:47:03 | im.approximatepsf(psf=fileroot+"/"+project+".quick.psf") |
| 28-Jul-2019 12:47:03 | if not telescopename in me.obslist(): |
| 28-Jul-2019 12:47:03 | casalog.filter() # set back to default level. |
| 28-Jul-2019 12:47:03 | @@ -1203,7 +1208,7 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | if not uvmode: #Single-dish |
| 28-Jul-2019 12:47:03 | msroot += ".sd" |
| 28-Jul-2019 12:47:03 | noisymsroot += ".sd" |
| 28-Jul-2019 12:47:03 | - |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | # Cosmic background radiation temperature in K. |
| 28-Jul-2019 12:47:03 | t_cmb = 2.725 |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | @@ -1243,7 +1248,7 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | # antenna efficiency |
| 28-Jul-2019 12:47:03 | eta_a = eta_p * eta_s * eta_b * eta_t |
| 28-Jul-2019 12:47:03 | - if verbose: |
| 28-Jul-2019 12:47:03 | + if verbose: |
| 28-Jul-2019 12:47:03 | msg('antenna efficiency = '+str(eta_a), origin="simobserve") |
| 28-Jul-2019 12:47:03 | msg('spillover efficiency = '+str(eta_s), origin="simobserve") |
| 28-Jul-2019 12:47:03 | msg('correlator efficiency = '+str(eta_q), origin="simobserve") |
| 28-Jul-2019 12:47:03 | @@ -1310,7 +1315,7 @@ def simobserve( |
| 28-Jul-2019 12:47:03 | if not uvmode: #Single-dish |
| 28-Jul-2019 12:47:03 | msg("Can't corrupt SD data with polarization leakage",priority="warn") |
| 28-Jul-2019 12:47:03 | if os.path.exists(msfile): |
| 28-Jul-2019 12:47:03 | - msg('copying '+msfile+' to ' + |
| 28-Jul-2019 12:47:03 | + msg('copying '+msfile+' to ' + |
| 28-Jul-2019 12:47:03 | noisymsroot+'.ms and adding polarization leakage', |
| 28-Jul-2019 12:47:03 | origin="noise",priority="warn") |
| 28-Jul-2019 12:47:03 | if os.path.exists(noisymsroot+".ms"): |
| 28-Jul-2019 12:47:03 | @@ -1395,7 +1400,7 @@ def plotpb(pb,axes,lims=None,color='k'): |
| 28-Jul-2019 12:47:03 | #ecy = lims[1][0] + bheight/2.*incy |
| 28-Jul-2019 12:47:03 | ccx = lims[0][0] + boxsize/2.*incx |
| 28-Jul-2019 12:47:03 | ccy = lims[1][0] + boxsize/2.*incy |
| 28-Jul-2019 12:47:03 | - |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | #box = Rectangle((lims[0][0],lims[1][0]),incx*bwidth,incy*bheight, |
| 28-Jul-2019 12:47:03 | box = Rectangle((lims[0][0],lims[1][0]),incx*boxsize,incy*boxsize, |
| 28-Jul-2019 12:47:03 | alpha=0.7,facecolor='w', |
| 28-Jul-2019 12:47:03 | diff --git a/src/tasks/task_tclean.py b/src/tasks/task_tclean.py |
| 28-Jul-2019 12:47:03 | index 7bf8e2f..f68d4f8 100644 |
| 28-Jul-2019 12:47:03 | --- a/src/tasks/task_tclean.py |
| 28-Jul-2019 12:47:03 | +++ b/src/tasks/task_tclean.py |
| 28-Jul-2019 12:47:03 | @@ -78,15 +78,16 @@ def tclean( |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | ### PB |
| 28-Jul-2019 12:47:03 | vptable, |
| 28-Jul-2019 12:47:03 | - usepointing, #=false |
| 28-Jul-2019 12:47:03 | mosweight, #=True |
| 28-Jul-2019 12:47:03 | aterm,#=True, |
| 28-Jul-2019 12:47:03 | psterm,#=True, |
| 28-Jul-2019 12:47:03 | wbawp ,#= True, |
| 28-Jul-2019 12:47:03 | conjbeams ,#= True, |
| 28-Jul-2019 12:47:03 | cfcache ,#= "", |
| 28-Jul-2019 12:47:03 | + usepointing, #=false |
| 28-Jul-2019 12:47:03 | computepastep ,#=360.0, |
| 28-Jul-2019 12:47:03 | rotatepastep ,#=360.0, |
| 28-Jul-2019 12:47:03 | + pointingoffsetsigdev ,#=0.0, |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | pblimit,#=0.01, |
| 28-Jul-2019 12:47:03 | normtype,#='flatnoise', |
| 28-Jul-2019 12:47:03 | @@ -216,6 +217,9 @@ def tclean( |
| 28-Jul-2019 12:47:03 | casalog.post(usemask+" is deprecated, will be removed in CASA 5.4. It is recommended to use auto-multithresh instead", "WARN") |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | #paramList.printParameters() |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + if pointingoffsetsigdev!=0.0 and usepointing==False: |
| 28-Jul-2019 12:47:03 | + casalog.post("pointingoffsetsigdev is only revelent when usepointing is True", "WARN") |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | pcube=False |
| 28-Jul-2019 12:47:03 | concattype='' |
| 28-Jul-2019 12:47:03 | diff --git a/tests/tasks/test_tclean.py b/tests/tasks/test_tclean.py |
| 28-Jul-2019 12:47:03 | index afa63ee..64ec8b4 100644 |
| 28-Jul-2019 12:47:03 | --- a/tests/tasks/test_tclean.py |
| 28-Jul-2019 12:47:03 | +++ b/tests/tasks/test_tclean.py |
| 28-Jul-2019 12:47:03 | @@ -97,7 +97,7 @@ import operator |
| 28-Jul-2019 12:47:03 | from casatasks.private.casa_transition import is_CASA6 |
| 28-Jul-2019 12:47:03 | if is_CASA6: |
| 28-Jul-2019 12:47:03 | from casatools import ctsys, quanta, measures, image, vpmanager, calibrater |
| 28-Jul-2019 12:47:03 | - from casatasks import casalog, delmod, imsubimage, tclean, uvsub, imhead, imsmooth, immath |
| 28-Jul-2019 12:47:03 | + from casatasks import casalog, delmod, imsubimage, tclean, uvsub, imhead, imsmooth, immath, widebandpbcor |
| 28-Jul-2019 12:47:03 | from casatasks.private.parallel.parallel_task_helper import ParallelTaskHelper |
| 28-Jul-2019 12:47:03 | from casatasks.private.imagerhelpers.parallel_imager_helper import PyParallelImagerHelper |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | @@ -134,7 +134,7 @@ else: |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | ## List to be run |
| 28-Jul-2019 12:47:03 | def suite(): |
| 28-Jul-2019 12:47:03 | - return [test_onefield, test_iterbot, test_multifield,test_stokes, test_modelvis, test_cube, test_mask, test_startmodel,test_widefield,test_pbcor] |
| 28-Jul-2019 12:47:03 | + return [test_onefield, test_iterbot, test_multifield,test_stokes, test_modelvis, test_cube, test_mask, test_startmodel,test_widefield,test_pbcor,test_mosaic_mtmfs,test_mosaic_cube] |
| 28-Jul-2019 12:47:03 | # return [test_onefield, test_iterbot, test_multifield,test_stokes,test_cube, test_widefield,test_mask, test_modelvis,test_startmodel,test_widefield_failing] |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | ## Base Test class with Utility functions |
| 28-Jul-2019 12:47:03 | @@ -2172,7 +2172,6 @@ class test_mask(testref_base): |
| 28-Jul-2019 12:47:03 | # report=self.th.checkall(imexist=[self.img+'.mask'], imval=[(self.img+'.mask',1.0,[50,50,0,0]),(self.img+'.mask',0.0,[50,85,0,0])]) |
| 28-Jul-2019 12:47:03 | # self.checkfinal(report) |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | - |
| 28-Jul-2019 12:47:03 | # the option, auto-thresh removed |
| 28-Jul-2019 12:47:03 | # def test_mask_autobox_pbmask(self): |
| 28-Jul-2019 12:47:03 | # """ [mask] test_mask_autobox_nmask : Autoboxing with pbmask""" |
| 28-Jul-2019 12:47:03 | @@ -2223,6 +2222,7 @@ class test_mask(testref_base): |
| 28-Jul-2019 12:47:03 | # report=self.th.checkall(imexist=[self.img+'.mask'], imval=[(self.img+'.mask',1.0,[50,50,0,0]),(self.img+'.mask',0.0,[50,85,0,0])]) |
| 28-Jul-2019 12:47:03 | # self.checkfinal(report) |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | + # AUTOMASK TESTS |
| 28-Jul-2019 12:47:03 | def test_mask_autobox_multithresh(self): |
| 28-Jul-2019 12:47:03 | """ [mask] test_mask__autobox_multithresh : multi-threshold Autobox (default)""" |
| 28-Jul-2019 12:47:03 | self.prepData('refim_twochan.ms') |
| 28-Jul-2019 12:47:03 | @@ -2279,6 +2279,19 @@ class test_mask(testref_base): |
| 28-Jul-2019 12:47:03 | report=self.th.checkall(imexist=[self.img+'.mask'], imval=[(self.img+'.mask',1.0,[50,50,0,0]),(self.img+'.mask',1.0,[60,40,0,0]),(self.img+'.mask',0.0,[65,50,0,0])]) |
| 28-Jul-2019 12:47:03 | self.checkfinal(report) |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | + def test_mask_autobox_multithresh_mfs_IQUV(self): |
| 28-Jul-2019 12:47:03 | + """ [mask] test_mask__autobox_multithresh_mtmfs_IQUV : multi-threshold Autobox (minbeamfrac=0.3) with cube full polarizaiton (IQUV) imaging """ |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_point_linRL.ms') |
| 28-Jul-2019 12:47:03 | + ret = tclean(vis=self.msfile,imagename=self.img,imsize=100,cell='8.0arcsec',niter=10, stokes='IQUV',interactive=0,specmode='mfs',interpolation='nearest',usemask="auto-multithresh", verbose=True, parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report=self.th.checkall(imexist=[self.img+'.mask'], imval=[(self.img+'.mask',1.0,[50,50,0,0]),(self.img+'.mask',1.0,[40,60,0,0]),(self.img+'.mask',0.0,[65,50,0,0]), (self.img+'.mask', 1.0,[40,60,3,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_mask_autobox_multithresh_cube_IQUV(self): |
| 28-Jul-2019 12:47:03 | + """ [mask] test_mask__autobox_multithresh_cube_IQUV : multi-threshold Autobox (minbeamfrac=0.05) with cube full polarizaiton (IQUV) imaging """ |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_point_linXY.ms') |
| 28-Jul-2019 12:47:03 | + ret = tclean(vis=self.msfile,imagename=self.img,imsize=100,cell='8.0arcsec',niter=10, stokes='IQUV',interactive=0,specmode='cube',interpolation='nearest',usemask="auto-multithresh", minbeamfrac=0.05, verbose=True, parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report=self.th.checkall(imexist=[self.img+'.mask'], imval=[(self.img+'.mask',1.0,[50,50,0,0]),(self.img+'.mask',1.0,[35,75,0,0]),(self.img+'.mask',0.0,[32,80,1,1]), (self.img+'.mask',1.0,[35,60,1,1]), (self.img+'.mask',1.0,[60,30,3,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | # def test_mask_outregion(self): |
| 28-Jul-2019 12:47:03 | # """ [mask] test_mask_outregion : Input mask has region that goes outside the image """ |
| 28-Jul-2019 12:47:03 | # self.prepData('refim_twochan.ms') |
| 28-Jul-2019 12:47:03 | @@ -3192,6 +3205,662 @@ class test_pbcor(testref_base): |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | self.checkfinal(report1+report2) |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | +##################################################### |
| 28-Jul-2019 12:47:03 | +##################################################### |
| 28-Jul-2019 12:47:03 | +##################################################### |
| 28-Jul-2019 12:47:03 | +##################################################### |
| 28-Jul-2019 12:47:03 | +class test_mosaic_mtmfs(testref_base): |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_standard_onefield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec', gridder='standard',field=field,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + widebandpbcor(vis=self.msfile, imagename=self.img,nterms=2,reffreq='1.5GHz',pbmin=0.1, field=field, spwlist=[0,1,2], chanlist=[0,0,0], weightlist=[1.0,1.0,1.0]) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.pbcor.image.tt0',1.021635,[512,596,0,0]),(self.img+'.pbcor.image.alpha', -0.4811849,[512,596,0,0]),(self.img+'.pbcor.workdirectory/' + self.img + '.pb.alpha',-1.526915,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec', gridder='standard',field=field,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + widebandpbcor(vis=self.msfile, imagename=self.img,nterms=2,reffreq='1.5GHz',pbmin=0.1, field=field, spwlist=[0,1,2], chanlist=[0,0,0], weightlist=[1.0,1.0,1.0]) |
| 28-Jul-2019 12:47:03 | + report2=self.th.checkall(imval=[(self.img+'.pbcor.image.tt0',1.02163529396,[512,596,0,0]),(self.img+'.pbcor.image.alpha', -0.481181353331 ,[512,596,0,0]),(self.img+'.pbcor.workdirectory/' + self.img + '.pb.alpha',-1.52691471577,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ########################### |
| 28-Jul-2019 12:47:03 | + #warning |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_standard_twofield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter ='J2000 19h59m28.5 +40d40m01.5' # pointing center of field0 |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec', gridder='standard',field=field,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.tt0',0.696767389774,[512,596,0,0]),(self.img+'.alpha', -1.22704184055,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + widebandpbcor(vis=self.msfile, imagename=self.img,nterms=2,reffreq='1.5GHz',pbmin=0.1, field=field, spwlist=[0,1,2], chanlist=[0,0,0], weightlist=[1.0,1.0,1.0]) |
| 28-Jul-2019 12:47:03 | + report2=self.th.checkall(imval=[(self.img+'.pbcor.image.tt0',0.7115,[512,596,0,0]),(self.img+'.pbcor.image.alpha',-1.2540572,[512,596,0,0]),(self.img+'.pbcor.workdirectory/' + self.img + '.pb.alpha',0.027015,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + os.system('rm -rf ' + self.img+'*') |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec', gridder='standard',field=field,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.tt0', 0.696766972542 ,[512,596,0,0]),(self.img+'.alpha',-1.2270386219,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + widebandpbcor(vis=self.msfile, imagename=self.img,nterms=2,reffreq='1.5GHz',pbmin=0.1, field=field, spwlist=[0,1,2], chanlist=[0,0,0], weightlist=[1.0,1.0,1.0]) |
| 28-Jul-2019 12:47:03 | + report4=self.th.checkall(imval=[(self.img+'.pbcor.image.tt0',0.7115,[512,596,0,0]),(self.img+'.pbcor.image.alpha',-1.2540572,[512,596,0,0]),(self.img+'.pbcor.workdirectory/' + self.img + '.pb.alpha',0.027015,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1 + report2 + report3 + report4 + '\n Warning: values must be theoretically validated') |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | +############################# |
| 28-Jul-2019 12:47:03 | + # |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_mosaic_cbFalse_onefield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.tt0',0.5302894,[512,596,0,0]),(self.img+'.pb.tt0',0.5144197,[512,596,0,0]),(self.img+'.alpha',-3.29580,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report2=self.th.checkall(imval=[(self.img+'.image.tt0',0.450807213783 ,[512,596,0,0]),(self.img+'.pb.tt0',0.514419734478,[512,596,0,0]),(self.img+'.alpha',-1.6693893671,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_mosaic_cbFalse_twofield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter ='J2000 19h59m28.5 +40d40m01.5' # pointing center of field0 |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.tt0', 0.97402,[512,596,0,0]),(self.img+'.pb.tt0', 0.982510,[512,596,0,0]),(self.img+'.alpha', -0.9880,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report2=self.th.checkall(imval=[(self.img+'.image.tt0', 0.982747793198,[512,596,0,0]),(self.img+'.pb.tt0', 0.982510685921,[512,596,0,0]),(self.img+'.alpha', -0.69006639719,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1 + report2 + '\n Warning: values must be theoretically validated') |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_mosaic_cbTrue_onefield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.tt0',0.492861807346,[512,596,0,0]),(self.img+'.pb.tt0', 0.492951124907,[512,596,0,0]),(self.img+'.alpha',-0.586531162262,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report2=self.th.checkall(imval=[(self.img+'.image.tt0',0.494806170464,[512,596,0,0]),(self.img+'.pb.tt0', 0.492951124907,[512,596,0,0]),(self.img+'.alpha',-0.548104763031,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1 + report2) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_mosaic_cbTrue_twofield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter ='J2000 19h59m28.5 +40d40m01.5' # pointing center of field0 |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.tt0',1.004476,[512,596,0,0]),(self.img+'.pb.tt0',0.98724,[512,596,0,0]),(self.img+'.alpha',-1.216693,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report2=self.th.checkall(imval=[(self.img+'.image.tt0',0.989150226116,[512,596,0,0]),(self.img+'.pb.tt0',0.987243115902,[512,596,0,0]),(self.img+'.alpha',-0.771613717079,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1 + report2 + '\n Warning: values must be theoretically validated') |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ################### |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_mosaic_cbTrue_onefield_use_standard_psf(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbTrue' |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Run the standard gridder imaging. |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename='std',niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='standard',field=field,cfcache=cfcache_path,conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Make the residual images using mosaic/awproject |
| 28-Jul-2019 12:47:03 | + os.system('rm -rf '+self.img+'.*'); |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field,cfcache=cfcache_path,conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Replace the PSFs with those from the standard gridder |
| 28-Jul-2019 12:47:03 | + os.system('rm -rf '+self.img+'.psf*') |
| 28-Jul-2019 12:47:03 | + os.system('cp -r std.psf.tt0 '+self.img+'.psf.tt0') |
| 28-Jul-2019 12:47:03 | + os.system('cp -r std.psf.tt1 '+self.img+'.psf.tt1') |
| 28-Jul-2019 12:47:03 | + os.system('cp -r std.psf.tt2 '+self.img+'.psf.tt2') |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Restart tclean with calcres=False and calcpsf=False |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field,cfcache=cfcache_path, conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,calcres=False, calcpsf=False,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + report=self.th.checkall(imval=[(self.img+'.image.tt0', 0.492902100086 ,[512,596,0,0]),(self.img+'.pb.tt0', 0.492951124907 ,[512,596,0,0]),(self.img+'.alpha', -0.591838240623,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_mosaic_cbTrue_twofield_use_standard_psf(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbTrue' |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Run the standard gridder imaging. |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename='std',niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='standard',field=field,cfcache=cfcache_path,conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Make the residual images using mosaic/awproject |
| 28-Jul-2019 12:47:03 | + os.system('rm -rf '+self.img+'.*'); |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field,cfcache=cfcache_path,conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Replace the PSFs with those from the standard gridder |
| 28-Jul-2019 12:47:03 | + os.system('rm -rf '+self.img+'.psf*') |
| 28-Jul-2019 12:47:03 | + os.system('cp -r std.psf.tt0 '+self.img+'.psf.tt0') |
| 28-Jul-2019 12:47:03 | + os.system('cp -r std.psf.tt1 '+self.img+'.psf.tt1') |
| 28-Jul-2019 12:47:03 | + os.system('cp -r std.psf.tt2 '+self.img+'.psf.tt2') |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Restart tclean with calcres=False and calcpsf=False |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field,cfcache=cfcache_path, conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,calcres=False, calcpsf=False,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report=self.th.checkall(imval=[(self.img+'.image.tt0', 0.982702195644,[512,596,0,0]),(self.img+'.pb.tt0', 0.987243115902,[512,596,0,0]),(self.img+'.alpha', -0.55044066906,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | +############################# |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_awproject_cbFalse_onefield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter ='' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbFalse' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field, cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.tt0',0.5269711,[512,596,0,0]),(self.img+'.pb.tt0',0.50752753,[512,596,0,0]),(self.img+'.alpha',-3.24132061,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field, cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report2=self.th.checkall(imval=[(self.img+'.image.tt0',0.451119929552,[512,596,0,0]),(self.img+'.pb.tt0',0.507527530193,[512,596,0,0]),(self.img+'.alpha',-1.65221953392,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1 + report2) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_awproject_cbFalse_twofield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter ='J2000 19h59m28.5 +40d40m01.5' # pointing center of field0 |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbFalse' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field, cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.tt0',0.97862583398,[512,596,0,0]),(self.img+'.pb.tt0',0.979142010212,[512,596,0,0]),(self.img+'.alpha',-1.24368548393,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject', field=field, cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report2=self.th.checkall(imval=[(self.img+'.image.tt0',0.975049376488,[512,596,0,0]),(self.img+'.pb.tt0',0.979141950607,[512,596,0,0]),(self.img+'.alpha',-0.80353230238,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1 + report2 + '\n Warning: values must be theoretically validated') |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_awproject_cbTrue_onefield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbTrue' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field, cfcache=cfcache_path, conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.tt0',0.477538466454,[512,596,0,0]),(self.img+'.pb.tt0', 0.479197412729,[512,596,0,0]),(self.img+'.alpha', -0.562356948853,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field, cfcache=cfcache_path, conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report2=self.th.checkall(imval=[(self.img+'.image.tt0',0.482177525759,[512,596,0,0]),(self.img+'.pb.tt0', 0.479197442532,[512,596,0,0]),(self.img+'.alpha', -0.568624258041,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1 + report2) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_awproject_cbTrue_twofield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter ='J2000 19h59m28.5 +40d40m01.5' # pointing center of field0 |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbTrue' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field, cfcache= cfcache_path, conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.tt0', 0.974482476711 ,[512,596,0,0]),(self.img+'.pb.tt0',0.979797422886,[512,596,0,0]),(self.img+'.alpha', -0.538577735424 ,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field, cfcache= cfcache_path, conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report2=self.th.checkall(imval=[(self.img+'.image.tt0', 0.97661459446 ,[512,596,0,0]),(self.img+'.pb.tt0',0.979797422886,[512,596,0,0]),(self.img+'.alpha', -0.538577854633 ,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2) |
| 28-Jul-2019 12:47:03 | +################### |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_awproject_cbTrue_onefield_use_standard_psf(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbTrue' |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Run the standard gridder imaging. |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename='std',niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='standard',field=field,cfcache=cfcache_path,conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Make the residual images using mosaic/awproject |
| 28-Jul-2019 12:47:03 | + os.system('rm -rf '+self.img+'.*'); |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path,conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Replace the PSFs with those from the standard gridder |
| 28-Jul-2019 12:47:03 | + os.system('rm -rf '+self.img+'.psf*') |
| 28-Jul-2019 12:47:03 | + os.system('cp -r std.psf.tt0 '+self.img+'.psf.tt0') |
| 28-Jul-2019 12:47:03 | + os.system('cp -r std.psf.tt1 '+self.img+'.psf.tt1') |
| 28-Jul-2019 12:47:03 | + os.system('cp -r std.psf.tt2 '+self.img+'.psf.tt2') |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Restart tclean with calcres=False and calcpsf=False |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,calcres=False, calcpsf=False,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + report=self.th.checkall(imval=[(self.img+'.image.tt0',0.477585822344,[512,596,0,0]),(self.img+'.pb.tt0', 0.479197412729 ,[512,596,0,0]),(self.img+'.alpha', -0.569523513317,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_awproject_cbTrue_twofield_use_standard_psf(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbTrue' |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Run the standard gridder imaging. |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename='std',niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='standard',field=field,cfcache=cfcache_path,conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Make the residual images using mosaic/awproject |
| 28-Jul-2019 12:47:03 | + os.system('rm -rf '+self.img+'.*'); |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path,conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Replace the PSFs with those from the standard gridder |
| 28-Jul-2019 12:47:03 | + os.system('rm -rf '+self.img+'.psf*') |
| 28-Jul-2019 12:47:03 | + os.system('cp -r std.psf.tt0 '+self.img+'.psf.tt0') |
| 28-Jul-2019 12:47:03 | + os.system('cp -r std.psf.tt1 '+self.img+'.psf.tt1') |
| 28-Jul-2019 12:47:03 | + os.system('cp -r std.psf.tt2 '+self.img+'.psf.tt2') |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + ## Restart tclean with calcres=False and calcpsf=False |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, conjbeams=True, wbawp=True, psterm=False,pblimit=0.1,calcres=False, calcpsf=False,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report=self.th.checkall(imval=[(self.img+'.image.tt0', 0.974559485912,[512,596,0,0]),(self.img+'.pb.tt0', 0.979797422886,[512,596,0,0]),(self.img+'.alpha', -0.542647540569,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_mosaic_cbFalse_onefield_upTrue(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, usepointing=True, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.tt0',0.5302894,[512,596,0,0]),(self.img+'.pb.tt0',0.5144197,[512,596,0,0]),(self.img+'.alpha',-3.29580,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, usepointing=True, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report2=self.th.checkall(imval=[(self.img+'.image.tt0',0.450807213783 ,[512,596,0,0]),(self.img+'.pb.tt0',0.514419734478,[512,596,0,0]),(self.img+'.alpha',-1.6693893671,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_mosaic_cbFalse_twofield_upTrue(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter ='J2000 19h59m28.5 +40d40m01.5' # pointing center of field0 |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, usepointing=True, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.tt0', 0.97402,[512,596,0,0]),(self.img+'.pb.tt0', 0.982510,[512,596,0,0]),(self.img+'.alpha', -0.9880,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, usepointing=True, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report2=self.th.checkall(imval=[(self.img+'.image.tt0', 0.982747793198,[512,596,0,0]),(self.img+'.pb.tt0', 0.982510685921,[512,596,0,0]),(self.img+'.alpha', -0.69006639719,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1 + report2 + '\n Warning: values must be theoretically validated') |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_awproject_cbFalse_onefield_upTrue(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter ='' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbFalse' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field, usepointing=True, cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.tt0',0.5269711,[512,596,0,0]),(self.img+'.pb.tt0',0.50752753,[512,596,0,0]),(self.img+'.alpha',-3.24132061,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field, usepointing=True, cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report2=self.th.checkall(imval=[(self.img+'.image.tt0',0.451119929552,[512,596,0,0]),(self.img+'.pb.tt0',0.507527530193,[512,596,0,0]),(self.img+'.alpha',-1.65221953392,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1 + report2) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_mtmfs_awproject_cbFalse_twofield_upTrue(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter ='J2000 19h59m28.5 +40d40m01.5' # pointing center of field0 |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbFalse' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field, usepointing=True, cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.tt0',0.97862583398,[512,596,0,0]),(self.img+'.pb.tt0',0.979142010212,[512,596,0,0]),(self.img+'.alpha',-1.24368548393,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='mfs',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject', field=field, usepointing=True, cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,deconvolver='mtmfs',nterms=2,reffreq='1.5GHz',pbcor=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report2=self.th.checkall(imval=[(self.img+'.image.tt0',0.975049376488,[512,596,0,0]),(self.img+'.pb.tt0',0.979141950607,[512,596,0,0]),(self.img+'.alpha',-0.80353230238,[512,596,0,0])]) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1 + report2 + '\n Warning: values must be theoretically validated') |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | +########################################################### |
| 28-Jul-2019 12:47:03 | +########################################################### |
| 28-Jul-2019 12:47:03 | +########################################################### |
| 28-Jul-2019 12:47:03 | +class test_mosaic_cube(testref_base): |
| 28-Jul-2019 12:47:03 | + def test_cube_standard_onefield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec', gridder='standard',field=field,pblimit=0.1,reffreq='1.5GHz',pbcor=True,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.pbcor', 1.10922563076,[512,596,0,0]),(self.img+'.image.pbcor',0.989520609379,[512,596,0,1]),(self.img+'.image.pbcor',0.90361648798,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report2 = self.th.checkval(spectral_index, -0.505622766021, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec', gridder='standard',field=field,pblimit=0.1,reffreq='1.5GHz',pbcor=True,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.pbcor', 1.10922467709,[512,596,0,0]),(self.img+'.image.pbcor', 0.989521086216,[512,596,0,1]),(self.img+'.image.pbcor',0.903617084026,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report4 = self.th.checkval(spectral_index, -0.505622766021, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2+report3+report4) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_cube_standard_twofield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter ='J2000 19h59m28.5 +40d40m01.5' # pointing center of field0 |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec', gridder='standard',field=field,pblimit=0.1,reffreq='1.5GHz',pbcor=True,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.pbcor',0.895364701748,[512,596,0,0]),(self.img+'.image.pbcor',0.701696813107,[512,596,0,1]),(self.img+'.image.pbcor',0.539412796497,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report2 = self.th.checkval(spectral_index, -1.249799253, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec', gridder='standard',field=field,pblimit=0.1,reffreq='1.5GHz',pbcor=True,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.pbcor', 0.895363807678,[512,596,0,0]),(self.img+'.image.pbcor',0.701696515083,[512,596,0,1]),(self.img+'.image.pbcor',0.539413094521,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report4 = self.th.checkval(spectral_index, -1.24979542771, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2+report3+report4) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | +########################################################### |
| 28-Jul-2019 12:47:03 | + def test_cube_mosaic_cbFalse_mwFalse_onefield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.pbcor',1.1021887064,[512,596,0,0]),(self.img+'.image.pbcor',0.983152985573,[512,596,0,1]),(self.img+'.image.pbcor',0.905956506729,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report2 = self.th.checkval(spectral_index, -0.483567460558, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.pbcor',1.11281049252,[512,596,0,0]),(self.img+'.image.pbcor',0.994237959385,[512,596,0,1]),(self.img+'.image.pbcor',0.912590324879,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report4 = self.th.checkval(spectral_index, -0.48920856272, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2+report3+report4) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_cube_mosaic_cbFalse_mwFalse_twofield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.pbcor',1.10354316235,[512,596,0,0]),(self.img+'.image.pbcor',0.981979846954,[512,596,0,1]),(self.img+'.image.pbcor', 0.895015060902,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report2 = self.th.checkval(spectral_index, -0.516544552, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.pbcor',1.11373543739,[512,596,0,0]),(self.img+'.image.pbcor',0.994349777699,[512,596,0,1]),(self.img+'.image.pbcor', 0.909038066864,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report4 = self.th.checkval(spectral_index, -0.500876470767, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2+report3+report4) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_cube_mosaic_cbFalse_mwTrue_onefield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=True,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.pbcor',1.1021887064,[512,596,0,0]),(self.img+'.image.pbcor',0.983152985573,[512,596,0,1]),(self.img+'.image.pbcor',0.905956506729,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report2 = self.th.checkval(spectral_index, -0.483548182038, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=True,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.pbcor',1.11281049252,[512,596,0,0]),(self.img+'.image.pbcor',0.994237959385,[512,596,0,1]),(self.img+'.image.pbcor',0.912590324879,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report4 = self.th.checkval(spectral_index, -0.48920856272, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2+report3+report4) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_cube_mosaic_cbFalse_mwTrue_twofield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=True,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.pbcor',1.10354316235,[512,596,0,0]),(self.img+'.image.pbcor', 0.981979727745,[512,596,0,1]),(self.img+'.image.pbcor', 0.895014822483,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report2 = self.th.checkval(spectral_index, -0.51654520997, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=True,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.pbcor',1.11373543739,[512,596,0,0]),(self.img+'.image.pbcor', 0.994349777699,[512,596,0,1]),(self.img+'.image.pbcor', 0.909038066864,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report4 = self.th.checkval(spectral_index, -0.50087647076, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2+report3+report4) |
| 28-Jul-2019 12:47:03 | +##################################################### |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_cube_awproject_cbFalse_mwFalse_onefield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbFalse' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.pbcor', 1.1262229681,[512,596,0,0]),(self.img+'.image.pbcor', 0.996681272984,[512,596,0,1]),(self.img+'.image.pbcor', 0.879481077194,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report2 = self.th.checkval(spectral_index, -0.641964, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.pbcor', 1.13023257256,[512,596,0,0]),(self.img+'.image.pbcor', 1.00326132774,[512,596,0,1]),(self.img+'.image.pbcor',0.8681204319,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report4 = self.th.checkval(spectral_index, -0.618663982179, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2+report3+report4) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_cube_awproject_cbFalse_mwFalse_twofield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbFalse' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.pbcor',1.11765027046,[512,596,0,0]),(self.img+'.image.pbcor',0.990656971931,[512,596,0,1]),(self.img+'.image.pbcor', 0.879846811295,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report2 = self.th.checkval(spectral_index, -0.590028509558, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.pbcor',1.12036144733,[512,596,0,0]),(self.img+'.image.pbcor',0.994982719421,[512,596,0,1]),(self.img+'.image.pbcor', 0.889532327652,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report4 = self.th.checkval(spectral_index, -0.569002802902, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2+report3+report4) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_cube_awproject_cbFalse_mwTrue_onefield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbFalse' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=True,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.pbcor',1.1262229681,[512,596,0,0]),(self.img+'.image.pbcor', 0.996681272984,[512,596,0,1]),(self.img+'.image.pbcor', 0.8681204319,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report2 = self.th.checkval(spectral_index, -0.641964870168, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=True,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.pbcor',1.13023257256,[512,596,0,0]),(self.img+'.image.pbcor', 1.00326132774,[512,596,0,1]),(self.img+'.image.pbcor', 0.879481077194,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report4 = self.th.checkval(spectral_index, -0.61866398217, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2+report3+report4) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_cube_awproject_cbFalse_mwTrue_twofield(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbFalse' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=True,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.pbcor',1.11765027046,[512,596,0,0]),(self.img+'.image.pbcor',0.990656971931,[512,596,0,1]),(self.img+'.image.pbcor', 0.879846811295,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report2 = self.th.checkval(spectral_index, -0.590028509558, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=True,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.pbcor',1.12036144733,[512,596,0,0]),(self.img+'.image.pbcor',0.994982719421,[512,596,0,1]),(self.img+'.image.pbcor', 0.889532327652,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report4 = self.th.checkval(spectral_index, -0.569002802902, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2+report3+report4) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | +##################################################### |
| 28-Jul-2019 12:47:03 | +##################################################### |
| 28-Jul-2019 12:47:03 | +##################################################### |
| 28-Jul-2019 12:47:03 | + def test_cube_mosaic_cbFalse_mwFalse_onefield_upTrue(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, usepointing = True, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.pbcor',1.1021887064,[512,596,0,0]),(self.img+'.image.pbcor',0.983152985573,[512,596,0,1]),(self.img+'.image.pbcor',0.905956506729,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report2 = self.th.checkval(spectral_index, -0.483567460558, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, usepointing = True, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.pbcor',1.11281049252,[512,596,0,0]),(self.img+'.image.pbcor',0.994237959385,[512,596,0,1]),(self.img+'.image.pbcor',0.912590324879,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report4 = self.th.checkval(spectral_index, -0.48920856272, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2+report3+report4) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_cube_mosaic_cbFalse_mwFalse_twofield_upTrue(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, usepointing = True, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.pbcor',1.10354316235,[512,596,0,0]),(self.img+'.image.pbcor',0.981979846954,[512,596,0,1]),(self.img+'.image.pbcor', 0.895015060902,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report2 = self.th.checkval(spectral_index, -0.516544552, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='mosaic',field=field, usepointing = True, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.pbcor',1.11373543739,[512,596,0,0]),(self.img+'.image.pbcor',0.994349777699,[512,596,0,1]),(self.img+'.image.pbcor', 0.909038066864,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report4 = self.th.checkval(spectral_index, -0.500876470767, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2+report3+report4) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_cube_awproject_cbFalse_mwFalse_onefield_upTrue(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbFalse' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, usepointing = True, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.pbcor', 1.1262229681,[512,596,0,0]),(self.img+'.image.pbcor', 0.996681272984,[512,596,0,1]),(self.img+'.image.pbcor', 0.879481077194,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report2 = self.th.checkval(spectral_index, -0.641964, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, usepointing = True, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.pbcor', 1.13023257256,[512,596,0,0]),(self.img+'.image.pbcor', 1.00326132774,[512,596,0,1]),(self.img+'.image.pbcor',0.8681204319,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report4 = self.th.checkval(spectral_index, -0.618663982179, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2+report3+report4) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + def test_cube_awproject_cbFalse_mwFalse_twofield_upTrue(self): |
| 28-Jul-2019 12:47:03 | + self.prepData('refim_oneshiftpoint.mosaic.ms') |
| 28-Jul-2019 12:47:03 | + phasecenter = '' |
| 28-Jul-2019 12:47:03 | + field='0,1' |
| 28-Jul-2019 12:47:03 | + cfcache_path = refdatapath + 'cfcache_oneshiftpoint_mosaic_cbFalse' |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=0,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, usepointing = True, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report1=self.th.checkall(imval=[(self.img+'.image.pbcor',1.11765027046,[512,596,0,0]),(self.img+'.image.pbcor',0.990656971931,[512,596,0,1]),(self.img+'.image.pbcor', 0.879846811295,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report2 = self.th.checkval(spectral_index, -0.590028509558, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + tclean(vis=self.msfile, imagename=self.img,niter=10,specmode='cube',spw='*',imsize=1024, phasecenter=phasecenter,cell='10.0arcsec',gridder='awproject',field=field,cfcache=cfcache_path, usepointing = True, conjbeams=False, wbawp=True, psterm=False,pblimit=0.1,reffreq='1.5GHz',pbcor=True,mosweight=False,parallel=self.parallel) |
| 28-Jul-2019 12:47:03 | + report3=self.th.checkall(imval=[(self.img+'.image.pbcor',1.12036144733,[512,596,0,0]),(self.img+'.image.pbcor',0.994982719421,[512,596,0,1]),(self.img+'.image.pbcor', 0.889532327652,[512,596,0,2])]) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | + source_flux_v0 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,0]) |
| 28-Jul-2019 12:47:03 | + source_flux_v2 = self.th.get_pix(self.img+'.image.pbcor',[512,596,0,2]) |
| 28-Jul-2019 12:47:03 | + v0 = 1.2 #In GHz |
| 28-Jul-2019 12:47:03 | + v2 = 1.8 #In GHz |
| 28-Jul-2019 12:47:03 | + spectral_index = np.log(source_flux_v0/source_flux_v2)/np.log(v0/v2) |
| 28-Jul-2019 12:47:03 | + report4 = self.th.checkval(spectral_index, -0.569002802902, valname='Spectral flux', exact=False) |
| 28-Jul-2019 12:47:03 | + self.checkfinal(report1+report2+report3+report4) |
| 28-Jul-2019 12:47:03 | + |
| 28-Jul-2019 12:47:03 | if is_CASA6: |
| 28-Jul-2019 12:47:03 | if __name__ == '__main__': |
| 28-Jul-2019 12:47:03 | unittest.main() |
| 28-Jul-2019 12:47:03 | Argument values: |
| 28-Jul-2019 12:47:03 | casatasks |
| 28-Jul-2019 12:47:03 | /export/home/cbt-el7-6/bamboohome/xml-data/build-dir/CASA-CTB-JOB1/casatasks |
| 28-Jul-2019 12:47:03 | master |
| 28-Jul-2019 12:47:03 | None |
| 28-Jul-2019 12:47:03 | Last Tag: 2019.118 |
| 28-Jul-2019 12:47:03 | New Tag: 2019.126 |
| 28-Jul-2019 12:47:03 | Creating a new tag |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | None |
| 28-Jul-2019 12:47:03 | Push enabled. Pushing new tag upstream. |
| 28-Jul-2019 12:47:03 | To ssh://git@open-bitbucket.nrao.edu:7999/casa/casatasks.git |
| 28-Jul-2019 12:47:03 | * [new tag] 2019.126 -> 2019.126 |
| 28-Jul-2019 12:47:03 | |
| 28-Jul-2019 12:47:03 | None |
| 28-Jul-2019 12:47:03 | Finished task 'Create wheel' with result: Success |
| 28-Jul-2019 12:47:03 | Running post build plugin 'Docker Container Cleanup' |
| 28-Jul-2019 12:47:03 | Beginning to execute external process for build 'CASA - Casa6 Tasks Wheel Build - Build and Package RHEL6 #123 (CASA-CTB-JOB1-123)' ... running command line: /bin/docker rm -f task-wheel-container ... in: /export/home/cbt-el7-6/bamboohome/xml-data/build-dir/CASA-CTB-JOB1 |
| 28-Jul-2019 12:47:04 | task-wheel-container |
| 28-Jul-2019 12:47:04 | Running post build plugin 'NCover Results Collector' |
| 28-Jul-2019 12:47:04 | Running post build plugin 'Clover Results Collector' |
| 28-Jul-2019 12:47:04 | Running post build plugin 'npm Cache Cleanup' |
| 28-Jul-2019 12:47:04 | Running post build plugin 'Artifact Copier' |
| 28-Jul-2019 12:47:04 | Publishing an artifact: EL6 wheel |
| 28-Jul-2019 12:47:04 | Finished publishing of artifact Non required shared artifact: [EL6 wheel], pattern: [*-any.whl] anchored at: [wheeldirectory/] in 138.9 ms |
| 28-Jul-2019 12:47:04 | Finalising the build... |
| 28-Jul-2019 12:47:04 | Stopping timer. |
| 28-Jul-2019 12:47:04 | Build CASA-CTB-JOB1-123 completed. |
| 28-Jul-2019 12:47:04 | Running on server: post build plugin 'NCover Results Collector' |
| 28-Jul-2019 12:47:04 | Running on server: post build plugin 'Build Hanging Detection Configuration' |
| 28-Jul-2019 12:47:04 | Running on server: post build plugin 'Clover Delta Calculator' |
| 28-Jul-2019 12:47:04 | Running on server: post build plugin 'Maven Dependencies Postprocessor' |
| 28-Jul-2019 12:47:05 | All post build plugins have finished |
| 28-Jul-2019 12:47:05 | Generating build results summary... |
| 28-Jul-2019 12:47:05 | Saving build results to disk... |
| 28-Jul-2019 12:47:05 | Logging substituted variables... |
| 28-Jul-2019 12:47:05 | Indexing build results... |
| 28-Jul-2019 12:47:05 | Finished building CASA-CTB-JOB1-123. |