(* ::Package:: *) (* what is given here are the contents of the brackets in eq. 56, 61, 64 and T_F*g_{q q'}^(2) in eq. 63*) matchingCoefficientFragmentingQuarkJetMeasureqOrder2 = (170*CA*DiracDelta[sDivmu2])/27 + 34*CF*DiracDelta[sDivmu2] - (8*CA*Pi^2*DiracDelta[sDivmu2])/9 - (5*CF*Pi^2*DiracDelta[sDivmu2])/2 - (88*TF*DiracDelta[sDivmu2])/9 - (76*nf*TF*DiracDelta[sDivmu2])/27 - (4*Pi^2*TF*DiracDelta[sDivmu2])/3 + (4*nf*Pi^2*TF*DiracDelta[sDivmu2])/ 9 + (5*CF*Pi^2*DiracDelta[sDivmu2])/(6*(-1 + z)) - (214*TF*DiracDelta[sDivmu2])/(27*z) + (638*CA*z*DiracDelta[sDivmu2])/ 27 - 34*CF*z*DiracDelta[sDivmu2] - (2*CA*Pi^2*z*DiracDelta[sDivmu2])/ 9 + (13*CF*Pi^2*z*DiracDelta[sDivmu2])/6 - (32*TF*z*DiracDelta[sDivmu2])/9 - (148*nf*TF*z*DiracDelta[sDivmu2])/27 - 2*Pi^2*TF*z*DiracDelta[sDivmu2] + (4*nf*Pi^2*TF*z*DiracDelta[sDivmu2])/ 9 + (574*TF*z^2*DiracDelta[sDivmu2])/27 - (17*CF*Pi^2*z^2*DiracDelta[sDivmu2])/(6*(-1 + z)) + (2428*CA*DiracDelta[sDivmu2]*DiracDelta[-1 + z])/81 - (67*CA*Pi^2*DiracDelta[sDivmu2]*DiracDelta[-1 + z])/18 - (CA*Pi^4*DiracDelta[sDivmu2]*DiracDelta[-1 + z])/9 + (7*CF*Pi^4*DiracDelta[sDivmu2]*DiracDelta[-1 + z])/30 - (656*nf*TF*DiracDelta[sDivmu2]*DiracDelta[-1 + z])/81 + (10*nf*Pi^2*TF*DiracDelta[sDivmu2]*DiracDelta[-1 + z])/9 + 2*CF*DiracDelta[sDivmu2]*HPL[{0}, z] - 40*TF*DiracDelta[sDivmu2]* HPL[{0}, z] + (4*Pi^2*TF*DiracDelta[sDivmu2]*HPL[{0}, z])/3 - (5*CF*Pi^2*DiracDelta[sDivmu2]*HPL[{0}, z])/(3*(1 - z)) - (214*CA*DiracDelta[sDivmu2]*HPL[{0}, z])/(9*(-1 + z)) + (16*CF*DiracDelta[sDivmu2]*HPL[{0}, z])/(-1 + z) + (2*CA*Pi^2*DiracDelta[sDivmu2]*HPL[{0}, z])/(3*(-1 + z)) + (CF*Pi^2*DiracDelta[sDivmu2]*HPL[{0}, z])/(3*(-1 + z)) + (44*nf*TF*DiracDelta[sDivmu2]*HPL[{0}, z])/(9*(-1 + z)) - (28*TF*DiracDelta[sDivmu2]*HPL[{0}, z])/(9*z) + 18*CF*z*DiracDelta[sDivmu2]*HPL[{0}, z] - 36*TF*z*DiracDelta[sDivmu2]* HPL[{0}, z] + (4*Pi^2*TF*z*DiracDelta[sDivmu2]*HPL[{0}, z])/3 + (86*CA*z*DiracDelta[sDivmu2]*HPL[{0}, z])/(3*(-1 + z)) - (12*CF*z*DiracDelta[sDivmu2]*HPL[{0}, z])/(-1 + z) - (16*nf*TF*z*DiracDelta[sDivmu2]*HPL[{0}, z])/(3*(-1 + z)) - (124*TF*z^2*DiracDelta[sDivmu2]*HPL[{0}, z])/9 - (7*CF*Pi^2*z^2*DiracDelta[sDivmu2]*HPL[{0}, z])/(3*(1 - z)) - (160*CA*z^2*DiracDelta[sDivmu2]*HPL[{0}, z])/(9*(-1 + z)) + (8*CF*z^2*DiracDelta[sDivmu2]*HPL[{0}, z])/(-1 + z) + (2*CA*Pi^2*z^2*DiracDelta[sDivmu2]*HPL[{0}, z])/(3*(-1 + z)) + (CF*Pi^2*z^2*DiracDelta[sDivmu2]*HPL[{0}, z])/(-1 + z) + (44*nf*TF*z^2*DiracDelta[sDivmu2]*HPL[{0}, z])/(9*(-1 + z)) - (22*CA*DiracDelta[sDivmu2]*HPL[{1}, z])/9 + 22*CF*DiracDelta[sDivmu2]* HPL[{1}, z] + (8*CF*Pi^2*DiracDelta[sDivmu2]*HPL[{1}, z])/3 + (80*TF*DiracDelta[sDivmu2]*HPL[{1}, z])/3 - (16*nf*TF*DiracDelta[sDivmu2]*HPL[{1}, z])/9 + (2*CA*Pi^2*DiracDelta[sDivmu2]*HPL[{1}, z])/(-1 + z) + (14*CF*Pi^2*DiracDelta[sDivmu2]*HPL[{1}, z])/(3*(-1 + z)) + (28*TF*DiracDelta[sDivmu2]*HPL[{1}, z])/(9*z) + (308*CA*z*DiracDelta[sDivmu2]*HPL[{1}, z])/9 - 24*CF*z*DiracDelta[sDivmu2]*HPL[{1}, z] + (8*CF*Pi^2*z*DiracDelta[sDivmu2]*HPL[{1}, z])/3 - (56*TF*z*DiracDelta[sDivmu2]*HPL[{1}, z])/3 - (64*nf*TF*z*DiracDelta[sDivmu2]*HPL[{1}, z])/9 - (100*TF*z^2*DiracDelta[sDivmu2]*HPL[{1}, z])/9 + (2*CA*Pi^2*z^2*DiracDelta[sDivmu2]*HPL[{1}, z])/(3*(-1 + z)) - (22*CF*Pi^2*z^2*DiracDelta[sDivmu2]*HPL[{1}, z])/(3*(-1 + z)) - 4*CA*DiracDelta[sDivmu2]*HPL[{2}, z] + 4*TF*DiracDelta[sDivmu2]* HPL[{2}, z] + (6*CF*DiracDelta[sDivmu2]*HPL[{2}, z])/(1 - z) - (6*CF*DiracDelta[sDivmu2]*HPL[{2}, z])/(-1 + z) - (16*TF*DiracDelta[sDivmu2]*HPL[{2}, z])/(3*z) - 4*CA*z*DiracDelta[sDivmu2]*HPL[{2}, z] + 16*TF*z*DiracDelta[sDivmu2]* HPL[{2}, z] - (4*CF*z*DiracDelta[sDivmu2]*HPL[{2}, z])/(1 - z) + (8*CF*z*DiracDelta[sDivmu2]*HPL[{2}, z])/(-1 + z) + (16*TF*z^2*DiracDelta[sDivmu2]*HPL[{2}, z])/3 + (10*CF*z^2*DiracDelta[sDivmu2]*HPL[{2}, z])/(1 - z) + (22*CF*z^2*DiracDelta[sDivmu2]*HPL[{2}, z])/(-1 + z) + 4*CF*DiracDelta[sDivmu2]*HPL[{3}, z] - 16*TF*DiracDelta[sDivmu2]* HPL[{3}, z] + 4*CF*z*DiracDelta[sDivmu2]*HPL[{3}, z] - 16*TF*z*DiracDelta[sDivmu2]*HPL[{3}, z] - (24*CF*z^2*DiracDelta[sDivmu2]*HPL[{3}, z])/(-1 + z) - 2*TF*DiracDelta[sDivmu2]*HPL[{0, 0}, z] + (2*CF*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(1 - z) - (35*CA*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(3*(-1 + z)) + (26*CF*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(-1 + z) + (4*nf*TF*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(3*(-1 + z)) + (16*TF*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(3*z) - 10*TF*z*DiracDelta[sDivmu2]*HPL[{0, 0}, z] - (12*CF*z*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(1 - z) + (4*CA*z*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(-1 + z) - (40*CF*z*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(-1 + z) - (2*CF*z^2*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(1 - z) + (CA*z^2*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(3*(-1 + z)) - (16*CF*z^2*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(-1 + z) + (4*nf*TF*z^2*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(3*(-1 + z)) - (10*CF*DiracDelta[sDivmu2]*HPL[{1, 0}, z])/(1 - z) + (2*CF*DiracDelta[sDivmu2]*HPL[{1, 0}, z])/(-1 + z) + (32*CF*z*DiracDelta[sDivmu2]*HPL[{1, 0}, z])/(1 - z) + (8*CF*z*DiracDelta[sDivmu2]*HPL[{1, 0}, z])/(-1 + z) - (10*CF*z^2*DiracDelta[sDivmu2]*HPL[{1, 0}, z])/(1 - z) + (2*CF*z^2*DiracDelta[sDivmu2]*HPL[{1, 0}, z])/(-1 + z) + (22*CA*DiracDelta[sDivmu2]*HPL[{1, 1}, z])/3 + 4*TF*DiracDelta[sDivmu2]*HPL[{1, 1}, z] - (8*nf*TF*DiracDelta[sDivmu2]*HPL[{1, 1}, z])/3 + (16*TF*DiracDelta[sDivmu2]*HPL[{1, 1}, z])/(3*z) + (22*CA*z*DiracDelta[sDivmu2]*HPL[{1, 1}, z])/3 - 4*TF*z*DiracDelta[sDivmu2]*HPL[{1, 1}, z] - (8*nf*TF*z*DiracDelta[sDivmu2]*HPL[{1, 1}, z])/3 - (16*TF*z^2*DiracDelta[sDivmu2]*HPL[{1, 1}, z])/3 - (8*CF*DiracDelta[sDivmu2]*HPL[{1, 2}, z])/(1 - z) - (4*CA*DiracDelta[sDivmu2]*HPL[{1, 2}, z])/(-1 + z) - (16*CF*DiracDelta[sDivmu2]*HPL[{1, 2}, z])/(-1 + z) - (8*CF*z^2*DiracDelta[sDivmu2]*HPL[{1, 2}, z])/(1 - z) - (4*CA*z^2*DiracDelta[sDivmu2]*HPL[{1, 2}, z])/(-1 + z) - (16*CF*z^2*DiracDelta[sDivmu2]*HPL[{1, 2}, z])/(-1 + z) - (12*CF*DiracDelta[sDivmu2]*HPL[{2, 0}, z])/(1 - z) + (4*CA*DiracDelta[sDivmu2]*HPL[{2, 0}, z])/(-1 + z) - (20*CF*z^2*DiracDelta[sDivmu2]*HPL[{2, 0}, z])/(1 - z) + (4*CA*z^2*DiracDelta[sDivmu2]*HPL[{2, 0}, z])/(-1 + z) - (8*CF*z^2*DiracDelta[sDivmu2]*HPL[{2, 0}, z])/(-1 + z) + 8*TF*DiracDelta[sDivmu2]*HPL[{2, 1}, z] - (12*CF*DiracDelta[sDivmu2]*HPL[{2, 1}, z])/(1 - z) - (28*CF*DiracDelta[sDivmu2]*HPL[{2, 1}, z])/(-1 + z) + 8*TF*z*DiracDelta[sDivmu2]*HPL[{2, 1}, z] - (4*CF*z^2*DiracDelta[sDivmu2]*HPL[{2, 1}, z])/(1 - z) - (12*CF*z^2*DiracDelta[sDivmu2]*HPL[{2, 1}, z])/(-1 + z) + 20*TF*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z] + (4*CF*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z])/(1 - z) - (10*CA*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z])/(-1 + z) + (34*CF*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z])/(-1 + z) + 20*TF*z*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z] + (12*CF*z^2*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z])/(1 - z) - (10*CA*z^2*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z])/(-1 + z) + (62*CF*z^2*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z])/(-1 + z) - (4*CF*DiracDelta[sDivmu2]*HPL[{1, 0, 0}, z])/(-1 + z) - (4*CF*z^2*DiracDelta[sDivmu2]*HPL[{1, 0, 0}, z])/(-1 + z) - (8*CF*DiracDelta[sDivmu2]*HPL[{1, 1, 0}, z])/(1 - z) + (4*CA*DiracDelta[sDivmu2]*HPL[{1, 1, 0}, z])/(-1 + z) - (24*CF*DiracDelta[sDivmu2]*HPL[{1, 1, 0}, z])/(-1 + z) - (8*CF*z^2*DiracDelta[sDivmu2]*HPL[{1, 1, 0}, z])/(1 - z) + (4*CA*z^2*DiracDelta[sDivmu2]*HPL[{1, 1, 0}, z])/(-1 + z) - (24*CF*z^2*DiracDelta[sDivmu2]*HPL[{1, 1, 0}, z])/(-1 + z) + 24*CF*DiracDelta[sDivmu2]*HPL[{1, 1, 1}, z] + 24*CF*z*DiracDelta[sDivmu2]*HPL[{1, 1, 1}, z] + (40*CA*PlusDistribution[sDivmu2^(-1)])/9 - 24*CF*PlusDistribution[sDivmu2^(-1)] + (2*CA*Pi^2*PlusDistribution[sDivmu2^(-1)])/3 + 4*CF*Pi^2*PlusDistribution[sDivmu2^(-1)] - (80*TF*PlusDistribution[sDivmu2^(-1)])/3 + (16*nf*TF*PlusDistribution[sDivmu2^(-1)])/9 + (4*Pi^2*TF*PlusDistribution[sDivmu2^(-1)])/3 - (28*TF*PlusDistribution[sDivmu2^(-1)])/(9*z) - (308*CA*z*PlusDistribution[sDivmu2^(-1)])/9 + 24*CF*z*PlusDistribution[sDivmu2^(-1)] + (2*CA*Pi^2*z*PlusDistribution[sDivmu2^(-1)])/3 + 4*CF*Pi^2*z*PlusDistribution[sDivmu2^(-1)] + (56*TF*z*PlusDistribution[sDivmu2^(-1)])/3 + (64*nf*TF*z*PlusDistribution[sDivmu2^(-1)])/9 + (4*Pi^2*TF*z*PlusDistribution[sDivmu2^(-1)])/3 + (100*TF*z^2*PlusDistribution[sDivmu2^(-1)])/9 - (808*CA*DiracDelta[-1 + z]*PlusDistribution[sDivmu2^(-1)])/27 + (22*CA*Pi^2*DiracDelta[-1 + z]*PlusDistribution[sDivmu2^(-1)])/9 + (224*nf*TF*DiracDelta[-1 + z]*PlusDistribution[sDivmu2^(-1)])/27 - (8*nf*Pi^2*TF*DiracDelta[-1 + z]*PlusDistribution[sDivmu2^(-1)])/9 + 4*CA*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)] - 4*TF*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)] + (2*CF*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/(1 - z) + (18*CF*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/(-1 + z) + (16*TF*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/(3*z) + 4*CA*z*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)] - 16*TF*z*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)] - (12*CF*z*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/(1 - z) - (32*CF*z*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/(-1 + z) - (16*TF*z^2*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/3 - (2*CF*z^2*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/(1 - z) - (10*CF*z^2*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/(-1 + z) - (22*CA*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)])/3 - 4*TF*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)] + (8*nf*TF*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)])/3 - (16*TF*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)])/(3*z) - (22*CA*z*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)])/3 + 4*TF*z*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)] + (8*nf*TF*z*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)])/3 + (16*TF*z^2*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)])/3 + 4*CF*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)] - 8*TF*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)] + (16*CF*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)])/(-1 + z) + 4*CF*z*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)] - 8*TF*z*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)] + 16*TF*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)] + (4*CF*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(1 - z) - (4*CA*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(-1 + z) + (20*CF*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(-1 + z) + 16*TF*z*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)] + (12*CF*z^2*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(1 - z) - (4*CA*z^2*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(-1 + z) + (44*CF*z^2*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(-1 + z) + (8*CF*HPL[{1, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(-1 + z) + (8*CF*z^2*HPL[{1, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(-1 + z) - 24*CF*HPL[{1, 1}, z]*PlusDistribution[sDivmu2^(-1)] - 24*CF*z*HPL[{1, 1}, z]*PlusDistribution[sDivmu2^(-1)] - (808*CA*DiracDelta[sDivmu2]*PlusDistribution[(1 - z)^(-1)])/27 + (22*CA*Pi^2*DiracDelta[sDivmu2]*PlusDistribution[(1 - z)^(-1)])/9 + (224*nf*TF*DiracDelta[sDivmu2]*PlusDistribution[(1 - z)^(-1)])/27 - (8*nf*Pi^2*TF*DiracDelta[sDivmu2]*PlusDistribution[(1 - z)^(-1)])/9 + (268*CA*PlusDistribution[sDivmu2^(-1)]*PlusDistribution[(1 - z)^(-1)])/ 9 - (4*CA*Pi^2*PlusDistribution[sDivmu2^(-1)]*PlusDistribution[ (1 - z)^(-1)])/3 - (20*CF*Pi^2*PlusDistribution[sDivmu2^(-1)]* PlusDistribution[(1 - z)^(-1)])/3 - (80*nf*TF*PlusDistribution[sDivmu2^(-1)]*PlusDistribution[(1 - z)^(-1)])/ 9 + (22*CA*PlusDistribution[Log[sDivmu2]/sDivmu2])/3 + 4*TF*PlusDistribution[Log[sDivmu2]/sDivmu2] - (8*nf*TF*PlusDistribution[Log[sDivmu2]/sDivmu2])/3 + (16*TF*PlusDistribution[Log[sDivmu2]/sDivmu2])/(3*z) + (22*CA*z*PlusDistribution[Log[sDivmu2]/sDivmu2])/3 - 4*TF*z*PlusDistribution[Log[sDivmu2]/sDivmu2] - (8*nf*TF*z*PlusDistribution[Log[sDivmu2]/sDivmu2])/3 - (16*TF*z^2*PlusDistribution[Log[sDivmu2]/sDivmu2])/3 + (268*CA*DiracDelta[-1 + z]*PlusDistribution[Log[sDivmu2]/sDivmu2])/9 - (4*CA*Pi^2*DiracDelta[-1 + z]*PlusDistribution[Log[sDivmu2]/sDivmu2])/ 3 - (20*CF*Pi^2*DiracDelta[-1 + z]*PlusDistribution[ Log[sDivmu2]/sDivmu2])/3 - (80*nf*TF*DiracDelta[-1 + z]* PlusDistribution[Log[sDivmu2]/sDivmu2])/9 + 8*TF*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] + (4*CF*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2])/(1 - z) + 8*TF*z*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] + (12*CF*z^2*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2])/(1 - z) + (16*CF*z^2*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2])/ (-1 + z) + 24*CF*HPL[{1}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] + 24*CF*z*HPL[{1}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] - (44*CA*PlusDistribution[(1 - z)^(-1)]*PlusDistribution[ Log[sDivmu2]/sDivmu2])/3 + (16*nf*TF*PlusDistribution[(1 - z)^(-1)]* PlusDistribution[Log[sDivmu2]/sDivmu2])/3 - 12*CF*PlusDistribution[Log[sDivmu2]^2/sDivmu2] - 12*CF*z*PlusDistribution[Log[sDivmu2]^2/sDivmu2] - (22*CA*DiracDelta[-1 + z]*PlusDistribution[Log[sDivmu2]^2/sDivmu2])/3 + (8*nf*TF*DiracDelta[-1 + z]*PlusDistribution[Log[sDivmu2]^2/sDivmu2])/ 3 + 24*CF*PlusDistribution[(1 - z)^(-1)]*PlusDistribution[ Log[sDivmu2]^2/sDivmu2] + 8*CF*DiracDelta[-1 + z]* PlusDistribution[Log[sDivmu2]^3/sDivmu2] + (268*CA*DiracDelta[sDivmu2]*PlusDistribution[Log[1 - z]/(1 - z)])/9 - (4*CA*Pi^2*DiracDelta[sDivmu2]*PlusDistribution[Log[1 - z]/(1 - z)])/3 - (20*CF*Pi^2*DiracDelta[sDivmu2]*PlusDistribution[Log[1 - z]/(1 - z)])/ 3 - (80*nf*TF*DiracDelta[sDivmu2]*PlusDistribution[Log[1 - z]/(1 - z)])/ 9 - (44*CA*PlusDistribution[sDivmu2^(-1)]*PlusDistribution[ Log[1 - z]/(1 - z)])/3 + (16*nf*TF*PlusDistribution[sDivmu2^(-1)]* PlusDistribution[Log[1 - z]/(1 - z)])/3 + 48*CF*PlusDistribution[Log[sDivmu2]/sDivmu2]*PlusDistribution[ Log[1 - z]/(1 - z)] - (22*CA*DiracDelta[sDivmu2]* PlusDistribution[Log[1 - z]^2/(1 - z)])/3 + (8*nf*TF*DiracDelta[sDivmu2]*PlusDistribution[Log[1 - z]^2/(1 - z)])/3 + 24*CF*PlusDistribution[sDivmu2^(-1)]*PlusDistribution[ Log[1 - z]^2/(1 - z)] + 8*CF*DiracDelta[sDivmu2]* PlusDistribution[Log[1 - z]^3/(1 - z)] + 8*TF*DiracDelta[sDivmu2]* Zeta[3] - (40*CF*DiracDelta[sDivmu2]*Zeta[3])/(3*(1 - z)) + (10*CA*DiracDelta[sDivmu2]*Zeta[3])/(-1 + z) + (152*CF*DiracDelta[sDivmu2]*Zeta[3])/(3*(-1 + z)) + 8*TF*z*DiracDelta[sDivmu2]*Zeta[3] - (152*CF*z^2*DiracDelta[sDivmu2]*Zeta[3])/(3*(1 - z)) - (18*CA*z^2*DiracDelta[sDivmu2]*Zeta[3])/(-1 + z) - (32*CF*z^2*DiracDelta[sDivmu2]*Zeta[3])/(3*(-1 + z)) - (110*CA*DiracDelta[sDivmu2]*DiracDelta[-1 + z]*Zeta[3])/9 + (40*nf*TF*DiracDelta[sDivmu2]*DiracDelta[-1 + z]*Zeta[3])/9 + 28*CA*DiracDelta[-1 + z]*PlusDistribution[sDivmu2^(-1)]*Zeta[3] + 32*CF*DiracDelta[-1 + z]*PlusDistribution[sDivmu2^(-1)]*Zeta[3] + 28*CA*DiracDelta[sDivmu2]*PlusDistribution[(1 - z)^(-1)]*Zeta[3] + (160*CF*DiracDelta[sDivmu2]*PlusDistribution[(1 - z)^(-1)]*Zeta[3])/3 - (64*CF*DiracDelta[sDivmu2]*PlusDistribution[(1 - z)^(-1)]*Zeta[3])/ (3*(1 - z)) + (64*CF*z*DiracDelta[sDivmu2]*PlusDistribution[ (1 - z)^(-1)]*Zeta[3])/(3*(1 - z)) PlusDistribution /: PlusDistribution::usage = "\!\(\*SubsuperscriptBox[\(\ \[Integral]\), \(0\), \(1\)]\)f[x] \ PlusDistribution[\!\(\*FractionBox[\(g[x]\), \ \(x\)]\)]\[DifferentialD]x\[Equal]\!\(\*SubsuperscriptBox[\(\[Integral]\), \ \(0\), \(1\)]\)\!\(\*FractionBox[\(\((f[x] - f[0])\)\\ g[x]\), \ \(x\)]\)\[DifferentialD]x" matchingCoefficientFragmentingQuarkJetMeasuregOrder2 = (-482*CA*DiracDelta[sDivmu2])/9 - 24*CF*DiracDelta[sDivmu2] - 4*CA*Pi^2*DiracDelta[sDivmu2] - (16*CF*Pi^2*DiracDelta[sDivmu2])/3 + (2338*CA*DiracDelta[sDivmu2])/(27*z) - (2*CA*Pi^2*DiracDelta[sDivmu2])/ z + (2*CF*Pi^2*DiracDelta[sDivmu2])/z - (136*CA*z*DiracDelta[sDivmu2])/ 9 + 19*CF*z*DiracDelta[sDivmu2] + (2*CA*Pi^2*z*DiracDelta[sDivmu2])/3 + (CF*Pi^2*z*DiracDelta[sDivmu2])/3 - (538*CA*z^2*DiracDelta[sDivmu2])/ 27 + (8*CA*Pi^2*DiracDelta[sDivmu2]*HPL[{-1}, z])/3 + (8*CA*Pi^2*DiracDelta[sDivmu2]*HPL[{-1}, z])/(3*z) + (4*CA*Pi^2*z*DiracDelta[sDivmu2]*HPL[{-1}, z])/3 + 90*CA*DiracDelta[sDivmu2]*HPL[{0}, z] + 3*CF*DiracDelta[sDivmu2]* HPL[{0}, z] - (8*CA*Pi^2*DiracDelta[sDivmu2]*HPL[{0}, z])/3 + (16*CF*Pi^2*DiracDelta[sDivmu2]*HPL[{0}, z])/3 + (304*CA*DiracDelta[sDivmu2]*HPL[{0}, z])/(9*z) + (28*CF*DiracDelta[sDivmu2]*HPL[{0}, z])/z - (8*CA*Pi^2*DiracDelta[sDivmu2]*HPL[{0}, z])/(3*z) - (4*CF*Pi^2*DiracDelta[sDivmu2]*HPL[{0}, z])/z + 38*CA*z*DiracDelta[sDivmu2]*HPL[{0}, z] - 13*CF*z*DiracDelta[sDivmu2]* HPL[{0}, z] - (8*CA*Pi^2*z*DiracDelta[sDivmu2]*HPL[{0}, z])/3 - (8*CF*Pi^2*z*DiracDelta[sDivmu2]*HPL[{0}, z])/3 + (124*CA*z^2*DiracDelta[sDivmu2]*HPL[{0}, z])/9 - (38*CA*DiracDelta[sDivmu2]*HPL[{1}, z])/3 + 28*CF*DiracDelta[sDivmu2]* HPL[{1}, z] - 8*CF*Pi^2*DiracDelta[sDivmu2]*HPL[{1}, z] - (52*CA*DiracDelta[sDivmu2]*HPL[{1}, z])/(9*z) - (28*CF*DiracDelta[sDivmu2]*HPL[{1}, z])/z + (8*CF*Pi^2*DiracDelta[sDivmu2]*HPL[{1}, z])/z + (32*CA*z*DiracDelta[sDivmu2]*HPL[{1}, z])/3 - 12*CF*z*DiracDelta[sDivmu2]*HPL[{1}, z] + 4*CF*Pi^2*z*DiracDelta[sDivmu2]*HPL[{1}, z] + (52*CA*z^2*DiracDelta[sDivmu2]*HPL[{1}, z])/9 - 8*CA*DiracDelta[sDivmu2]*HPL[{2}, z] + (160*CA*DiracDelta[sDivmu2]*HPL[{2}, z])/(3*z) + (12*CF*DiracDelta[sDivmu2]*HPL[{2}, z])/z - 4*CA*z*DiracDelta[sDivmu2]* HPL[{2}, z] - 8*CF*z*DiracDelta[sDivmu2]*HPL[{2}, z] - (16*CA*z^2*DiracDelta[sDivmu2]*HPL[{2}, z])/3 + 32*CA*DiracDelta[sDivmu2]*HPL[{3}, z] + (48*CA*DiracDelta[sDivmu2]*HPL[{3}, z])/z - (16*CF*DiracDelta[sDivmu2]*HPL[{3}, z])/z + 40*CA*z*DiracDelta[sDivmu2]* HPL[{3}, z] + 16*CA*DiracDelta[sDivmu2]*HPL[{-2, 0}, z] + (16*CA*DiracDelta[sDivmu2]*HPL[{-2, 0}, z])/z + 8*CA*z*DiracDelta[sDivmu2]*HPL[{-2, 0}, z] + 4*CA*z*DiracDelta[sDivmu2]* HPL[{-1, 0}, z] - 16*CA*DiracDelta[sDivmu2]*HPL[{-1, 2}, z] - (16*CA*DiracDelta[sDivmu2]*HPL[{-1, 2}, z])/z - 8*CA*z*DiracDelta[sDivmu2]*HPL[{-1, 2}, z] - 32*CA*DiracDelta[sDivmu2]* HPL[{0, 0}, z] - 16*CF*DiracDelta[sDivmu2]*HPL[{0, 0}, z] - (160*CA*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(3*z) - 10*CA*z*DiracDelta[sDivmu2]*HPL[{0, 0}, z] + 13*CF*z*DiracDelta[sDivmu2]* HPL[{0, 0}, z] - 12*CA*z*DiracDelta[sDivmu2]*HPL[{1, 0}, z] + 32*CA*DiracDelta[sDivmu2]*HPL[{1, 1}, z] + 32*CF*DiracDelta[sDivmu2]* HPL[{1, 1}, z] - (124*CA*DiracDelta[sDivmu2]*HPL[{1, 1}, z])/(3*z) - (24*CF*DiracDelta[sDivmu2]*HPL[{1, 1}, z])/z - 2*CF*z*DiracDelta[sDivmu2]*HPL[{1, 1}, z] + (16*CA*z^2*DiracDelta[sDivmu2]*HPL[{1, 1}, z])/3 - 16*CA*DiracDelta[sDivmu2]*HPL[{1, 2}, z] - 8*CF*DiracDelta[sDivmu2]* HPL[{1, 2}, z] + (16*CA*DiracDelta[sDivmu2]*HPL[{1, 2}, z])/z + (8*CF*DiracDelta[sDivmu2]*HPL[{1, 2}, z])/z + 8*CA*z*DiracDelta[sDivmu2]* HPL[{1, 2}, z] + 4*CF*z*DiracDelta[sDivmu2]*HPL[{1, 2}, z] + 16*CA*DiracDelta[sDivmu2]*HPL[{2, 0}, z] + 16*CF*DiracDelta[sDivmu2]* HPL[{2, 0}, z] - (16*CA*DiracDelta[sDivmu2]*HPL[{2, 0}, z])/z - (16*CF*DiracDelta[sDivmu2]*HPL[{2, 0}, z])/z - 8*CA*z*DiracDelta[sDivmu2]*HPL[{2, 0}, z] - 8*CF*z*DiracDelta[sDivmu2]* HPL[{2, 0}, z] - 8*CA*DiracDelta[sDivmu2]*HPL[{2, 1}, z] - 40*CF*DiracDelta[sDivmu2]*HPL[{2, 1}, z] - (24*CA*DiracDelta[sDivmu2]*HPL[{2, 1}, z])/z + (48*CF*DiracDelta[sDivmu2]*HPL[{2, 1}, z])/z - 20*CA*z*DiracDelta[sDivmu2]*HPL[{2, 1}, z] + 20*CF*z*DiracDelta[sDivmu2]* HPL[{2, 1}, z] + 24*CA*DiracDelta[sDivmu2]*HPL[{-1, 0, 0}, z] + (24*CA*DiracDelta[sDivmu2]*HPL[{-1, 0, 0}, z])/z + 12*CA*z*DiracDelta[sDivmu2]*HPL[{-1, 0, 0}, z] - 40*CA*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z] + 20*CF*DiracDelta[sDivmu2]* HPL[{0, 0, 0}, z] - (80*CA*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z])/z - 60*CA*z*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z] - 10*CF*z*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z] + 24*CA*DiracDelta[sDivmu2]*HPL[{1, 0, 0}, z] - 32*CF*DiracDelta[sDivmu2]* HPL[{1, 0, 0}, z] - (24*CA*DiracDelta[sDivmu2]*HPL[{1, 0, 0}, z])/z + (32*CF*DiracDelta[sDivmu2]*HPL[{1, 0, 0}, z])/z - 12*CA*z*DiracDelta[sDivmu2]*HPL[{1, 0, 0}, z] + 16*CF*z*DiracDelta[sDivmu2]*HPL[{1, 0, 0}, z] - 32*CA*DiracDelta[sDivmu2]*HPL[{1, 1, 0}, z] + 8*CF*DiracDelta[sDivmu2]* HPL[{1, 1, 0}, z] + (32*CA*DiracDelta[sDivmu2]*HPL[{1, 1, 0}, z])/z - (8*CF*DiracDelta[sDivmu2]*HPL[{1, 1, 0}, z])/z + 16*CA*z*DiracDelta[sDivmu2]*HPL[{1, 1, 0}, z] - 4*CF*z*DiracDelta[sDivmu2]*HPL[{1, 1, 0}, z] + 8*CA*DiracDelta[sDivmu2]*HPL[{1, 1, 1}, z] + 40*CF*DiracDelta[sDivmu2]* HPL[{1, 1, 1}, z] - (8*CA*DiracDelta[sDivmu2]*HPL[{1, 1, 1}, z])/z - (40*CF*DiracDelta[sDivmu2]*HPL[{1, 1, 1}, z])/z - 4*CA*z*DiracDelta[sDivmu2]*HPL[{1, 1, 1}, z] - 20*CF*z*DiracDelta[sDivmu2]*HPL[{1, 1, 1}, z] + (44*CA*PlusDistribution[sDivmu2^(-1)])/3 - 30*CF*PlusDistribution[sDivmu2^(-1)] - (8*CA*Pi^2*PlusDistribution[sDivmu2^(-1)])/3 + 8*CF*Pi^2*PlusDistribution[sDivmu2^(-1)] + (52*CA*PlusDistribution[sDivmu2^(-1)])/(9*z) + (28*CF*PlusDistribution[sDivmu2^(-1)])/z - (20*CF*Pi^2*PlusDistribution[sDivmu2^(-1)])/(3*z) - (32*CA*z*PlusDistribution[sDivmu2^(-1)])/3 + 12*CF*z*PlusDistribution[sDivmu2^(-1)] - (4*CA*Pi^2*z*PlusDistribution[sDivmu2^(-1)])/3 - 4*CF*Pi^2*z*PlusDistribution[sDivmu2^(-1)] - (52*CA*z^2*PlusDistribution[sDivmu2^(-1)])/9 + 8*CA*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)] - (160*CA*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/(3*z) - (12*CF*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/z + 8*CA*z*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)] + (16*CA*z^2*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/3 - 32*CA*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)] - 32*CF*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)] + (124*CA*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)])/(3*z) + (24*CF*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)])/z - 4*CA*z*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)] + 6*CF*z*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)] - (16*CA*z^2*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)])/3 + 16*CA*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)] + 24*CF*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)] + (16*CA*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)])/z - (32*CF*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)])/z + 16*CA*z*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)] - 12*CF*z*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)] + 16*CA*HPL[{-1, 0}, z]*PlusDistribution[sDivmu2^(-1)] + (16*CA*HPL[{-1, 0}, z]*PlusDistribution[sDivmu2^(-1)])/z + 8*CA*z*HPL[{-1, 0}, z]*PlusDistribution[sDivmu2^(-1)] - 32*CA*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)] + 16*CF*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)] - (48*CA*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)])/z - 40*CA*z*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)] - 8*CF*z*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)] + 32*CA*HPL[{1, 0}, z]*PlusDistribution[sDivmu2^(-1)] - 16*CF*HPL[{1, 0}, z]*PlusDistribution[sDivmu2^(-1)] - (32*CA*HPL[{1, 0}, z]*PlusDistribution[sDivmu2^(-1)])/z + (16*CF*HPL[{1, 0}, z]*PlusDistribution[sDivmu2^(-1)])/z - 16*CA*z*HPL[{1, 0}, z]*PlusDistribution[sDivmu2^(-1)] + 8*CF*z*HPL[{1, 0}, z]*PlusDistribution[sDivmu2^(-1)] - 16*CA*HPL[{1, 1}, z]*PlusDistribution[sDivmu2^(-1)] - 32*CF*HPL[{1, 1}, z]*PlusDistribution[sDivmu2^(-1)] + (16*CA*HPL[{1, 1}, z]*PlusDistribution[sDivmu2^(-1)])/z + (32*CF*HPL[{1, 1}, z]*PlusDistribution[sDivmu2^(-1)])/z + 8*CA*z*HPL[{1, 1}, z]*PlusDistribution[sDivmu2^(-1)] + 16*CF*z*HPL[{1, 1}, z]*PlusDistribution[sDivmu2^(-1)] + 32*CA*PlusDistribution[Log[sDivmu2]/sDivmu2] + 32*CF*PlusDistribution[Log[sDivmu2]/sDivmu2] - (124*CA*PlusDistribution[Log[sDivmu2]/sDivmu2])/(3*z) - (24*CF*PlusDistribution[Log[sDivmu2]/sDivmu2])/z + 4*CA*z*PlusDistribution[Log[sDivmu2]/sDivmu2] - 6*CF*z*PlusDistribution[Log[sDivmu2]/sDivmu2] + (16*CA*z^2*PlusDistribution[Log[sDivmu2]/sDivmu2])/3 - 16*CA*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] - 8*CF*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] - (16*CA*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2])/z + (16*CF*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2])/z - 16*CA*z*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] + 4*CF*z*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] + 16*CA*HPL[{1}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] + 32*CF*HPL[{1}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] - (16*CA*HPL[{1}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2])/z - (32*CF*HPL[{1}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2])/z - 8*CA*z*HPL[{1}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] - 16*CF*z*HPL[{1}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] - 24*CF*PlusDistribution[Log[sDivmu2]^2/sDivmu2] + (24*CF*PlusDistribution[Log[sDivmu2]^2/sDivmu2])/z + 12*CF*z*PlusDistribution[Log[sDivmu2]^2/sDivmu2] + 8*CA*DiracDelta[sDivmu2]*Zeta[3] - 8*CF*DiracDelta[sDivmu2]*Zeta[3] - (64*CA*DiracDelta[sDivmu2]*Zeta[3])/z + (16*CF*DiracDelta[sDivmu2]*Zeta[3])/z - 40*CA*z*DiracDelta[sDivmu2]* Zeta[3] + 4*CF*z*DiracDelta[sDivmu2]*Zeta[3] matchingCoefficientFragmentingQuarkJetMeasureqbarprimeOrder2 = (-88*TF*DiracDelta[sDivmu2])/9 - (4*Pi^2*TF*DiracDelta[sDivmu2])/3 - (214*TF*DiracDelta[sDivmu2])/(27*z) - (32*TF*z*DiracDelta[sDivmu2])/9 - 2*Pi^2*TF*z*DiracDelta[sDivmu2] + (574*TF*z^2*DiracDelta[sDivmu2])/27 - 40*TF*DiracDelta[sDivmu2]*HPL[{0}, z] + (4*Pi^2*TF*DiracDelta[sDivmu2]*HPL[{0}, z])/3 - (28*TF*DiracDelta[sDivmu2]*HPL[{0}, z])/(9*z) - 36*TF*z*DiracDelta[sDivmu2]*HPL[{0}, z] + (4*Pi^2*TF*z*DiracDelta[sDivmu2]*HPL[{0}, z])/3 - (124*TF*z^2*DiracDelta[sDivmu2]*HPL[{0}, z])/9 + (80*TF*DiracDelta[sDivmu2]*HPL[{1}, z])/3 + (28*TF*DiracDelta[sDivmu2]*HPL[{1}, z])/(9*z) - (56*TF*z*DiracDelta[sDivmu2]*HPL[{1}, z])/3 - (100*TF*z^2*DiracDelta[sDivmu2]*HPL[{1}, z])/9 + 4*TF*DiracDelta[sDivmu2]*HPL[{2}, z] - (16*TF*DiracDelta[sDivmu2]*HPL[{2}, z])/(3*z) + 16*TF*z*DiracDelta[sDivmu2]*HPL[{2}, z] + (16*TF*z^2*DiracDelta[sDivmu2]*HPL[{2}, z])/3 - 16*TF*DiracDelta[sDivmu2]*HPL[{3}, z] - 16*TF*z*DiracDelta[sDivmu2]* HPL[{3}, z] - 2*TF*DiracDelta[sDivmu2]*HPL[{0, 0}, z] + (16*TF*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(3*z) - 10*TF*z*DiracDelta[sDivmu2]*HPL[{0, 0}, z] + 4*TF*DiracDelta[sDivmu2]* HPL[{1, 1}, z] + (16*TF*DiracDelta[sDivmu2]*HPL[{1, 1}, z])/(3*z) - 4*TF*z*DiracDelta[sDivmu2]*HPL[{1, 1}, z] - (16*TF*z^2*DiracDelta[sDivmu2]*HPL[{1, 1}, z])/3 + 8*TF*DiracDelta[sDivmu2]*HPL[{2, 1}, z] + 8*TF*z*DiracDelta[sDivmu2]* HPL[{2, 1}, z] + 20*TF*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z] + 20*TF*z*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z] - (80*TF*PlusDistribution[sDivmu2^(-1)])/3 + (4*Pi^2*TF*PlusDistribution[sDivmu2^(-1)])/3 - (28*TF*PlusDistribution[sDivmu2^(-1)])/(9*z) + (56*TF*z*PlusDistribution[sDivmu2^(-1)])/3 + (4*Pi^2*TF*z*PlusDistribution[sDivmu2^(-1)])/3 + (100*TF*z^2*PlusDistribution[sDivmu2^(-1)])/9 - 4*TF*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)] + (16*TF*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/(3*z) - 16*TF*z*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)] - (16*TF*z^2*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/3 - 4*TF*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)] - (16*TF*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)])/(3*z) + 4*TF*z*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)] + (16*TF*z^2*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)])/3 - 8*TF*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)] - 8*TF*z*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)] + 16*TF*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)] + 16*TF*z*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)] + 4*TF*PlusDistribution[Log[sDivmu2]/sDivmu2] + (16*TF*PlusDistribution[Log[sDivmu2]/sDivmu2])/(3*z) - 4*TF*z*PlusDistribution[Log[sDivmu2]/sDivmu2] - (16*TF*z^2*PlusDistribution[Log[sDivmu2]/sDivmu2])/3 + 8*TF*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] + 8*TF*z*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] + 8*TF*DiracDelta[sDivmu2]*Zeta[3] + 8*TF*z*DiracDelta[sDivmu2]*Zeta[3] matchingCoefficientFragmentingQuarkJetMeasureqbarOrder2 = CA*DiracDelta[sDivmu2] - 2*CF*DiracDelta[sDivmu2] - (CA*Pi^2*DiracDelta[sDivmu2])/3 + (2*CF*Pi^2*DiracDelta[sDivmu2])/3 - (88*TF*DiracDelta[sDivmu2])/9 - (4*Pi^2*TF*DiracDelta[sDivmu2])/3 - (214*TF*DiracDelta[sDivmu2])/(27*z) - CA*z*DiracDelta[sDivmu2] + 2*CF*z*DiracDelta[sDivmu2] - (CA*Pi^2*z*DiracDelta[sDivmu2])/3 + (2*CF*Pi^2*z*DiracDelta[sDivmu2])/3 - (32*TF*z*DiracDelta[sDivmu2])/9 - 2*Pi^2*TF*z*DiracDelta[sDivmu2] + (574*TF*z^2*DiracDelta[sDivmu2])/27 + (4*CA*Pi^2*DiracDelta[sDivmu2]*HPL[{-1}, z])/(3*(1 + z)) - (8*CF*Pi^2*DiracDelta[sDivmu2]*HPL[{-1}, z])/(3*(1 + z)) + (4*CA*Pi^2*z^2*DiracDelta[sDivmu2]*HPL[{-1}, z])/(3*(1 + z)) - (8*CF*Pi^2*z^2*DiracDelta[sDivmu2]*HPL[{-1}, z])/(3*(1 + z)) - 7*CA*DiracDelta[sDivmu2]*HPL[{0}, z] + 14*CF*DiracDelta[sDivmu2]* HPL[{0}, z] - 40*TF*DiracDelta[sDivmu2]*HPL[{0}, z] + (4*Pi^2*TF*DiracDelta[sDivmu2]*HPL[{0}, z])/3 - (28*TF*DiracDelta[sDivmu2]*HPL[{0}, z])/(9*z) + CA*z*DiracDelta[sDivmu2]*HPL[{0}, z] - 2*CF*z*DiracDelta[sDivmu2]* HPL[{0}, z] - 36*TF*z*DiracDelta[sDivmu2]*HPL[{0}, z] + (4*Pi^2*TF*z*DiracDelta[sDivmu2]*HPL[{0}, z])/3 - (124*TF*z^2*DiracDelta[sDivmu2]*HPL[{0}, z])/9 + 8*CA*DiracDelta[sDivmu2]*HPL[{1}, z] - 16*CF*DiracDelta[sDivmu2]* HPL[{1}, z] + (80*TF*DiracDelta[sDivmu2]*HPL[{1}, z])/3 + (28*TF*DiracDelta[sDivmu2]*HPL[{1}, z])/(9*z) - 8*CA*z*DiracDelta[sDivmu2]*HPL[{1}, z] + 16*CF*z*DiracDelta[sDivmu2]* HPL[{1}, z] - (56*TF*z*DiracDelta[sDivmu2]*HPL[{1}, z])/3 - (100*TF*z^2*DiracDelta[sDivmu2]*HPL[{1}, z])/9 + 4*CA*DiracDelta[sDivmu2]*HPL[{2}, z] - 8*CF*DiracDelta[sDivmu2]* HPL[{2}, z] + 4*TF*DiracDelta[sDivmu2]*HPL[{2}, z] - (16*TF*DiracDelta[sDivmu2]*HPL[{2}, z])/(3*z) + 4*CA*z*DiracDelta[sDivmu2]*HPL[{2}, z] - 8*CF*z*DiracDelta[sDivmu2]* HPL[{2}, z] + 16*TF*z*DiracDelta[sDivmu2]*HPL[{2}, z] + (16*TF*z^2*DiracDelta[sDivmu2]*HPL[{2}, z])/3 - 16*TF*DiracDelta[sDivmu2]*HPL[{3}, z] - 16*TF*z*DiracDelta[sDivmu2]* HPL[{3}, z] + (4*CA*DiracDelta[sDivmu2]*HPL[{3}, z])/(1 + z) - (8*CF*DiracDelta[sDivmu2]*HPL[{3}, z])/(1 + z) + (4*CA*z^2*DiracDelta[sDivmu2]*HPL[{3}, z])/(1 + z) - (8*CF*z^2*DiracDelta[sDivmu2]*HPL[{3}, z])/(1 + z) + (8*CA*DiracDelta[sDivmu2]*HPL[{-2, 0}, z])/(1 + z) - (16*CF*DiracDelta[sDivmu2]*HPL[{-2, 0}, z])/(1 + z) + (8*CA*z^2*DiracDelta[sDivmu2]*HPL[{-2, 0}, z])/(1 + z) - (16*CF*z^2*DiracDelta[sDivmu2]*HPL[{-2, 0}, z])/(1 + z) + 4*CA*DiracDelta[sDivmu2]*HPL[{-1, 0}, z] - 8*CF*DiracDelta[sDivmu2]* HPL[{-1, 0}, z] + 4*CA*z*DiracDelta[sDivmu2]*HPL[{-1, 0}, z] - 8*CF*z*DiracDelta[sDivmu2]*HPL[{-1, 0}, z] - (8*CA*DiracDelta[sDivmu2]*HPL[{-1, 2}, z])/(1 + z) + (16*CF*DiracDelta[sDivmu2]*HPL[{-1, 2}, z])/(1 + z) - (8*CA*z^2*DiracDelta[sDivmu2]*HPL[{-1, 2}, z])/(1 + z) + (16*CF*z^2*DiracDelta[sDivmu2]*HPL[{-1, 2}, z])/(1 + z) - 8*CA*DiracDelta[sDivmu2]*HPL[{0, 0}, z] + 16*CF*DiracDelta[sDivmu2]* HPL[{0, 0}, z] - 2*TF*DiracDelta[sDivmu2]*HPL[{0, 0}, z] + (16*TF*DiracDelta[sDivmu2]*HPL[{0, 0}, z])/(3*z) - 8*CA*z*DiracDelta[sDivmu2]*HPL[{0, 0}, z] + 16*CF*z*DiracDelta[sDivmu2]* HPL[{0, 0}, z] - 10*TF*z*DiracDelta[sDivmu2]*HPL[{0, 0}, z] + 4*TF*DiracDelta[sDivmu2]*HPL[{1, 1}, z] + (16*TF*DiracDelta[sDivmu2]*HPL[{1, 1}, z])/(3*z) - 4*TF*z*DiracDelta[sDivmu2]*HPL[{1, 1}, z] - (16*TF*z^2*DiracDelta[sDivmu2]*HPL[{1, 1}, z])/3 + 8*TF*DiracDelta[sDivmu2]*HPL[{2, 1}, z] + 8*TF*z*DiracDelta[sDivmu2]* HPL[{2, 1}, z] + (12*CA*DiracDelta[sDivmu2]*HPL[{-1, 0, 0}, z])/ (1 + z) - (24*CF*DiracDelta[sDivmu2]*HPL[{-1, 0, 0}, z])/(1 + z) + (12*CA*z^2*DiracDelta[sDivmu2]*HPL[{-1, 0, 0}, z])/(1 + z) - (24*CF*z^2*DiracDelta[sDivmu2]*HPL[{-1, 0, 0}, z])/(1 + z) + 20*TF*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z] + 20*TF*z*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z] - (10*CA*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z])/(1 + z) + (20*CF*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z])/(1 + z) - (10*CA*z^2*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z])/(1 + z) + (20*CF*z^2*DiracDelta[sDivmu2]*HPL[{0, 0, 0}, z])/(1 + z) - 8*CA*PlusDistribution[sDivmu2^(-1)] + 16*CF*PlusDistribution[sDivmu2^(-1)] - (80*TF*PlusDistribution[sDivmu2^(-1)])/3 + (4*Pi^2*TF*PlusDistribution[sDivmu2^(-1)])/3 - (28*TF*PlusDistribution[sDivmu2^(-1)])/(9*z) + 8*CA*z*PlusDistribution[sDivmu2^(-1)] - 16*CF*z*PlusDistribution[sDivmu2^(-1)] + (56*TF*z*PlusDistribution[sDivmu2^(-1)])/3 + (4*Pi^2*TF*z*PlusDistribution[sDivmu2^(-1)])/3 + (100*TF*z^2*PlusDistribution[sDivmu2^(-1)])/9 + (2*CA*Pi^2*PlusDistribution[sDivmu2^(-1)])/(3*(1 + z)) - (4*CF*Pi^2*PlusDistribution[sDivmu2^(-1)])/(3*(1 + z)) + (2*CA*Pi^2*z^2*PlusDistribution[sDivmu2^(-1)])/(3*(1 + z)) - (4*CF*Pi^2*z^2*PlusDistribution[sDivmu2^(-1)])/(3*(1 + z)) - 4*CA*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)] + 8*CF*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)] - 4*TF*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)] + (16*TF*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/(3*z) - 4*CA*z*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)] + 8*CF*z*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)] - 16*TF*z*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)] - (16*TF*z^2*HPL[{0}, z]*PlusDistribution[sDivmu2^(-1)])/3 - 4*TF*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)] - (16*TF*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)])/(3*z) + 4*TF*z*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)] + (16*TF*z^2*HPL[{1}, z]*PlusDistribution[sDivmu2^(-1)])/3 - 8*TF*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)] - 8*TF*z*HPL[{2}, z]*PlusDistribution[sDivmu2^(-1)] + (8*CA*HPL[{-1, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(1 + z) - (16*CF*HPL[{-1, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(1 + z) + (8*CA*z^2*HPL[{-1, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(1 + z) - (16*CF*z^2*HPL[{-1, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(1 + z) + 16*TF*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)] + 16*TF*z*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)] - (4*CA*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(1 + z) + (8*CF*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(1 + z) - (4*CA*z^2*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(1 + z) + (8*CF*z^2*HPL[{0, 0}, z]*PlusDistribution[sDivmu2^(-1)])/(1 + z) + 4*TF*PlusDistribution[Log[sDivmu2]/sDivmu2] + (16*TF*PlusDistribution[Log[sDivmu2]/sDivmu2])/(3*z) - 4*TF*z*PlusDistribution[Log[sDivmu2]/sDivmu2] - (16*TF*z^2*PlusDistribution[Log[sDivmu2]/sDivmu2])/3 + 8*TF*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] + 8*TF*z*HPL[{0}, z]*PlusDistribution[Log[sDivmu2]/sDivmu2] + 8*TF*DiracDelta[sDivmu2]*Zeta[3] + 8*TF*z*DiracDelta[sDivmu2]*Zeta[3] - (6*CA*DiracDelta[sDivmu2]*Zeta[3])/(1 + z) + (12*CF*DiracDelta[sDivmu2]*Zeta[3])/(1 + z) - (6*CA*z^2*DiracDelta[sDivmu2]*Zeta[3])/(1 + z) + (12*CF*z^2*DiracDelta[sDivmu2]*Zeta[3])/(1 + z)