Hyperbolic surfaces with sublinearly many systoles that fill

For any $\varepsilon>0$, we construct a closed hyperbolic surface of genus $g=g(\varepsilon)$ with a set of at most $\varepsilon g$ systoles that fill, meaning that each component of the complement of their union is contractible. This surface is also a critical point of index at most $\varepsilon g$ for the systole function, disproving the lower bound of $2g-1$ posited by Schmutz Schaller.