Date - 2001.12.06 Tester - S.T. Myers (AOC) Platform - Linux (sandrock) 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: Sort out polarization calibration problems in K and Q band data. 1. Date - 2001.12.06 Version - Weekly (1.7 build 033) Directory - /home/sandrock/smyers/Testing/2001-12-06/ Data - 20010204AX.UVF Script - 2001-12-06-myers.g Make script using 2001-03-20-myers.g - include '2001-12-06-myers.g'; Change script to set up record par with needed parameters: par.pwdir := '/home/sandrock/smyers/Testing/2001-12-06/'; par.prefix := '20010204.xband.'; # Prefix for files generated par.fitsf := '20010204AX.UVF'; # Name of UVFITS file par.msf := '20010204AX.ms'; # Name of ms par.abscal := '1331+305'; # Name of primary flux cal par.absid := 6; # FIELD_ID of primary flux cal par.angcal := '1331+305'; # Polarization angle calibrator par.absang := 66.00; # R-L Phase Difference of PCal par.polcal := ['1159+292']; # Name of d-term calibrator par.polid :=1; # FIELD_ID of d-term calibrator par.refsrc := ['0713+438', '0854+201', '0927+390', '1159+292', '1229+020', '1256-057', '1310+323']; # Names of secondary sources par.nfld := 8; # Number of targets par.fieldlist := ['0713+438', '0854+201', '0927+390', '1159+292', '1229+020', '1256-057', '1310+323', '1331+305']; # Names of targets par.gsolint := 60.0; # solution interval (secs) for gain par.grefant := 8; # reference antenna for gain sols par.clniter := 200; # Number of clean iterations desired par.clcell := '0.4arcsec'; # image cell size (Xband) - makemsf(par) Note that fitstoms no longer prints summary as default - add m.summary() to script. MeasurementSet Name: 20010204AX.ms MS Version 2 Observer: AC524 Project: Observation: VLA(29 antennas) Data records: 127406 Total integration time = 8080 seconds Observed from 04-Feb-2001/08:22:50 to 04-Feb-2001/10:37:30 Fields: 8 ID Name Right Ascension Declination Epoch 1 1159+292 11:59:31.83 +29.14.43.83 J2000 2 0713+438 07:13:38.16 +43.49.17.20 J2000 3 0927+390 09:27:03.01 +39.02.20.85 J2000 4 0854+201 08:54:48.87 +20.06.30.64 J2000 5 1310+323 13:10:28.66 +32.20.43.78 J2000 6 1331+305 13:31:08.29 +30.30.32.96 J2000 7 1229+020 12:29:06.70 +02.03.08.60 J2000 8 1256-057 12:56:11.17 -05.47.21.52 J2000 Data descriptions: 2 (2 spectral windows and 1 polarization setups) ID Ref.Freq #Chans Resolution TotalBW Correlations 1 8435.1 MHz 1 50000 kHz 50000 kHz RR RL LR LL 2 8485.1 MHz 1 50000 kHz 50000 kHz RR RL LR LL Antennas: 26 ID= 1-7: 1=VLA:N72, 2=VLA:N16, 3=VLA:N24, 4=VLA:N56, 6=VLA:E8, 7=VLA:E16, ID= 8-12: 8=VLA:N8, 9=VLA:E28, 10=VLA:W24, 11=VLA:W32, 12=VLA:E20, ID= 13-18: 13=VLA:W36, 14=VLA:E4, 16=VLA:W12, 17=VLA:E36, 18=VLA:N40, ID= 19-23: 19=VLA:N48, 20=VLA:E32, 21=VLA:W4, 22=VLA:W8, 23=VLA:W16, ID= 24-26: 24=VLA:W20, 25=VLA:N32, 26=VLA:W28 Tables(rows): (-1 = table absent) MAIN(127406) ANTENNA(29) DATA_DESCRIPTION(2) DOPPLER(-1) FEED(29) FIELD(8) FLAG_CMD(0) FREQ_OFFSET(-1) HISTORY(523) OBSERVATION(1) POINTING(348) POLARIZATION(1) PROCESSOR(0) SOURCE(-1) SPECTRAL_WINDOW(2) STATE(0) SYSCAL(-1) WEATHER(-1) - automap(par) 1331+305 spwid= 1 [I=5.216, Q=0, U=0, V=0] Jy, (Perley-Taylor 95) 1331+305 spwid= 2 [I=5.195, Q=0, U=0, V=0] Jy, (Perley-Taylor 95) Flux density for 0713+438 (spw=1) is: 1.161 +/- 0.002 Jy Flux density for 0713+438 (spw=2) is: 1.158 +/- 0.002 Jy Flux density for 0854+201 (spw=1) is: 2.344 +/- 0.011 Jy Flux density for 0854+201 (spw=2) is: 2.361 +/- 0.011 Jy Flux density for 0927+390 (spw=1) is: 11.764 +/- 0.053 Jy Flux density for 0927+390 (spw=2) is: 11.807 +/- 0.054 Jy Flux density for 1159+292 (spw=1) is: 1.546 +/- 0.001 Jy Flux density for 1159+292 (spw=2) is: 1.553 +/- 0.001 Jy Flux density for 1229+020 (spw=1) is: 23.653 +/- 0.026 Jy Flux density for 1229+020 (spw=2) is: 23.718 +/- 0.027 Jy Flux density for 1256-057 (spw=1) is: 22.520 +/- 0.024 Jy Flux density for 1256-057 (spw=2) is: 22.613 +/- 0.025 Jy Flux density for 1310+323 (spw=1) is: 1.605 +/- 0.001 Jy Flux density for 1310+323 (spw=2) is: 1.610 +/- 0.001 Jy 0713+438 IF1 I= 1.159 P= 0.011 chi= -167.67 0713+438 IF2 I= 1.156 P= 0.018 chi= 140.53 0854+201 IF1 I= 2.343 P= 0.071 chi= -144.01 0854+201 IF2 I= 2.360 P= 0.088 chi= -153.43 0927+390 IF1 I= 11.744 P= 0.055 chi= 21.99 0927+390 IF2 I= 11.782 P= 0.148 chi= 60.72 1159+292 IF1 I= 1.543 P= 0.032 chi= -104.59 1159+292 IF2 I= 1.551 P= 0.027 chi= -100.86 1229+020 IF1 I= 23.445 P= 1.316 chi= -171.03 1229+020 IF2 I= 23.571 P= 1.177 chi= -167.24 1256-057 IF1 I= 22.483 P= 2.328 chi= 24.43 1256-057 IF2 I= 22.577 P= 2.436 chi= 17.20 1310+323 IF1 I= 1.604 P= 0.058 chi= 178.10 1310+323 IF2 I= 1.609 P= 0.059 chi= 163.04 1331+305 IF1 I= 5.224 P= 0.614 chi= 66.00 1331+305 IF2 I= 5.205 P= 0.655 chi= 66.00 Timing: 5m20s (expected speedup on new sandrock versus old cpu) Compare to (from 2001-03-10-myers.txt): 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 solved fluxes for 1331+305 are higher (is this really Perley-Taylor 95 or later?) but the fluxes seem closer to what AIPS gave before. Check with the aips++ group on what flux scale is used. Note that the AIPS run did use 1995.2 = 5.192,5.170 so I'll bet aips++ uses 1999.2. >>>BUG: Asked Athol about this - indeed the default is 1999.2 despite what is reported. For consistency with AIPS, use the 1995.2 values if possible. How is this done (nothing in Users Guide for imager.setjy). >>>BUG: Athols says that this is in the development version (as per a defect of Greg's. More worrisome is that the RLPDs seem broke now. Asked Athol, and he said that the sign of U has been fixed and -U is no longer needed. 2. Fix and rerun - include '2001-12-06-myers.g'; - makemsf(par) - automap(par) 0713+438 IF1 I= 1.159 P= 0.003 chi= 138.52 0713+438 IF2 I= 1.156 P= 0.003 chi= 7.87 0854+201 IF1 I= 2.343 P= 0.060 chi= -76.89 0854+201 IF2 I= 2.360 P= 0.063 chi= -74.77 0927+390 IF1 I= 11.747 P= 0.127 chi= 130.43 0927+390 IF2 I= 11.782 P= 0.135 chi= 121.35 1159+292 IF1 I= 1.544 P= 0.035 chi= -131.86 1159+292 IF2 I= 1.551 P= 0.033 chi= -132.64 1229+020 IF1 I= 23.458 P= 1.458 chi= -64.01 1229+020 IF2 I= 23.596 P= 1.431 chi= -62.58 1256-057 IF1 I= 22.471 P= 2.177 chi= 105.14 1256-057 IF2 I= 22.560 P= 2.213 chi= 108.66 1310+323 IF1 I= 1.604 P= 0.047 chi= -52.68 1310+323 IF2 I= 1.609 P= 0.047 chi= -47.92 1331+305 IF1 I= 5.225 P= 0.616 chi= 66.00 1331+305 IF2 I= 5.202 P= 0.619 chi= 66.00 This now reproduces the 2001-03-20-myers.g results. The script is 2001-12-06-myers.g 3. This script seems to work well, and the automap function may be useful to general users who can fill in the par record. What is needed to use aips++ for the polarization calibration program? o Full editing capability in msplot, with differences from average as in TVPLT especially. o Noise estimates in script o Output of results to table (probably via an output record from automap) o Output of d-term solutions and gains to table (probably via record) Try to incorporate the last of these next, and test the new msplot for editing. ------------------------------------------------------------------------------ Testing logs archived at http://www.aoc.nrao.edu/~smyers/aips++/testlogs/ This file can be found as 2001-12-06-myers.txt. The scripts can be found as 2001-12-06-myers.g