  We present a calculation of Pi, D and B mesons production at RHIC and LHC
energies based upon the KKT model of gluon saturation. We discuss dependence of
the nuclear modification factor on rapidity and transverse momentum.
