Up to now, we have calculated the inclusive cross sections for heavy quarks production at hadron colliders. These calculations were performed using approach based on the unintegrated parton distributions functions. I have tested some of the models in photoproduction and in hadroproduction. For the ccbar and bbar production at high-energies the gluon-gluon fusion is assumed to be the dominant mechanism. This process was calculated in the NLO collinear as well as in the kt-factorisation approaches. Now, I study production of the charm and bottom quarks for following subprocesses: gg->Q anti-Q, γg->Q anti-Q, gγ->Q anti-Q, γγ->Q anti-Q. I used MRST_QED (Martin, Roberts, Stirling, Thorne) parton distributions. I would like to show detail analyzed that other processes ignored so far should be carefully evaluated. The number of potential contributions is not small.