Results are presented from four CDF analyses involving heavy quark production in proton-antiproton collisions at center of mass energy 1.96 TeV. The shapes of b-jets are found to be broader than inclusive predictions and broader than both PYTHIA and HERWIG defaults. A measurement of the production cross section for (psi)(2S) is consistent with Run 1 results and with theoretical predictions associated with parton distribution function energy dependence. The inclusive b-jet production cross section is also consistent with theoretical predictions over six orders of magnitude. The b(bar b) differential production cross section is compared to several theoretical models and found to be best described by MC(at)NLO + JIMMY.