These processes represent the production of a Standard Model Higgs boson that decays into a pair of bottom quarks, in association with a further bottom quark. The initial state at lowest order is a bottom quark and a gluon. The calculation may be performed at NLO, although radiation from the bottom quarks in the Higgs decay is not included. For more details on this calculation, please see Ref. [1].
For this process, the matrix elements are divided up into a number of different sub-processes, so the user must sum over these after performing more runs than usual. At lowest order one can proceed as normal, using nproc=136. For a NLO calculation, the sequence of runs is as follows:
Run nproc=136 with part=virt and part=real (or, both at the same time using part=tota);
Run nproc=137 with part=real.
The sum of these yields the cross-section with one identified b-quark in the final state. To calculate the contribution with two b-quarks in the final state, one should use nproc=138 with part=real.
When removebr is .true., the Higgs boson does not decay.
nplotter_auto.f is the default plotting routine.
[1] J.M. Campbell, R.K. Ellis, F. Maltoni and S. Willenbrock, Higgs-Boson Production in Association with a Single Bottom Quark, Phys. Rev. D 67 (2003) 095002 [hep-ph/0204093].