MC@NLO

Since recently the Monte Carlo generator program MC@NLO exists [130,131,132]. It provides a calculation of DGLAP evolution equations with NLO Matrix Elements which is matched with the generation of parton showers to describe higher orders which are not contained in the NLO calculations. It has been shown (see fig.7) to successfully describe the beauty production measurements in $p\bar{p}$ collisions with data collected in Run-II at the Tevatron [133] but is yet to be adapted for the processes in $ep$ collisions at HERA.

Figure: CDF $J/\psi$ spectrum from $B$ decays. Two MC@NLO predictions are shown as histograms (see [132] for details). The theory bands represent the FONLL prediction and systematic uncertainties [91].