In the framework of the NRQCD factorization formalism, we calculate the decay rate for the process Υ(1S) → ccgg to the next-to-leading order (NLO) in the relative velocity v of the b quark in the bottomonium rest frame. We also study the momentum distributions of the charm quark and the charmed-hadron in the decay. The momentum distribution of the charmed-hadron is obtained by convolving the charm quark momentum distribution with a fragmentation function of the charm quark into the hadron. In addition, we fit the nonperturbative NRQCD matrix element v 2 Υ through comparing the theoretical prediction with the measurement from the BaBar collaboration for the decay rate of Υ(1S) → D * + X. In return, taking this matrix element as an input parameter, we predict the decay rates as well as the momentum distributions for a collection of charmed-hadrons in the process Υ(1S) → ccgg → hX.