Presented is the search for anomalous Higgs boson decays into two axion-like particles (ALPs) using the full Run 2 data set of 140 fb$^{-1}$ of proton-proton collisions at a centre-of-mass energy of 13 TeV recorded by the ATLAS experiment. The ALPs are assumed to decay into two photons, providing sensitivity to recently proposed models that could explain the $(g-2)_\mu$ discrepancy. This analysis covers an ALP mass range from 100 MeV to 62 GeV and ALP-photon couplings in the range $10^{-5}\, \text{TeV}^{-1}
Cut-flow and number of events in a region with $110 \text{ GeV} <m_{aa}<130 \text{ GeV}$ for selected signal samples
Overview of relative systematic uncertainties on the signal normalization for selected couplings and mass points in the dominant categories. All photon related uncertainties are summarized under "standard photon", while all customised photon related uncertainties (e.g. displaced photons, NN-based photon IDs) are summarized under "custom photon".
Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow 4\gamma)$ at 95% CL as a function of the axion mass and for ALP-photon coupling $C_{a\gamma\gamma}=1$.