Date - 2001.03.10 Tester - S.T. Myers (AOC) Platform - Linux (kernow) Note: bugs are denoted by lines with prefix: >>>BUG: and queries are denoted by lines with prefix: >>>QUERY: while comments are prefixed by: >>>COMMENT: Some enhancement requests are given after >>>REQUEST: ------------------------------------------------------------------------------- GOAL: Test Athol's PCAL emulation on different test data, now at X band 1. Date - 2001.03.10 Version - Stable (1.5 build 200) Data - /home/sandrock/smyers/VLA/Pol/20010204A/20010204AX.UVF Script - /home/sandrock/smyers/Testing/Mar01/10Mar01/2001-03-10-myers.g include '2001-03-10-myers.g' automap(fitsf='20010204AX.UVF',prefix='test4cm.10mar01.'); Results: 1331+305 spwid= 1 [I=5.191, Q=0, U=0, V=0] Jy, (Perley-Taylor 95) 1331+305 spwid= 2 [I=5.17, Q=0, U=0, V=0] Jy, (Perley-Taylor 95) Flux density for 0713+438 (spw=1) is: 1.240 +/- 0.005 Jy Flux density for 0713+438 (spw=2) is: 1.239 +/- 0.004 Jy Flux density for 0854+201 (spw=1) is: 2.515 +/- 0.013 Jy Flux density for 0854+201 (spw=2) is: 2.499 +/- 0.012 Jy Flux density for 0927+390 (spw=1) is: 12.626 +/- 0.063 Jy Flux density for 0927+390 (spw=2) is: 12.498 +/- 0.060 Jy Flux density for 1159+292 (spw=1) is: 1.649 +/- 0.005 Jy Flux density for 1159+292 (spw=2) is: 1.634 +/- 0.005 Jy Flux density for 1229+020 (spw=1) is: 24.757 +/- 0.089 Jy Flux density for 1229+020 (spw=2) is: 24.372 +/- 0.076 Jy Flux density for 1256-057 (spw=1) is: 24.234 +/- 0.087 Jy Flux density for 1256-057 (spw=2) is: 23.990 +/- 0.077 Jy Flux density for 1310+323 (spw=1) is: 1.733 +/- 0.006 Jy Flux density for 1310+323 (spw=2) is: 1.713 +/- 0.005 Jy This was a clean dataset, and so at first disable the flagging. 0713+438 IF1 I= 1.237 P= 0.003 chi= 137.07 0713+438 IF2 I= 1.236 P= 0.003 chi= 4.59 0854+201 IF1 I= 2.513 P= 0.064 chi= -76.56 0854+201 IF2 I= 2.495 P= 0.067 chi= -74.87 0927+390 IF1 I= 12.596 P= 0.136 chi= 130.56 0927+390 IF2 I= 12.459 P= 0.143 chi= 121.32 1159+292 IF1 I= 1.645 P= 0.037 chi= -131.49 1159+292 IF2 I= 1.631 P= 0.035 chi= -132.81 1229+020 IF1 I= 24.318 P= 1.510 chi= -63.69 1229+020 IF2 I= 23.967 P= 1.451 chi= -62.36 1256-057 IF1 I= 24.176 P= 2.339 chi= 105.41 1256-057 IF2 I= 23.921 P= 2.344 chi= 108.58 1310+323 IF1 I= 1.730 P= 0.050 chi= -52.22 1310+323 IF2 I= 1.711 P= 0.049 chi= -47.96 1331+305 IF1 I= 5.137 P= 0.607 chi= 66.00 1331+305 IF2 I= 5.116 P= 0.610 chi= 66.00 The source 0713+438 is too weakly polarized (0.1%) at X-band to give a reliable angle measurement. We had from the AIPS processing: I Peak mJy P Peak mJy Pol % RLPD deg -------------- -------------- ----------- -------------- 0713+438 1131.51 ± 0.08 1.09 ± 0.63 0.10 ± 0.06 -70.84 ± 43.64 0854+201 2298.93 ± 8.20 60.85 ± 1.37 2.65 ± 0.05 -76.06 ± 0.99 0927+390 11558.10 ± 26.01 121.92 ± 1.88 1.05 ± 0.01 123.33 ± 2.28 1159+292 1512.09 ± 3.36 33.20 ± 0.64 2.20 ± 0.05 -130.08 ± 0.18 1229+020 23620.53 ± 13.21 1422.99 ± 16.25 6.02 ± 0.07 -62.93 ± 1.25 1256-057 22226.38 ± 46.70 2190.59 ± 19.76 9.86 ± 0.07 106.89 ± 1.97 1310+323 1574.11 ± 2.12 46.30 ± 0.50 2.94 ± 0.03 -50.66 ± 1.20 1331+305 5071.61 ± 0.09 605.44 ± 1.81 11.94 ± 0.04 66.00 ± 0.00 Note that the quoted ± values on the RL Phase Difference (2xEVPA) are 1/2 the difference between the phase measured in the 2 IFs. Looks like this worked fine. 2. Date - 2001.03.10 Version - Weekly (1.5 build 256) Data - /home/sandrock/smyers/VLA/Pol/20010204A/20010204AK.UVF Now try the K-band data. My log from AIPS reduction says that AN13 had low and variable gains and had to be deleted. Use fitstoms -> 20010204AK.ms, print summary: Data descriptions: 2 (2 spectral windows and 1 polarization setups) ID Ref.Freq #Chans Resolution TotalBW Correlations 1 22485.1MHz 1 50000 kHz 50000 kHz RR RL LR LL 2 22435.1MHz 1 50000 kHz 50000 kHz RR RL LR LL Antennas: 29 ID= 1-6: 1=VLA:N72, 2=VLA:N16, 3=VLA:N24, 4=VLA:N56, 5=VLA:OUT, 6=VLA:E8, ID= 7-11: 7=VLA:E16, 8=VLA:N8, 9=VLA:E28, 10=VLA:W24, 11=VLA:W32, ID= 12-16: 12=VLA:E20, 13=VLA:W36, 14=VLA:E4, 15=VLA:E24, 16=VLA:W12, ID= 17-21: 17=VLA:E36, 18=VLA:N40, 19=VLA:N48, 20=VLA:E32, 21=VLA:W4, ID= 22-26: 22=VLA:W8, 23=VLA:W16, 24=VLA:W20, 25=VLA:N32, 26=VLA:W28, ID= 27-29: 27=VLA:N64, 28=VLA:E12, 29=VLA:VPT AN8 is indeed in the N:8 pad, so this is the same as the other freqs. Done with ms, and fire up msplot: Want X-Y plot, data vs. time for source 1159+292 Got Field_Id from MS nicely using plot. Try to select RR LL, and Go >>>BUG: got into some hangup 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. CNTL-C exit glish (y/n)? y Crap. Still hung. Have to kill-9 from outside... I don't want to mess with this right now :-( Lets try without editing and see what happens... 3. Date - 2001.03.10 Version - Stable (1.5 build 200) Data - /home/sandrock/smyers/VLA/Pol/20010204A/20010204AK.UVF Script - /home/sandrock/smyers/Testing/Mar01/10Mar01/2001-03-10-myers.g include '2001-03-10-myers.g' automap(fitsf='20010204AK.UVF',prefix='test2cm.10mar01.'); Looks OK: 1331+305 spwid= 1 [I=2.496, Q=0, U=0, V=0] Jy, (Perley-Taylor 95) 1331+305 spwid= 2 [I=2.5, Q=0, U=0, V=0] Jy, (Perley-Taylor 95) Flux density for 0713+438 (spw=1) is: 0.503 +/- 0.002 Jy Flux density for 0713+438 (spw=2) is: 0.502 +/- 0.002 Jy Flux density for 0854+201 (spw=1) is: 5.197 +/- 0.218 Jy Flux density for 0854+201 (spw=2) is: 5.187 +/- 0.218 Jy Flux density for 0927+390 (spw=1) is: 15.752 +/- 0.665 Jy Flux density for 0927+390 (spw=2) is: 15.958 +/- 0.684 Jy Flux density for 1159+292 (spw=1) is: 1.337 +/- 0.001 Jy Flux density for 1159+292 (spw=2) is: 1.338 +/- 0.001 Jy Flux density for 1229+020 (spw=1) is: 14.809 +/- 0.040 Jy Flux density for 1229+020 (spw=2) is: 14.801 +/- 0.037 Jy Flux density for 1256-057 (spw=1) is: 22.794 +/- 0.069 Jy Flux density for 1256-057 (spw=2) is: 22.741 +/- 0.066 Jy Flux density for 1310+323 (spw=1) is: 1.202 +/- 0.001 Jy Flux density for 1310+323 (spw=2) is: 1.195 +/- 0.001 Jy 0713+438 IF1 I= 0.453 P= 0.452 chi= 3.75 0713+438 IF2 I= 0.485 P= 0.016 chi= -83.33 0854+201 IF1 I= 5.140 P= 0.216 chi= -94.06 0854+201 IF2 I= 5.119 P= 0.206 chi= -94.32 0927+390 IF1 I= 15.523 P= 0.227 chi= 144.90 0927+390 IF2 I= 15.788 P= 0.272 chi= 125.47 1159+292 IF1 I= 1.322 P= 0.021 chi= -89.29 1159+292 IF2 I= 1.326 P= 0.019 chi= -87.91 1229+020 IF1 I= 14.743 P= 0.979 chi= -88.77 1229+020 IF2 I= 14.754 P= 1.000 chi= -85.41 1256-057 IF1 I= 22.658 P= 2.519 chi= 104.28 1256-057 IF2 I= 22.636 P= 2.479 chi= 102.60 1310+323 IF1 I= 1.191 P= 0.030 chi= -51.70 1310+323 IF2 I= 1.182 P= 0.026 chi= -23.53 1331+305 IF1 I= 2.492 P= 0.314 chi= 66.00 1331+305 IF2 I= 2.496 P= 0.314 chi= 66.00 Versus AIPS results: I Peak mJy P Peak mJy Pol % RLPD deg -------------- -------------- ----------- -------------- 0713+438 511.29 ± 0.90 1.41 ± 0.23 0.28 ± 0.04 31.76 ± 7.78 0854+201 2668.92 ± 0.65 115.23 ± 1.90 4.32 ± 0.07 -88.22 ± 0.92 0927+390 8181.94 ± 2.52 142.57 ± 7.00 1.74 ± 0.08 124.92 ± 1.45 1159+292 1377.02 ± 1.44 15.70 ± 0.64 1.14 ± 0.05 -100.76 ± 1.93 1229+020 14218.67 ± 12.01 858.19 ± 12.05 6.04 ± 0.08 -86.88 ± 1.81 1256-057 21928.05 ± 34.34 2544.29 ± 1.50 11.60 ± 0.01 97.89 ± 0.21 1310+323 1186.02 ± 0.02 23.72 ± 0.42 2.00 ± 0.04 -46.76 ± 0.94 1331+305 2408.60 ± 0.45 324.17 ± 1.42 13.46 ± 0.06 66.00 ± 0.00 There are significant phase differences between the two. Could be due to editing. Looks like 6 degrees for most sources, and higher on others (eg. 0927+390 IF1, 1310+323 IF2). There are clear problems with the flux density solutions for 0854+201 and 0927+390 also. 4. Date - 2001.03.13 Version - Stable (1.5 build 200) Data - /home/sandrock/smyers/VLA/Pol/20010204A/20010204AK.UVF Script - /home/sandrock/smyers/Testing/Mar01/10Mar01/2001-03-13-myers.g Q-Band definitely needs some editing. Known problems: AN4 has intermittent RR polarization AN13,18 have abnormally low gains and are variable AN1 has low gain and seems fairly stable (probably OK to keep) Need to separate out mstofits and subsequent processing. Change script 2001-03-13-myers.g include '2001-03-13-myers.g' makemsf(fitsf='20010204AQ.UVF',msf='20010204AQ.ms') Try without any editing first automap(msf='20010204AQ.ms',prefix='test1cm.13mar01.') Output: 1331+305 spwid= 1 [I=1.453, Q=0, U=0, V=0] Jy, (Perley-Taylor 95) 1331+305 spwid= 2 [I=1.452, Q=0, U=0, V=0] Jy, (Perley-Taylor 95) Flux density for 0713+438 (spw=1) is: 1.392 +/- 0.086 Jy Flux density for 0713+438 (spw=2) is: 1.536 +/- 0.094 Jy Flux density for 0854+201 (spw=1) is: 291.219 +/- 24.120 Jy Flux density for 0854+201 (spw=2) is: 339.420 +/- 28.305 Jy Flux density for 0927+390 (spw=1) is: 619.032 +/- 51.154 Jy Flux density for 0927+390 (spw=2) is: 685.939 +/- 57.909 Jy Flux density for 1159+292 (spw=1) is: 2.995 +/- 0.159 Jy Flux density for 1159+292 (spw=2) is: 3.280 +/- 0.186 Jy Flux density for 1229+020 (spw=1) is: 13.849 +/- 0.134 Jy Flux density for 1229+020 (spw=2) is: 13.388 +/- 0.093 Jy Flux density for 1256-057 (spw=1) is: 21.750 +/- 0.149 Jy Flux density for 1256-057 (spw=2) is: 21.310 +/- 0.111 Jy Flux density for 1310+323 (spw=1) is: 1.007 +/- 0.004 Jy Flux density for 1310+323 (spw=2) is: 0.988 +/- 0.003 Jy 0713+438 IF1 I= 144.458 P= 84.643 chi= -168.47 0713+438 IF2 I= 144.501 P= 112.562 chi= -151.89 0854+201 IF1 I= -2867.000 P= 451.535 chi= 88.78 0854+201 IF2 I= 5498.365 P= 460.777 chi= -48.48 0927+390 IF1 I= 991.997 P= 29.219 chi= 171.60 0927+390 IF2 I= 654.728 P= 17.418 chi= 152.37 1159+292 IF1 I= 2.943 P= 0.081 chi= -140.59 1159+292 IF2 I= 3.231 P= 0.085 chi= -124.67 1229+020 IF1 I= 13.470 P= 5.706 chi= -93.20 1229+020 IF2 I= 13.276 P= 3.048 chi= 89.32 1256-057 IF1 I= 21.505 P= 13.403 chi= -90.93 1256-057 IF2 I= 46.216 P= 73.004 chi= -39.21 1310+323 IF1 I= 1.039 P= 0.516 chi= -54.95 1310+323 IF2 I= 3.315 P= 9.798 chi= -127.54 1331+305 IF1 I= 1.291 P= 0.333 chi= 66.00 1331+305 IF2 I= 1.312 P= 0.330 chi= 66.00 0713+438 IF1 I= 144.458 P= 84.643 chi= -168.47 0713+438 IF2 I= 144.501 P= 112.562 chi= -151.89 0854+201 IF1 I= -2867.000 P= 451.535 chi= 88.78 0854+201 IF2 I= 5498.365 P= 460.777 chi= -48.48 0927+390 IF1 I= 991.997 P= 29.219 chi= 171.60 0927+390 IF2 I= 654.728 P= 17.418 chi= 152.37 1159+292 IF1 I= 2.943 P= 0.081 chi= -140.59 1159+292 IF2 I= 3.231 P= 0.085 chi= -124.67 1229+020 IF1 I= 13.470 P= 5.706 chi= -93.20 1229+020 IF2 I= 13.276 P= 3.048 chi= 89.32 1256-057 IF1 I= 21.505 P= 13.403 chi= -90.93 1256-057 IF2 I= 46.216 P= 73.004 chi= -39.21 1310+323 IF1 I= 1.039 P= 0.516 chi= -54.95 1310+323 IF2 I= 3.315 P= 9.798 chi= -127.54 1331+305 IF1 I= 1.291 P= 0.333 chi= 66.00 1331+305 IF2 I= 1.312 P= 0.330 chi= 66.00 Versus AIPS results: I Peak mJy P Peak mJy Pol % RLPD deg -------------- -------------- ----------- -------------- 0713+438 266.94 ± 0.99 0.43 ± 0.63 0.16 ± 0.24 -123.72 ± 37.49 0854+201 2641.82 ± 1.66 110.76 ± 2.44 4.19 ± 0.09 -67.15 ± 2.37 0927+390 5636.36 ± 4.86 91.46 ± 3.10 1.62 ± 0.05 134.59 ± 0.74 1159+292 1189.96 ± 2.27 27.79 ± 0.31 2.34 ± 0.02 -99.14 ± 0.04 1229+020 11706.89 ± 7.95 369.76 ± 7.44 3.16 ± 0.06 -100.59 ± 0.56 1256-057 19792.85 ± 3.55 1788.23 ± 7.06 9.03 ± 0.04 99.77 ± 0.54 1310+323 965.03 ± 2.07 10.25 ± 1.02 1.06 ± 0.10 -69.85 ± 3.14 1331+305 1401.21 ± 1.06 188.42 ± 2.65 13.45 ± 0.20 66.00 ± 0.00 There are some severe problems (eg. 0854+201, 0927+390, 1256-057). Get a load of the 300 and 600 Jy fluxes on 0854+201, 0927+390! Some problems due to resolution effects on angle calibrator 1331+305? Will need to pursue editing and flagging. 4. Try to fire up msplot. I thought it was in the ms Module but it seems to be a separate item under the General package (in weekly #2 above it was in ms). I can't even Create an msplot tool. Have to use weekly. Lets wait until next testing. 5. Date - 2001.03.19 Version - Stable (1.5 build 200) Data - /home/sandrock/smyers/VLA/Pol/20010204A/20010204AX.UVF Script - /home/sandrock/smyers/Testing/Mar01/10Mar01/2001-03-13-myers.g Do some timing with X-Band data. In AIPS on sandrock: 7:55 include '2001-03-13-myers.g' makemsf(fitsf='20010204AX.UVF',msf='20010204AX.ms') time = 1:04 automap(msf='20010204AX.ms',prefix='test4cm.19mar01.') time = 17:22 The aips++ was over a factor of 2 slower than the equivalent AIPS. However, I notice that much of the overhead seemed to be in the imaging which did niter=5000 which is overkill. In the AIPS version, NITER=200 is sufficient for these nearly point sources. Lets try niter=500 for aips++ (add clniter variable to script preamble): include '2001-03-19-myers.g' automap(msf='20010204AX.ms',prefix='test4cm.19mar01.') time = 8:59 Thats better. The two operations take essentially the same time (note niter=200 takes 8:29, even closer). Results: 0713+438 IF1 I= 1.238 P= 0.003 chi= 137.01 0713+438 IF2 I= 1.236 P= 0.003 chi= 4.59 0854+201 IF1 I= 2.513 P= 0.064 chi= -76.56 0854+201 IF2 I= 2.495 P= 0.067 chi= -74.87 0927+390 IF1 I= 12.596 P= 0.136 chi= 130.56 0927+390 IF2 I= 12.459 P= 0.143 chi= 121.32 1159+292 IF1 I= 1.645 P= 0.037 chi= -131.49 1159+292 IF2 I= 1.631 P= 0.035 chi= -132.81 1229+020 IF1 I= 24.318 P= 1.510 chi= -63.69 1229+020 IF2 I= 23.967 P= 1.451 chi= -62.36 1256-057 IF1 I= 24.176 P= 2.339 chi= 105.41 1256-057 IF2 I= 23.921 P= 2.344 chi= 108.58 1310+323 IF1 I= 1.730 P= 0.050 chi= -52.22 1310+323 IF2 I= 1.711 P= 0.049 chi= -47.96 1331+305 IF1 I= 5.137 P= 0.607 chi= 66.00 1331+305 IF2 I= 5.116 P= 0.610 chi= 66.00 0713+438 IF2 I= 1.237 P= 0.003 chi= 4.57 0854+201 IF1 I= 2.513 P= 0.064 chi= -76.80 0854+201 IF2 I= 2.496 P= 0.067 chi= -74.90 0927+390 IF1 I= 12.598 P= 0.137 chi= 130.67 0927+390 IF2 I= 12.464 P= 0.143 chi= 121.08 1159+292 IF1 I= 1.646 P= 0.037 chi= -131.60 1159+292 IF2 I= 1.631 P= 0.035 chi= -132.86 1229+020 IF1 I= 24.354 P= 1.512 chi= -63.69 1229+020 IF2 I= 23.990 P= 1.453 chi= -62.40 1256-057 IF1 I= 24.178 P= 2.340 chi= 105.33 1256-057 IF2 I= 23.922 P= 2.344 chi= 108.60 1310+323 IF1 I= 1.731 P= 0.050 chi= -52.54 1310+323 IF2 I= 1.711 P= 0.050 chi= -47.98 1331+305 IF1 I= 5.136 P= 0.607 chi= 66.00 1331+305 IF2 I= 5.116 P= 0.610 chi= 66.00 These are virtually indistinguishable with the niter=500 results (as expected). ------------------------------------------------------------------------------ This file can be found as 2001-03-10-myers.txt. The scripts can be found as 2001-03-10-myers.g, 2001-03-13-myers.g and 2001-03-19-myers.g