A novel method for calculating Bose-Einstein correlation functions with Coulomb final-state interaction

TL;DR

This paper introduces a new computational method to accurately calculate Bose-Einstein correlation functions, including Coulomb interactions, especially for Lévystable sources, enhancing analysis capabilities in femtoscopy.

Contribution

The paper presents a novel, exact integral calculation method for correlation functions with Coulomb interaction, applicable to any source with a known Fourier transform, notably improving analysis of Lévystable sources.

Findings

The new method is more robust than previous techniques.

It allows wider parameter range investigation.

A ready-to-use software package is provided.

Abstract

Measurements of Bose-Einstein correlations played a crucial role in the discovery and the subsequent detailed exploration of the Quark-Gluon-Plasma (QGP) created in high-energy collisions of heavy nuclei. Such measurements gave rise to femtoscopy, a flourishing sub-field of high-energy physics, and there are important new directions to explore and discoveries to be made in the near future. One of these important current topics is the precise investigation of the shape of the correlation functions utilizing L\'evy-stable sources. In this paper, we present a novel method of calculating the shape of the two-particle correlation functions, including the Coulomb final-state interaction. This method relies on an exact calculation of a large part of the necessary integrals of the Coulomb wave function and can be utilized to calculate the correlation function for any source function with an easily accessible Fourier transform. In this way, it is eminently applicable to L\'evy-stable source functions. In this particular case, we find that the new method is more robust and allows the investigation of a wider parameter range than the previously utilized techniques. We present an easily applicable software package that is ready to use in experimental studies.