$\eta$ and $\eta'$ masses and decay constants from lattice QCD with $N_f=2+1+1$ quark flavours