Fast Jet Finding in Julia

Abstract: Jet reconstruction remains a critical task in the analysis of data from HEP colliders. We describe in this paper a new, highly performant, Julia package for jet reconstruction, JetReconstruction.jl, which integrates into the growing ecosystem of Julia packages for HEP. With this package users can run sequential reconstruction algorithms for jets. In particular, for LHC events, the Anti-${k}<em>\text{T}$, Cambridge/Aachen and Inclusive-${k}</em>\text{T}$ algorithms can be used. For FCCee studies the use of alternative algorithms such as the Generalised ${k}_\text{T}$ for $e<sup>+e<sup>-$ and Durham are also supported. The performance of the core algorithms is better than Fastjet's C++ implementation, for typical LHC and FCCee events, thanks to the Julia compiler's exploitation of single-instruction-multiple-data (SIMD), as well as ergonomic compact data layouts. The full reconstruction history is made available, allowing inclusive and exclusive jets to be retrieved. The package also provides the means to visualise the reconstruction. Substructure algorithms have been added that allow advanced analysis techniques to be employed. The package can read event data from EDM4hep files and reconstruct jets from these directly, opening the door to FCCee and other future collider studies in Julia.