Coulomb interaction for L\'evy sources

During the study of Bose-Einstein correlations in heavy ion collisions one has to take into account the final state interactions, amongst them the Coulomb interaction playing a prominent role for charged particles. In some cases measurements have shown that the correlation function can be best described by L\'evy sources, and three dimensional measurements have indicated the possibility of deviation from spherical symmetry. Therefore, one would like to study the Coulomb interaction for non-spherical L\'evy sources. We resort to numerical methods which are most commonly used in order to take into account the Coulomb interaction such measurements. Here, we utilize the Metropolis-Hastings algorithm. The symmetric L\'evy distribution that describes the source can be characterized by three L\'evy scale parameters and the L\'evy exponent. We investigate the roles of these parameters in the correlation function. We show the results for the Bose-Einstein correlation functions for ellipsoidal L\'evy sources with Coulomb interaction. We also compare our results with previous ways to treat the Coulomb interaction in the presence of L\'evy sources.