From gtaylor@cv3.cv.nrao.edu Wed Mar 28 16:36:35 2001 Date: Wed, 28 Mar 2001 16:30:38 -0700 (MST) From: gtaylor To: Steven T. Myers Cc: AIPS++ Core Testing Group , dshepher@zia.aoc.NRAO.EDU, efomalon@cv3.cv.nrao.edu, gmoellen@zia.aoc.NRAO.EDU, gtaylor@zia.aoc.NRAO.EDU, gvanmoor@zia.aoc.NRAO.EDU, jbenson@zia.aoc.NRAO.EDU, jhibbard@cv3.cv.nrao.edu, julvesta@zia.aoc.NRAO.EDU, tcornwel@zia.aoc.NRAO.EDU Subject: Testing log (March 22-28) Testing AIPS++ Detailed Log Mar 22-28, 2001, session: note: TTT = trouble BBB = bug QQQ = question Athol and George report that bugs in calibrator have been found and fixed so I am back at work on feb19LL data. Calibration appears nominal through to fluxscale which reports: aips++ Flux density for 0319+415 (spw=1) is: 24.409 +/- 0.251 Jy Flux density for 0319+415 (spw=2) is: 23.177 +/- 0.139 Jy Flux density for 0521+166 (spw=1) is: 9.741 +/- 0.099 Jy Flux density for 0521+166 (spw=2) is: 8.777 +/- 0.049 Jy Flux density for 0713+438 (spw=1) is: 2.154 +/- 0.019 Jy Flux density for 0713+438 (spw=2) is: 2.038 +/- 0.009 Jy These values are a bit high. Lets flag some antennas on 3C286 and try again. Went into msplot and flagged antennas 24, 25 and 26 at the end of the North arm. Data won't stay flagged!!!!! Trying again on March 28 in weekly again. Tried to plot the data. BBB: - LocalExec::SetStatus: abnormal child termination for /home/cluster/aips++/weekly/linux_gnu/bin/ms Server 'ms' has failed unexpectedly! You will need to create the relevant tool again. If that causes unexpected behavior, please restart AIPS++ Please submit a bug-report using bug() if you can reproduce the problem. and again. So I'm bugging Athol. His diagnosis is that I had a bad save/get file. To get rid of them: rm aips++.inputs* rm *.gc.* Msplot works and the data appear to be flagged. Reran my calibration script and the results are much worse: Set the flux density scale using reference calibrators Flux density for 0319+415 (spw=1) is: 22.074 +/- 0.070 Jy Flux density for 0319+415 (spw=2) is: 33.205 +/- 1.154 Jy Flux density for 0521+166 (spw=1) is: 8.818 +/- 0.029 Jy Flux density for 0521+166 (spw=2) is: 12.695 +/- 0.451 Jy Flux density for 0713+438 (spw=1) is: 1.971 +/- 0.005 Jy Flux density for 0713+438 (spw=2) is: 2.800 +/- 0.087 Jy What is the deal? This is a lot worse than I had last week? George thinks that it is the antenna numbering that is confused so I should fill the data and start from scratch. Started the filler and it crashed aips++. Restarted and it seems to work. Finally got a reasonable calibration (24, 25, and 26 left in for now): My simple script is: include 'flagger.g' fg:=flagger(msfile="/export/home/pegasus/aips++/newfeb19LL.ms"); ok:=fg.auto(); ok:=fg.done(); include 'imager.g' myimager:=imager(filename="/export/home/pegasus/aips++/newfeb19LL.ms" ); ok:=myimager.setjy(fieldid=2, spwid=-1, fluxdensity=1.0); ok:=myimager.setjy(fieldid=1, spwid=-1, fluxdensity=1.0); ok:=myimager.setjy(fieldid=3, spwid=-1, fluxdensity=1.0); ok:=myimager.setjy(fieldid=4, spwid=-1, fluxdensity=-1.0); include 'calibrater.g'; mycalibrater:=calibrater(filename="/export/home/pegasus/aips++/newfeb19LL.ms" ); ok:=mycalibrater.setdata(mode="none" , nchan=1, start=1, step=1, mstart=[value=0.0, unit="km/s" ], mstep=[value=0.0, unit="km/s" ], uvrange=[0, 40] , msselect='(field_id == 2 ) && (antenna1 != antenna2)'); ok:=mycalibrater.setsolve(type="G" , t=300.0, preavg=0.0, phaseonly=F, refant=27, table="newfeb19SN3" , append=T); ok:=mycalibrater.solve(); ok:=mycalibrater.setdata(mode="none" , nchan=1, start=1, step=1, mstart=[value=0.0, unit="km/s" ], mstep=[value=0.0, unit="km/s" ], uvrange=[0, 18] , msselect='(field_id == 4 ) && (antenna1 != antenna2)'); ok:=mycalibrater.setsolve(type="G" , t=300.0, preavg=0.0, phaseonly=F, refant=27, table="newfeb19SN3" , append=T); ok:=mycalibrater.solve(); ok:=mycalibrater.setdata(mode="none" , nchan=1, start=1, step=1, mstart=[value=0.0, unit="km/s" ], mstep=[value=0.0, unit="km/s" ], uvrange=0, msselect='(field_id == 1 || field_id == 3 ) && (antenna1 != antenna2)'); ok:=mycalibrater.setsolve(type="G" , t=300.0, preavg=0.0, phaseonly=F, refant=27, table="newfeb19SN3" , append=T); ok:=mycalibrater.solve(); ok:=mycalibrater.plotcal(plottype="AMP" , tablename="/export/home/pegasus/aips++/newfeb19SN3" , antennas=[], polarization=1, spwids=2, timeslot=1, multiplot=F); ok:=mycalibrater.fluxscale(tablein="/export/home/pegasus/aips++/newfeb19SN3" , tableout="/export/home/pegasus/aips++/newfeb19CL3" , reference="1331+305" , transfer=['0521+166', '0319+415 ', '0713+438'] ); Result: Set the flux density scale using reference calibrators Flux density for 0319+415 (spw=1) is: 21.515 +/- 0.030 Jy Flux density for 0319+415 (spw=2) is: 22.437 +/- 0.065 Jy Flux density for 0521+166 (spw=1) is: 8.575 +/- 0.012 Jy Flux density for 0521+166 (spw=2) is: 8.525 +/- 0.023 Jy Flux density for 0713+438 (spw=1) is: 1.982 +/- 0.003 Jy Flux density for 0713+438 (spw=2) is: 2.048 +/- 0.006 Jy Using George's script to punt some antennas off the end of the arms as is common practice: # (watch the filenames that I've used...) include 'imager.g'; myimager:=imager(filename="newfeb19LL.ms" ); ok:=myimager.setjy(fieldid=-1, spwid=-1, fluxdensity=1.0); ok:=myimager.setjy(fieldid=4, spwid=-1, fluxdensity=-1.0); include 'calibrater.g' mycalibrater:=calibrater(filename="newfeb19LL.ms" ); ok:=mycalibrater.setsolve(type="G" , t=50.0, preavg=0.0, phaseonly=F, refant=13, table="newfeb19LLG.cal" , append=F); # limited 3C286: ok:=mycalibrater.setdata(mode="none" , uvrange=[0, 18], msselect='antenna1!=19 && antenna1!=18 && antenna1!=26 && antenna1!=24 && antenna1!=25 && antenna2!=19 && antenna2!=18 && antenna2!=26 && antenna2!=24 && antenna2!=25 && antenna1!=8 && antenna2!=8 && antenna1!=16 && antenna2!=16 && field_id==4') ; ok:=mycalibrater.solve(); # field 2 w/ uvrange ok:=mycalibrater.setdata(mode="none" , uvrange=[0, 40] , msselect='(field_id == 2 ) && (antenna1 != antenna2) && antenna1!=19 && antenna1!=18 && antenna1!=26 && antenna1!=24 && antenna1!=25 && antenna2!=19 && antenna2!=18 && antenna2!=26 && antenna2!=24 && antenna2!=25 && antenna1!=8 && antenna2!=8 && antenna1!=16 && antenna2!=16'); ok:=mycalibrater.setsolve(type="G" , t=50.0, preavg=0.0, phaseonly=F, refant=13, table="newfeb19LLG.cal" , append=T); ok:=mycalibrater.solve(); # field 1 an 3 with no uvrange: ok:=mycalibrater.setdata(mode="none" , uvrange=[0, 0] , msselect='(field_id == 1 || field_id == 3 ) && (antenna1 != antenna2)'); ok:=mycalibrater.setsolve(type="G" , t=50.0, preavg=0.0, phaseonly=F, refant=13, table="newfeb19LLG.cal" , append=T); ok:=mycalibrater.solve(); ok:=mycalibrater.fluxscale(tablein="newfeb19LLG.cal" , tableout="newfeb19LLG.flx" , reference="1331+305" , transfer=['0521+166', '0319+415 ', '0713+438'] ); ok:=mycalibrater.setapply(type="G" , t=0.0, table="newfeb19LLG.flx" , select=''); ok:=mycalibrater.correct(); Results in: Set the flux density scale using reference calibrators Flux density for 0319+415 (spw=1) is: 21.522 +/- 0.010 Jy Flux density for 0319+415 (spw=2) is: 22.026 +/- 0.018 Jy Flux density for 0521+166 (spw=1) is: 8.466 +/- 0.002 Jy Flux density for 0521+166 (spw=2) is: 8.275 +/- 0.002 Jy Flux density for 0713+438 (spw=1) is: 1.964 +/- 0.001 Jy Flux density for 0713+438 (spw=2) is: 1.991 +/- 0.001 Jy Starting calibrater::correct These are indeed getting pretty close. Better than what George reported to me, but I don't know why: George AIPS aips++ diff (%) 0319+415 21.477 21.469 0.037 IF( 2) 22.080 21.974 0.480 0521+166 8.415 8.444 0.343 IF( 2) 8.234 8.259 0.304 0713+438 1.958 1.959 0.051 IF( 2) 1.993 1.987 0.301 1331+305 15.007 IF( 2) 14.654 So in the end calibration appears to work on this data set in the current weekly version. But note the rather messy antenna specification needed in the calibrator. We really need a GUI to allow for easy antenna selection. --- end of testing. Frustration meter reads: |********** | low high