  We extend the recent argument by Ding and Zhuang from nearest-neighbor to
long-range interactions and prove the phase transition in the class of
ferromagnetic random field Ising models. Our proof combines a generalization of
Fr\"ohlich-Spencer contours to the multidimensional setting, proposed by two of
us, with the coarse-graining procedure introduced by Fisher, Fr\"ohlich and
Spencer. The result shows that the Ding-Zhuang strategy is also useful for
interactions $J_{xy}=|x-y|^{- \alpha}$ when $\alpha > d$ in dimension $d\geq 3$
if we have a suitable system of contours. We can consider i.i.d. random fields
with Gaussian or Bernoulli distributions. Our main result is an alternative
proof that does not use the Renormalization Group Method (RGM), since Bricmont
and Kupiainen claimed that the RGM should also work on this generality.
