def beamy(self, feed, freq):
bpat = visibility.cylinder_beam(self._angpos, self.zenith,
self.illumination_y * self.cylinder_width / self.wavelengths[freq])
bm = np.zeros_like(self._angpos)
if self.ortho_pol:
bm[:, 0] = bpat
else:
thatz, phatz = coord.thetaphi_plane_cart(self.zenith)
thatp, phatp = coord.thetaphi_plane_cart(self._angpos)
bm[:, 0] = np.dot(thatp, thatz) * bpat
bm[:, 1] = np.dot(phatp, thatz) * bpat
return bm
def beamy(self, feed, freq):
bpat = visibility.cylinder_beam(
self._angpos, self.zenith,
self.illumination_y * self.cylinder_width / self.wavelengths[freq])
bm = np.zeros_like(self._angpos)
if self.ortho_pol:
bm[:, 0] = bpat
else:
thatz, phatz = coord.thetaphi_plane_cart(self.zenith)
thatp, phatp = coord.thetaphi_plane_cart(self._angpos)
bm[:, 0] = np.dot(thatp, thatz) * bpat
bm[:, 1] = np.dot(phatp, thatz) * bpat
return bm
def beam(self, feed, freq):
"""Beam for a particular feed.
Parameters
----------
feed : integer
Index for the feed.
freq : integer
Index for the frequency.
Returns
-------
beam : np.ndarray
A Healpix map (of size self._nside) of the beam. Potentially
complex.
"""
return visibility.cylinder_beam(self._angpos, self.zenith,
self.cylinder_width / self.wavelengths[freq])
def beam(self, feed, freq):
"""Beam for a particular feed.
Parameters
----------
feed : integer
Index for the feed.
freq : integer
Index for the frequency.
Returns
-------
beam : np.ndarray
A Healpix map (of size self._nside) of the beam. Potentially
complex.
"""
return visibility.cylinder_beam(
self._angpos, self.zenith,
self.cylinder_width / self.wavelengths[freq])