We calculate cross section for production of χ c pairs in proton-proton collisions. The cross section for the gg → χ cJ 1 χ cJ 2 is considerably smaller (especially for χ c1 χ c1 final state) than that obtained recently in the k T -factorization approach. We calculate therefore next-to-leading order contributions with χ c pair and one extra associated (mini-)jet. We find these contributions to be much larger than those for the 2 → 2 contribution. Especially the emission of a leading gluon (carrying a large momentum fraction of one of the incoming gluons) are important. These emissions in the k T -factorization approach are absorbed into the initial state unintegrated gluon distributions. A smaller contribution to the cross section comes from the production of central gluons emitted with rapidities between the χ c -mesons. They do lead, however, to an enhancement of the χ c -pair production at large rapidity distance between the mesons. Our present study explains the size of the cross section for the χ c pair production obtained previously in the k T -factorization approach. Several differential distributions are presented.