  The advent of ultra short high intensity lasers has paved the way to new and
promising, yet challenging, areas of research in the laser-plasma interaction
physics. The success of constructing petawatt femtosecond lasers, for instance
the Apollon laser in France, will help understanding and designing future
particle accelerators and next generation of light sources. Achieving this goal
intrinsically relies on the combination between experiments and massively
parallel simulations. So far, Particle-In-Cell (PIC) codes have been the
ultimate tool to accurately describe the laser-plasma interaction especially in
the field of Laser WakeField Acceleration (LWFA) . Nevertheless, the numerical
modelling of laser plasma accelerators in 3D can be a very challenging task.
This is due to the large dispersity between the scales involved in this
process. In order to make such simulations feasible with a significant speed
up, we need to use reduced numerical models which simplify the problem while
retaining a high fidelity. Among these models, Fourier field decomposition in
azimuthal modes for the cylindrical geometry is a promising reduced model
especially for physical problems that have close to cylindrical symmetry which
is the case in LWFA. This geometry has been implemented in the open-source code
Smilei in Finite Difference Time Domain (FDTD) discretization scheme for the
Maxwell solver. In this paper we will study the case of a realistic laser
measurement from Apollon facility, the ability of this method to describe it
correctly and the determination of the necessary number of modes for this
purpose. We will also show the importance of higher modes inclusion in the case
of realistic laser profiles to insure fidelity in simulation.
