  Valuable theoretical predictions of nuclear dipole excitations in the whole
chart are of great interest for different nuclear applications, including in
particular nuclear astrophysics. Here we present large-scale calculations of
the $E1$ $\gamma$-ray strength function obtained in the framework of the
axially-symmetric deformed QRPA based on the finite-range Gogny force. This
approach is applied to even-even nuclei, the strength function for odd nuclei
being derived by interpolation. The convergence with respect to the adopted
number of harmonic oscillator shells and the cut-off energy introduced in the
2-quasiparticle (2-$qp$) excitation space is analyzed. The calculations
performed with two different Gogny interactions, namely D1S and D1M, are
compared. A systematic energy shift of the $E1$ strength is found for D1M
relative to D1S, leading to a lower energy centroid and a smaller
energy-weighted sum rule for D1M. When comparing with experimental
photoabsorption data, the Gogny-QRPA predictions are found to overestimate the
giant dipole energy by typically $\sim$2 MeV. Despite the microscopic nature of
our self-consistent Hartree-Fock-Bogoliubov plus QRPA calculation, some
phenomenological corrections need to be included to take into account the
effects beyond the standard 2-$qp$ QRPA excitations and the coupling between
the single-particle and low-lying collective phonon degrees of freedom. For
this purpose, three prescriptions of folding procedure are considered and
adjusted to reproduce experimental photoabsorption data at best. All of them
are shown to lead to rather similar predictions of the $E1$ strength, both at
low energies and for exotic neutron-rich nuclei. Predictions of $\gamma$-ray
strength functions and Maxwellian-averaged neutron capture rates for the whole
Sn isotopic chain are also discussed and compared with previous theoretical
calculations.
