We develop an analytical and numerical framework based on the disentanglement approach to study the ground states of many-body quantum spins systems. In this approach, observables are expressed as functional integrals over scalar fields, where the relevant measure is the Wiener measure. We identify the leading contribution to these integrals, given by the saddle point field configuration. Analytically, this can be used to develop an exact field-theoretical expansion of the functional integrals, performed by means of appropriate Feynman rules. The expansion can be truncated to the desired order to obtain approximate analytical results for ground state expectation values. Numerically, the saddle point configuration can be used to compute physical observables by means of an exact importance sampling scheme. We illustrate our methods by considering the quantum Ising model in 1, 2 and 3 spatial dimensions. Our analytical and numerical results are applicable to a broad class of many-body quantum spin systems, bridging concepts from quantum lattice models, continuum field theory, and classical stochastic processes.