Theoretical beam shapes and ALMA measured beams

There are two sets of model beam images for DA and DV antennas. TICRA supplied the electric field simulation data for the ones containing EFP in the filename. According to Dirk, the ones containing VP in the filename are "for a baseline", and are apparently derived from products of the EFP files. See Dirk's presentation at the 2012 CALIM workshop for more information. Unfortunately, all the results below inadvertently used the VP files. However, the VP files differ from the EFP noticeably only at 84 and 116 GHz.

The following table was produced by running getfwhm, and (for comparison) CASA's imfit on the TICRA model beams for the DA and DV antennas using the attached script (after having converted the complex images to amplitude with immath using this script). Note that the TICRA files for the PM antennas are numerically identical to the DV antennas, as they were not independently modeled. So I did not bother fitting those. The FWHM factor and effective taper were then computed with au.primaryBeamArcsec and au.effectiveTaper which account for the ~0.5% effect of the central obstruction using au.centralObstructionFactor. The cubic formula from Baars (2007) was used. The CASA image files are stored here: /lustre/cv/thunter/alma/beamMaps/ticra/*.square.

TICRA Filename
(azmin_aznominal_azmax _elmin_elnominal_elmax _freqmin_freqnominal_freqmax)
Stokes FWHM (arcsec) from au.getfwhm FWHM factor inferred taper (dB) Geometric mean (arcsec) from imfit FWHM factor inferred taper (dB) Major axis (arcsec) Minor axis (arcsec) Position angle (deg) Major axis uncertainty (arcsec) Minor axis uncertainty (arcsec) P.A. uncertainty (deg)
Median over all bands, polarizations and antenna types: 1.14258 11.12   1.0936 6.54  
DA__0_0_360_0_45_90_84_84_92_GHz XX 72.79 1.187 16.11 69.50 1.133 10.18 69.59 69.41 10.34 0.83 0.83 6.37
DA__0_0_360_0_45_90_84_84_92_GHz YY 72.79 1.187 16.11 69.48 1.133 10.16 69.72 69.25 9.92 0.84 0.84 2.33
DA__0_0_360_0_45_90_92_100_108_GHz XX 60.76 1.179 15.13 57.97 1.125 9.41 58.07 57.86 9.63 0.71 0.71 4.57
DA__0_0_360_0_45_90_92_100_108_GHz YY 60.76 1.179 15.13 57.95 1.125 9.39 58.14 57.77 9.66 0.71 0.71 2.56
DA__0_0_360_0_45_90_108_116_116_GHz XX 50.16 1.129 9.81 47.68 1.073 4.64 47.77 47.59 10.02 0.72 0.72 5.32
DA__0_0_360_0_45_90_108_116_116_GHz YY 50.16 1.129 9.81 47.65 1.073 4.57 48.62 46.70 99.69 0.72 0.72 0.53
DA__0_0_360_0_45_90_125_125_134.5_GHz XX 48.02 1.165 13.44 45.33 1.100 7.10 45.42 45.25 80.47 0.59 0.59 4.46
DA__0_0_360_0_45_90_125_125_134.5_GHz YY 48.02 1.165 13.44 45.33 1.100 7.09 45.42 45.25 80.47 0.59 0.59 4.49
DA__0_0_360_0_45_90_134.5_144_153.5_GHz XX 40.58 1.134 10.30 39.28 1.098 6.91 39.33 39.22 80.39 0.51 0.51 6.03
DA__0_0_360_0_45_90_134.5_144_153.5_GHz YY 40.58 1.134 10.30 39.28 1.098 6.91 39.33 39.22 80.40 0.51 0.51 6.06
DA__0_0_360_0_45_90_153.5_163_163_GHz XX 35.84 1.134 10.26 34.66 1.096 6.78 34.70 34.61 79.95 0.45 0.45 6.45
DA__0_0_360_0_45_90_153.5_163_163_GHz YY 35.84 1.134 10.26 34.66 1.096 6.78 34.70 34.61 79.94 0.45 0.45 6.51
DA__0_0_360_0_45_90_211_211_227_GHz XX 27.83 1.140 10.84 26.69 1.093 6.46 26.76 26.61 135.00 0.35 0.35 3.20
DA__0_0_360_0_45_90_211_211_227_GHz YY 27.83 1.140 10.84 26.69 1.093 6.46 26.76 26.61 135.00 0.35 0.35 3.18
DA__0_0_360_0_45_90_227_243_259_GHz XX 23.84 1.124 9.34 23.12 1.090 6.23 23.21 23.03 135.00 0.30 0.30 2.16
DA__0_0_360_0_45_90_227_243_259_GHz YY 23.84 1.124 9.34 23.12 1.090 6.23 23.21 23.03 135.00 0.30 0.30 2.15
DA__0_0_360_0_45_90_259_275_275_GHz XX 21.58 1.152 12.02 20.37 1.087 5.94 20.48 20.26 135.00 0.26 0.26 1.59
DA__0_0_360_0_45_90_259_275_275_GHz YY 21.58 1.152 12.02 20.37 1.087 5.94 20.48 20.26 135.00 0.26 0.26 1.59
DA__0_0_360_0_45_90_275_275_299.5_GHz XX 21.41 1.142 11.09 20.45 1.092 6.35 20.48 20.43 91.71 0.27 0.27 6.66
DA__0_0_360_0_45_90_275_275_299.5_GHz YY 21.41 1.142 11.09 20.44 1.091 6.26 20.47 20.40 94.05 0.27 0.27 5.05
DA__0_0_360_0_45_90_299.5_324_348.5_GHz XX 18.18 1.143 11.15 17.28 1.086 5.87 17.31 17.25 92.34 0.23 0.23 4.71
DA__0_0_360_0_45_90_299.5_324_348.5_GHz YY 18.18 1.143 11.15 17.27 1.086 5.81 17.30 17.23 93.21 0.23 0.23 4.14
DA__0_0_360_0_45_90_348.5_373_373_GHz XX 15.87 1.149 11.72 14.96 1.083 5.53 15.00 14.91 90.07 0.20 0.20 2.88
DA__0_0_360_0_45_90_348.5_373_373_GHz YY 15.87 1.149 11.72 14.95 1.082 5.48 15.00 14.90 90.73 0.20 0.20 2.63
DA__0_0_360_0_45_90_602_602_631.5_GHz XX 9.99 1.167 13.69 9.37 1.094 6.59 9.39 9.34 0.85 0.12 0.12 3.55
DA__0_0_360_0_45_90_602_602_631.5_GHz YY 9.99 1.167 13.69 9.36 1.094 6.58 9.39 9.34 1.27 0.12 0.12 3.51
DA__0_0_360_0_45_90_631.5_661_690.5_GHz XX 8.89 1.141 10.93 8.51 1.091 6.33 8.53 8.49 0.55 0.11 0.11 3.36
DA__0_0_360_0_45_90_631.5_661_690.5_GHz YY 8.89 1.141 10.93 8.51 1.091 6.32 8.53 8.49 1.03 0.11 0.11 3.36
DA__0_0_360_0_45_90_690.5_720_720_GHz XX 8.17 1.142 11.03 7.80 1.090 6.20 7.82 7.78 178.30 0.10 0.10 2.94
DA__0_0_360_0_45_90_690.5_720_720_GHz YY 8.17 1.142 11.03 7.80 1.090 6.20 7.82 7.78 178.63 0.10 0.10 2.95
-----------------------------------------------------------------------------------------------------------------------
DV__0_0_360_0_45_90_84_84_92_GHz XX 73.03 1.190 16.64 69.74 1.137 10.57 69.84 69.65 10.96 0.82 0.82 5.88
DV__0_0_360_0_45_90_84_84_92_GHz YY 73.03 1.190 16.64 69.73 1.137 10.54 69.97 69.49 10.59 0.82 0.82 2.27
DV__0_0_360_0_45_90_92_100_108_GHz XX 60.76 1.179 15.13 58.13 1.128 9.70 58.22 58.03 9.98 0.70 0.70 4.68
DV__0_0_360_0_45_90_92_100_108_GHz YY 60.76 1.179 15.13 58.11 1.128 9.68 58.29 57.93 10.12 0.70 0.70 2.60
DV__0_0_360_0_45_90_108_116_116_GHz XX 50.16 1.129 9.81 47.86 1.077 5.03 47.95 47.78 10.76 0.70 0.70 5.33
DV__0_0_360_0_45_90_108_116_116_GHz YY 50.16 1.129 9.81 47.83 1.077 4.96 48.78 46.90 100.28 0.70 0.70 0.51
DV__0_0_360_0_45_90_125_125_134.5_GHz XX 48.02 1.165 13.44 45.49 1.103 7.44 45.58 45.40 79.39 0.58 0.58 4.23
DV__0_0_360_0_45_90_125_125_134.5_GHz YY 48.02 1.165 13.44 45.49 1.103 7.44 45.58 45.40 79.39 0.58 0.58 4.26
DV__0_0_360_0_45_90_134.5_144_153.5_GHz XX 40.58 1.134 10.30 39.36 1.100 7.12 39.41 39.31 79.78 0.50 0.50 6.37
DV__0_0_360_0_45_90_134.5_144_153.5_GHz YY 40.58 1.134 10.30 39.36 1.100 7.12 39.41 39.31 79.80 0.50 0.50 6.37
DV__0_0_360_0_45_90_153.5_163_163_GHz XX 36.75 1.163 13.20 34.74 1.099 7.03 34.79 34.70 79.93 0.44 0.44 6.64
DV__0_0_360_0_45_90_153.5_163_163_GHz YY 36.75 1.163 13.20 34.74 1.099 7.03 34.79 34.70 79.93 0.44 0.44 6.65
DV__0_0_360_0_45_90_211_211_227_GHz XX 27.83 1.140 10.84 26.80 1.097 6.89 26.89 26.72 134.95 0.34 0.34 2.64
DV__0_0_360_0_45_90_211_211_227_GHz YY 27.83 1.140 10.84 26.80 1.097 6.89 26.89 26.72 134.95 0.34 0.34 2.63
DV__0_0_360_0_45_90_227_243_259_GHz XX 23.84 1.124 9.34 23.19 1.093 6.53 23.29 23.08 135.37 0.29 0.29 1.84
DV__0_0_360_0_45_90_227_243_259_GHz YY 23.84 1.124 9.34 23.19 1.093 6.53 23.29 23.08 135.37 0.29 0.29 1.84
DV__0_0_360_0_45_90_259_275_275_GHz XX 21.58 1.152 12.02 20.42 1.090 6.16 20.53 20.30 136.56 0.25 0.25 1.46
DV__0_0_360_0_45_90_259_275_275_GHz YY 21.58 1.152 12.02 20.42 1.090 6.16 20.53 20.30 136.57 0.25 0.25 1.46
DV__0_0_360_0_45_90_275_275_299.5_GHz XX 21.41 1.142 11.09 20.51 1.095 6.64 20.54 20.48 91.47 0.27 0.27 5.88
DV__0_0_360_0_45_90_275_275_299.5_GHz YY 21.41 1.142 11.09 20.49 1.094 6.55 20.53 20.46 93.79 0.27 0.27 4.59
DV__0_0_360_0_45_90_299.5_324_348.5_GHz XX 18.23 1.146 11.50 17.34 1.090 6.25 17.38 17.30 91.89 0.22 0.22 3.82
DV__0_0_360_0_45_90_299.5_324_348.5_GHz YY 18.23 1.146 11.50 17.33 1.090 6.19 17.38 17.29 92.78 0.22 0.22 3.48
DV__0_0_360_0_45_90_348.5_373_373_GHz XX 15.87 1.149 11.72 15.01 1.086 5.87 15.05 14.97 90.10 0.19 0.19 3.09
DV__0_0_360_0_45_90_348.5_373_373_GHz YY 15.87 1.149 11.72 15.00 1.086 5.82 15.05 14.95 90.85 0.19 0.19 2.81
DV__0_0_360_0_45_90_602_602_631.5_GHz XX 9.99 1.167 13.69 9.41 1.099 7.05 9.44 9.38 0.36 0.12 0.12 2.63
DV__0_0_360_0_45_90_602_602_631.5_GHz YY 9.99 1.167 13.69 9.41 1.099 7.04 9.44 9.38 0.70 0.12 0.12 2.61
DV__0_0_360_0_45_90_631.5_661_690.5_GHz XX 8.89 1.141 10.93 8.53 1.095 6.65 8.57 8.50 179.96 0.11 0.11 2.33
DV__0_0_360_0_45_90_631.5_661_690.5_GHz YY 8.89 1.141 10.93 8.53 1.095 6.64 8.56 8.50 0.30 0.11 0.11 2.33
DV__0_0_360_0_45_90_690.5_720_720_GHz XX 8.17 1.142 11.03 7.82 1.093 6.46 7.85 7.79 179.89 0.10 0.10 2.12
DV__0_0_360_0_45_90_690.5_720_720_GHz YY 8.17 1.142 11.03 7.82 1.093 6.46 7.85 7.79 0.14 0.10 0.10 2.12


Generic theoretical illumination profiles and their resulting farfield beam profiles

The following plots were made with scipy.fftpack and scipy.signal. One dimensional version:

beams1d.png

Two dimensional version: beams2d.png

ALMA measured beams (data obtained from Stuartt back in April 2012)

The following are measured beams in Band 3 on one of each of the 12m antenna models: almabeams2d.png

Now compared to theory (on the same color scale): almavstheory.png

-- ToddHunter - 2012-10-23

Topic attachments
I Attachment Action Size Date Who Comment
DA.pdfpdf DA.pdf manage 1 MB 2014-02-21 - 14:13 ToddHunter overlay of 1D cuts through the VP and EFP (squared) images of DA antennas
DV.pdfpdf DV.pdf manage 1 MB 2014-02-21 - 14:13 ToddHunter overlay of 1D cuts through the VP and EFP (squared) images of DV antennas
all_square.py.txttxt all_square.py.txt manage 4 K 2013-08-26 - 18:01 ToddHunter version with au.getfwhm and median calculation
all_vp.py.txttxt all_vp.py.txt manage 2 K 2013-08-02 - 17:33 ToddHunter  
alma-imaging-dpetry-calim2012.pdfpdf alma-imaging-dpetry-calim2012.pdf manage 1 MB 2014-03-07 - 10:31 ToddHunter  
almabeams2d.pngpng almabeams2d.png manage 126 K 2012-10-23 - 10:32 ToddHunter  
almabeams2d.py.txttxt almabeams2d.py.txt manage 1 K 2012-10-23 - 10:45 ToddHunter  
almavstheory.pngpng almavstheory.png manage 175 K 2012-10-23 - 10:44 ToddHunter  
almavstheory.py.txttxt almavstheory.py.txt manage 3 K 2012-10-23 - 10:44 ToddHunter  
beams1d.pngpng beams1d.png manage 86 K 2012-10-23 - 10:31 ToddHunter  
beams1d.py.txttxt beams1d.py.txt manage 3 K 2012-10-23 - 10:44 ToddHunter  
beams2d.pngpng beams2d.png manage 176 K 2012-10-23 - 10:31 ToddHunter  
beams2d.py.txttxt beams2d.py.txt manage 4 K 2012-10-23 - 10:45 ToddHunter  
make_amplitude.py.txttxt make_amplitude.py.txt manage 3 K 2013-08-02 - 14:57 ToddHunter converts TICRA complex images to amplitude using immath
make_square.py.txttxt make_square.py.txt manage 3 K 2013-08-02 - 22:38 ToddHunter  
uid___A002_X2d0501_X1a2-DA42-H-USB.fitsfits uid___A002_X2d0501_X1a2-DA42-H-USB.fits manage 11 K 2012-10-23 - 10:45 ToddHunter  
uid___A002_X2d0501_X1a2-DV03-H-USB.fitsfits uid___A002_X2d0501_X1a2-DV03-H-USB.fits manage 11 K 2012-10-23 - 10:45 ToddHunter  
uid___A002_X2d0501_X1a2-PM04-H-USB.fitsfits uid___A002_X2d0501_X1a2-PM04-H-USB.fits manage 11 K 2012-10-23 - 10:46 ToddHunter  
vpefp.py.txttxt vpefp.py.txt manage 1 K 2014-02-21 - 14:12 ToddHunter script to overlay cuts of the VP and EFP images
This topic: Main > TWikiUsers > ToddHunter > HunterCASA > ALMABeamsVsTheoretical
Topic revision: 2014-03-07, ToddHunter
This site is powered by FoswikiCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding NRAO Public Wiki? Send feedback