We examine carefully bottomonia hadroproduction in proton colliders, especially focusing on the LHC, as a way of probing the gluon density in protons. To this end we develop some previous work, getting quantitative predictions and concluding that our proposal can be useful to perform consistency checks of the parameterization sets of different parton distribution functions.