  In this paper, the soft gluon resummation effect in double heavy quarkonium
production at the LHC is studied. By applying the transverse momentum dependent
factorization formalism, the large logarithms introduced by the small total
transverse momentum of heavy quarkonium pair final state system, are resummed
to all orders in the expansion of the strong interaction coupling at the
Next-to-Leading Logarithm accuracy. We also compare our result with the LHC
data. We find that the distribution shape predicted by resummation calculation
is consistent with experimental data very well. Since this process is mainly
initiated by gluon fusion, it supplies us an important channel to study the
gluon parton property in hadrons.
