The pure-quantum self-consistent harmonic approximation (PQSCHA) permits to study a quantum system by means of an effective classical Hamiltonian. In this work the PQSCHA is reformulated in terms of the holomorphic variables connected to a set of bosonic operators. The holomorphic formulation, based on the olomorphic path integral for the Weyl symbol of the density matrix, makes it possible to directly approach general Hamiltonians given in terms of bosonic creation and annihilation operators.