The Fourier expansion of Eisenstein series $E_k$ $(k \ge 4)$, which are modular forms, as well as the quasimodular $E_2$, involves powers-of-divisors $\sigma_{k-1}(n) = \sum_{d|n} d^{k-1}$.
Is there some generalization of a modular form (quasimodular, mock modular, etc.) where you can find some similar Fourier expansion of the form $$f(z) = 1 + C \sum_{n=0}^{\infty} \sigma_0(n) q^n, \; q = e^{2\pi i z}$$ involving just the number-of-divisors function $\sigma_0(n)$? I haven't been able to find this.