rsigc["LO", "gg"] = PlusDistribution[-1, 1 - x] + 
     ep*PlusDistribution[-1, 1 - x] + ep^2*PlusDistribution[-1, 1 - x] + 
     ep^3*PlusDistribution[-1, 1 - x]
 
rsigc["NLO", "gg"] = (11*(-1 + x)^3)/2 + (6*(1 - x + x^2)^2*HPL[{0}, x])/
      (-1 + x) + 12*x*(2 - x + x^2)*HPL[{1}, x] + 
     Pi^2*PlusDistribution[-1, 1 - x] + HPL[{0}, mu^2/Mh^2]*
      (6*x*(2 - x + x^2) - 6*PlusDistribution[0, 1 - x]) + 
     12*PlusDistribution[1, 1 - x] + ep*((-3*Pi^2*x*(2 - x + x^2))/2 + 
       (2*(-1 + x)*(25 - 41*x + 25*x^2))/3 + 
       ((23 - 68*x + 102*x^2 - 68*x^3 + 23*x^4)*HPL[{0}, x])/(2*(-1 + x)) + 
       (-11 + 57*x - 45*x^2 + 23*x^3)*HPL[{1}, x] + 
       (12*(1 - x + x^2)^2*HPL[{2}, x])/(-1 + x) + 
       (6*(1 - x + x^2)^2*HPL[{0, 0}, x])/(-1 + x) + 
       (12*(1 - x + x^2)^2*HPL[{1, 0}, x])/(-1 + x) + 
       24*x*(2 - x + x^2)*HPL[{1, 1}, x] + (-3 + Pi^2)*PlusDistribution[-1, 
         1 - x] + HPL[{0, 0}, mu^2/Mh^2]*(6*x*(2 - x + x^2) - 
         6*PlusDistribution[0, 1 - x]) + (3*Pi^2*PlusDistribution[0, 1 - x])/
        2 + 12*PlusDistribution[1, 1 - x] + HPL[{0}, mu^2/Mh^2]*
        ((-11 + 57*x - 45*x^2 + 23*x^3)/2 + (6*(1 - x + x^2)^2*HPL[{0}, x])/
          (-1 + x) + 12*x*(2 - x + x^2)*HPL[{1}, x] + 
         Pi^2*PlusDistribution[-1, 1 - x] - 6*PlusDistribution[0, 1 - x] + 
         12*PlusDistribution[1, 1 - x]) - 12*PlusDistribution[2, 1 - x]) + 
     ep^2*((8*(-1 + x)*(44 - 61*x + 44*x^2))/9 + 
       (Pi^2*(11 - 57*x + 45*x^2 - 23*x^3))/8 + 
       ((23 - 68*x + 102*x^2 - 68*x^3 + 23*x^4)*HPL[{2}, x])/(-1 + x) + 
       12*x*(2 - x + x^2)*HPL[{3}, x] + HPL[{1}, x]*
        (-3*Pi^2*x*(2 - x + x^2) + (4*(-25 + 84*x - 75*x^2 + 34*x^3))/3 + 
         12*x*(2 - x + x^2)*HPL[{0, 0}, mu^2/Mh^2]) + 
       ((23 - 68*x + 102*x^2 - 68*x^3 + 23*x^4)*HPL[{1, 0}, x])/(-1 + x) + 
       2*(-11 + 57*x - 45*x^2 + 23*x^3)*HPL[{1, 1}, x] + 
       24*x*(2 - x + x^2)*HPL[{1, 2}, x] + 12*x*(2 - x + x^2)*
        HPL[{2, 0}, x] + 24*x*(2 - x + x^2)*HPL[{2, 1}, x] + 
       12*x*(2 - x + x^2)*HPL[{1, 0, 0}, x] + 24*x*(2 - x + x^2)*
        HPL[{1, 1, 0}, x] + 48*x*(2 - x + x^2)*HPL[{1, 1, 1}, x] + 
       (-12 + Pi^2 - (11*Pi^4)/60)*PlusDistribution[-1, 1 - x] + 
       HPL[{0, 0, 0}, mu^2/Mh^2]*(6*x*(2 - x + x^2) - 
         6*PlusDistribution[0, 1 - x]) + HPL[{0, 0, 0}, x]*
        (6*x*(2 - x + x^2) - 6*PlusDistribution[0, 1 - x]) + 
       (12 - 3*Pi^2)*PlusDistribution[1, 1 - x] + HPL[{0, 0}, x]*
        ((23 - 68*x + 102*x^2 - 68*x^3 + 23*x^4)/(2*(-1 + x)) + 
         12*PlusDistribution[1, 1 - x]) + HPL[{0, 0}, mu^2/Mh^2]*
        ((-11 + 57*x - 45*x^2 + 23*x^3)/2 + Pi^2*PlusDistribution[-1, 
           1 - x] - 6*PlusDistribution[0, 1 - x] + 
         12*PlusDistribution[1, 1 - x]) + HPL[{0}, x]*
        ((-3*Pi^2*(1 - x + x^2)^2)/(2*(-1 + x)) + 
         (2*(34 - 109*x + 159*x^2 - 109*x^3 + 34*x^4))/(3*(-1 + x)) + 
         (6*(1 - x + x^2)^2*HPL[{0, 0}, mu^2/Mh^2])/(-1 + x) - 
         12*PlusDistribution[2, 1 - x]) + HPL[{0}, mu^2/Mh^2]*
        ((-3*Pi^2*x*(2 - x + x^2))/2 + (2*(-25 + 84*x - 75*x^2 + 34*x^3))/3 + 
         ((23 - 68*x + 102*x^2 - 68*x^3 + 23*x^4)*HPL[{0}, x])/(2*(-1 + x)) + 
         (-11 + 57*x - 45*x^2 + 23*x^3)*HPL[{1}, x] + 
         (12*(1 - x + x^2)^2*HPL[{2}, x])/(-1 + x) + 
         (6*(1 - x + x^2)^2*HPL[{0, 0}, x])/(-1 + x) + 
         (12*(1 - x + x^2)^2*HPL[{1, 0}, x])/(-1 + x) + 24*x*(2 - x + x^2)*
          HPL[{1, 1}, x] + (-3 + Pi^2)*PlusDistribution[-1, 1 - x] + 
         (-6 + (3*Pi^2)/2)*PlusDistribution[0, 1 - x] + 
         12*PlusDistribution[1, 1 - x] - 12*PlusDistribution[2, 1 - x]) - 
       12*PlusDistribution[2, 1 - x] + 8*PlusDistribution[3, 1 - x] - 
       14*x*(2 - x + x^2)*Zeta[3] + PlusDistribution[0, 1 - x]*
        ((3*Pi^2)/2 + 14*Zeta[3]))
 
rsigc["NLO", "qb"] = (-32*(-1 + x)^3)/27 + 
     ep*((-64*(-1 + x)^3)/81 - (32*(-1 + x)^3*HPL[{0}, mu^2/Mh^2])/27 - 
       (32*(-1 + x)^3*HPL[{0}, x])/27 - (64*(-1 + x)^3*HPL[{1}, x])/27) + 
     ep^2*((-416*(-1 + x)^3)/243 + (8*Pi^2*(-1 + x)^3)/27 - 
       (64*(-1 + x)^3*HPL[{0}, x])/81 - (128*(-1 + x)^3*HPL[{1}, x])/81 + 
       HPL[{0}, mu^2/Mh^2]*((-64*(-1 + x)^3)/81 - (32*(-1 + x)^3*HPL[{0}, x])/
          27 - (64*(-1 + x)^3*HPL[{1}, x])/27) - (64*(-1 + x)^3*HPL[{2}, x])/
        27 - (32*(-1 + x)^3*HPL[{0, 0}, mu^2/Mh^2])/27 - 
       (32*(-1 + x)^3*HPL[{0, 0}, x])/27 - (64*(-1 + x)^3*HPL[{1, 0}, x])/
        27 - (128*(-1 + x)^3*HPL[{1, 1}, x])/27)
 
rsigc["NLO", "qg"] = (-3 + 6*x - x^2)/3 - 
     (2*(2 - 2*x + x^2)*HPL[{0}, mu^2/Mh^2])/3 - 
     (2*(2 - 2*x + x^2)*HPL[{0}, x])/3 - (4*(2 - 2*x + x^2)*HPL[{1}, x])/3 + 
     ep*((Pi^2*(2 - 2*x + x^2))/6 - (2*(5 - 8*x + 2*x^2))/3 - 
       ((-1 + x)*(-7 + 3*x)*HPL[{0}, x])/3 - 
       (2*(-1 + x)*(-7 + 3*x)*HPL[{1}, x])/3 + HPL[{0}, mu^2/Mh^2]*
        (-((-1 + x)*(-7 + 3*x))/3 - (2*(2 - 2*x + x^2)*HPL[{0}, x])/3 - 
         (4*(2 - 2*x + x^2)*HPL[{1}, x])/3) - (4*(2 - 2*x + x^2)*HPL[{2}, x])/
        3 - (2*(2 - 2*x + x^2)*HPL[{0, 0}, mu^2/Mh^2])/3 - 
       (2*(2 - 2*x + x^2)*HPL[{0, 0}, x])/3 - 
       (4*(2 - 2*x + x^2)*HPL[{1, 0}, x])/3 - 
       (8*(2 - 2*x + x^2)*HPL[{1, 1}, x])/3) + 
     ep^2*((Pi^2*(-1 + x)*(-7 + 3*x))/12 - (2*(12 - 18*x + 5*x^2))/3 - 
       (2*(-1 + x)*(-7 + 3*x)*HPL[{2}, x])/3 - 
       (4*(2 - 2*x + x^2)*HPL[{3}, x])/3 - 
       ((-1 + x)*(-7 + 3*x)*HPL[{0, 0}, mu^2/Mh^2])/3 + 
       HPL[{1}, x]*((-4*(-1 + x)*(-7 + 3*x))/3 + (Pi^2*(2 - 2*x + x^2))/3 - 
         (4*(2 - 2*x + x^2)*HPL[{0, 0}, mu^2/Mh^2])/3) + 
       HPL[{0}, x]*((-2*(-1 + x)*(-7 + 3*x))/3 + (Pi^2*(2 - 2*x + x^2))/6 - 
         (2*(2 - 2*x + x^2)*HPL[{0, 0}, mu^2/Mh^2])/3) - 
       ((-1 + x)*(-7 + 3*x)*HPL[{0, 0}, x])/3 - 
       (2*(-1 + x)*(-7 + 3*x)*HPL[{1, 0}, x])/3 - 
       (4*(-1 + x)*(-7 + 3*x)*HPL[{1, 1}, x])/3 + HPL[{0}, mu^2/Mh^2]*
        ((-2*(-1 + x)*(-7 + 3*x))/3 + (Pi^2*(2 - 2*x + x^2))/6 - 
         ((-1 + x)*(-7 + 3*x)*HPL[{0}, x])/3 - (2*(-1 + x)*(-7 + 3*x)*
           HPL[{1}, x])/3 - (4*(2 - 2*x + x^2)*HPL[{2}, x])/3 - 
         (2*(2 - 2*x + x^2)*HPL[{0, 0}, x])/3 - 
         (4*(2 - 2*x + x^2)*HPL[{1, 0}, x])/3 - 
         (8*(2 - 2*x + x^2)*HPL[{1, 1}, x])/3) - 
       (8*(2 - 2*x + x^2)*HPL[{1, 2}, x])/3 - 
       (4*(2 - 2*x + x^2)*HPL[{2, 0}, x])/3 - 
       (8*(2 - 2*x + x^2)*HPL[{2, 1}, x])/3 - 
       (2*(2 - 2*x + x^2)*HPL[{0, 0, 0}, mu^2/Mh^2])/3 - 
       (2*(2 - 2*x + x^2)*HPL[{0, 0, 0}, x])/3 - 
       (4*(2 - 2*x + x^2)*HPL[{1, 0, 0}, x])/3 - 
       (8*(2 - 2*x + x^2)*HPL[{1, 1, 0}, x])/3 - 
       (16*(2 - 2*x + x^2)*HPL[{1, 1, 1}, x])/3 + 
       (14*(2 - 2*x + x^2)*Zeta[3])/9)
 
rsigc["NNLO", "gg"] = (-451035 + 12923*nl + 500121*x - 2982*nl*x - 
       533277*x^2 + 3801*nl*x^2 + 571455*x^3 - 17774*nl*x^3)/1296 + 
     (Pi^2*(-432 - 4*nl + 5400*x + 142*nl*x - 15228*x^2 + 99*nl*x^2 + 
        23922*x^3 - 387*nl*x^3 - 13365*x^4 + 132*nl*x^4))/(216*(-1 + x)) + 
     (3*Pi^2*(11 + 22*x + 33*x^2 + 16*x^3 + 8*x^4)*HPL[{-1}, x])/
      (4*(1 + x)) - ((-2376 + 34*nl + 5697*x + 14*nl*x - 8424*x^2 + 
        42*nl*x^2 + 9855*x^3 - 122*nl*x^3 - 5049*x^4 + 50*nl*x^4)*
       HPL[{2}, x])/(9*(-1 + x)) + 
     (2*(189 + 81*x - 8*nl*x + 513*x^2 - 8*nl*x^2 - 81*x^3 + 8*nl*x^3 - 
        513*x^4 + 8*nl*x^4 + 243*x^5)*HPL[{3}, x])/(3*(-1 + x)*(1 + x)) - 
     (9*(2 + 4*x + 6*x^2 + 4*x^3 + 3*x^4)*HPL[{-2, 0}, x])/(1 + x) + 
     (3*(-14 - 12*x + 21*x^2 + 11*x^3)*HPL[{-1, 0}, x])/2 - 
     (36*(1 + x + x^2)^2*HPL[{-1, 2}, x])/(1 + x) + 
     HPL[{1}, x]*((-54729 + 922*nl + 58293*x - 294*nl*x - 51138*x^2 + 
         249*nl*x^2 + 61965*x^3 - 1570*nl*x^3)/108 + 
       (Pi^2*(-162 - 2*nl + 2484*x + 4*nl*x - 3726*x^2 - 3*nl*x^2 + 
          2592*x^3 + nl*x^3 - 1242*x^4))/(72*(-1 + x)) + 
       72*x*(2 - x + x^2)*HPL[{0, 0}, mu^2/Mh^2]) + 
     HPL[{0}, x]*(-(-55863 + 1282*nl + 137241*x - 598*nl*x - 152685*x^2 + 
          333*nl*x^2 + 173529*x^3 - 3041*nl*x^3 - 109458*x^4 + 2384*nl*x^4)/
        (216*(-1 + x)) - (Pi^2*(18 + 54*x - 4*nl*x + 162*x^2 - 4*nl*x^2 - 
          54*x^3 + 4*nl*x^3 - 171*x^4 + 4*nl*x^4 + 90*x^5))/
        (6*(-1 + x)*(1 + x)) + (4*(27 - nl*x + 81*x^2 - 108*x^3 + nl*x^3 + 
          27*x^4)*HPL[{0, 0}, mu^2/Mh^2])/(3*(-1 + x))) - 
     ((-4158 + 68*nl + 10503*x + 52*nl*x - 18873*x^2 + 117*nl*x^2 + 
        22329*x^3 - 351*nl*x^3 - 10098*x^4 + 132*nl*x^4)*HPL[{0, 0}, x])/
      (36*(-1 + x)) - ((-5769 + 68*nl + 11142*x - 82*nl*x - 11043*x^2 + 
        63*nl*x^2 + 10620*x^3 - 87*nl*x^3 - 5445*x^4 + 68*nl*x^4)*
       HPL[{1, 0}, x])/(12*(-1 + x)) - 
     (2*(2673 - 16*nl - 3429*x + 24*nl*x + 3132*x^2 - 6*nl*x^2 - 2970*x^3 + 
        34*nl*x^3)*HPL[{1, 1}, x])/9 + (279*(1 - x + x^2)^2*HPL[{1, 2}, x])/
      (-1 + x) + ((360 + 108*x - 8*nl*x + 792*x^2 - 8*nl*x^2 - 108*x^3 + 
        8*nl*x^3 - 783*x^4 + 8*nl*x^4 + 333*x^5)*HPL[{2, 0}, x])/
      (2*(-1 + x)*(1 + x)) + ((837 - 810*x - 16*nl*x + 2511*x^2 - 2538*x^3 + 
        16*nl*x^3 + 837*x^4)*HPL[{2, 1}, x])/(3*(-1 + x)) + 
     (27*(1 + 2*x + 3*x^2)*HPL[{-1, -1, 0}, x])/(1 + x) - 
     (9*(9 + 18*x + 27*x^2 + 8*x^3 + 4*x^4)*HPL[{-1, 0, 0}, x])/(2*(1 + x)) + 
     ((108 + 189*x - 10*nl*x + 486*x^2 - 10*nl*x^2 - 189*x^3 + 10*nl*x^3 - 
        459*x^4 + 10*nl*x^4 + 243*x^5)*HPL[{0, 0, 0}, x])/
      (3*(-1 + x)*(1 + x)) + ((1296 - 2*nl - 2592*x + 4*nl*x + 3888*x^2 - 
        3*nl*x^2 - 2376*x^3 + nl*x^3 + 1296*x^4)*HPL[{1, 0, 0}, x])/
      (12*(-1 + x)) + ((3186 - 2*nl - 6372*x + 4*nl*x + 9558*x^2 - 3*nl*x^2 - 
        6264*x^3 + nl*x^3 + 3186*x^4)*HPL[{1, 1, 0}, x])/(12*(-1 + x)) + 
     432*x*(2 - x + x^2)*HPL[{1, 1, 1}, x] + ((201 - 10*nl)/3 - 15*Pi^2)*
      PlusDistribution[1, 1 - x] + HPL[{0, 0}, mu^2/Mh^2]*
      ((-2376 + 16*nl + 2538*x - 24*nl*x - 2241*x^2 + 6*nl*x^2 + 2673*x^3 - 
         34*nl*x^3)/18 - 6*Pi^2*PlusDistribution[-1, 1 - x] + 
       ((-33 + 2*nl)*PlusDistribution[0, 1 - x])/2 + 
       72*PlusDistribution[1, 1 - x]) + (-33 + 2*nl)*
      PlusDistribution[2, 1 - x] + 72*PlusDistribution[3, 1 - x] + 
     HPL[{0}, mu^2/Mh^2]*((-51111 + 922*nl + 48087*x - 294*nl*x - 40851*x^2 + 
         222*nl*x^2 + 58347*x^3 - 1570*nl*x^3)/216 - 
       (Pi^2*(27 - 108*x + 8*nl*x - 648*x^2 + 16*nl*x^2 - 378*x^3 + 
          8*nl*x^3 + 162*x^4))/(18*(1 + x)) - 
       ((-2079 + 34*nl + 4833*x + 14*nl*x - 7533*x^2 + 42*nl*x^2 + 9531*x^3 - 
          122*nl*x^3 - 5049*x^4 + 50*nl*x^4)*HPL[{0}, x])/(18*(-1 + x)) + 
       ((-2673 + 16*nl + 3429*x - 24*nl*x - 3132*x^2 + 6*nl*x^2 + 2970*x^3 - 
          34*nl*x^3)*HPL[{1}, x])/9 + (2*(189 - 162*x - 4*nl*x + 567*x^2 - 
          594*x^3 + 4*nl*x^3 + 189*x^4)*HPL[{2}, x])/(3*(-1 + x)) - 
       (18*(1 + x + x^2)^2*HPL[{-1, 0}, x])/(1 + x) + 
       (2*(54 + 81*x - 4*nl*x + 216*x^2 - 4*nl*x^2 - 81*x^3 + 4*nl*x^3 - 
          216*x^4 + 4*nl*x^4 + 81*x^5)*HPL[{0, 0}, x])/(3*(-1 + x)*(1 + x)) + 
       (126*(1 - x + x^2)^2*HPL[{1, 0}, x])/(-1 + x) + 
       216*x*(2 - x + x^2)*HPL[{1, 1}, x] + ((-201 + 10*nl)/6 + (15*Pi^2)/2)*
        PlusDistribution[0, 1 - x] + (33 - 2*nl)*PlusDistribution[1, 1 - x] - 
       108*PlusDistribution[2, 1 - x] + PlusDistribution[-1, 1 - x]*
        ((81 - 11*nl)/6 + ((33 - 2*nl)*Pi^2)/12 - (171*Zeta[3])/2)) - 
     ((-972 - 2*nl - 1404*x - 14*nl*x + 513*x^2 - 15*nl*x^2 + 351*x^3 + 
        14*nl*x^3 - 405*x^4 + 17*nl*x^4 + 891*x^5)*Zeta[3])/
      (6*(-1 + x)*(1 + x)) + PlusDistribution[0, 1 - x]*
      ((-303 + 14*nl)/9 + ((33 - 2*nl)*Pi^2)/6 + (351*Zeta[3])/2) + 
     PlusDistribution[-1, 1 - x]*((7533 - 988*nl)/144 + 
       ((201 - 10*nl)*Pi^2)/36 - Pi^4/80 + (5*(-99 + 2*nl)*Zeta[3])/12) + 
     ep*((-5268285 + 150313*nl + 4986639*x - 48234*nl*x - 5569551*x^2 + 
         37431*nl*x^2 + 6375357*x^3 - 163318*nl*x^3)/2592 + 
       (Pi^2*(-218673 + 4942*nl + 879606*x - 1192*nl*x - 1293894*x^2 + 
          2940*nl*x^2 + 1666494*x^3 - 23992*nl*x^3 - 1022733*x^4 + 
          16726*nl*x^4))/(2592*(-1 + x)) + 
       (Pi^4*(4293 + 30*nl + 783*x - 134*nl*x + 4455*x^2 - 119*nl*x^2 + 
          1647*x^3 + 134*nl*x^3 - 6237*x^4 + 89*nl*x^4 + 3375*x^5))/
        (1080*(-1 + x)*(1 + x)) + (3*Pi^2*(22 + 44*x + 66*x^2 + 44*x^3 + 
          37*x^4)*HPL[{-2}, x])/(4*(1 + x)) - 
       ((-9135 + 156*nl + 9738*x + 216*nl*x - 15633*x^2 + 195*nl*x^2 + 
          3645*x^3 - 308*nl*x^3 + 22167*x^4 - 261*nl*x^4 - 16308*x^5 + 
          182*nl*x^5)*HPL[{3}, x])/(9*(-1 + x)*(1 + x)) + 
       (2*(-153 - 243*x - 8*nl*x + 504*x^2 - 8*nl*x^2 - 90*x^3 + 8*nl*x^3 - 
          486*x^4 + 8*nl*x^4 + 270*x^5)*HPL[{4}, x])/((-1 + x)*(1 + x)) - 
       (9*(4 + 8*x + 12*x^2 + 8*x^3 + 7*x^4)*HPL[{-3, 0}, x])/(1 + x) + 
       (3*(-40 - 82*x - 27*x^2 + 48*x^3 + 15*x^4)*HPL[{-2, 0}, x])/
        (2*(1 + x)) - (36*(2 + 4*x + 6*x^2 + 4*x^3 + 3*x^4)*HPL[{-2, 2}, x])/
        (1 + x) + (3*Pi^2*(-37 - 74*x - 111*x^2 + 4*x^3 + 2*x^4)*
         HPL[{-1, -1}, x])/(4*(1 + x)) + 
       (6*(-20 - 38*x - 9*x^2 + 20*x^3 + 5*x^4)*HPL[{-1, 2}, x])/(1 + x) - 
       (18*(9 + 18*x + 27*x^2 + 8*x^3 + 4*x^4)*HPL[{-1, 3}, x])/(1 + x) + 
       HPL[{-1, 0}, x]*((-206 - 222*x + 303*x^2 + 234*x^3)/2 + 
         (3*Pi^2*(41 + 82*x + 123*x^2 + 20*x^3 + 10*x^4))/(4*(1 + x)) - 
         (36*(1 + x + x^2)^2*HPL[{0, 0}, mu^2/Mh^2])/(1 + x)) + 
       HPL[{2}, x]*(-(-66663 + 1174*nl + 161874*x - 660*nl*x - 195264*x^2 + 
            556*nl*x^2 + 217044*x^3 - 2656*nl*x^3 - 123003*x^4 + 1898*nl*x^4)/
          (36*(-1 + x)) - (Pi^2*(4752 + 2*nl + 216*x - 90*nl*x + 9504*x^2 - 
            89*nl*x^2 - 216*x^3 + 90*nl*x^3 - 9450*x^4 + 87*nl*x^4 + 
            4590*x^5))/(36*(-1 + x)*(1 + x)) + 
         (8*(45 - 36*x - nl*x + 135*x^2 - 144*x^3 + nl*x^3 + 45*x^4)*
           HPL[{0, 0}, mu^2/Mh^2])/(-1 + x)) + 
       (-(Pi^2*(1998 - 2*nl - 3996*x + 4*nl*x + 5994*x^2 - 3*nl*x^2 - 
             3816*x^3 + nl*x^3 + 1998*x^4))/(24*(-1 + x)) - 
         (-302535 + 4642*nl + 556335*x - 5242*nl*x - 515403*x^2 + 
           2820*nl*x^2 + 526878*x^3 - 6070*nl*x^3 - 286011*x^4 + 4930*nl*x^4)/
          (108*(-1 + x)))*HPL[{1, 0}, x] + 
       ((-139284 + 1724*nl + 151920*x - 852*nl*x - 140085*x^2 + 483*nl*x^2 + 
           151308*x^3 - 2540*nl*x^3)/36 - 
         (6*Pi^2*(1 - 64*x + 96*x^2 - 64*x^3 + 32*x^4))/(-1 + x))*
        HPL[{1, 1}, x] - ((-28692 + 380*nl + 55404*x - 326*nl*x - 60912*x^2 + 
          249*nl*x^2 + 53838*x^3 - 341*nl*x^3 - 25974*x^4 + 272*nl*x^4)*
         HPL[{1, 2}, x])/(9*(-1 + x)) + 
       ((-1134 - 2*nl - 4860*x + 4*nl*x + 7290*x^2 - 3*nl*x^2 - 4644*x^3 + 
          nl*x^3 + 2430*x^4)*HPL[{1, 3}, x])/(3*(-1 + x)) - 
       ((-48150 + 760*nl + 40734*x + 798*nl*x - 52785*x^2 + 575*nl*x^2 + 
          24381*x^3 - 1112*nl*x^3 + 86346*x^4 - 849*nl*x^4 - 79056*x^5 + 
          800*nl*x^5)*HPL[{2, 0}, x])/(36*(-1 + x)*(1 + x)) - 
       ((-6369 + 104*nl + 12654*x + 8*nl*x - 20196*x^2 + 68*nl*x^2 + 
          22338*x^3 - 208*nl*x^3 - 10440*x^4 + 100*nl*x^4)*HPL[{2, 1}, x])/
        (3*(-1 + x)) + (2*(-108 - 1620*x - 32*nl*x + 3402*x^2 - 32*nl*x^2 - 
          162*x^3 + 32*nl*x^3 - 3375*x^4 + 32*nl*x^4 + 1593*x^5)*
         HPL[{2, 2}, x])/(3*(-1 + x)*(1 + x)) + 
       ((-972 - 2808*x - 104*nl*x + 7560*x^2 - 104*nl*x^2 - 1188*x^3 + 
          104*nl*x^3 - 7479*x^4 + 104*nl*x^4 + 3321*x^5)*HPL[{3, 0}, x])/
        (6*(-1 + x)*(1 + x)) + (2*(-513 - 1539*x - 40*nl*x + 2997*x^2 - 
          40*nl*x^2 - 243*x^3 + 40*nl*x^3 - 2997*x^4 + 40*nl*x^4 + 1485*x^5)*
         HPL[{3, 1}, x])/(3*(-1 + x)*(1 + x)) + 
       (9*(6 + 12*x + 18*x^2 + 12*x^3 + 13*x^4)*HPL[{-2, -1, 0}, x])/
        (1 + x) - (9*(18 + 36*x + 54*x^2 + 36*x^3 + 35*x^4)*
         HPL[{-2, 0, 0}, x])/(2*(1 + x)) + 
       (63*(1 + 2*x + 3*x^2)*HPL[{-1, -2, 0}, x])/(1 + x) - 
       (3*(-74 - 158*x - 45*x^2 + 152*x^3 + 77*x^4)*HPL[{-1, -1, 0}, x])/
        (2*(1 + x)) + (108*(1 + 2*x + 3*x^2)*HPL[{-1, -1, 2}, x])/(1 + x) + 
       (3*(-222 - 450*x - 99*x^2 + 376*x^3 + 163*x^4)*HPL[{-1, 0, 0}, x])/
        (4*(1 + x)) - (18*(4 + 8*x + 12*x^2 + 10*x^3 + 5*x^4)*
         HPL[{-1, 2, 0}, x])/(1 + x) - (144*(1 + x + x^2)^2*
         HPL[{-1, 2, 1}, x])/(1 + x) - 
       ((-20217 + 280*nl + 37062*x - 258*nl*x - 37611*x^2 + 200*nl*x^2 + 
          34236*x^3 - 276*nl*x^3 - 17433*x^4 + 204*nl*x^4)*HPL[{1, 0, 0}, x])/
        (12*(-1 + x)) - ((-112338 + 1532*nl + 215136*x - 1316*nl*x - 
          236844*x^2 + 993*nl*x^2 + 211680*x^3 - 1325*nl*x^3 - 103410*x^4 + 
          1088*nl*x^4)*HPL[{1, 1, 0}, x])/(36*(-1 + x)) - 
       (4*(2772 - 16*nl - 4374*x + 24*nl*x + 3753*x^2 - 6*nl*x^2 - 3393*x^3 + 
          34*nl*x^3)*HPL[{1, 1, 1}, x])/3 + 
       ((-1512 - 2*nl - 10368*x + 4*nl*x + 15552*x^2 - 3*nl*x^2 - 10260*x^3 + 
          nl*x^3 + 5184*x^4)*HPL[{1, 1, 2}, x])/(3*(-1 + x)) + 
       ((-1782 - 2*nl - 24948*x + 4*nl*x + 37422*x^2 - 3*nl*x^2 - 24624*x^3 + 
          nl*x^3 + 12474*x^4)*HPL[{1, 2, 0}, x])/(12*(-1 + x)) + 
       (9*(-51 - 394*x + 591*x^2 - 394*x^3 + 197*x^4)*HPL[{1, 2, 1}, x])/
        (-1 + x) + ((-972 - 2*nl - 2592*x - 70*nl*x + 6912*x^2 - 71*nl*x^2 - 
          1404*x^3 + 70*nl*x^3 - 6777*x^4 + 73*nl*x^4 + 2619*x^5)*
         HPL[{2, 0, 0}, x])/(6*(-1 + x)*(1 + x)) + 
       ((-810 - 2*nl - 6318*x - 126*nl*x + 13230*x^2 - 127*nl*x^2 - 810*x^3 + 
          126*nl*x^3 - 13176*x^4 + 129*nl*x^4 + 6156*x^5)*HPL[{2, 1, 0}, x])/
        (6*(-1 + x)*(1 + x)) + ((-459 - 1818*x - 32*nl*x + 5319*x^2 - 
          5274*x^3 + 32*nl*x^3 + 1773*x^4)*HPL[{2, 1, 1}, x])/(-1 + x) + 
       (9*(-13 - 26*x - 39*x^2 + 4*x^3 + 2*x^4)*HPL[{-1, -1, -1, 0}, x])/
        (1 + x) - (9*(-35 - 70*x - 105*x^2 + 12*x^3 + 6*x^4)*
         HPL[{-1, -1, 0, 0}, x])/(2*(1 + x)) - 
       (9*(13 + 26*x + 39*x^2 + 8*x^3 + 4*x^4)*HPL[{-1, 0, 0, 0}, x])/
        (1 + x) + ((-2295 - 4*nl - 3402*x + 8*nl*x + 5103*x^2 - 6*nl*x^2 - 
          2808*x^3 + 2*nl*x^3 + 1701*x^4)*HPL[{1, 0, 0, 0}, x])/
        (6*(-1 + x)) + ((-5454 - 14*nl - 17604*x + 28*nl*x + 26406*x^2 - 
          21*nl*x^2 - 16632*x^3 + 7*nl*x^3 + 8802*x^4)*HPL[{1, 1, 0, 0}, x])/
        (12*(-1 + x)) + ((-1620 - 2*nl - 10152*x + 4*nl*x + 15228*x^2 - 
          3*nl*x^2 - 10044*x^3 + nl*x^3 + 5076*x^4)*HPL[{1, 1, 1, 0}, x])/
        (3*(-1 + x)) + (9*(-781 + 1280*x - 640*x^2 + 640*x^3)*
         HPL[{1, 1, 1, 1}, x])/2 + HPL[{0, 0, 0, 0}, x]*
        ((-891 - 702*x - 26*nl*x + 1350*x^2 - 26*nl*x^2 - 513*x^3 + 
           26*nl*x^3 - 1269*x^4 + 26*nl*x^4 + 702*x^5)/(3*(-1 + x)*(1 + x)) - 
         405*PlusDistribution[0, 1 - x]) + HPL[{0, 0, 0}, mu^2/Mh^2]*
        ((-2376 + 16*nl + 2538*x - 24*nl*x - 2241*x^2 + 6*nl*x^2 + 2673*x^3 - 
           34*nl*x^3)/6 - 18*Pi^2*PlusDistribution[-1, 1 - x] + 
         (3*(-33 + 2*nl)*PlusDistribution[0, 1 - x])/2 + 
         216*PlusDistribution[1, 1 - x]) + HPL[{0, 0, 0}, x]*
        (-(-5172 + 104*nl + 4449*x + 158*nl*x - 10008*x^2 + 158*nl*x^2 + 
            3780*x^3 - 242*nl*x^3 + 13860*x^4 - 214*nl*x^4 - 9981*x^5 + 
            132*nl*x^5)/(12*(-1 + x)*(1 + x)) + 
         ((97 - 6*nl)*PlusDistribution[0, 1 - x])/2 + 
         666*PlusDistribution[1, 1 - x]) + HPL[{0, 0}, x]*
        (-(-181953 + 3522*nl + 483975*x - 1722*nl*x - 597213*x^2 + 
            2151*nl*x^2 + 650619*x^3 - 10067*nl*x^3 - 368064*x^4 + 
            6764*nl*x^4)/(216*(-1 + x)) - 
         (Pi^2*(378 + 405*x - 32*nl*x + 1620*x^2 - 32*nl*x^2 - 405*x^3 + 
            32*nl*x^3 - 1593*x^4 + 32*nl*x^4 + 1107*x^5))/
          (18*(-1 + x)*(1 + x)) + (-97 + 6*nl)*PlusDistribution[1, 1 - x] - 
         594*PlusDistribution[2, 1 - x]) + ((-501 + 26*nl)/3 + 93*Pi^2)*
        PlusDistribution[2, 1 - x] - 2*(-69 + 2*nl)*PlusDistribution[3, 
         1 - x] - (3*(-141 + 781*x)*PlusDistribution[4, 1 - x])/16 + 
       HPL[{0, 0}, mu^2/Mh^2]*((-63027 + 862*nl + 57789*x - 426*nl*x - 
           51777*x^2 + 210*nl*x^2 + 69039*x^3 - 1270*nl*x^3)/72 - 
         (Pi^2*(9 - 72*x + 4*nl*x - 342*x^2 + 8*nl*x^2 - 198*x^3 + 4*nl*x^3 + 
            72*x^4))/(3*(1 + x)) + (4*(81 + 108*x - 5*nl*x + 297*x^2 - 
            5*nl*x^2 - 108*x^3 + 5*nl*x^3 - 297*x^4 + 5*nl*x^4 + 108*x^5)*
           HPL[{0, 0}, x])/(3*(-1 + x)*(1 + x)) + 
         (360*(1 - x + x^2)^2*HPL[{1, 0}, x])/(-1 + x) + 
         648*x*(2 - x + x^2)*HPL[{1, 1}, x] + ((-501 + 26*nl)/6 + 21*Pi^2)*
          PlusDistribution[0, 1 - x] - 3*(-57 + 2*nl)*PlusDistribution[1, 
           1 - x] - 324*PlusDistribution[2, 1 - x] + 
         PlusDistribution[-1, 1 - x]*((81 - 11*nl)/3 + ((9 - 2*nl)*Pi^2)/4 - 
           243*Zeta[3])) + PlusDistribution[1, 1 - x]*((1491 - 86*nl)/9 + 
         ((-543 + 22*nl)*Pi^2)/12 - 990*Zeta[3]) + 
       ((-1008 - 20*nl + 16020*x + 26*nl*x - 15534*x^2 + 167*nl*x^2 - 
          2232*x^3 - 224*nl*x^3 + 19359*x^4 - 213*nl*x^4 - 12375*x^5 + 
          132*nl*x^5)*Zeta[3])/(12*(-1 + x)*(1 + x)) + 
       PlusDistribution[0, 1 - x]*((2*(-1365 + 62*nl))/27 + 
         ((1803 - 94*nl)*Pi^2)/72 + (9*Pi^4)/4 + ((615 - 16*nl)*Zeta[3])/2) + 
       HPL[{-1}, x]*((-3*Pi^2*(-78 - 154*x - 39*x^2 + 104*x^3 + 39*x^4))/
          (8*(1 + x)) + (3*Pi^2*(11 + 22*x + 33*x^2 + 16*x^3 + 8*x^4)*
           HPL[{0}, mu^2/Mh^2])/(2*(1 + x)) + 
         (9*(49 + 98*x + 147*x^2 + 24*x^3 + 12*x^4)*Zeta[3])/(2*(1 + x))) + 
       HPL[{1}, x]*((-325782 + 6442*nl + 323901*x - 2028*nl*x - 329958*x^2 + 
           1647*nl*x^2 + 367245*x^3 - 7972*nl*x^3)/108 + 
         (Pi^2*(-60048 + 340*nl + 142452*x - 868*nl*x - 153738*x^2 + 
            663*nl*x^2 + 139644*x^3 - 919*nl*x^3 - 67878*x^4 + 748*nl*x^4))/
          (216*(-1 + x)) + ((-2673 + 16*nl + 3861*x - 24*nl*x - 3348*x^2 + 
            6*nl*x^2 + 3186*x^3 - 34*nl*x^3)*HPL[{0, 0}, mu^2/Mh^2])/3 + 
         216*x*(2 - x + x^2)*HPL[{0, 0, 0}, mu^2/Mh^2] - 
         ((-6318 - 10*nl - 11124*x + 20*nl*x + 16686*x^2 - 15*nl*x^2 - 
            10368*x^3 + 5*nl*x^3 + 5562*x^4)*Zeta[3])/(12*(-1 + x))) + 
       HPL[{0}, x]*(-(-619614 + 14786*nl + 1538658*x - 14778*nl*x - 
            1691739*x^2 + 5093*nl*x^2 + 1982322*x^3 - 27275*nl*x^3 - 
            1245087*x^4 + 24130*nl*x^4)/(432*(-1 + x)) + 
         (Pi^2*(-13689 + 238*nl + 29052*x + 474*nl*x - 53568*x^2 + 
            605*nl*x^2 + 8964*x^3 - 880*nl*x^3 + 65745*x^4 - 771*nl*x^4 - 
            42444*x^5 + 478*nl*x^5))/(216*(-1 + x)*(1 + x)) - 
         ((-2295 + 34*nl + 4608*x + 4*nl*x - 7857*x^2 + 34*nl*x^2 + 
            9900*x^3 - 104*nl*x^3 - 4869*x^4 + 50*nl*x^4)*
           HPL[{0, 0}, mu^2/Mh^2])/(6*(-1 + x)) + 
         (4*(27 - nl*x + 81*x^2 - 108*x^3 + nl*x^3 + 27*x^4)*
           HPL[{0, 0, 0}, mu^2/Mh^2])/(-1 + x) + (97 - 6*nl)*
          PlusDistribution[2, 1 - x] + 372*PlusDistribution[3, 1 - x] - 
         ((243 - 2*nl + 81*x - 26*nl*x + 1323*x^2 - 27*nl*x^2 - 81*x^3 + 
            26*nl*x^3 - 1350*x^4 + 29*nl*x^4 + 1269*x^5)*Zeta[3])/
          (3*(-1 + x)*(1 + x))) + HPL[{0}, mu^2/Mh^2]*
        ((-305433 + 6361*nl + 276057*x - 1974*nl*x - 283761*x^2 + 
           1467*nl*x^2 + 348921*x^3 - 7918*nl*x^3)/216 + 
         (Pi^2*(-6480 + 40*nl + 18846*x + 300*nl*x - 38259*x^2 + 500*nl*x^2 + 
            24057*x^3 - 766*nl*x^3 + 45657*x^4 - 576*nl*x^4 - 42633*x^5 + 
            430*nl*x^5))/(216*(-1 + x)*(1 + x)) + 
         (-(-60651 + 1174*nl + 144846*x - 660*nl*x - 174960*x^2 + 
              538*nl*x^2 + 206550*x^3 - 2638*nl*x^3 - 121797*x^4 + 
              1898*nl*x^4)/(72*(-1 + x)) - (Pi^2*(63 + 63*x - 5*nl*x + 
              261*x^2 - 5*nl*x^2 - 63*x^3 + 5*nl*x^3 - 270*x^4 + 5*nl*x^4 + 
              135*x^5))/(3*(-1 + x)*(1 + x)))*HPL[{0}, x] + 
         ((-69039 + 862*nl + 74097*x - 426*nl*x - 68139*x^2 + 228*nl*x^2 + 
             75051*x^3 - 1270*nl*x^3)/36 + (Pi^2*(-162 - 2*nl + 5508*x + 
              4*nl*x - 8262*x^2 - 3*nl*x^2 + 5616*x^3 + nl*x^3 - 2754*x^4))/
            (36*(-1 + x)))*HPL[{1}, x] - ((-2754 + 34*nl + 5904*x + 4*nl*x - 
            9531*x^2 + 34*nl*x^2 + 10836*x^3 - 104*nl*x^3 - 5130*x^4 + 
            50*nl*x^4)*HPL[{2}, x])/(3*(-1 + x)) + 
         (8*(135 + 54*x - 5*nl*x + 351*x^2 - 5*nl*x^2 - 54*x^3 + 5*nl*x^3 - 
            351*x^4 + 5*nl*x^4 + 162*x^5)*HPL[{3}, x])/(3*(-1 + x)*(1 + x)) - 
         (18*(2 + 4*x + 6*x^2 + 4*x^3 + 3*x^4)*HPL[{-2, 0}, x])/(1 + x) + 
         (3*(-20 - 38*x - 9*x^2 + 20*x^3 + 5*x^4)*HPL[{-1, 0}, x])/(1 + x) - 
         (72*(1 + x + x^2)^2*HPL[{-1, 2}, x])/(1 + x) - 
         ((-6885 + 102*nl + 7452*x + 162*nl*x - 13365*x^2 + 195*nl*x^2 + 
            5913*x^3 - 308*nl*x^3 + 19008*x^4 - 261*nl*x^4 - 14931*x^5 + 
            182*nl*x^5)*HPL[{0, 0}, x])/(18*(-1 + x)*(1 + x)) - 
         ((-26109 + 272*nl + 52542*x - 326*nl*x - 56619*x^2 + 249*nl*x^2 + 
            50976*x^3 - 341*nl*x^3 - 25137*x^4 + 272*nl*x^4)*HPL[{1, 0}, x])/
          (18*(-1 + x)) - (2*(2772 - 16*nl - 4374*x + 24*nl*x + 3753*x^2 - 
            6*nl*x^2 - 3393*x^3 + 34*nl*x^3)*HPL[{1, 1}, x])/3 + 
         (846*(1 - x + x^2)^2*HPL[{1, 2}, x])/(-1 + x) + 
         ((1512 + 324*x - 32*nl*x + 3240*x^2 - 32*nl*x^2 - 324*x^3 + 
            32*nl*x^3 - 3213*x^4 + 32*nl*x^4 + 1431*x^5)*HPL[{2, 0}, x])/
          (3*(-1 + x)*(1 + x)) + (2*(423 - 414*x - 8*nl*x + 1269*x^2 - 
            1278*x^3 + 8*nl*x^3 + 423*x^4)*HPL[{2, 1}, x])/(-1 + x) + 
         (54*(1 + 2*x + 3*x^2)*HPL[{-1, -1, 0}, x])/(1 + x) - 
         (9*(9 + 18*x + 27*x^2 + 8*x^3 + 4*x^4)*HPL[{-1, 0, 0}, x])/(1 + x) + 
         (2*(54 + 81*x - 4*nl*x + 216*x^2 - 4*nl*x^2 - 81*x^3 + 4*nl*x^3 - 
            207*x^4 + 4*nl*x^4 + 99*x^5)*HPL[{0, 0, 0}, x])/
          ((-1 + x)*(1 + x)) + ((1944 - 2*nl - 3888*x + 4*nl*x + 5832*x^2 - 
            3*nl*x^2 - 3672*x^3 + nl*x^3 + 1944*x^4)*HPL[{1, 0, 0}, x])/
          (6*(-1 + x)) + ((4914 - 2*nl - 9828*x + 4*nl*x + 14742*x^2 - 
            3*nl*x^2 - 9720*x^3 + nl*x^3 + 4914*x^4)*HPL[{1, 1, 0}, x])/
          (6*(-1 + x)) + 1440*x*(2 - x + x^2)*HPL[{1, 1, 1}, x] + 
         ((501 - 26*nl)/3 - 72*Pi^2)*PlusDistribution[1, 1 - x] + 
         3*(-69 + 2*nl)*PlusDistribution[2, 1 - x] + 
         240*PlusDistribution[3, 1 - x] - ((-1188 - 2*nl - 1620*x - 22*nl*x + 
            1161*x^2 - 23*nl*x^2 + 135*x^3 + 22*nl*x^3 - 1053*x^4 + 
            25*nl*x^4 + 1107*x^5)*Zeta[3])/(3*(-1 + x)*(1 + x)) + 
         PlusDistribution[0, 1 - x]*((-1491 + 86*nl)/18 + 
           ((543 - 22*nl)*Pi^2)/24 + 495*Zeta[3]) + 
         PlusDistribution[-1, 1 - x]*((7911 - 1084*nl)/72 + 
           ((501 - 26*nl)*Pi^2)/36 - Pi^4/8 + ((-504 + 5*nl)*Zeta[3])/3)) + 
       PlusDistribution[-1, 1 - x]*((226719 - 32054*nl)/864 + 
         ((-51 + 64*nl)*Pi^2)/54 + ((-19683 + 988*nl)*Pi^4)/6480 + 
         ((-6309 + 67*nl)/36 + (609*Pi^2)/8)*Zeta[3] - (1161*Zeta[5])/2))
 
rsigc["NNLO", "qb"] = (4*(-1 + x)*(-723 + 46*nl + 4983*x - 176*nl*x - 
        1953*x^2 + 82*nl*x^2))/243 + (8*Pi^2*(5 + 75*x - 102*x^2 + 47*x^3))/
      243 - (4*Pi^2*(2 + 2*x + x^2)*HPL[{-1}, x])/81 + 
     (4*(-1 + x)*(75 + 8*nl + 967*x - 16*nl*x - 384*x^2 + 8*nl*x^2)*
       HPL[{1}, x])/81 - (8*(76 + 72*x - 189*x^2 + 64*x^3)*HPL[{2}, x])/81 - 
     (16*(2 + x)^2*HPL[{3}, x])/9 + (16*x*(2 + x + 6*x^2)*HPL[{-1, 0}, x])/
      27 + (16*(-1 + x)*(3 + x)*HPL[{0, 0}, mu^2/Mh^2])/9 + 
     HPL[{0}, x]*((4*Pi^2*(2 + x)^2)/27 + 
       (2*(-135 - 24*nl - 1448*x + 120*nl*x + 2111*x^2 - 96*nl*x^2 - 
          768*x^3 + 32*nl*x^3))/81 - (8*(2 + x)^2*HPL[{0, 0}, mu^2/Mh^2])/
        9) - (16*(27 + 39*x - 81*x^2 + 44*x^3)*HPL[{0, 0}, x])/81 + 
     HPL[{0}, mu^2/Mh^2]*((8*Pi^2*(2 + x)^2)/27 + 
       (2*(-1 + x)*(63 + 24*nl + 1055*x - 48*nl*x - 396*x^2 + 24*nl*x^2))/
        81 - (4*(108 + 24*x - 141*x^2 + 64*x^3)*HPL[{0}, x])/81 - 
       (32*(-1 + x)*(-19 - 25*x + 8*x^2)*HPL[{1}, x])/81 - 
       (16*(2 + x)^2*HPL[{2}, x])/9 - (8*(2 + x)^2*HPL[{0, 0}, x])/9) - 
     (16*(-1 + x)*(-44 - 101*x + 37*x^2)*HPL[{1, 0}, x])/81 - 
     (64*(-1 + x)*(-14 - 35*x + 13*x^2)*HPL[{1, 1}, x])/81 - 
     (8*(2 + x)^2*HPL[{2, 0}, x])/3 - (32*(2 + x)^2*HPL[{2, 1}, x])/9 - 
     (16*(2 + 2*x + x^2)*HPL[{-1, -1, 0}, x])/27 + 
     (8*(2 + 2*x + x^2)*HPL[{-1, 0, 0}, x])/9 - 
     (8*(2 + x)^2*HPL[{0, 0, 0}, x])/9 - (16*(2 + 2*x + x^2)*Zeta[3])/27 + 
     ep*((-2*Pi^4*(16 + 16*x + 9*x^2))/405 + 
       (Pi^2*(1145 + 168*nl + 12224*x - 888*nl*x - 19241*x^2 + 696*nl*x^2 + 
          6388*x^3 - 232*nl*x^3))/486 - (8*Pi^2*(2 + 2*x + x^2)*HPL[{-2}, x])/
        81 - (8*(300 + 696*x - 1059*x^2 + 512*x^3)*HPL[{3}, x])/81 - 
       (16*(2 + x)^2*HPL[{4}, x])/3 + (16*x*(2 + x + 6*x^2)*HPL[{-2, 0}, x])/
        27 + (44*Pi^2*(2 + 2*x + x^2)*HPL[{-1, -1}, x])/81 + 
       ((-52*Pi^2*(2 + 2*x + x^2))/81 + (32*x*(1 - x + 2*x^2))/27)*
        HPL[{-1, 0}, x] + (64*x*(2 + x + 6*x^2)*HPL[{-1, 2}, x])/27 + 
       (32*(2 + 2*x + x^2)*HPL[{-1, 3}, x])/9 + HPL[{2}, x]*
        ((44*Pi^2*(2 + x)^2)/27 + (4*(-1591 - 120*nl - 10968*x + 648*nl*x + 
            18276*x^2 - 504*nl*x^2 - 6032*x^3 + 168*nl*x^3))/243 - 
         (16*(2 + x)^2*HPL[{0, 0}, mu^2/Mh^2])/3) + 
       HPL[{1}, x]*((8*Pi^2*(-1 + x)*(-49 - 595*x + 248*x^2))/243 + 
         (2*(-1 + x)*(247 + 80*nl + 14955*x - 384*nl*x - 5816*x^2 + 
            176*nl*x^2))/81 - (32*(-1 + x)*(-19 - 25*x + 8*x^2)*
           HPL[{0, 0}, mu^2/Mh^2])/27) + ((14*Pi^2*(2 + x)^2)/27 + 
         (2*(-1647 - 216*nl - 11640*x + 1080*nl*x + 21639*x^2 - 864*nl*x^2 - 
            7772*x^3 + 288*nl*x^3))/243)*HPL[{0, 0}, x] + 
       HPL[{0, 0}, mu^2/Mh^2]*((8*Pi^2*(2 + x)^2)/9 + 
         (2*(-1 + x)*(135 + 24*nl + 1079*x - 48*nl*x - 396*x^2 + 24*nl*x^2))/
          27 - (8*(2 + x)^2*HPL[{0, 0}, x])/3) + 
       (8*(-1 + x)*(731 + 72*nl + 6374*x - 144*nl*x - 2509*x^2 + 72*nl*x^2)*
         HPL[{1, 0}, x])/243 + (8*(-1 + x)*(363 + 24*nl + 3161*x - 48*nl*x - 
          1230*x^2 + 24*nl*x^2)*HPL[{1, 1}, x])/81 - 
       (64*(-1 + x)*(-37 - 178*x + 71*x^2)*HPL[{1, 2}, x])/81 - 
       (8*(328 + 648*x - 981*x^2 + 340*x^3)*HPL[{2, 0}, x])/81 - 
       (16*(228 + 528*x - 789*x^2 + 280*x^3)*HPL[{2, 1}, x])/81 - 
       (128*(2 + x)^2*HPL[{2, 2}, x])/9 - (64*(2 + x)^2*HPL[{3, 0}, x])/9 - 
       (32*(2 + x)^2*HPL[{3, 1}, x])/3 - 
       (32*(2 + 2*x + x^2)*HPL[{-2, -1, 0}, x])/27 + 
       (16*(2 + 2*x + x^2)*HPL[{-2, 0, 0}, x])/9 - 
       (16*(2 + 2*x + x^2)*HPL[{-1, -2, 0}, x])/27 - 
       (32*x*(1 + 3*x^2)*HPL[{-1, -1, 0}, x])/27 - 
       (64*(2 + 2*x + x^2)*HPL[{-1, -1, 2}, x])/27 + 
       (16*x*(7 + 2*x + 21*x^2)*HPL[{-1, 0, 0}, x])/27 + 
       (16*(-1 + x)*(3 + x)*HPL[{0, 0, 0}, mu^2/Mh^2])/3 - 
       (16*(33 + 46*x - 79*x^2 + 47*x^3)*HPL[{0, 0, 0}, x])/27 - 
       (16*(-1 + x)*(-94 - 379*x + 149*x^2)*HPL[{1, 0, 0}, x])/81 - 
       (256*(-1 + x)*(-1 - 5*x + 2*x^2)*HPL[{1, 1, 0}, x])/9 - 
       (64*(-1 + x)*(-53 - 272*x + 109*x^2)*HPL[{1, 1, 1}, x])/81 - 
       8*(2 + x)^2*HPL[{2, 0, 0}, x] - (128*(2 + x)^2*HPL[{2, 1, 0}, x])/9 - 
       (64*(2 + x)^2*HPL[{2, 1, 1}, x])/3 + 
       (16*(2 + 2*x + x^2)*HPL[{-1, -1, -1, 0}, x])/9 - 
       (104*(2 + 2*x + x^2)*HPL[{-1, -1, 0, 0}, x])/27 + 
       (80*(2 + 2*x + x^2)*HPL[{-1, 0, 0, 0}, x])/27 - 
       (8*(2 + x)^2*HPL[{0, 0, 0, 0}, x])/3 + 
       ((-1 + x)*(-13131 + 940*nl + 243261*x - 7088*nl*x - 86094*x^2 + 
          3172*nl*x^2 + 648*PolyGamma[2, 1] + 216*x*PolyGamma[2, 1]))/729 + 
       (16*(62 + 243*x - 354*x^2 + 229*x^3)*Zeta[3])/81 + 
       HPL[{-1}, x]*((-8*Pi^2*x*(9 + 4*x + 27*x^2))/81 - 
         (8*Pi^2*(2 + 2*x + x^2)*HPL[{0}, mu^2/Mh^2])/81 - 
         (88*(2 + 2*x + x^2)*Zeta[3])/27) + HPL[{0}, x]*
        ((Pi^2*(140 + 1328*x - 1665*x^2 + 888*x^3))/81 - 
         (4*(132 + 48*x - 135*x^2 + 64*x^3)*HPL[{0, 0}, mu^2/Mh^2])/27 - 
         (8*(2 + x)^2*HPL[{0, 0, 0}, mu^2/Mh^2])/3 + 
         (-9 - 464*nl - 46160*x + 2640*nl*x + 81857*x^2 - 2904*nl*x^2 - 
           27492*x^3 + 952*nl*x^3 - 144*PolyGamma[2, 1] - 
           144*x*PolyGamma[2, 1] - 36*x^2*PolyGamma[2, 1])/243 + 
         (32*(6 + 6*x + x^2)*Zeta[3])/27) + HPL[{0}, mu^2/Mh^2]*
        (((-1 + x)*(441 + 464*nl + 47237*x - 1600*nl*x - 17688*x^2 + 
            752*nl*x^2))/243 + (4*Pi^2*(141 + 510*x - 618*x^2 + 284*x^3))/
          243 + ((14*Pi^2*(2 + x)^2)/27 + (2*(-1647 - 216*nl - 11328*x + 
              936*nl*x + 17988*x^2 - 792*nl*x^2 - 6052*x^3 + 264*nl*x^3))/
            243)*HPL[{0}, x] + (4*(-1 + x)*(1159 + 120*nl + 9439*x - 
            240*nl*x - 3620*x^2 + 120*nl*x^2)*HPL[{1}, x])/243 - 
         (8*(300 + 336*x - 597*x^2 + 224*x^3)*HPL[{2}, x])/81 - 
         (16*(2 + x)^2*HPL[{3}, x])/3 + (32*x*(2 + x + 6*x^2)*
           HPL[{-1, 0}, x])/27 - (4*(396 + 408*x - 771*x^2 + 416*x^3)*
           HPL[{0, 0}, x])/81 - (32*(-1 + x)*(-55 - 142*x + 53*x^2)*
           HPL[{1, 0}, x])/81 - (64*(-1 + x)*(-13 - 37*x + 14*x^2)*
           HPL[{1, 1}, x])/27 - (64*(2 + x)^2*HPL[{2, 0}, x])/9 - 
         (32*(2 + x)^2*HPL[{2, 1}, x])/3 - (32*(2 + 2*x + x^2)*
           HPL[{-1, -1, 0}, x])/27 + (16*(2 + 2*x + x^2)*HPL[{-1, 0, 0}, x])/
          9 - (8*(2 + x)^2*HPL[{0, 0, 0}, x])/3 + 
         (16*(8 + 8*x + x^2)*Zeta[3])/27))
 
rsigc["NNLO", "qg"] = (Pi^2*(221 - 738*x + 144*x^2 - 292*x^3))/162 + 
     (-199423 + 3180*nl + 176196*x - 5016*nl*x + 10743*x^2 + 2148*nl*x^2 + 
       6148*x^3)/1944 + Pi^2*(2 + 2*x + x^2)*HPL[{-1}, x] + 
     ((-1665 + 12*nl + 2040*x - 12*nl*x - 174*x^2 + 6*nl*x^2 + 400*x^3)*
       HPL[{2}, x])/27 - (8*(38 + 21*x + 39*x^2)*HPL[{3}, x])/9 - 
     2*(2 + 2*x + x^2)*HPL[{-2, 0}, x] + ((-166 - 222*x - 33*x^2 - 4*x^3)*
       HPL[{-1, 0}, x])/27 - 4*(2 + 2*x + x^2)*HPL[{-1, 2}, x] + 
     ((-345 + 4*nl + 290*x - 4*nl*x - 8*x^2 + 2*nl*x^2 + 36*x^3)*
       HPL[{0, 0}, mu^2/Mh^2])/9 + HPL[{1}, x]*
      ((59*Pi^2*(2 - 2*x + x^2))/54 + (-21295 + 234*nl + 14796*x - 288*nl*x + 
         3834*x^2 + 162*nl*x^2 + 784*x^3)/162 - 
       (62*(2 - 2*x + x^2)*HPL[{0, 0}, mu^2/Mh^2])/9) + 
     HPL[{0}, x]*((2*Pi^2*(13 + 35*x + 32*x^2))/27 + 
       (-22042 + 348*nl + 26424*x - 456*nl*x + 7143*x^2 + 228*nl*x^2 + 
         2464*x^3)/324 - (4*(27 + 25*x + 28*x^2)*HPL[{0, 0}, mu^2/Mh^2])/9) + 
     ((-837 + 12*nl + 1464*x - 12*nl*x - 234*x^2 + 6*nl*x^2 + 226*x^3)*
       HPL[{0, 0}, x])/27 + ((-5935 + 24*nl + 4956*x - 24*nl*x + 363*x^2 + 
        12*nl*x^2 + 652*x^3)*HPL[{1, 0}, x])/54 + 
     ((-2592 + 4*nl + 2278*x - 4*nl*x + 111*x^2 + 2*nl*x^2 + 288*x^3)*
       HPL[{1, 1}, x])/18 + HPL[{0}, mu^2/Mh^2]*
      ((Pi^2*(62 + 128*x + 89*x^2))/27 + (-11111 + 174*nl + 7974*x - 
         228*nl*x + 1701*x^2 + 114*nl*x^2 + 392*x^3)/162 + 
       ((-837 + 12*nl + 1032*x - 12*nl*x - 132*x^2 + 6*nl*x^2 + 200*x^3)*
         HPL[{0}, x])/27 + ((-687 + 4*nl + 616*x - 4*nl*x + 2*x^2 + 
          2*nl*x^2 + 72*x^3)*HPL[{1}, x])/9 - 
       (4*(76 + 28*x + 67*x^2)*HPL[{2}, x])/9 - 2*(2 + 2*x + x^2)*
        HPL[{-1, 0}, x] - (2*(54 + 64*x + 67*x^2)*HPL[{0, 0}, x])/9 - 
       (106*(2 - 2*x + x^2)*HPL[{1, 0}, x])/9 - 
       (62*(2 - 2*x + x^2)*HPL[{1, 1}, x])/3) - 
     (239*(2 - 2*x + x^2)*HPL[{1, 2}, x])/9 - 
     (4*(107 + 58*x + 97*x^2)*HPL[{2, 0}, x])/9 + 
     ((-642 - 190*x - 553*x^2)*HPL[{2, 1}, x])/9 + 
     4*(2 + 2*x + x^2)*HPL[{-1, -1, 0}, x] - 6*(2 + 2*x + x^2)*
      HPL[{-1, 0, 0}, x] + ((-108 - 142*x - 145*x^2)*HPL[{0, 0, 0}, x])/9 - 
     (107*(2 - 2*x + x^2)*HPL[{1, 0, 0}, x])/9 - 
     (230*(2 - 2*x + x^2)*HPL[{1, 1, 0}, x])/9 - 
     (367*(2 - 2*x + x^2)*HPL[{1, 1, 1}, x])/9 + 
     (2*(46 - 81*x + 45*x^2)*Zeta[3])/9 + 
     ep*((Pi^4*(95 + 66*x - 15*x^2))/405 + 
       (Pi^2*(50125 - 522*nl - 90462*x + 684*nl*x - 31407*x^2 - 342*nl*x^2 - 
          15672*x^3))/1944 + (-2246949 + 30868*nl + 1831952*x - 49192*nl*x + 
         289297*x^2 + 20228*nl*x^2 + 79956*x^3)/3888 + 
       2*Pi^2*(2 + 2*x + x^2)*HPL[{-2}, x] + 
       ((-6033 + 48*nl + 7266*x - 48*nl*x - 2415*x^2 + 24*nl*x^2 + 1160*x^3)*
         HPL[{3}, x])/27 - (2*(146 + 90*x + 153*x^2)*HPL[{4}, x])/3 - 
       4*(2 + 2*x + x^2)*HPL[{-3, 0}, x] - (8*(1 + x)*(55 + 14*x + x^2)*
         HPL[{-2, 0}, x])/27 - 8*(2 + 2*x + x^2)*HPL[{-2, 2}, x] - 
       4*Pi^2*(2 + 2*x + x^2)*HPL[{-1, -1}, x] - 
       (16*(1 + x)*(55 + 14*x + x^2)*HPL[{-1, 2}, x])/27 - 
       24*(2 + 2*x + x^2)*HPL[{-1, 3}, x] + HPL[{-1, 0}, x]*
        ((9*Pi^2*(2 + 2*x + x^2))/2 + (-2494 - 3342*x - 813*x^2 - 46*x^3)/
          81 - 4*(2 + 2*x + x^2)*HPL[{0, 0}, mu^2/Mh^2]) + 
       HPL[{2}, x]*((2*Pi^2*(423 + 167*x + 371*x^2))/27 + 
         (2*(-17850 + 162*nl + 20115*x - 216*nl*x + 4377*x^2 + 99*nl*x^2 + 
            2564*x^3))/81 - (2*(146 + 62*x + 131*x^2)*HPL[{0, 0}, mu^2/Mh^2])/
          3) + ((Pi^2*(318 + 392*x + 461*x^2))/54 + 
         (-71670 + 948*nl + 107544*x - 1272*nl*x + 25431*x^2 + 564*nl*x^2 + 
           12724*x^3)/324)*HPL[{0, 0}, x] + ((391*Pi^2*(2 - 2*x + x^2))/54 + 
         (-51802 + 324*nl + 36831*x - 432*nl*x + 9807*x^2 + 198*nl*x^2 + 
           3751*x^3)/81)*HPL[{1, 0}, x] + ((887*Pi^2*(2 - 2*x + x^2))/54 + 
         (-50225 + 282*nl + 36690*x - 372*nl*x + 8649*x^2 + 180*nl*x^2 + 
           3488*x^3)/54)*HPL[{1, 1}, x] + HPL[{0, 0}, mu^2/Mh^2]*
        ((Pi^2*(62 + 134*x + 89*x^2))/9 + (-13127 + 158*nl + 9768*x - 
           212*nl*x + 1791*x^2 + 94*nl*x^2 + 872*x^3)/54 - 
         (4*(81 + 89*x + 95*x^2)*HPL[{0, 0}, x])/9 - 
         (100*(2 - 2*x + x^2)*HPL[{1, 0}, x])/3 - 
         (548*(2 - 2*x + x^2)*HPL[{1, 1}, x])/9) + 
       ((-17678 + 60*nl + 15672*x - 60*nl*x - 117*x^2 + 30*nl*x^2 + 1736*x^3)*
         HPL[{1, 2}, x])/27 - (242*(2 - 2*x + x^2)*HPL[{1, 3}, x])/3 + 
       ((-7939 + 48*nl + 6060*x - 48*nl*x - 1926*x^2 + 24*nl*x^2 + 1384*x^3)*
         HPL[{2, 0}, x])/27 + (2*(-1989 + 10*nl + 1794*x - 10*nl*x - 
          520*x^2 + 5*nl*x^2 + 352*x^3)*HPL[{2, 1}, x])/9 - 
       (16*(158 + 59*x + 137*x^2)*HPL[{2, 2}, x])/9 - 
       (14*(84 + 50*x + 83*x^2)*HPL[{3, 0}, x])/9 - 
       (2*(962 + 398*x + 917*x^2)*HPL[{3, 1}, x])/9 + 
       8*(2 + 2*x + x^2)*HPL[{-2, -1, 0}, x] - 12*(2 + 2*x + x^2)*
        HPL[{-2, 0, 0}, x] + 8*(2 + 2*x + x^2)*HPL[{-1, -2, 0}, x] + 
       (16*(1 + x)*(55 + 14*x + x^2)*HPL[{-1, -1, 0}, x])/27 + 
       16*(2 + 2*x + x^2)*HPL[{-1, -1, 2}, x] - 
       (8*(1 + x)*(55 + 14*x + x^2)*HPL[{-1, 0, 0}, x])/9 - 
       8*(2 + 2*x + x^2)*HPL[{-1, 2, 0}, x] - 16*(2 + 2*x + x^2)*
        HPL[{-1, 2, 1}, x] + ((-345 + 4*nl + 290*x - 4*nl*x - 8*x^2 + 
          2*nl*x^2 + 36*x^3)*HPL[{0, 0, 0}, mu^2/Mh^2])/3 + 
       ((-2835 + 36*nl + 3726*x - 36*nl*x - 1176*x^2 + 18*nl*x^2 + 610*x^3)*
         HPL[{0, 0, 0}, x])/27 + ((-19655 + 96*nl + 16596*x - 96*nl*x + 
          483*x^2 + 48*nl*x^2 + 1964*x^3)*HPL[{1, 0, 0}, x])/54 + 
       ((-17296 + 60*nl + 15180*x - 60*nl*x - 3*x^2 + 30*nl*x^2 + 1732*x^3)*
         HPL[{1, 1, 0}, x])/27 + ((-17452 + 44*nl + 15790*x - 44*nl*x - 
          297*x^2 + 22*nl*x^2 + 1728*x^3)*HPL[{1, 1, 1}, x])/18 - 
       (1436*(2 - 2*x + x^2)*HPL[{1, 1, 2}, x])/9 - 
       (854*(2 - 2*x + x^2)*HPL[{1, 2, 0}, x])/9 - 
       (1453*(2 - 2*x + x^2)*HPL[{1, 2, 1}, x])/9 - 
       (2*(644 + 382*x + 583*x^2)*HPL[{2, 0, 0}, x])/9 - 
       (2*(1210 + 490*x + 1069*x^2)*HPL[{2, 1, 0}, x])/9 + 
       ((-3910 - 1082*x - 3347*x^2)*HPL[{2, 1, 1}, x])/9 - 
       16*(2 + 2*x + x^2)*HPL[{-1, -1, -1, 0}, x] + 24*(2 + 2*x + x^2)*
        HPL[{-1, -1, 0, 0}, x] - 18*(2 + 2*x + x^2)*HPL[{-1, 0, 0, 0}, x] + 
       ((-324 - 398*x - 413*x^2)*HPL[{0, 0, 0, 0}, x])/9 - 
       (295*(2 - 2*x + x^2)*HPL[{1, 0, 0, 0}, x])/9 - 
       (664*(2 - 2*x + x^2)*HPL[{1, 1, 0, 0}, x])/9 - 
       (460*(2 - 2*x + x^2)*HPL[{1, 1, 1, 0}, x])/3 - 
       (775*(2 - 2*x + x^2)*HPL[{1, 1, 1, 1}, x])/3 + 
       ((2555 - 24*nl - 6738*x + 24*nl*x + 1506*x^2 - 12*nl*x^2 - 616*x^3)*
         Zeta[3])/27 + HPL[{1}, x]*
        ((Pi^2*(21983 - 36*nl - 18816*x + 36*nl*x - 522*x^2 - 18*nl*x^2 - 
            2384*x^3))/324 + (-746693 + 6930*nl + 566334*x - 10332*nl*x + 
           123777*x^2 + 4518*nl*x^2 + 38924*x^3)/972 + 
         ((-2227 + 16*nl + 2014*x - 16*nl*x - 87*x^2 + 8*nl*x^2 + 216*x^3)*
           HPL[{0, 0}, mu^2/Mh^2])/9 - (62*(2 - 2*x + x^2)*
           HPL[{0, 0, 0}, mu^2/Mh^2])/3 + (151*(2 - 2*x + x^2)*Zeta[3])/3) + 
       HPL[{-1}, x]*((4*Pi^2*(1 + x)*(55 + 14*x + x^2))/27 + 
         2*Pi^2*(2 + 2*x + x^2)*HPL[{0}, mu^2/Mh^2] + 30*(2 + 2*x + x^2)*
          Zeta[3]) + HPL[{0}, x]*((Pi^2*(4733 - 36*nl - 10548*x + 36*nl*x + 
            3612*x^2 - 18*nl*x^2 - 1472*x^3))/324 + 
         (-734282 + 9960*nl + 954696*x - 14928*nl*x + 215703*x^2 + 
           6312*nl*x^2 + 87900*x^3)/1944 + 
         ((-945 + 12*nl + 832*x - 12*nl*x - 275*x^2 + 6*nl*x^2 + 176*x^3)*
           HPL[{0, 0}, mu^2/Mh^2])/9 - (4*(27 + 25*x + 28*x^2)*
           HPL[{0, 0, 0}, mu^2/Mh^2])/3 + (4*(95 + 17*x + 113*x^2)*Zeta[3])/
          9) + HPL[{0}, mu^2/Mh^2]*((Pi^2*(4175 - 36*nl - 3942*x + 36*nl*x + 
            2598*x^2 - 18*nl*x^2 - 1572*x^3))/324 + 
         (-382603 + 4980*nl + 298458*x - 7464*nl*x + 54027*x^2 + 
           3156*nl*x^2 + 19462*x^3)/972 + ((Pi^2*(53 + 63*x + 75*x^2))/9 + 
           (-35835 + 474*nl + 41796*x - 636*nl*x + 7854*x^2 + 282*nl*x^2 + 
             5128*x^3)/162)*HPL[{0}, x] + ((35*Pi^2*(2 - 2*x + x^2))/6 + 
           (2*(-6541 + 54*nl + 4941*x - 72*nl*x + 936*x^2 + 33*nl*x^2 + 
              436*x^3))/27)*HPL[{1}, x] + 
         ((-2011 + 16*nl + 1814*x - 16*nl*x - 547*x^2 + 8*nl*x^2 + 352*x^3)*
           HPL[{2}, x])/9 - (2*(438 + 242*x + 437*x^2)*HPL[{3}, x])/9 - 
         4*(2 + 2*x + x^2)*HPL[{-2, 0}, x] - (8*(1 + x)*(55 + 14*x + x^2)*
           HPL[{-1, 0}, x])/27 - 8*(2 + 2*x + x^2)*HPL[{-1, 2}, x] + 
         ((-2835 + 36*nl + 3276*x - 36*nl*x - 1095*x^2 + 18*nl*x^2 + 580*x^3)*
           HPL[{0, 0}, x])/27 + ((-8875 + 48*nl + 7800*x - 48*nl*x - 81*x^2 + 
            24*nl*x^2 + 868*x^3)*HPL[{1, 0}, x])/27 + 
         ((-4500 + 20*nl + 4122*x - 20*nl*x - 159*x^2 + 10*nl*x^2 + 432*x^3)*
           HPL[{1, 1}, x])/9 - (706*(2 - 2*x + x^2)*HPL[{1, 2}, x])/9 - 
         (4*(294 + 140*x + 263*x^2)*HPL[{2, 0}, x])/9 - 
         (2*(962 + 286*x + 829*x^2)*HPL[{2, 1}, x])/9 + 
         8*(2 + 2*x + x^2)*HPL[{-1, -1, 0}, x] - 12*(2 + 2*x + x^2)*
          HPL[{-1, 0, 0}, x] - (2*(54 + 64*x + 67*x^2)*HPL[{0, 0, 0}, x])/3 - 
         (302*(2 - 2*x + x^2)*HPL[{1, 0, 0}, x])/9 - 
         (688*(2 - 2*x + x^2)*HPL[{1, 1, 0}, x])/9 - 
         (1190*(2 - 2*x + x^2)*HPL[{1, 1, 1}, x])/9 + 
         (2*(190 - 52*x + 197*x^2)*Zeta[3])/9))
 
rsigc["NNLO", "qp"] = (2*(-1 + x)*(105 + 11*x))/9 - 
     (4*Pi^2*(-6 + 4*x + 3*x^2))/27 + (4*(-1 + x)*(17 + 5*x)*HPL[{1}, x])/3 + 
     (8*(-12 + 8*x + 5*x^2)*HPL[{2}, x])/9 - (16*(2 + x)^2*HPL[{3}, x])/9 + 
     (16*(-1 + x)*(3 + x)*HPL[{0, 0}, mu^2/Mh^2])/9 + 
     HPL[{0}, x]*((4*Pi^2*(2 + x)^2)/27 + (2*(-59 + 44*x + 29*x^2))/9 - 
       (8*(2 + x)^2*HPL[{0, 0}, mu^2/Mh^2])/9) + 
     (16*(-3 + 4*x + x^2)*HPL[{0, 0}, x])/9 + HPL[{0}, mu^2/Mh^2]*
      ((8*Pi^2*(2 + x)^2)/27 + (2*(-1 + x)*(17 + 5*x))/3 + 
       (4*(-12 + 8*x + 5*x^2)*HPL[{0}, x])/9 + 
       (32*(-1 + x)*(3 + x)*HPL[{1}, x])/9 - (16*(2 + x)^2*HPL[{2}, x])/9 - 
       (8*(2 + x)^2*HPL[{0, 0}, x])/9) + (16*(-1 + x)*(3 + x)*HPL[{1, 0}, x])/
      3 + (64*(-1 + x)*(3 + x)*HPL[{1, 1}, x])/9 - 
     (8*(2 + x)^2*HPL[{2, 0}, x])/3 - (32*(2 + x)^2*HPL[{2, 1}, x])/9 - 
     (8*(2 + x)^2*HPL[{0, 0, 0}, x])/9 + 
     ep*((2*Pi^4*(2 + x)^2)/405 + 2*(-1 + x)*(65 + 8*x) + 
       (Pi^2*(469 - 348*x - 245*x^2))/54 + (8*(-44 + 32*x + 11*x^2)*
         HPL[{3}, x])/9 - (16*(2 + x)^2*HPL[{4}, x])/3 + 
       HPL[{2}, x]*((44*Pi^2*(2 + x)^2)/27 + (4*(-193 + 140*x + 84*x^2))/9 - 
         (16*(2 + x)^2*HPL[{0, 0}, mu^2/Mh^2])/3) + 
       HPL[{1}, x]*((-88*Pi^2*(-1 + x)*(3 + x))/27 + 
         (2*(-1 + x)*(727 + 131*x))/9 + (32*(-1 + x)*(3 + x)*
           HPL[{0, 0}, mu^2/Mh^2])/3) + ((14*Pi^2*(2 + x)^2)/27 + 
         (2*(-193 + 172*x + 99*x^2))/9)*HPL[{0, 0}, x] + 
       HPL[{0, 0}, mu^2/Mh^2]*((8*Pi^2*(2 + x)^2)/9 + 
         (2*(-1 + x)*(177 + 53*x))/9 - (8*(2 + x)^2*HPL[{0, 0}, x])/3) + 
       (32*(-1 + x)*(10 + 3*x)*HPL[{1, 0}, x])/3 + 
       (8*(-1 + x)*(177 + 53*x)*HPL[{1, 1}, x])/9 + 
       (256*(-1 + x)*(3 + x)*HPL[{1, 2}, x])/9 + 
       (8*(-60 + 20*x + 17*x^2)*HPL[{2, 0}, x])/9 + 
       (16*(-44 + 16*x + 13*x^2)*HPL[{2, 1}, x])/9 - 
       (128*(2 + x)^2*HPL[{2, 2}, x])/9 - (64*(2 + x)^2*HPL[{3, 0}, x])/9 - 
       (32*(2 + x)^2*HPL[{3, 1}, x])/3 + (16*(-1 + x)*(3 + x)*
         HPL[{0, 0, 0}, mu^2/Mh^2])/3 + (8*(-22 + 20*x + 5*x^2)*
         HPL[{0, 0, 0}, x])/9 + 16*(-1 + x)*(3 + x)*HPL[{1, 0, 0}, x] + 
       (256*(-1 + x)*(3 + x)*HPL[{1, 1, 0}, x])/9 + 
       (128*(-1 + x)*(3 + x)*HPL[{1, 1, 1}, x])/3 - 
       8*(2 + x)^2*HPL[{2, 0, 0}, x] - (128*(2 + x)^2*HPL[{2, 1, 0}, x])/9 - 
       (64*(2 + x)^2*HPL[{2, 1, 1}, x])/3 - 
       (8*(2 + x)^2*HPL[{0, 0, 0, 0}, x])/3 - 
       (8*(-18 + 28*x + 5*x^2)*Zeta[3])/9 + HPL[{0}, mu^2/Mh^2]*
        (((-1 + x)*(727 + 131*x))/9 - (4*Pi^2*(-29 + 6*x + 8*x^2))/27 + 
         ((14*Pi^2*(2 + x)^2)/27 + (2*(-193 + 140*x + 84*x^2))/9)*
          HPL[{0}, x] + (4*(-1 + x)*(177 + 53*x)*HPL[{1}, x])/9 + 
         (8*(-44 + 16*x + 13*x^2)*HPL[{2}, x])/9 - (16*(2 + x)^2*HPL[{3}, x])/
          3 + (4*(-44 + 32*x + 11*x^2)*HPL[{0, 0}, x])/9 + 
         (128*(-1 + x)*(3 + x)*HPL[{1, 0}, x])/9 + 
         (64*(-1 + x)*(3 + x)*HPL[{1, 1}, x])/3 - 
         (64*(2 + x)^2*HPL[{2, 0}, x])/9 - (32*(2 + x)^2*HPL[{2, 1}, x])/3 - 
         (8*(2 + x)^2*HPL[{0, 0, 0}, x])/3 + (16*(2 + x)^2*Zeta[3])/9) + 
       HPL[{0}, x]*((Pi^2*(100 - 104*x - 23*x^2))/27 + 
         (-711 + 880*x + 229*x^2)/9 + (4*(-44 + 16*x + 13*x^2)*
           HPL[{0, 0}, mu^2/Mh^2])/9 - (8*(2 + x)^2*HPL[{0, 0, 0}, 
            mu^2/Mh^2])/3 + (8*(2 + x)^2*Zeta[3])/3))
 
rsigc["NNLO", "qq"] = (4*(-1 + x)*(160 + 27*x))/27 - 
     (4*Pi^2*(-6 + 4*x + 3*x^2))/27 + (4*(-1 + x)*(17 + 5*x)*HPL[{1}, x])/3 + 
     (8*(-12 + 8*x + 5*x^2)*HPL[{2}, x])/9 - (16*(2 + x)^2*HPL[{3}, x])/9 + 
     (16*(-1 + x)*(3 + x)*HPL[{0, 0}, mu^2/Mh^2])/9 + 
     HPL[{0}, x]*((4*Pi^2*(2 + x)^2)/27 + (2*(-177 + 124*x + 69*x^2))/27 - 
       (8*(2 + x)^2*HPL[{0, 0}, mu^2/Mh^2])/9) + 
     (8*(-18 + 26*x + 9*x^2)*HPL[{0, 0}, x])/27 + HPL[{0}, mu^2/Mh^2]*
      ((8*Pi^2*(2 + x)^2)/27 + (2*(-1 + x)*(17 + 5*x))/3 + 
       (4*(-12 + 8*x + 5*x^2)*HPL[{0}, x])/9 + 
       (32*(-1 + x)*(3 + x)*HPL[{1}, x])/9 - (16*(2 + x)^2*HPL[{2}, x])/9 - 
       (8*(2 + x)^2*HPL[{0, 0}, x])/9) + (16*(-1 + x)*(3 + x)*HPL[{1, 0}, x])/
      3 + (64*(-1 + x)*(3 + x)*HPL[{1, 1}, x])/9 - 
     (8*(2 + x)^2*HPL[{2, 0}, x])/3 - (32*(2 + x)^2*HPL[{2, 1}, x])/9 - 
     (8*(2 + x)^2*HPL[{0, 0, 0}, x])/9 + 
     (8*(2 - 2*x + x^2)*HPL[{1, 0, 0}, x])/27 - (8*(2 - 2*x + x^2)*Zeta[3])/
      27 + ep*((-4*Pi^4*(-3 + x)*(-1 + 2*x))/405 + ((-1 + x)*(3563 + 669*x))/
        27 + (Pi^2*(1407 - 980*x - 591*x^2))/162 + 
       (8*(-132 + 104*x + 45*x^2)*HPL[{3}, x])/27 - 
       (16*(2 + x)^2*HPL[{4}, x])/3 + HPL[{2}, x]*((44*Pi^2*(2 + x)^2)/27 + 
         (4*(-579 + 404*x + 216*x^2))/27 - 
         (16*(2 + x)^2*HPL[{0, 0}, mu^2/Mh^2])/3) + 
       ((14*Pi^2*(2 + x)^2)/27 + (2*(-579 + 492*x + 251*x^2))/27)*
        HPL[{0, 0}, x] + HPL[{0, 0}, mu^2/Mh^2]*((8*Pi^2*(2 + x)^2)/9 + 
         (2*(-1 + x)*(177 + 53*x))/9 - (8*(2 + x)^2*HPL[{0, 0}, x])/3) + 
       ((32*(-1 + x)*(10 + 3*x))/3 - (16*Pi^2*(2 - 2*x + x^2))/81)*
        HPL[{1, 0}, x] + (8*(-1 + x)*(177 + 53*x)*HPL[{1, 1}, x])/9 + 
       (256*(-1 + x)*(3 + x)*HPL[{1, 2}, x])/9 + 
       (32*(2 - 2*x + x^2)*HPL[{1, 3}, x])/27 + 
       (8*(-60 + 20*x + 17*x^2)*HPL[{2, 0}, x])/9 + 
       (16*(-44 + 16*x + 13*x^2)*HPL[{2, 1}, x])/9 - 
       (128*(2 + x)^2*HPL[{2, 2}, x])/9 - (64*(2 + x)^2*HPL[{3, 0}, x])/9 - 
       (32*(2 + x)^2*HPL[{3, 1}, x])/3 + (16*(-1 + x)*(3 + x)*
         HPL[{0, 0, 0}, mu^2/Mh^2])/3 + (16*(-11 + 11*x + 4*x^2)*
         HPL[{0, 0, 0}, x])/9 + (8*(-162 + 108*x + 53*x^2)*HPL[{1, 0, 0}, x])/
        27 + (256*(-1 + x)*(3 + x)*HPL[{1, 1, 0}, x])/9 + 
       (128*(-1 + x)*(3 + x)*HPL[{1, 1, 1}, x])/3 - 
       (8*(104 + 112*x + 25*x^2)*HPL[{2, 0, 0}, x])/27 - 
       (128*(2 + x)^2*HPL[{2, 1, 0}, x])/9 - (64*(2 + x)^2*HPL[{2, 1, 1}, x])/
        3 - (8*(2 + x)^2*HPL[{0, 0, 0, 0}, x])/3 + 
       (8*(2 - 2*x + x^2)*HPL[{1, 0, 0, 0}, x])/9 + 
       (16*(2 - 2*x + x^2)*HPL[{1, 1, 0, 0}, x])/27 - 
       (16*(-27 + 46*x + 13*x^2)*Zeta[3])/27 + HPL[{1}, x]*
        ((-88*Pi^2*(-1 + x)*(3 + x))/27 + (2*(-1 + x)*(2201 + 477*x))/27 + 
         (32*(-1 + x)*(3 + x)*HPL[{0, 0}, mu^2/Mh^2])/3 - 
         (16*(2 - 2*x + x^2)*Zeta[3])/9) + HPL[{0}, mu^2/Mh^2]*
        (((-1 + x)*(2201 + 477*x))/27 - (4*Pi^2*(-29 + 6*x + 8*x^2))/27 + 
         ((14*Pi^2*(2 + x)^2)/27 + (2*(-579 + 404*x + 216*x^2))/27)*
          HPL[{0}, x] + (4*(-1 + x)*(177 + 53*x)*HPL[{1}, x])/9 + 
         (8*(-44 + 16*x + 13*x^2)*HPL[{2}, x])/9 - (16*(2 + x)^2*HPL[{3}, x])/
          3 + (4*(-132 + 104*x + 45*x^2)*HPL[{0, 0}, x])/27 + 
         (128*(-1 + x)*(3 + x)*HPL[{1, 0}, x])/9 + 
         (64*(-1 + x)*(3 + x)*HPL[{1, 1}, x])/3 - 
         (64*(2 + x)^2*HPL[{2, 0}, x])/9 - (32*(2 + x)^2*HPL[{2, 1}, x])/3 - 
         (8*(2 + x)^2*HPL[{0, 0, 0}, x])/3 + (16*(2 - 2*x + x^2)*
           HPL[{1, 0, 0}, x])/27 + (32*(5 + 7*x + x^2)*Zeta[3])/27) + 
       HPL[{0}, x]*((Pi^2*(300 - 344*x - 117*x^2))/81 + 
         (-2153 + 2544*x + 721*x^2)/27 + (4*(-44 + 16*x + 13*x^2)*
           HPL[{0, 0}, mu^2/Mh^2])/9 - (8*(2 + x)^2*HPL[{0, 0, 0}, 
            mu^2/Mh^2])/3 + (8*(32 + 40*x + 7*x^2)*Zeta[3])/27))
 
Attributes[PlusDistribution] = {NHoldAll}
 
PlusDistribution /: PlusDistribution::usage = "PlusDistribution[i,1-x] is \
generalized function (Log[1-x]^i/(1-x))_+, regularized at x = 1. By \
definition, PlusDistribution[-1,1-x] = delta(1-x)."
