Get beam responses
The beam
subcommand can be used to obtain beam responses from any of the
supported beam types. The output format is
tab-separated values (tsv
).
The responses are calculated by moving the zenith angle from 0 to the --max-za
in steps of --step
, then for each of these zenith angles, moving from 0 to
\( 2 \pi \) in steps of --step
for the azimuth. Using a smaller --step
will generate many more responses, so be aware that it might take a while.
If CUDA or HIP is available to you, the --gpu
flag will generate the beam
responses on the GPU, vastly decreasing the time taken.
#!/usr/bin/env python3
import numpy as np
import matplotlib.pyplot as plt
data = np.genfromtxt(fname="beam_responses.tsv", delimiter="\t", skip_header=0)
fig, ax = plt.subplots(1, 2, subplot_kw=dict(projection="polar"))
p = ax[0].scatter(data[:, 0], data[:, 1], c=data[:, 2])
plt.colorbar(p)
p = ax[1].scatter(data[:, 0], data[:, 1], c=np.log10(data[:, 2]))
plt.colorbar(p)
plt.show()