The Monte Carlo generator CASCADE [104] contains an implementation of the factorization approach (see section 2.5) using the CCFM evolution equation [65]. In CASCADE the direct heavy quark production processes and are implemented using off-shell matrix elements convoluted with -unintegrated parton distributions in the proton. In CASCADE, higher order parton emissions based on the CCFM [65] evolution equations are matched to matrix elements in which the incoming parton can be off-shell. This implies that the parton kinematics are treated correctly. For the final state parton showers PYTHIA [108] is used. The unintegrated gluon density in the proton has been determined in fits to the inclusive structure function data, in the range and GeV. With this input CASCADE has been able to correctly reproduce the quark production cross-sections in -collisions at the Tevatron [128].