  We present an analysis of the most precise set of HERA data within the color
dipole formalism, by using an analytical gluon density, based on the
double-logarithm approximation of the DGLAP equations in the asymptotic limit
of the scaling variable, $\sigma=\log{(1/x)}\log{(\log{(Q^2/Q_
0^2)})}\rightarrow \infty$. Fits to data, including charm and bottom quarks are
performed and demonstrate the efficiency of the model in describing the reduced
cross section, $\sigma_{r}$, in the wide range $Q^2:(1.5,500)$ GeV$^2$ for two
dipole models including parton saturation effects. We also give predictions to
$F_{2}^{c\bar{c}}$ , $F_{2}^{b\bar{b}}$ and $F_{L}$, all describing the data
reasonably well in the range $Q^2:(2.5,120)$ GeV$^2$. Total cross sections of
exclusive photoproduction of $J/\psi$ and $\rho$ are also calculated and
successfully compared to HERA data and recent measurements at LHCb.
