Hierarchical Markovian algorithm in QCD evolution

A general formalism was introduced for reorganization of QCD evolution equations and derivation of hierarchical solution to DGLAP equation. The hierarchical solution separates two types of parton emissions: the flavour changing emissions and the diagonal ones (bremsstrahlung). It has been shown that both of the achieved processes are of the Markovian type, what gives a nice possibility of Monte Carlo implementation. Hierarchical algorithm has been implemented and crosschecked with standard Markovian algorithm implemented in program EvolFMC.