  Double parton distributions can be computed from the perturbative splitting
of one parton into two if the distance between the two observed partons is
small. We develop schemes to take into account quark mass effects in this
computation, and we study these schemes numerically at leading order in the
strong coupling. Furthermore, we investigate in detail the structure of the
next-to-leading order corrections to the splitting kernels that include quark
mass effects.
