The extraction of the spectrum of glueballs in full QCD at physical quark masses is still an open question which requires a coordinated effort from the Lattice community. Given that a number of experiments such as PANDA [1]} and BESIII [2]} are currently looking for such states, the extraction of the associated spectrum and the understanding of how light dynamical quarks affect the spectrum of glueballs is a timely question. Recent reviews on searches for glueballs can be found in the Lattice 2022 plenary presentation by Davide Vadacchino [3]} as well as in Ref. [4]}. Glueball states have large fluctuations and possibly a broad resonance width and so the precise determination of their spectrum requires high statistics of the order of \({{\cal O}} (10-100{\rm K})\)  configurations. Hadron physics performed using current Lattice methodology and machinery, on the other hand, makes use of sets of typically a few hundred configurations. Nevertheless, for a number of measurements on the lattice such as the extraction of renormalisation factors, large sets of configurations are produced and can be, thus, used for exploratory investigations of the glueball spectrum. This work results from the availability of accessing such large sets of configurations.
