  Recently we constructed a renormalizable field theory up to two loops for the
quasi-static depinning of elastic manifolds in a disordered environment. Here
we explore further properties of the theory. We show how higher correlation
functions of the displacement field can be computed. Drastic simplifications
occur, unveiling much simpler diagrammatic rules than anticipated. This is
applied to the universal scaled width-distribution. The expansion in
d=4-epsilon predicts that the scaled distribution coincides to the lowest
orders with the one for a Gaussian theory with propagator G(q)=1/q^(d+2 \zeta),
zeta being the roughness exponent. The deviations from this Gaussian result are
small and involve higher correlation functions, which are computed here for
different boundary conditions. Other universal quantities are defined and
evaluated: We perform a general analysis of the stability of the fixed point.
We find that the correction-to-scaling exponent is omega=-epsilon and not
-epsilon/3 as used in the analysis of some simulations. A more detailed study
of the upper critical dimension is given, where the roughness of interfaces
grows as a power of a logarithm instead of a pure power.
