Kernel Density Estimators for Axisymmetric Particle Beams

Bright beams are commonly represented by sampled data in the numerical algorithms used to simulate their properties. However, in these calculations and the analyses of their outputs, the beam's density is sometimes required and must be calculated from the samples. Axisymmetric beams, which possess a rotational symmetry and are naturally expressed in polar coordinates, pose a particular challenge to density estimators. The area element in polar coordinates shrinks as the radius becomes small, and weighting the samples to account for their reduced frequency may cause unwelcome artifacts. In this work, we derive analytical expressions for two kernel density estimators, which solve these problems in the spatial coordinates and in the transverse phase space. We show how the kernels can be found by averaging the Gaussian kernel in Cartesian coordinates over the polar angle and demonstrate their use on test problems. These results show that particle beam symmetries can be taken advantage of in density estimation while avoiding artifacts.