We employ the recently proposed formalism of the "horizon wave-function" to investigate the emergence of a horizon in models of black holes as Bose-Einstein condensates of gravitons. We start from the Klein-Gordon equation for a massless scalar (toy graviton) field coupled to a static matter current. The (spherically symmetric) classical field reproduces the Newtonian potential generated by the matter source, and the corresponding quantum state is given by a coherent superposition of scalar modes with continuous occupation number. Assuming an attractive self-interaction that allows for bound states, one finds that (approximately) only one mode is allowed, and the system can be confined in a region of the size of the Schwarzschild radius. This radius is then shown to correspond to a proper horizon, by means of the horizon wave-function of the quantum system, with an uncertainty in size naturally related to the expected typical energy of Hawking modes. In particular, this uncertainty decreases for larger black hole mass (with larger number of light scalar quanta), in agreement with semiclassical expectations, a result which does not hold for a single very massive particle. We finally speculate that a phase transition should occur during the gravitational collapse of a star, ideally represented by a static matter current and Newtonian potential, that leads to a black hole, again ideally represented by the condensate of toy gravitons, and suggest an effective order parameter that could be used to investigate this transition. *