rsigc["LO", "qb"] = TF*x*coeff["Cii[qi,bj]*(vi^2+ai^2)"]*
     (2*PlusDistribution[-1, 1 - x] - 2*ep*PlusDistribution[-1, 1 - x])
 
rsigc["NLO", "qb"] = CF*TF*x*coeff["Cii[qi,bj]*(vi^2+ai^2)*alphas/Pi"]*
     ((2*(1 + x^2)*HPL[{0}, x])/(-1 + x) + 4*(1 + x)*HPL[{1}, x] + 
      (-8 + (2*Pi^2)/3)*PlusDistribution[-1, 1 - x] + 
      ep^2*Lv^3*((1 + x)/3 - PlusDistribution[-1, 1 - x]/2 - 
        (2*PlusDistribution[0, 1 - x])/3) + 8*PlusDistribution[1, 1 - x] + 
      Lv^2*(ep*(1 + x - (3*PlusDistribution[-1, 1 - x])/2 - 
          2*PlusDistribution[0, 1 - x]) + 
        ep^2*(-1 - x + ((1 + x^2)*HPL[{0}, x])/(-1 + x) + 
          2*(1 + x)*HPL[{1}, x] + (-5/2 + Pi^2/3)*PlusDistribution[-1, 
            1 - x] + 2*PlusDistribution[0, 1 - x] + 
          4*PlusDistribution[1, 1 - x])) + 
      Lv*(2*(1 + x) - 3*PlusDistribution[-1, 1 - x] - 
        4*PlusDistribution[0, 1 - x] + ep*(-2*(1 + x) + 
          (2*(1 + x^2)*HPL[{0}, x])/(-1 + x) + 4*(1 + x)*HPL[{1}, x] + 
          (-5 + (2*Pi^2)/3)*PlusDistribution[-1, 1 - x] + 
          4*PlusDistribution[0, 1 - x] + 8*PlusDistribution[1, 1 - x]) + 
        ep^2*((Pi^2*(-1 - x))/2 + 4*(-1 + x) - (2*(1 + x^2)*HPL[{0}, x])/
           (-1 + x) - 4*(1 + x)*HPL[{1}, x] + (4*(1 + x^2)*HPL[{2}, x])/
           (-1 + x) + (2*(1 + x^2)*HPL[{0, 0}, x])/(-1 + x) + 
          (4*(1 + x^2)*HPL[{1, 0}, x])/(-1 + x) + 8*(1 + x)*HPL[{1, 1}, x] + 
          (-8 + (13*Pi^2)/12)*PlusDistribution[-1, 1 - x] + 
          Pi^2*PlusDistribution[0, 1 - x] - 8*PlusDistribution[1, 1 - x] - 
          8*PlusDistribution[2, 1 - x])) + 
      ep*((Pi^2*(-1 - x))/2 + 4*(-1 + x) - (2*(1 + x^2)*HPL[{0}, x])/
         (-1 + x) - 4*(1 + x)*HPL[{1}, x] + (4*(1 + x^2)*HPL[{2}, x])/
         (-1 + x) + (2*(1 + x^2)*HPL[{0, 0}, x])/(-1 + x) + 
        (4*(1 + x^2)*HPL[{1, 0}, x])/(-1 + x) + 8*(1 + x)*HPL[{1, 1}, x] + 
        (-8 + (13*Pi^2)/12)*PlusDistribution[-1, 1 - x] + 
        Pi^2*PlusDistribution[0, 1 - x] - 8*PlusDistribution[1, 1 - x] - 
        8*PlusDistribution[2, 1 - x]) + ep^2*(4*(-1 + x) + (Pi^2*(1 + x))/2 + 
        (4*(-1 + x) - (Pi^2*(1 + x^2))/(2*(-1 + x)))*HPL[{0}, x] + 
        (Pi^2*(-1 - x) + 8*(-1 + x))*HPL[{1}, x] - (4*(1 + x^2)*HPL[{2}, x])/
         (-1 + x) + (4*(1 + x^2)*HPL[{3}, x])/(-1 + x) - 
        (2*(1 + x^2)*HPL[{0, 0}, x])/(-1 + x) - (4*(1 + x^2)*HPL[{1, 0}, x])/
         (-1 + x) - 8*(1 + x)*HPL[{1, 1}, x] + (8*(1 + x^2)*HPL[{1, 2}, x])/
         (-1 + x) + (4*(1 + x^2)*HPL[{2, 0}, x])/(-1 + x) + 
        (8*(1 + x^2)*HPL[{2, 1}, x])/(-1 + x) + 
        (2*(1 + x^2)*HPL[{0, 0, 0}, x])/(-1 + x) + 
        (4*(1 + x^2)*HPL[{1, 0, 0}, x])/(-1 + x) + 
        (8*(1 + x^2)*HPL[{1, 1, 0}, x])/(-1 + x) + 
        16*(1 + x)*HPL[{1, 1, 1}, x] - 2*Pi^2*PlusDistribution[1, 1 - x] + 
        8*PlusDistribution[2, 1 - x] + (16*PlusDistribution[3, 1 - x])/3 - 
        (14*(1 + x)*Zeta[3])/3 + PlusDistribution[-1, 1 - x]*
         (-16 + (35*Pi^2)/12 - (11*Pi^4)/90 + 7*Zeta[3]) + 
        PlusDistribution[0, 1 - x]*(-Pi^2 + (28*Zeta[3])/3)))
 
rsigc["NLO", "qg"] = TF*x*coeff["Sum[k,QB]*Cif[qi,qk]*(vi^2+ai^2)*alphas/Pi"]*
     (-((-1 + x)*(1 + 7*x))/4 + (ep^2*Lv^3*(-1 + 2*x - 2*x^2))/12 + 
      ((-1 + 2*x - 2*x^2)*HPL[{0}, x])/2 + (-1 + 2*x - 2*x^2)*HPL[{1}, x] + 
      Lv^2*((ep*(-1 + 2*x - 2*x^2))/4 + ep^2*((3 + 2*x - 3*x^2)/8 + 
          ((-1 + 2*x - 2*x^2)*HPL[{0}, x])/4 + ((-1 + 2*x - 2*x^2)*
            HPL[{1}, x])/2)) + ep*(-((-1 + x)*(5 + 7*x))/4 + 
        (Pi^2*(1 - 2*x + 2*x^2))/8 + ((3 + 2*x - 3*x^2)*HPL[{0}, x])/4 + 
        ((3 + 2*x - 3*x^2)*HPL[{1}, x])/2 + (-1 + 2*x - 2*x^2)*HPL[{2}, x] + 
        ((-1 + 2*x - 2*x^2)*HPL[{0, 0}, x])/2 + (-1 + 2*x - 2*x^2)*
         HPL[{1, 0}, x] - 2*(1 - 2*x + 2*x^2)*HPL[{1, 1}, x]) + 
      Lv*((-1 + 2*x - 2*x^2)/2 + ep*((3 + 2*x - 3*x^2)/4 + 
          ((-1 + 2*x - 2*x^2)*HPL[{0}, x])/2 + (-1 + 2*x - 2*x^2)*
           HPL[{1}, x]) + ep^2*(-((-1 + x)*(5 + 7*x))/4 + 
          (Pi^2*(1 - 2*x + 2*x^2))/8 + ((3 + 2*x - 3*x^2)*HPL[{0}, x])/4 + 
          ((3 + 2*x - 3*x^2)*HPL[{1}, x])/2 + (-1 + 2*x - 2*x^2)*
           HPL[{2}, x] + ((-1 + 2*x - 2*x^2)*HPL[{0, 0}, x])/2 + 
          (-1 + 2*x - 2*x^2)*HPL[{1, 0}, x] - 2*(1 - 2*x + 2*x^2)*
           HPL[{1, 1}, x])) + ep^2*(-((-1 + x)*(5 + 7*x))/2 + 
        (Pi^2*(-3 - 2*x + 3*x^2))/16 + (-((-1 + x)*(5 + 7*x))/4 + 
          (Pi^2*(1 - 2*x + 2*x^2))/8)*HPL[{0}, x] + 
        (-((-1 + x)*(5 + 7*x))/2 + (Pi^2*(1 - 2*x + 2*x^2))/4)*HPL[{1}, x] + 
        ((3 + 2*x - 3*x^2)*HPL[{2}, x])/2 + (-1 + 2*x - 2*x^2)*HPL[{3}, x] + 
        ((3 + 2*x - 3*x^2)*HPL[{0, 0}, x])/4 + 
        ((3 + 2*x - 3*x^2)*HPL[{1, 0}, x])/2 + (3 + 2*x - 3*x^2)*
         HPL[{1, 1}, x] - 2*(1 - 2*x + 2*x^2)*HPL[{1, 2}, x] + 
        (-1 + 2*x - 2*x^2)*HPL[{2, 0}, x] - 2*(1 - 2*x + 2*x^2)*
         HPL[{2, 1}, x] + ((-1 + 2*x - 2*x^2)*HPL[{0, 0, 0}, x])/2 + 
        (-1 + 2*x - 2*x^2)*HPL[{1, 0, 0}, x] - 2*(1 - 2*x + 2*x^2)*
         HPL[{1, 1, 0}, x] - 4*(1 - 2*x + 2*x^2)*HPL[{1, 1, 1}, x] + 
        (7*(1 - 2*x + 2*x^2)*Zeta[3])/6))
 
rsigc["NNLO", "gg"] = 
    x*coeff["Sum[k,l,Q]*Cff[qk,bl]*(vk^2+ak^2)*alphas^2/Pi^2*TF/CA/CF"]*
     (CA^2*((Pi^2*(1 + x)^2)/36 + ((-1 + x)*(47 + 191*x))/48 + 
        (Pi^2*(1 + x)^2*HPL[{-1}, x])/12 + ((-6 - 38*x - 75*x^2)*HPL[{0}, x])/
         24 + ((1 + x)^2*HPL[{-1, 0}, x])/3 + 
        ((-2 + 2*x + 25*x^2)*HPL[{0, 0}, x])/12 + 
        (1 + x)^2*HPL[{-1, -1, 0}, x] - (3*(1 + x)^2*HPL[{-1, 0, 0}, x])/2 + 
        ((-1 + x)^2*HPL[{1, 0, 0}, x])/2 + ((1 + 6*x + x^2)*Zeta[3])/2 + 
        ep*Lv*((Pi^2*(1 + x)^2)/18 + ((-1 + x)*(47 + 191*x))/24 + 
          (Pi^2*(1 + x)^2*HPL[{-1}, x])/6 + ((-6 - 38*x - 75*x^2)*
            HPL[{0}, x])/12 + (2*(1 + x)^2*HPL[{-1, 0}, x])/3 + 
          ((-2 + 2*x + 25*x^2)*HPL[{0, 0}, x])/6 + 2*(1 + x)^2*
           HPL[{-1, -1, 0}, x] - 3*(1 + x)^2*HPL[{-1, 0, 0}, x] + 
          (-1 + x)^2*HPL[{1, 0, 0}, x] + (1 + 6*x + x^2)*Zeta[3]) + 
        ep*(((-1 + x)*(2021 + 8141*x))/288 + (Pi^4*(1 + 6*x + x^2))/24 + 
          (Pi^2*(46 + 281*x + 469*x^2))/216 + (Pi^2*(1 + x)^2*HPL[{-2}, x])/
           6 + ((-6 - 38*x - 75*x^2)*HPL[{2}, x])/6 + 
          ((-2 + 2*x + 25*x^2)*HPL[{3}, x])/3 + ((1 + x)^2*HPL[{-2, 0}, x])/
           3 - (11*Pi^2*(1 + x)^2*HPL[{-1, -1}, x])/12 + 
          ((13*Pi^2*(1 + x)^2)/12 + (10 + 53*x + 19*x^2)/18)*
           HPL[{-1, 0}, x] + (4*(1 + x)^2*HPL[{-1, 2}, x])/3 - 
          6*(1 + x)^2*HPL[{-1, 3}, x] + ((-74 - 322*x - 227*x^2)*
            HPL[{0, 0}, x])/72 - (Pi^2*(-1 + x)^2*HPL[{1, 0}, x])/3 + 
          2*(-1 + x)^2*HPL[{1, 3}, x] + 2*(1 + x)^2*HPL[{-2, -1, 0}, x] - 
          3*(1 + x)^2*HPL[{-2, 0, 0}, x] + (1 + x)^2*HPL[{-1, -2, 0}, x] + 
          ((2 + 7*x + 2*x^2)*HPL[{-1, -1, 0}, x])/3 + 4*(1 + x)^2*
           HPL[{-1, -1, 2}, x] + ((-2 - 13*x - 2*x^2)*HPL[{-1, 0, 0}, x])/6 + 
          ((-8 + 2*x + 73*x^2)*HPL[{0, 0, 0}, x])/12 + 
          ((1 - 3*x + x^2)*HPL[{1, 0, 0}, x])/2 + (-1 + x)^2*
           HPL[{2, 0, 0}, x] - 3*(1 + x)^2*HPL[{-1, -1, -1, 0}, x] + 
          (13*(1 + x)^2*HPL[{-1, -1, 0, 0}, x])/2 - 5*(1 + x)^2*
           HPL[{-1, 0, 0, 0}, x] + (3*(-1 + x)^2*HPL[{1, 0, 0, 0}, x])/2 + 
          (-1 + x)^2*HPL[{1, 1, 0, 0}, x] - ((1 + 4*x)*(-10 + 11*x)*Zeta[3])/
           6 + HPL[{1}, x]*(((-1 + x)*(47 + 191*x))/12 - 3*(-1 + x)^2*
             Zeta[3]) + HPL[{-1}, x]*((Pi^2*(-2 - 3*x - 2*x^2))/12 + 
            (11*(1 + x)^2*Zeta[3])/2) + HPL[{0}, x]*
           ((-396 - 1544*x - 285*x^2)/144 + (Pi^2*(5 - 2*x - 49*x^2))/36 + 
            (1 + 6*x + x^2)*Zeta[3]))) + CA*CF*((Pi^2*x*(5 + 2*x))/12 + 
        ((-1 + x)*(16 + 49*x))/4 - (Pi^2*(1 + x)^2*HPL[{-1}, x])/6 + 
        ((Pi^2*(1 + 3*x + 3*x^2))/6 + (-23 - 64*x + 105*x^2)/8)*HPL[{0}, x] + 
        ep*Lv^3*(((-1 + x)*(1 + 3*x))/2 - ((1 + 2*x)^2*HPL[{0}, x])/4) + 
        ((-1 + x)*(7 + 67*x)*HPL[{1}, x])/4 + 
        ((-1 - 8*x + 4*x^2)*HPL[{2}, x])/2 - (1 + 2*x)^2*HPL[{3}, x] + 
        (1 + 2*x + 2*x^2)*HPL[{-2, 0}, x] + (1 + x)*HPL[{-1, 0}, x] - 
        ((1 + x)*(3 + 4*x)*HPL[{0, 0}, x])/2 + 
        Lv^2*(((-1 + x)*(1 + 3*x))/2 - ((1 + 2*x)^2*HPL[{0}, x])/4 + 
          ep*((Pi^2*(1 + 2*x)^2)/4 + ((-1 + x)*(13 + 177*x))/16 + 
            ((-1 - 16*x + 20*x^2)*HPL[{0}, x])/8 + 3*(-1 + x)*(1 + 3*x)*
             HPL[{1}, x] - (3*(1 + 2*x)^2*HPL[{2}, x])/2 - 
            (3*(1 + 2*x)^2*HPL[{0, 0}, x])/4)) + 3*(-1 + x)*(1 + 3*x)*
         HPL[{1, 0}, x] + 4*(-1 + x)*(1 + 3*x)*HPL[{1, 1}, x] - 
        (3*(1 + 2*x)^2*HPL[{2, 0}, x])/2 - 2*(1 + 2*x)^2*HPL[{2, 1}, x] - 
        2*(1 + x)^2*HPL[{-1, -1, 0}, x] + 3*(1 + x)^2*HPL[{-1, 0, 0}, x] + 
        ((-3 - 8*x - 8*x^2)*HPL[{0, 0, 0}, x])/2 - 
        (-1 + x)^2*HPL[{1, 0, 0}, x] + ((1 - 6*x + 4*x^2)*Zeta[3])/2 + 
        ep*(((-1 + x)*(104 + 359*x))/8 + (Pi^2*(149 + 372*x - 693*x^2))/96 + 
          (Pi^4*(-21 - 160*x - 10*x^2))/360 + (Pi^2*(-13 - 26*x - 22*x^2)*
            HPL[{-2}, x])/12 + ((11*Pi^2*(1 + 2*x)^2)/12 + 
            (-36 - 156*x + 245*x^2)/4)*HPL[{2}, x] + 
          ((-11 - 28*x - 4*x^2)*HPL[{3}, x])/2 + (-7 - 20*x - 20*x^2)*
           HPL[{4}, x] + (1 + 2*x + 2*x^2)*HPL[{-3, 0}, x] + 
          (1 + 5*x + 2*x^2)*HPL[{-2, 0}, x] + 4*(1 + 2*x + 2*x^2)*
           HPL[{-2, 2}, x] + (11*Pi^2*(1 + x)^2*HPL[{-1, -1}, x])/6 + 
          (3 + 2*x - (13*Pi^2*(1 + x)^2)/6)*HPL[{-1, 0}, x] + 
          4*(1 + x)*HPL[{-1, 2}, x] + 12*(1 + x)^2*HPL[{-1, 3}, x] + 
          ((Pi^2*(25 + 64*x + 64*x^2))/24 + (-63 - 176*x + 211*x^2)/8)*
           HPL[{0, 0}, x] + ((2*Pi^2*(-1 + x)^2)/3 + 2*(-1 + x)*(2 + 29*x))*
           HPL[{1, 0}, x] + ((-1 + x)*(13 + 177*x)*HPL[{1, 1}, x])/2 + 
          16*(-1 + x)*(1 + 3*x)*HPL[{1, 2}, x] - 4*(-1 + x)^2*
           HPL[{1, 3}, x] + ((-1 - 20*x + 28*x^2)*HPL[{2, 0}, x])/2 + 
          (-1 - 16*x + 20*x^2)*HPL[{2, 1}, x] - 8*(1 + 2*x)^2*
           HPL[{2, 2}, x] - 4*(1 + 2*x)^2*HPL[{3, 0}, x] - 
          6*(1 + 2*x)^2*HPL[{3, 1}, x] + (-5 - 10*x - 6*x^2)*
           HPL[{-2, -1, 0}, x] + ((19 + 38*x + 26*x^2)*HPL[{-2, 0, 0}, x])/
           2 - 2*(1 + x)^2*HPL[{-1, -2, 0}, x] - (3 + x)*(1 + 2*x)*
           HPL[{-1, -1, 0}, x] - 8*(1 + x)^2*HPL[{-1, -1, 2}, x] + 
          ((13 + 25*x + 6*x^2)*HPL[{-1, 0, 0}, x])/2 + 
          ((-9 - 26*x - 12*x^2)*HPL[{0, 0, 0}, x])/2 + (-10 - 15*x + 26*x^2)*
           HPL[{1, 0, 0}, x] + 16*(-1 + x)*(1 + 3*x)*HPL[{1, 1, 0}, x] + 
          24*(-1 + x)*(1 + 3*x)*HPL[{1, 1, 1}, x] + 
          ((-13 - 28*x - 40*x^2)*HPL[{2, 0, 0}, x])/2 - 
          8*(1 + 2*x)^2*HPL[{2, 1, 0}, x] - 12*(1 + 2*x)^2*
           HPL[{2, 1, 1}, x] + 6*(1 + x)^2*HPL[{-1, -1, -1, 0}, x] - 
          13*(1 + x)^2*HPL[{-1, -1, 0, 0}, x] + 10*(1 + x)^2*
           HPL[{-1, 0, 0, 0}, x] + (-5 - 13*x - 13*x^2)*HPL[{0, 0, 0, 0}, 
            x] - 3*(-1 + x)^2*HPL[{1, 0, 0, 0}, x] - 2*(-1 + x)^2*
           HPL[{1, 1, 0, 0}, x] + ((14 + 21*x + 18*x^2)*Zeta[3])/2 + 
          HPL[{1}, x]*((-11*Pi^2*(-1 + x)*(1 + 3*x))/6 + 
            (3*(-1 + x)*(27 + 191*x))/8 + 6*(-1 + x)^2*Zeta[3]) + 
          HPL[{-1}, x]*((Pi^2*(-11 - 15*x - 2*x^2))/12 - 11*(1 + x)^2*
             Zeta[3]) + HPL[{0}, x]*((Pi^2*(47 + 108*x + 60*x^2))/48 + 
            (-169 - 448*x + 999*x^2)/16 + 5*(1 + x + 3*x^2)*Zeta[3])) + 
        Lv*((Pi^2*(1 + 2*x)^2)/6 + ((-1 + x)*(7 + 67*x))/8 + 
          ((-1 - 8*x + 4*x^2)*HPL[{0}, x])/4 + 2*(-1 + x)*(1 + 3*x)*
           HPL[{1}, x] - (1 + 2*x)^2*HPL[{2}, x] - 
          ((1 + 2*x)^2*HPL[{0, 0}, x])/2 + ep*((3*(-1 + x)*(27 + 191*x))/16 + 
            (Pi^2*(-2 + 8*x - 5*x^2))/12 - (Pi^2*(1 + x)^2*HPL[{-1}, x])/3 + 
            ((Pi^2*(11 + 36*x + 36*x^2))/24 + (-36 - 156*x + 245*x^2)/8)*
             HPL[{0}, x] + ((-1 + x)*(13 + 177*x)*HPL[{1}, x])/4 + 
            ((-1 - 16*x + 20*x^2)*HPL[{2}, x])/2 - 3*(1 + 2*x)^2*
             HPL[{3}, x] + 2*(1 + 2*x + 2*x^2)*HPL[{-2, 0}, x] + 
            2*(1 + x)*HPL[{-1, 0}, x] + ((-11 - 28*x - 4*x^2)*HPL[{0, 0}, x])/
             4 + 8*(-1 + x)*(1 + 3*x)*HPL[{1, 0}, x] + 12*(-1 + x)*(1 + 3*x)*
             HPL[{1, 1}, x] - 4*(1 + 2*x)^2*HPL[{2, 0}, x] - 
            6*(1 + 2*x)^2*HPL[{2, 1}, x] - 4*(1 + x)^2*HPL[{-1, -1, 0}, x] + 
            6*(1 + x)^2*HPL[{-1, 0, 0}, x] + ((-7 - 20*x - 20*x^2)*
              HPL[{0, 0, 0}, x])/2 - 2*(-1 + x)^2*HPL[{1, 0, 0}, x] + 
            2*(1 - x + 4*x^2)*Zeta[3]))))
 
rsigc["NNLO", "qb"] = 
    CF*TF*x*coeff["delta[i,j]*Sum[k,l,Q]*Cff[qk,bl]*(vk^2+ak^2)*alphas^2/Pi^2\
"]*((-5*(-1 + x)*(1 + x))/3 - (Pi^2*(1 + x)^2)/9 + 
       ((3 + 4*x + 3*x^2)*HPL[{0}, x])/3 - (4*(1 + x)^2*HPL[{-1, 0}, x])/3 + 
       (2*(1 + x)^2*HPL[{0, 0}, x])/3 + ep*Lv*((-10*(-1 + x)*(1 + x))/3 - 
         (2*Pi^2*(1 + x)^2)/9 + (2*(3 + 4*x + 3*x^2)*HPL[{0}, x])/3 - 
         (8*(1 + x)^2*HPL[{-1, 0}, x])/3 + (4*(1 + x)^2*HPL[{0, 0}, x])/3) + 
       ep*((-79*(-1 + x)*(1 + x))/9 + (Pi^2*(-20 - 31*x - 20*x^2))/27 + 
         Pi^2*(1 + x)^2*HPL[{-1}, x] + ((-5*Pi^2*(1 + x)^2)/9 + 
           (99 + 28*x - 51*x^2)/18)*HPL[{0}, x] - 
         (20*(-1 + x)*(1 + x)*HPL[{1}, x])/3 + 
         (4*(3 + 4*x + 3*x^2)*HPL[{2}, x])/3 + (8*(1 + x)^2*HPL[{3}, x])/3 - 
         (4*(1 + x)^2*HPL[{-2, 0}, x])/3 - (4*(2 + 7*x + 2*x^2)*
           HPL[{-1, 0}, x])/9 - (16*(1 + x)^2*HPL[{-1, 2}, x])/3 + 
         ((31 + 50*x + 31*x^2)*HPL[{0, 0}, x])/9 + 
         (4*(1 + x)^2*HPL[{-1, -1, 0}, x])/3 - 
         (14*(1 + x)^2*HPL[{-1, 0, 0}, x])/3 + 
         (8*(1 + x)^2*HPL[{0, 0, 0}, x])/3 - (14*(1 + x)^2*Zeta[3])/3)) + 
     x*coeff[
       "delta[i,j]*Sum[k,QB]*(Cif[qi,bk]+Cif[bi,qk])*(vi^2+ai^2)*alphas^2/Pi^\
2"]*(CA*CF*TF*((Pi^2*(-2 + x)*(1 + x))/4 + ((-1 + x)*(27 + 13*x))/8 + 
         (Pi^2*(1 + x)^2*HPL[{-1}, x])/4 + ((-9 - 11*x)/4 + 
           (Pi^2*(-3 - 8*x - 3*x^2))/12)*HPL[{0}, x] - 
         (1 + x)^2*HPL[{-2, 0}, x] - (3*(1 + x)^2*HPL[{-1, 0}, x])/2 + 
         ((-6 + 8*x + 15*x^2)*HPL[{0, 0}, x])/4 + 
         (9*(-1 + x)*(1 + x)*HPL[{1, 0}, x])/4 + (-1 - 3*x - x^2)*
          HPL[{2, 0}, x] + 3*(1 + x)^2*HPL[{-1, -1, 0}, x] - 
         (5*(1 + x)^2*HPL[{-1, 0, 0}, x])/2 + 
         ((-1 - 4*x - x^2)*HPL[{0, 0, 0}, x])/2 - 2*(1 + 3*x + x^2)*Zeta[3] + 
         ep*Lv*((Pi^2*(-2 + x)*(1 + x))/2 + ((-1 + x)*(27 + 13*x))/4 + 
           (Pi^2*(1 + x)^2*HPL[{-1}, x])/2 + ((-9 - 11*x)/2 + 
             (Pi^2*(-3 - 8*x - 3*x^2))/6)*HPL[{0}, x] - 
           2*(1 + x)^2*HPL[{-2, 0}, x] - 3*(1 + x)^2*HPL[{-1, 0}, x] + 
           ((-6 + 8*x + 15*x^2)*HPL[{0, 0}, x])/2 + 
           (9*(-1 + x)*(1 + x)*HPL[{1, 0}, x])/2 - 2*(1 + 3*x + x^2)*
            HPL[{2, 0}, x] + 6*(1 + x)^2*HPL[{-1, -1, 0}, x] - 
           5*(1 + x)^2*HPL[{-1, 0, 0}, x] + (-1 - 4*x - x^2)*
            HPL[{0, 0, 0}, x] - 4*(1 + 3*x + x^2)*Zeta[3]) + 
         ep*(((-1 + x)*(121 + 43*x))/8 + (Pi^2*(-1 + 13*x + 3*x^2))/24 + 
           (Pi^4*(63 + 172*x + 63*x^2))/360 + (5*Pi^2*(1 + x)^2*HPL[{-2}, x])/
            4 + ((-3*Pi^2*(-1 + x)*(1 + x))/4 + ((-1 + x)*(27 + 13*x))/2)*
            HPL[{1}, x] + (-9 - 11*x + (Pi^2*(1 + 3*x + x^2))/3)*
            HPL[{2}, x] + (-6 + 8*x + 15*x^2)*HPL[{3}, x] - 
           2*(1 + 4*x + x^2)*HPL[{4}, x] - 3*(1 + x)^2*HPL[{-3, 0}, x] + 
           ((-7 - 16*x - 7*x^2)*HPL[{-2, 0}, x])/2 - 4*(1 + x)^2*
            HPL[{-2, 2}, x] - (13*Pi^2*(1 + x)^2*HPL[{-1, -1}, x])/4 + 
           ((31*Pi^2*(1 + x)^2)/12 + (-17 - 31*x - 17*x^2)/2)*
            HPL[{-1, 0}, x] - 6*(1 + x)^2*HPL[{-1, 2}, x] - 
           10*(1 + x)^2*HPL[{-1, 3}, x] + ((Pi^2*(-1 - x^2))/4 + 
             (-15 + 7*x + 27*x^2)/2)*HPL[{0, 0}, x] + 5*(-1 + x)*(1 + x)*
            HPL[{1, 0}, x] + 9*(-1 + x)*(1 + x)*HPL[{1, 2}, x] + 
           ((-1 + x)*(4 + 5*x)*HPL[{2, 0}, x])/2 - 4*(1 + 3*x + x^2)*
            HPL[{2, 2}, x] - (2 + x)*(1 + 2*x)*HPL[{3, 0}, x] + 
           7*(1 + x)^2*HPL[{-2, -1, 0}, x] - (17*(1 + x)^2*HPL[{-2, 0, 0}, 
              x])/2 + 7*(1 + x)^2*HPL[{-1, -2, 0}, x] + 
           (3*(5 + 12*x + 5*x^2)*HPL[{-1, -1, 0}, x])/2 + 
           12*(1 + x)^2*HPL[{-1, -1, 2}, x] + ((-41 - 92*x - 41*x^2)*
             HPL[{-1, 0, 0}, x])/4 + 2*(1 + x)^2*HPL[{-1, 2, 0}, x] + 
           ((-19 + 38*x + 62*x^2)*HPL[{0, 0, 0}, x])/4 + 
           (45*(-1 + x)*(1 + x)*HPL[{1, 0, 0}, x])/4 + 
           (9*(-1 + x)*(1 + x)*HPL[{1, 1, 0}, x])/2 - 5*(1 + 3*x + x^2)*
            HPL[{2, 0, 0}, x] - 2*(1 + 3*x + x^2)*HPL[{2, 1, 0}, x] - 
           15*(1 + x)^2*HPL[{-1, -1, -1, 0}, x] + 
           (41*(1 + x)^2*HPL[{-1, -1, 0, 0}, x])/2 - 9*(1 + x)^2*
            HPL[{-1, 0, 0, 0}, x] - 2*(1 + 4*x + x^2)*HPL[{0, 0, 0, 0}, x] + 
           ((-5 - 27*x - 20*x^2)*Zeta[3])/2 + HPL[{-1}, x]*
            ((Pi^2*(13 + 28*x + 13*x^2))/8 + (37*(1 + x)^2*Zeta[3])/2) + 
           HPL[{0}, x]*((Pi^2*(9 - 50*x - 57*x^2))/24 + (-22 + 8*x + 13*x^2)/
              2 + ((-9 - 16*x - 9*x^2)*Zeta[3])/2))) + 
       CF^2*TF*(-(Pi^2*(-2 + x)*(1 + x))/2 - ((-1 + x)*(27 + 13*x))/4 - 
         (Pi^2*(1 + x)^2*HPL[{-1}, x])/2 + 
         ((9 + 11*x)/2 + (Pi^2*(3 + 8*x + 3*x^2))/6)*HPL[{0}, x] + 
         2*(1 + x)^2*HPL[{-2, 0}, x] + 3*(1 + x)^2*HPL[{-1, 0}, x] + 
         ((6 - 8*x - 15*x^2)*HPL[{0, 0}, x])/2 - 
         (9*(-1 + x)*(1 + x)*HPL[{1, 0}, x])/2 + 2*(1 + 3*x + x^2)*
          HPL[{2, 0}, x] - 6*(1 + x)^2*HPL[{-1, -1, 0}, x] + 
         5*(1 + x)^2*HPL[{-1, 0, 0}, x] + (1 + 4*x + x^2)*HPL[{0, 0, 0}, x] + 
         4*(1 + 3*x + x^2)*Zeta[3] + ep*Lv*(-(Pi^2*(-2 + x)*(1 + x)) - 
           ((-1 + x)*(27 + 13*x))/2 - Pi^2*(1 + x)^2*HPL[{-1}, x] + 
           (9 + 11*x + (Pi^2*(3 + 8*x + 3*x^2))/3)*HPL[{0}, x] + 
           4*(1 + x)^2*HPL[{-2, 0}, x] + 6*(1 + x)^2*HPL[{-1, 0}, x] + 
           (6 - 8*x - 15*x^2)*HPL[{0, 0}, x] - 9*(-1 + x)*(1 + x)*
            HPL[{1, 0}, x] + 4*(1 + 3*x + x^2)*HPL[{2, 0}, x] - 
           12*(1 + x)^2*HPL[{-1, -1, 0}, x] + 10*(1 + x)^2*
            HPL[{-1, 0, 0}, x] + 2*(1 + 4*x + x^2)*HPL[{0, 0, 0}, x] + 
           8*(1 + 3*x + x^2)*Zeta[3]) + ep*(-((-1 + x)*(121 + 43*x))/4 + 
           (Pi^4*(-63 - 172*x - 63*x^2))/180 + (Pi^2*(1 - 13*x - 3*x^2))/12 - 
           (5*Pi^2*(1 + x)^2*HPL[{-2}, x])/2 + ((3*Pi^2*(-1 + x)*(1 + x))/2 - 
             (-1 + x)*(27 + 13*x))*HPL[{1}, x] + 
           (2*(9 + 11*x) - (2*Pi^2*(1 + 3*x + x^2))/3)*HPL[{2}, x] - 
           2*(-6 + 8*x + 15*x^2)*HPL[{3}, x] + 4*(1 + 4*x + x^2)*
            HPL[{4}, x] + 6*(1 + x)^2*HPL[{-3, 0}, x] + (7 + 16*x + 7*x^2)*
            HPL[{-2, 0}, x] + 8*(1 + x)^2*HPL[{-2, 2}, x] + 
           (13*Pi^2*(1 + x)^2*HPL[{-1, -1}, x])/2 + 
           (17 + 31*x + 17*x^2 - (31*Pi^2*(1 + x)^2)/6)*HPL[{-1, 0}, x] + 
           12*(1 + x)^2*HPL[{-1, 2}, x] + 20*(1 + x)^2*HPL[{-1, 3}, x] + 
           (15 - 7*x - 27*x^2 + (Pi^2*(1 + x^2))/2)*HPL[{0, 0}, x] - 
           10*(-1 + x)*(1 + x)*HPL[{1, 0}, x] - 18*(-1 + x)*(1 + x)*
            HPL[{1, 2}, x] - (-1 + x)*(4 + 5*x)*HPL[{2, 0}, x] + 
           8*(1 + 3*x + x^2)*HPL[{2, 2}, x] + 2*(2 + x)*(1 + 2*x)*
            HPL[{3, 0}, x] - 14*(1 + x)^2*HPL[{-2, -1, 0}, x] + 
           17*(1 + x)^2*HPL[{-2, 0, 0}, x] - 14*(1 + x)^2*HPL[{-1, -2, 0}, 
             x] - 3*(5 + 12*x + 5*x^2)*HPL[{-1, -1, 0}, x] - 
           24*(1 + x)^2*HPL[{-1, -1, 2}, x] + ((41 + 92*x + 41*x^2)*
             HPL[{-1, 0, 0}, x])/2 - 4*(1 + x)^2*HPL[{-1, 2, 0}, x] + 
           ((19 - 38*x - 62*x^2)*HPL[{0, 0, 0}, x])/2 - 
           (45*(-1 + x)*(1 + x)*HPL[{1, 0, 0}, x])/2 - 9*(-1 + x)*(1 + x)*
            HPL[{1, 1, 0}, x] + 10*(1 + 3*x + x^2)*HPL[{2, 0, 0}, x] + 
           4*(1 + 3*x + x^2)*HPL[{2, 1, 0}, x] + 30*(1 + x)^2*
            HPL[{-1, -1, -1, 0}, x] - 41*(1 + x)^2*HPL[{-1, -1, 0, 0}, x] + 
           18*(1 + x)^2*HPL[{-1, 0, 0, 0}, x] + 4*(1 + 4*x + x^2)*
            HPL[{0, 0, 0, 0}, x] + (5 + 27*x + 20*x^2)*Zeta[3] + 
           HPL[{-1}, x]*((Pi^2*(-13 - 28*x - 13*x^2))/4 - 37*(1 + x)^2*
              Zeta[3]) + HPL[{0}, x]*(22 - 8*x - 13*x^2 + 
             (Pi^2*(-9 + 50*x + 57*x^2))/12 + (9 + 16*x + 9*x^2)*
              Zeta[3])))) + x*coeff["Cii[qi,bj]*(vi^2+ai^2)*alphas^2/Pi^2"]*
      (CF*nl*(TF*((5*(-1 + x)*(1 + x))/6 + (Pi^2*(1 + x)^2)/18 + 
           ((-3 - 4*x - 3*x^2)*HPL[{0}, x])/6 + (2*(1 + x)^2*HPL[{-1, 0}, x])/
            3 - ((1 + x)^2*HPL[{0, 0}, x])/3 + 
           ep*Lv*((5*(-1 + x)*(1 + x))/3 + (Pi^2*(1 + x)^2)/9 + 
             ((-3 - 4*x - 3*x^2)*HPL[{0}, x])/3 + (4*(1 + x)^2*HPL[{-1, 0}, 
                x])/3 - (2*(1 + x)^2*HPL[{0, 0}, x])/3) + 
           ep*((79*(-1 + x)*(1 + x))/18 + (Pi^2*(20 + 31*x + 20*x^2))/54 - 
             (Pi^2*(1 + x)^2*HPL[{-1}, x])/2 + ((5*Pi^2*(1 + x)^2)/18 + 
               (-99 - 28*x + 51*x^2)/36)*HPL[{0}, x] + 
             (10*(-1 + x)*(1 + x)*HPL[{1}, x])/3 - (2*(3 + 4*x + 3*x^2)*HPL[
                {2}, x])/3 - (4*(1 + x)^2*HPL[{3}, x])/3 + 
             (2*(1 + x)^2*HPL[{-2, 0}, x])/3 + (2*(2 + 7*x + 2*x^2)*HPL[
                {-1, 0}, x])/9 + (8*(1 + x)^2*HPL[{-1, 2}, x])/3 + 
             ((-31 - 50*x - 31*x^2)*HPL[{0, 0}, x])/18 - 
             (2*(1 + x)^2*HPL[{-1, -1, 0}, x])/3 + (7*(1 + x)^2*HPL[
                {-1, 0, 0}, x])/3 - (4*(1 + x)^2*HPL[{0, 0, 0}, x])/3 + 
             (7*(1 + x)^2*Zeta[3])/3)) + TF^2*((92 - 103*x - 45*x^2)/27 - 
           (Pi^2*(1 - x - 2*x^2 + x^3))/(9*(-1 + x)) + 
           ((-12 + 9*x - 10*x^2 + 3*x^3)*HPL[{0}, x])/(3*(-1 + x)) - 
           (4*(-1 + 11*x)*HPL[{1}, x])/9 - (8*(1 + x^2)*HPL[{2}, x])/
            (3*(-1 + x)) - (4*(1 + x)^2*HPL[{-1, 0}, x])/3 + 
           ((-9 - 2*x - 3*x^2 + 2*x^3)*HPL[{0, 0}, x])/(3*(-1 + x)) - 
           (2*(4 + 3*x^2)*HPL[{1, 0}, x])/(3*(-1 + x)) - 
           (8*(1 + x)*HPL[{1, 1}, x])/3 + ep*Lv^3*((-1 - x)/3 + 
             PlusDistribution[-1, 1 - x]/2 + (2*PlusDistribution[0, 1 - x])/
              3) + (56/27 - (4*Pi^2)/9)*PlusDistribution[0, 1 - x] + 
           Lv^2*((-1 - x)/3 + PlusDistribution[-1, 1 - x]/2 + 
             (2*PlusDistribution[0, 1 - x])/3 + ep*((5 - 19*x)/9 - 
               (5*(1 + x^2)*HPL[{0}, x])/(3*(-1 + x)) - 2*(1 + x)*
                HPL[{1}, x] + (11/3 - Pi^2/9)*PlusDistribution[-1, 1 - x] + 
               (14*PlusDistribution[0, 1 - x])/9 - 4*PlusDistribution[1, 
                 1 - x])) - (40*PlusDistribution[1, 1 - x])/9 + 
           (8*PlusDistribution[2, 1 - x])/3 + PlusDistribution[-1, 1 - x]*
            (127/24 - (14*Pi^2)/27 + 2*Zeta[3]) + 
           ep*((1519 - 968*x - 711*x^2)/81 - (Pi^2*(-47 + 68*x - 65*x^2 + 
                40*x^3))/(54*(-1 + x)) + Pi^2*(1 + x)^2*HPL[{-1}, x] + 
             (-(Pi^2*(-16 - 5*x - 3*x^2 + 5*x^3))/(9*(-1 + x)) - 
               (231 - 339*x + 109*x^2 + 51*x^3)/(18*(-1 + x)))*HPL[{0}, x] + 
             ((Pi^2*(-9 + 11*x^2))/(9*(-1 + x)) - (4*(-107 + 88*x + 45*x^2))/
                27)*HPL[{1}, x] + (4*(-2 + x)*(5 - 2*x + 3*x^2)*HPL[{2}, x])/
              (3*(-1 + x)) + (8*(-5 - x - 2*x^2 + x^3)*HPL[{3}, x])/
              (3*(-1 + x)) - (4*(1 + x)^2*HPL[{-2, 0}, x])/3 - 
             (4*(2 + 7*x + 2*x^2)*HPL[{-1, 0}, x])/9 - 
             (16*(1 + x)^2*HPL[{-1, 2}, x])/3 + ((-79 + 41*x - 37*x^2 + 
                31*x^3)*HPL[{0, 0}, x])/(9*(-1 + x)) - 
             (2*(48 - 57*x + 43*x^2)*HPL[{1, 0}, x])/(9*(-1 + x)) - 
             (8*(-5 + 19*x)*HPL[{1, 1}, x])/9 - (8*(5 + 4*x^2)*HPL[{1, 2}, 
                x])/(3*(-1 + x)) - (2*(14 + 13*x^2)*HPL[{2, 0}, x])/
              (3*(-1 + x)) - (40*(1 + x^2)*HPL[{2, 1}, x])/(3*(-1 + x)) + 
             (4*(1 + x)^2*HPL[{-1, -1, 0}, x])/3 - (14*(1 + x)^2*HPL[
                {-1, 0, 0}, x])/3 + ((-29 - 8*x - 5*x^2 + 8*x^3)*HPL[
                {0, 0, 0}, x])/(3*(-1 + x)) - (2*(14 + 9*x^2)*HPL[{1, 0, 0}, 
                x])/(3*(-1 + x)) - (4*(9 + 8*x^2)*HPL[{1, 1, 0}, x])/
              (3*(-1 + x)) - 16*(1 + x)*HPL[{1, 1, 1}, x] + 
             (-104/27 + (22*Pi^2)/9)*PlusDistribution[1, 1 - x] + 
             (56*PlusDistribution[2, 1 - x])/9 - (16*PlusDistribution[3, 
                1 - x])/3 + PlusDistribution[0, 1 - x]*(160/81 - (23*Pi^2)/
                27 - (32*Zeta[3])/3) + PlusDistribution[-1, 1 - x]*
              (991/48 - (2147*Pi^2)/648 - (2*Zeta[3])/3) - 
             (2*(1 - 7*x + 7*x^3)*Zeta[3])/(3*(-1 + x))) + 
           Lv*((-2*(-1 + 11*x))/9 - (4*(1 + x^2)*HPL[{0}, x])/(3*(-1 + x)) - 
             (4*(1 + x)*HPL[{1}, x])/3 + (17*PlusDistribution[-1, 1 - x])/6 + 
             (20*PlusDistribution[0, 1 - x])/9 - (8*PlusDistribution[1, 
                1 - x])/3 + ep*((-2*(-107 + 88*x + 45*x^2))/27 - 
               (Pi^2*(7 - 4*x - 11*x^2 + 4*x^3))/(18*(-1 + x)) + (2*(-2 + x)*
                 (5 - 2*x + 3*x^2)*HPL[{0}, x])/(3*(-1 + x)) - (4*(-5 + 19*x)*
                 HPL[{1}, x])/9 - (20*(1 + x^2)*HPL[{2}, x])/(3*(-1 + x)) - 
               (8*(1 + x)^2*HPL[{-1, 0}, x])/3 + (4*(-5 - x - 2*x^2 + x^3)*
                 HPL[{0, 0}, x])/(3*(-1 + x)) - (4*(5 + 4*x^2)*HPL[{1, 0}, 
                  x])/(3*(-1 + x)) - 8*(1 + x)*HPL[{1, 1}, x] + (52/27 - 
                 (11*Pi^2)/9)*PlusDistribution[0, 1 - x] - 
               (56*PlusDistribution[1, 1 - x])/9 + 8*PlusDistribution[2, 
                 1 - x] + PlusDistribution[-1, 1 - x]*(157/12 - (175*Pi^2)/
                  108 + 4*Zeta[3]))))) + 
       CF*(TF*(-((-1 + x)*(-116 + 265*x + 703*x^2))/(216*x) - 
           (Pi^2*(8 + 63*x + 24*x^2 + 16*x^3))/(36*x) + 
           (Pi^2*(2 + 2*x + x^2)*HPL[{-1}, x])/(2*x) + 
           ((Pi^2*(-9 + x))/6 + (-165 + 192*x - 20*x^2)/36)*HPL[{0}, x] + 
           ep*Lv^3*(((-1 + x)*(4 + 7*x + 4*x^2))/(12*x) + 
             ((-1 - x)*HPL[{0}, x])/2) + ((Pi^2*(4 - 6*x + 3*x^2))/(12*x) - 
             ((-1 + x)*(22 - 17*x + 22*x^2))/(9*x))*HPL[{1}, x] + 
           (3 + 6*x + 4*x^2)*HPL[{2}, x] - 4*(1 + x)*HPL[{3}, x] + 
           (-10 - x)*HPL[{-2, 0}, x] - 5*(1 + x)*HPL[{-1, 0}, x] + 
           ((15 + 153*x + 40*x^2)*HPL[{0, 0}, x])/12 + 
           Lv^2*(((-1 + x)*(4 + 7*x + 4*x^2))/(12*x) + ((-1 - x)*HPL[{0}, x])/
              2 + ep*((Pi^2*(1 + x))/2 - ((-1 + x)*(29 - 16*x + 11*x^2))/
                (18*x) + ((30 + 45*x + 28*x^2)*HPL[{0}, x])/12 + ((-1 + x)*
                 (4 + 7*x + 4*x^2)*HPL[{1}, x])/(2*x) - 3*(1 + x)*
                HPL[{2}, x] - (5*(1 + x)*HPL[{0, 0}, x])/2)) + 
           ((-16 - 36*x + 33*x^2 + 16*x^3)*HPL[{1, 0}, x])/(6*x) + 
           (2*(-1 + x)*(4 + 7*x + 4*x^2)*HPL[{1, 1}, x])/(3*x) + 
           ((-20 - 9*x)*HPL[{2, 0}, x])/2 - 4*(1 + x)*HPL[{2, 1}, x] + 
           (6*(2 + 2*x + x^2)*HPL[{-1, -1, 0}, x])/x - 
           (5*(2 + 2*x + x^2)*HPL[{-1, 0, 0}, x])/x + 
           ((-9 - 5*x)*HPL[{0, 0, 0}, x])/2 + ((4 - 6*x + 3*x^2)*
             HPL[{1, 0, 0}, x])/x + ((4 - 6*x + 3*x^2)*HPL[{1, 1, 0}, x])/
            (2*x) - 4*(3 + x)*Zeta[3] + Lv*((Pi^2*(1 + x))/3 - 
             ((-1 + x)*(22 - 17*x + 22*x^2))/(18*x) + 
             ((3 + 6*x + 4*x^2)*HPL[{0}, x])/2 + ((-1 + x)*(4 + 7*x + 4*x^2)*
               HPL[{1}, x])/(3*x) - 2*(1 + x)*HPL[{2}, x] - 
             2*(1 + x)*HPL[{0, 0}, x] + ep*(-((-1 + x)*(-164 + 1135*x + 
                   187*x^2))/(216*x) - (Pi^2*(36 + 303*x + 129*x^2 + 76*x^3))/
                (72*x) + (Pi^2*(2 + 2*x + x^2)*HPL[{-1}, x])/x + 
               ((Pi^2*(-33 + 7*x))/12 + (-291 + 282*x - 44*x^2)/36)*
                HPL[{0}, x] + ((Pi^2*(4 - 6*x + 3*x^2))/(6*x) - 
                 (2*(-1 + x)*(29 - 16*x + 11*x^2))/(9*x))*HPL[{1}, x] + 
               ((30 + 45*x + 28*x^2)*HPL[{2}, x])/3 - 10*(1 + x)*
                HPL[{3}, x] - 2*(10 + x)*HPL[{-2, 0}, x] - 10*(1 + x)*
                HPL[{-1, 0}, x] + ((33 + 168*x + 44*x^2)*HPL[{0, 0}, x])/6 + 
               ((-20 - 39*x + 36*x^2 + 20*x^3)*HPL[{1, 0}, x])/(3*x) + 
               (2*(-1 + x)*(4 + 7*x + 4*x^2)*HPL[{1, 1}, x])/x - 11*(2 + x)*
                HPL[{2, 0}, x] - 12*(1 + x)*HPL[{2, 1}, x] + 
               (12*(2 + 2*x + x^2)*HPL[{-1, -1, 0}, x])/x - (10*(2 + 2*x + 
                  x^2)*HPL[{-1, 0, 0}, x])/x - 2*(5 + 3*x)*HPL[{0, 0, 0}, 
                 x] + (2*(4 - 6*x + 3*x^2)*HPL[{1, 0, 0}, x])/x + 
               ((4 - 6*x + 3*x^2)*HPL[{1, 1, 0}, x])/x - 2*(11 + 3*x)*
                Zeta[3])) + ep*((Pi^4*(209 + 17*x))/180 - 
             ((-1 + x)*(-2236 + 16697*x + 4901*x^2))/(1296*x) + 
             (Pi^2*(166 - 339*x - 489*x^2 + 30*x^3))/(216*x) + 
             (Pi^2*(42 + 5*x)*HPL[{-2}, x])/4 + ((Pi^2*(25 + 14*x))/6 + 
               (-291 + 282*x - 44*x^2)/18)*HPL[{2}, x] + 
             ((33 + 168*x + 44*x^2)*HPL[{3}, x])/3 - 4*(5 + 3*x)*
              HPL[{4}, x] + (-22 - 3*x)*HPL[{-3, 0}, x] + 
             ((-22 - 21*x)*HPL[{-2, 0}, x])/2 - 4*(10 + x)*HPL[{-2, 2}, x] - 
             (13*Pi^2*(2 + 2*x + x^2)*HPL[{-1, -1}, x])/(2*x) + 
             ((-40 - 41*x)/2 + (31*Pi^2*(2 + 2*x + x^2))/(6*x))*
              HPL[{-1, 0}, x] - 20*(1 + x)*HPL[{-1, 2}, x] - 
             (20*(2 + 2*x + x^2)*HPL[{-1, 3}, x])/x + 
             ((2*Pi^2*(-2 + x))/3 + (-510 + 1491*x - 80*x^2)/36)*
              HPL[{0, 0}, x] + ((-5*Pi^2*(4 - 6*x + 3*x^2))/(12*x) - 
               (-376 + 882*x - 621*x^2 + 124*x^3)/(36*x))*HPL[{1, 0}, x] - 
             (4*(-1 + x)*(29 - 16*x + 11*x^2)*HPL[{1, 1}, x])/(9*x) + 
             (2*(-20 - 39*x + 36*x^2 + 20*x^3)*HPL[{1, 2}, x])/(3*x) + 
             (4*(4 - 6*x + 3*x^2)*HPL[{1, 3}, x])/x + 
             ((132 + 297*x + 160*x^2)*HPL[{2, 0}, x])/12 + 
             (2*(30 + 45*x + 28*x^2)*HPL[{2, 1}, x])/3 - 22*(2 + x)*
              HPL[{2, 2}, x] + ((-52 - 31*x)*HPL[{3, 0}, x])/2 - 
             20*(1 + x)*HPL[{3, 1}, x] + (46 + 7*x)*HPL[{-2, -1, 0}, x] + 
             ((-130 - 17*x)*HPL[{-2, 0, 0}, x])/2 + (14*(2 + 2*x + x^2)*HPL[
                {-1, -2, 0}, x])/x + (23 + 17*x)*HPL[{-1, -1, 0}, x] + 
             (24*(2 + 2*x + x^2)*HPL[{-1, -1, 2}, x])/x - 
             (5*(13 + 11*x)*HPL[{-1, 0, 0}, x])/2 + (4*(2 + 2*x + x^2)*HPL[
                {-1, 2, 0}, x])/x + (5*(7 + 39*x + 8*x^2)*HPL[{0, 0, 0}, x])/
              4 + ((-28 - 81*x + 69*x^2 + 28*x^3)*HPL[{1, 0, 0}, x])/(3*x) + 
             ((-48 - 68*x + 61*x^2 + 48*x^3)*HPL[{1, 1, 0}, x])/(4*x) + 
             (4*(-1 + x)*(4 + 7*x + 4*x^2)*HPL[{1, 1, 1}, x])/x + 
             (2*(4 - 6*x + 3*x^2)*HPL[{1, 1, 2}, x])/x + 
             (3*(4 - 6*x + 3*x^2)*HPL[{1, 2, 0}, x])/(2*x) - 
             (11*(8 + 3*x)*HPL[{2, 0, 0}, x])/2 + (-30 - 19*x)*
              HPL[{2, 1, 0}, x] - 24*(1 + x)*HPL[{2, 1, 1}, x] - 
             (30*(2 + 2*x + x^2)*HPL[{-1, -1, -1, 0}, x])/x + 
             (41*(2 + 2*x + x^2)*HPL[{-1, -1, 0, 0}, x])/x - 
             (18*(2 + 2*x + x^2)*HPL[{-1, 0, 0, 0}, x])/x + 
             ((-29 - 15*x)*HPL[{0, 0, 0, 0}, x])/2 + (11*(4 - 6*x + 3*x^2)*
               HPL[{1, 0, 0, 0}, x])/(2*x) + (9*(4 - 6*x + 3*x^2)*HPL[
                {1, 1, 0, 0}, x])/(2*x) + (2*(4 - 6*x + 3*x^2)*HPL[
                {1, 1, 1, 0}, x])/x + ((-123 - 243*x - 40*x^2)*Zeta[3])/6 + 
             HPL[{0}, x]*((Pi^2*(-177 - 600*x - 116*x^2))/72 + (858 - 606*x + 
                 173*x^2)/216 + (-33 - x)*Zeta[3]) + HPL[{-1}, x]*
              ((Pi^2*(21 + 19*x))/4 + (37*(2 + 2*x + x^2)*Zeta[3])/x) + 
             HPL[{1}, x]*(-((-1 + x)*(-164 + 1135*x + 187*x^2))/(108*x) - 
               (Pi^2*(-104 - 174*x + 171*x^2 + 104*x^3))/(72*x) - 
               (7*(4 - 6*x + 3*x^2)*Zeta[3])/(2*x)))) + 
         TF^2*(((-1 + x)*(-116 + 265*x + 703*x^2))/(108*x) + 
           (Pi^2*(8 + 63*x + 24*x^2 + 16*x^3))/(18*x) - 
           (Pi^2*(2 + 2*x + x^2)*HPL[{-1}, x])/x + 
           ((Pi^2*(9 - x))/3 + (165 - 192*x + 20*x^2)/18)*HPL[{0}, x] + 
           ep*Lv^3*(-((-1 + x)*(4 + 7*x + 4*x^2))/(6*x) + 
             (1 + x)*HPL[{0}, x]) + (-(Pi^2*(4 - 6*x + 3*x^2))/(6*x) + 
             (2*(-1 + x)*(22 - 17*x + 22*x^2))/(9*x))*HPL[{1}, x] - 
           2*(3 + 6*x + 4*x^2)*HPL[{2}, x] + 8*(1 + x)*HPL[{3}, x] + 
           2*(10 + x)*HPL[{-2, 0}, x] + 10*(1 + x)*HPL[{-1, 0}, x] + 
           ((-15 - 153*x - 40*x^2)*HPL[{0, 0}, x])/6 + 
           Lv^2*(-((-1 + x)*(4 + 7*x + 4*x^2))/(6*x) + (1 + x)*HPL[{0}, x] + 
             ep*(Pi^2*(-1 - x) + ((-1 + x)*(29 - 16*x + 11*x^2))/(9*x) + 
               ((-30 - 45*x - 28*x^2)*HPL[{0}, x])/6 - ((-1 + x)*(4 + 7*x + 
                  4*x^2)*HPL[{1}, x])/x + 6*(1 + x)*HPL[{2}, x] + 5*(1 + x)*
                HPL[{0, 0}, x])) - ((-16 - 36*x + 33*x^2 + 16*x^3)*
             HPL[{1, 0}, x])/(3*x) - (4*(-1 + x)*(4 + 7*x + 4*x^2)*
             HPL[{1, 1}, x])/(3*x) + (20 + 9*x)*HPL[{2, 0}, x] + 
           8*(1 + x)*HPL[{2, 1}, x] - (12*(2 + 2*x + x^2)*HPL[{-1, -1, 0}, 
              x])/x + (10*(2 + 2*x + x^2)*HPL[{-1, 0, 0}, x])/x + 
           (9 + 5*x)*HPL[{0, 0, 0}, x] - (2*(4 - 6*x + 3*x^2)*
             HPL[{1, 0, 0}, x])/x - ((4 - 6*x + 3*x^2)*HPL[{1, 1, 0}, x])/x + 
           8*(3 + x)*Zeta[3] + Lv*((-2*Pi^2*(1 + x))/3 + 
             ((-1 + x)*(22 - 17*x + 22*x^2))/(9*x) + (-3 - 6*x - 4*x^2)*
              HPL[{0}, x] - (2*(-1 + x)*(4 + 7*x + 4*x^2)*HPL[{1}, x])/
              (3*x) + 4*(1 + x)*HPL[{2}, x] + 4*(1 + x)*HPL[{0, 0}, x] + 
             ep*(((-1 + x)*(-164 + 1135*x + 187*x^2))/(108*x) + 
               (Pi^2*(36 + 303*x + 129*x^2 + 76*x^3))/(36*x) - (2*Pi^2*
                 (2 + 2*x + x^2)*HPL[{-1}, x])/x + ((Pi^2*(33 - 7*x))/6 + 
                 (291 - 282*x + 44*x^2)/18)*HPL[{0}, x] + 
               (-(Pi^2*(4 - 6*x + 3*x^2))/(3*x) + (4*(-1 + x)*(29 - 16*x + 
                    11*x^2))/(9*x))*HPL[{1}, x] - (2*(30 + 45*x + 28*x^2)*
                 HPL[{2}, x])/3 + 20*(1 + x)*HPL[{3}, x] + 4*(10 + x)*
                HPL[{-2, 0}, x] + 20*(1 + x)*HPL[{-1, 0}, x] + 
               ((-33 - 168*x - 44*x^2)*HPL[{0, 0}, x])/3 - (2*(-20 - 39*x + 
                  36*x^2 + 20*x^3)*HPL[{1, 0}, x])/(3*x) - (4*(-1 + x)*
                 (4 + 7*x + 4*x^2)*HPL[{1, 1}, x])/x + 22*(2 + x)*
                HPL[{2, 0}, x] + 24*(1 + x)*HPL[{2, 1}, x] - 
               (24*(2 + 2*x + x^2)*HPL[{-1, -1, 0}, x])/x + (20*(2 + 2*x + 
                  x^2)*HPL[{-1, 0, 0}, x])/x + 4*(5 + 3*x)*HPL[{0, 0, 0}, 
                 x] - (4*(4 - 6*x + 3*x^2)*HPL[{1, 0, 0}, x])/x - 
               (2*(4 - 6*x + 3*x^2)*HPL[{1, 1, 0}, x])/x + 4*(11 + 3*x)*
                Zeta[3])) + ep*((Pi^4*(-209 - 17*x))/90 + 
             ((-1 + x)*(-2236 + 16697*x + 4901*x^2))/(648*x) - 
             (Pi^2*(166 - 339*x - 489*x^2 + 30*x^3))/(108*x) + 
             (Pi^2*(-42 - 5*x)*HPL[{-2}, x])/2 + ((Pi^2*(-25 - 14*x))/3 + 
               (291 - 282*x + 44*x^2)/9)*HPL[{2}, x] - 
             (2*(33 + 168*x + 44*x^2)*HPL[{3}, x])/3 + 8*(5 + 3*x)*
              HPL[{4}, x] + 2*(22 + 3*x)*HPL[{-3, 0}, x] + 
             (22 + 21*x)*HPL[{-2, 0}, x] + 8*(10 + x)*HPL[{-2, 2}, x] + 
             (13*Pi^2*(2 + 2*x + x^2)*HPL[{-1, -1}, x])/x + 
             (40 + 41*x - (31*Pi^2*(2 + 2*x + x^2))/(3*x))*HPL[{-1, 0}, x] + 
             40*(1 + x)*HPL[{-1, 2}, x] + (40*(2 + 2*x + x^2)*HPL[{-1, 3}, 
                x])/x + ((-4*Pi^2*(-2 + x))/3 + (510 - 1491*x + 80*x^2)/18)*
              HPL[{0, 0}, x] + ((5*Pi^2*(4 - 6*x + 3*x^2))/(6*x) + 
               (-376 + 882*x - 621*x^2 + 124*x^3)/(18*x))*HPL[{1, 0}, x] + 
             (8*(-1 + x)*(29 - 16*x + 11*x^2)*HPL[{1, 1}, x])/(9*x) - 
             (4*(-20 - 39*x + 36*x^2 + 20*x^3)*HPL[{1, 2}, x])/(3*x) - 
             (8*(4 - 6*x + 3*x^2)*HPL[{1, 3}, x])/x + 
             ((-132 - 297*x - 160*x^2)*HPL[{2, 0}, x])/6 - 
             (4*(30 + 45*x + 28*x^2)*HPL[{2, 1}, x])/3 + 44*(2 + x)*
              HPL[{2, 2}, x] + (52 + 31*x)*HPL[{3, 0}, x] + 
             40*(1 + x)*HPL[{3, 1}, x] - 2*(46 + 7*x)*HPL[{-2, -1, 0}, x] + 
             (130 + 17*x)*HPL[{-2, 0, 0}, x] - (28*(2 + 2*x + x^2)*HPL[
                {-1, -2, 0}, x])/x - 2*(23 + 17*x)*HPL[{-1, -1, 0}, x] - 
             (48*(2 + 2*x + x^2)*HPL[{-1, -1, 2}, x])/x + 5*(13 + 11*x)*
              HPL[{-1, 0, 0}, x] - (8*(2 + 2*x + x^2)*HPL[{-1, 2, 0}, x])/x - 
             (5*(7 + 39*x + 8*x^2)*HPL[{0, 0, 0}, x])/2 - 
             (2*(-28 - 81*x + 69*x^2 + 28*x^3)*HPL[{1, 0, 0}, x])/(3*x) - 
             ((-48 - 68*x + 61*x^2 + 48*x^3)*HPL[{1, 1, 0}, x])/(2*x) - 
             (8*(-1 + x)*(4 + 7*x + 4*x^2)*HPL[{1, 1, 1}, x])/x - 
             (4*(4 - 6*x + 3*x^2)*HPL[{1, 1, 2}, x])/x - 
             (3*(4 - 6*x + 3*x^2)*HPL[{1, 2, 0}, x])/x + 11*(8 + 3*x)*
              HPL[{2, 0, 0}, x] + 2*(30 + 19*x)*HPL[{2, 1, 0}, x] + 
             48*(1 + x)*HPL[{2, 1, 1}, x] + (60*(2 + 2*x + x^2)*HPL[
                {-1, -1, -1, 0}, x])/x - (82*(2 + 2*x + x^2)*HPL[
                {-1, -1, 0, 0}, x])/x + (36*(2 + 2*x + x^2)*HPL[{-1, 0, 0, 
                 0}, x])/x + (29 + 15*x)*HPL[{0, 0, 0, 0}, x] - 
             (11*(4 - 6*x + 3*x^2)*HPL[{1, 0, 0, 0}, x])/x - 
             (9*(4 - 6*x + 3*x^2)*HPL[{1, 1, 0, 0}, x])/x - 
             (4*(4 - 6*x + 3*x^2)*HPL[{1, 1, 1, 0}, x])/x + 
             ((123 + 243*x + 40*x^2)*Zeta[3])/3 + HPL[{0}, x]*
              ((-858 + 606*x - 173*x^2)/108 + (Pi^2*(177 + 600*x + 116*x^2))/
                36 + 2*(33 + x)*Zeta[3]) + HPL[{-1}, x]*((Pi^2*(-21 - 19*x))/
                2 - (74*(2 + 2*x + x^2)*Zeta[3])/x) + HPL[{1}, x]*
              (((-1 + x)*(-164 + 1135*x + 187*x^2))/(54*x) + (Pi^2*
                 (-104 - 174*x + 171*x^2 + 104*x^3))/(36*x) + 
               (7*(4 - 6*x + 3*x^2)*Zeta[3])/x))) + 
         CA*TF*((-373 + 548*x)/27 - (Pi^2*(-1 - 12*x + 24*x^2))/
            (36*(-1 + x)) + (-(Pi^2*(3 + 4*x^2))/(6*(-1 + x)) + 
             (81 - 88*x + 77*x^2)/(6*(-1 + x)))*HPL[{0}, x] + 
           ((-5*Pi^2*x^2)/(6*(-1 + x)) + (-106 + 365*x)/18)*HPL[{1}, x] + 
           (4*(7 + 4*x^2)*HPL[{2}, x])/(3*(-1 + x)) + 
           (2*(1 + x^2)*HPL[{3}, x])/(-1 + x) + 
           ((89 + 43*x^2)*HPL[{0, 0}, x])/(12*(-1 + x)) + 
           ((38 + 6*x + 33*x^2)*HPL[{1, 0}, x])/(6*(-1 + x)) + 
           (22*(1 + x)*HPL[{1, 1}, x])/3 + (2*HPL[{1, 2}, x])/(-1 + x) - 
           (x^2*HPL[{2, 0}, x])/(-1 + x) + ((5 + 3*x^2)*HPL[{0, 0, 0}, x])/
            (2*(-1 + x)) - (2*(-2 + x^2)*HPL[{1, 0, 0}, x])/(-1 + x) - 
           (3*x^2*HPL[{1, 1, 0}, x])/(-1 + x) + ep*Lv^3*((11*(1 + x))/12 - 
             (11*PlusDistribution[-1, 1 - x])/8 - (11*PlusDistribution[0, 
                1 - x])/6) + (134/9 - (2*Pi^2)/3)*PlusDistribution[1, 
             1 - x] - (22*PlusDistribution[2, 1 - x])/3 - 
           ((-3 + 5*x^2)*Zeta[3])/(2*(-1 + x)) + PlusDistribution[-1, 1 - x]*
            (-1535/96 + (37*Pi^2)/27 - Pi^4/120 + (7*Zeta[3])/2) + 
           PlusDistribution[0, 1 - x]*(-202/27 + (11*Pi^2)/9 + 7*Zeta[3]) + 
           Lv^2*((11*(1 + x))/12 - (11*PlusDistribution[-1, 1 - x])/8 - 
             (11*PlusDistribution[0, 1 - x])/6 + ep*((Pi^2*(-1 - x))/6 + 
               (-139 + 341*x)/36 + ((67 + 43*x^2)*HPL[{0}, x])/
                (12*(-1 + x)) + (11*(1 + x)*HPL[{1}, x])/2 + ((1 + x^2)*
                 HPL[{0, 0}, x])/(-1 + x) + (-101/18 + Pi^2/3)*
                PlusDistribution[0, 1 - x] + 11*PlusDistribution[1, 1 - x] + 
               PlusDistribution[-1, 1 - x]*(-31/3 + (11*Pi^2)/36 + 
                 3*Zeta[3]))) + Lv*((Pi^2*(-1 - x))/6 + (-53 + 187*x)/18 + 
             (2*(7 + 4*x^2)*HPL[{0}, x])/(3*(-1 + x)) + 
             (11*(1 + x)*HPL[{1}, x])/3 + ((1 + x^2)*HPL[{0, 0}, x])/
              (-1 + x) + (-67/9 + Pi^2/3)*PlusDistribution[0, 1 - x] + 
             (22*PlusDistribution[1, 1 - x])/3 + PlusDistribution[-1, 1 - x]*
              (-193/24 + 3*Zeta[3]) + ep*((-1531 + 1829*x)/54 - 
               (Pi^2*(-25 - 48*x + 117*x^2))/(72*(-1 + x)) + 
               (-(Pi^2*(3 + 4*x^2))/(3*(-1 + x)) + (67 - 88*x + 69*x^2)/
                  (3*(-1 + x)))*HPL[{0}, x] + ((-5*Pi^2*x^2)/(3*(-1 + x)) + 
                 (-139 + 332*x)/9)*HPL[{1}, x] + ((67 + 43*x^2)*HPL[{2}, x])/
                (3*(-1 + x)) + (4*(1 + x^2)*HPL[{3}, x])/(-1 + x) + 
               ((47 + 24*x^2)*HPL[{0, 0}, x])/(3*(-1 + x)) + 
               ((49 + 6*x + 44*x^2)*HPL[{1, 0}, x])/(3*(-1 + x)) + 22*(1 + x)*
                HPL[{1, 1}, x] + (4*HPL[{1, 2}, x])/(-1 + x) - (2*x^2*
                 HPL[{2, 0}, x])/(-1 + x) + ((5 + 3*x^2)*HPL[{0, 0, 0}, x])/
                (-1 + x) - (4*(-2 + x^2)*HPL[{1, 0, 0}, x])/(-1 + x) - 
               (6*x^2*HPL[{1, 1, 0}, x])/(-1 + x) + (202/9 - (4*Pi^2)/3)*
                PlusDistribution[1, 1 - x] - 22*PlusDistribution[2, 1 - x] - 
               ((-3 + 5*x^2)*Zeta[3])/(-1 + x) + PlusDistribution[-1, 1 - x]*
                (-1853/48 + (1877*Pi^2)/432 - Pi^4/60 + 4*Zeta[3]) + 
               PlusDistribution[0, 1 - x]*(-203/27 + (109*Pi^2)/36 + 
                 14*Zeta[3]))) + ep*((5*(-964 + 1037*x))/81 + 
             (Pi^4*(-25 + 18*x^2))/(90*(-1 + x)) - (Pi^2*(695 - 2088*x + 
                1569*x^2))/(216*(-1 + x)) + ((-4*Pi^2*(1 + x^2))/
                (3*(-1 + x)) + (2*(67 - 82*x + 66*x^2))/(3*(-1 + x)))*
              HPL[{2}, x] + (2*(47 + 24*x^2)*HPL[{3}, x])/(3*(-1 + x)) + 
             (2*(5 + 3*x^2)*HPL[{4}, x])/(-1 + x) + (-(Pi^2*(23 + 17*x^2))/
                (12*(-1 + x)) + (969 - 822*x + 691*x^2)/(36*(-1 + x)))*
              HPL[{0, 0}, x] + (-(Pi^2*(8 + 5*x^2))/(2*(-1 + x)) + 
               (453 - 570*x + 379*x^2)/(9*(-1 + x)))*HPL[{1, 0}, x] + 
             ((5*(-122 + 271*x))/18 - (Pi^2*(1 + 3*x^2))/(-1 + x))*
              HPL[{1, 1}, x] + (2*(46 + 9*x + 41*x^2)*HPL[{1, 2}, x])/
              (3*(-1 + x)) - (2*(-11 + x^2)*HPL[{1, 3}, x])/(-1 + x) + 
             ((178 + 12*x + 113*x^2)*HPL[{2, 0}, x])/(6*(-1 + x)) + 
             (2*(67 + 43*x^2)*HPL[{2, 1}, x])/(3*(-1 + x)) - 
             4*(1 + x)*HPL[{2, 2}, x] + ((4 + 3*x^2)*HPL[{3, 0}, x])/
              (-1 + x) + (8*(1 + x^2)*HPL[{3, 1}, x])/(-1 + x) + 
             ((225 + 101*x^2)*HPL[{0, 0, 0}, x])/(12*(-1 + x)) + 
             ((88 + 90*x + 63*x^2)*HPL[{1, 0, 0}, x])/(6*(-1 + x)) + 
             ((84 + 33*x + 79*x^2)*HPL[{1, 1, 0}, x])/(3*(-1 + x)) + 
             44*(1 + x)*HPL[{1, 1, 1}, x] - 8*(1 + x)*HPL[{1, 1, 2}, x] - 
             ((-4 + 7*x^2)*HPL[{1, 2, 0}, x])/(-1 + x) + 
             (2*(6 + x^2)*HPL[{1, 2, 1}, x])/(-1 + x) - 
             ((-8 + 9*x^2)*HPL[{2, 0, 0}, x])/(-1 + x) - 
             (8*x^2*HPL[{2, 1, 0}, x])/(-1 + x) + ((15 + 7*x^2)*HPL[
                {0, 0, 0, 0}, x])/(2*(-1 + x)) - ((-16 + 5*x^2)*HPL[
                {1, 0, 0, 0}, x])/(-1 + x) - ((-8 + 19*x^2)*HPL[{1, 1, 0, 0}, 
                x])/(-1 + x) - (6*(1 + 3*x^2)*HPL[{1, 1, 1, 0}, x])/
              (-1 + x) + (-202/9 + (4*Pi^2)/3)*PlusDistribution[2, 1 - x] + 
             (44*PlusDistribution[3, 1 - x])/3 + PlusDistribution[1, 1 - x]*
              (406/27 - (109*Pi^2)/18 - 28*Zeta[3]) - 
             ((-58 + 18*x + 23*x^2)*Zeta[3])/(6*(-1 + x)) + 
             PlusDistribution[0, 1 - x]*(-608/81 + (337*Pi^2)/108 - Pi^4/18 + 
               (67*Zeta[3])/3) + HPL[{0}, x]*(-(Pi^2*(130 - 24*x + 61*x^2))/
                (36*(-1 + x)) + (1455 - 3052*x + 1949*x^2)/(36*(-1 + x)) - 
               ((17 + 6*x^2)*Zeta[3])/(-1 + x)) + HPL[{1}, x]*
              ((-3413 + 3874*x)/54 - (Pi^2*(-81 - 54*x + 127*x^2))/
                (36*(-1 + x)) - ((-2 + 7*x^2)*Zeta[3])/(-1 + x)) + 
             PlusDistribution[-1, 1 - x]*(-11783/192 + (22747*Pi^2)/2592 + 
               (7*Pi^4)/40 + (535/12 - (59*Pi^2)/12)*Zeta[3] - (43*Zeta[5])/
                4)))) + CF^2*TF*((29 - 27*x)/2 - (Pi^2*(-3 - 17*x + 17*x^2))/
          (6*(-1 + x)) + (-(30 + 7*x + 3*x^2)/(2*(-1 + x)) - 
           (Pi^2*(-1 + 5*x^2))/(6*(-1 + x)))*HPL[{0}, x] + 
         ((-32 - 31*x)/2 - (Pi^2*(-9 + 5*x^2))/(6*(-1 + x)))*HPL[{1}, x] + 
         (2*(-1 - 11*x + 9*x^2)*HPL[{2}, x])/(-1 + x) + 
         (4*(2 + 5*x^2)*HPL[{3}, x])/(-1 + x) + 
         ((-15 - 16*x + 16*x^2)*HPL[{0, 0}, x])/(2*(-1 + x)) + 
         (3*(3 - 7*x + 3*x^2)*HPL[{1, 0}, x])/(-1 + x) + 
         16*(-1 + x)*HPL[{1, 1}, x] + (4*(7 + 8*x^2)*HPL[{1, 2}, x])/
          (-1 + x) + ((13 + 27*x^2)*HPL[{2, 0}, x])/(-1 + x) + 
         ((23 + 39*x^2)*HPL[{2, 1}, x])/(-1 + x) + 
         ((1 + 19*x^2)*HPL[{0, 0, 0}, x])/(2*(-1 + x)) + 
         ((7 + 19*x^2)*HPL[{1, 0, 0}, x])/(-1 + x) + 
         ((29 + 35*x^2)*HPL[{1, 1, 0}, x])/(-1 + x) + 
         48*(1 + x)*HPL[{1, 1, 1}, x] + (-32 - (8*Pi^2)/3)*
          PlusDistribution[1, 1 - x] + ep*Lv^3*(-5 - x + 
           ((1 + 3*x^2)*HPL[{0}, x])/(-1 + x) + 4*(1 + x)*HPL[{1}, x] + 
           (9/4 - (2*Pi^2)/3)*PlusDistribution[-1, 1 - x] + 
           6*PlusDistribution[0, 1 - x] + 8*PlusDistribution[1, 1 - x]) + 
         16*PlusDistribution[3, 1 - x] + 
         Lv^2*(-5 - x + ((1 + 3*x^2)*HPL[{0}, x])/(-1 + x) + 
           4*(1 + x)*HPL[{1}, x] + (9/4 - (2*Pi^2)/3)*PlusDistribution[-1, 
             1 - x] + 6*PlusDistribution[0, 1 - x] + 
           8*PlusDistribution[1, 1 - x] + ep*(-2*Pi^2*(1 + x) - 2*(4 + 5*x) + 
             ((-11 - 28*x + 7*x^2)*HPL[{0}, x])/(2*(-1 + x)) + 
             (-25 - x)*HPL[{1}, x] + (2*(7 + 13*x^2)*HPL[{2}, x])/(-1 + x) + 
             (2*(1 + 5*x^2)*HPL[{0, 0}, x])/(-1 + x) + 
             (20*(1 + x^2)*HPL[{1, 0}, x])/(-1 + x) + 36*(1 + x)*
              HPL[{1, 1}, x] + (18 + 2*Pi^2)*PlusDistribution[0, 1 - x] - 
             26*PlusDistribution[1, 1 - x] - 36*PlusDistribution[2, 1 - x] + 
             PlusDistribution[-1, 1 - x]*(123/8 - Pi^2/3 - 30*Zeta[3]))) + 
         PlusDistribution[-1, 1 - x]*(511/32 - (35*Pi^2)/24 + Pi^4/180 - 
           (15*Zeta[3])/2) - ((-31 + 11*x^2)*Zeta[3])/(-1 + x) + 
         32*PlusDistribution[0, 1 - x]*Zeta[3] + 
         Lv*(-7 - 9*x - (4*Pi^2*(1 + x))/3 + ((-4 - 10*x + 5*x^2)*
             HPL[{0}, x])/(-1 + x) + 2*(-7 + x)*HPL[{1}, x] + 
           (2*(5 + 9*x^2)*HPL[{2}, x])/(-1 + x) + 
           ((1 + 7*x^2)*HPL[{0, 0}, x])/(-1 + x) + 
           (14*(1 + x^2)*HPL[{1, 0}, x])/(-1 + x) + 24*(1 + x)*
            HPL[{1, 1}, x] + (16 + (4*Pi^2)/3)*PlusDistribution[0, 1 - x] - 
           12*PlusDistribution[1, 1 - x] - 24*PlusDistribution[2, 1 - x] + 
           PlusDistribution[-1, 1 - x]*(93/8 - Pi^2/2 - 22*Zeta[3]) + 
           ep*(-14*(-1 + 2*x) - (Pi^2*(19 - 46*x + 21*x^2))/(6*(-1 + x)) + 
             ((-20 - 13*x + 2*x^2)/(-1 + x) - (Pi^2*(5 + 23*x^2))/
                (6*(-1 + x)))*HPL[{0}, x] + (-17*(2 + x) - (Pi^2*
                 (-23 + 19*x^2))/(3*(-1 + x)))*HPL[{1}, x] + 
             (4*(-4 - 15*x + 6*x^2)*HPL[{2}, x])/(-1 + x) + 
             (8*(3 + 7*x^2)*HPL[{3}, x])/(-1 + x) + ((-17 - 24*x + 12*x^2)*
               HPL[{0, 0}, x])/(-1 + x) + (2*(3 - 29*x + 3*x^2)*HPL[{1, 0}, 
                x])/(-1 + x) + 12*(-7 + x)*HPL[{1, 1}, x] + 
             (8*(11 + 12*x^2)*HPL[{1, 2}, x])/(-1 + x) + 
             (2*(19 + 37*x^2)*HPL[{2, 0}, x])/(-1 + x) + 
             (2*(35 + 59*x^2)*HPL[{2, 1}, x])/(-1 + x) + 
             ((3 + 25*x^2)*HPL[{0, 0, 0}, x])/(-1 + x) + 
             (2*(13 + 25*x^2)*HPL[{1, 0, 0}, x])/(-1 + x) + 
             (6*(15 + 17*x^2)*HPL[{1, 1, 0}, x])/(-1 + x) + 
             160*(1 + x)*HPL[{1, 1, 1}, x] + (-52 - (44*Pi^2)/3)*
              PlusDistribution[1, 1 - x] + 36*PlusDistribution[2, 1 - x] + 
             (160*PlusDistribution[3, 1 - x])/3 - (2*(-45 + 17*x^2)*Zeta[3])/
              (-1 + x) + PlusDistribution[-1, 1 - x]*(709/16 - (121*Pi^2)/
                24 - Pi^4/90 + 7*Zeta[3]) + PlusDistribution[0, 1 - x]*
              (16 - (19*Pi^2)/3 + 96*Zeta[3]))) + 
         ep*((127 - 121*x)/2 - (Pi^4*(-33 + 5*x^2))/(30*(-1 + x)) - 
           (Pi^2*(-5 - 6*x + 5*x^2))/(4*(-1 + x)) + 
           ((-20 - 55*x + 29*x^2)/(-1 + x) - (Pi^2*(57 + 101*x^2))/
              (6*(-1 + x)))*HPL[{2}, x] + (2*(-14 - 30*x + 15*x^2)*
             HPL[{3}, x])/(-1 + x) + (22*(1 + 3*x^2)*HPL[{4}, x])/(-1 + x) + 
           ((-54 - 57*x + 46*x^2)/(2*(-1 + x)) - (Pi^2*(7 + 65*x^2))/
              (12*(-1 + x)))*HPL[{0, 0}, x] + 
           (-((16 + 17*x + 4*x^2)/(-1 + x)) - (Pi^2*(16 + 25*x^2))/
              (3*(-1 + x)))*HPL[{1, 0}, x] + ((-100 - 21*x)/2 - 
             (8*Pi^2*(-8 + 7*x^2))/(3*(-1 + x)))*HPL[{1, 1}, x] + 
           (4*(7 - 32*x + 7*x^2)*HPL[{1, 2}, x])/(-1 + x) + 
           (16*(4 + 7*x^2)*HPL[{1, 3}, x])/(-1 + x) + 
           ((-19 - 71*x + 32*x^2)*HPL[{2, 0}, x])/(-1 + x) + 
           ((-21 - 126*x + 61*x^2)*HPL[{2, 1}, x])/(-1 + x) + 
           (2*(41 + 81*x^2)*HPL[{2, 2}, x])/(-1 + x) + 
           (2*(17 + 39*x^2)*HPL[{3, 0}, x])/(-1 + x) + 
           (2*(31 + 63*x^2)*HPL[{3, 1}, x])/(-1 + x) + 
           ((-48 - 30*x + 13*x^2)*HPL[{0, 0, 0}, x])/(2*(-1 + x)) + 
           ((16 - 69*x + 12*x^2)*HPL[{1, 0, 0}, x])/(-1 + x) + 
           (2*(14 - 62*x + 13*x^2)*HPL[{1, 1, 0}, x])/(-1 + x) + 
           48*(-3 + x)*HPL[{1, 1, 1}, x] + (4*(45 + 53*x^2)*
             HPL[{1, 1, 2}, x])/(-1 + x) + ((109 + 131*x^2)*
             HPL[{1, 2, 0}, x])/(-1 + x) + (20*(9 + 10*x^2)*
             HPL[{1, 2, 1}, x])/(-1 + x) + ((19 + 89*x^2)*HPL[{2, 0, 0}, x])/
            (-1 + x) + (80*(1 + 2*x^2)*HPL[{2, 1, 0}, x])/(-1 + x) + 
           ((149 + 245*x^2)*HPL[{2, 1, 1}, x])/(-1 + x) + 
           (3*(1 + 19*x^2)*HPL[{0, 0, 0, 0}, x])/(2*(-1 + x)) + 
           (2*(8 + 29*x^2)*HPL[{1, 0, 0, 0}, x])/(-1 + x) + 
           (3*(23 + 41*x^2)*HPL[{1, 1, 0, 0}, x])/(-1 + x) + 
           (24*(8 + 9*x^2)*HPL[{1, 1, 1, 0}, x])/(-1 + x) + 
           320*(1 + x)*HPL[{1, 1, 1, 1}, x] + (32 + (58*Pi^2)/3)*
            PlusDistribution[2, 1 - x] - 16*PlusDistribution[3, 1 - x] - 
           (80*PlusDistribution[4, 1 - x])/3 + PlusDistribution[1, 1 - x]*
            (-32 + (29*Pi^2)/3 - 192*Zeta[3]) + PlusDistribution[0, 1 - x]*
            (-4*Pi^2 + (5*Pi^4)/9 - 32*Zeta[3]) - 
           (5*(5 - 9*x + 5*x^2)*Zeta[3])/(-1 + x) + HPL[{0}, x]*
            (-(Pi^2*(-34 - 96*x + 35*x^2))/(12*(-1 + x)) - 
             (83 - 93*x + 42*x^2)/(2*(-1 + x)) - ((-5 + 37*x^2)*Zeta[3])/
              (-1 + x)) + HPL[{1}, x]*((41 - 88*x)/2 - 
             (Pi^2*(73 - 92*x + 17*x^2))/(6*(-1 + x)) - 
             ((-169 + 61*x^2)*Zeta[3])/(-1 + x)) + 
           PlusDistribution[-1, 1 - x]*(3183/64 - (571*Pi^2)/96 - 
             (13*Pi^4)/180 + (-107/2 + (131*Pi^2)/6)*Zeta[3] - 116*Zeta[5]))))
 
rsigc["NNLO", "qg"] = 
    x*coeff["Sum[k,QB]*Cif[qi,qk]*(vi^2+ai^2)*alphas^2/Pi^2"]*
     (CA*TF*((Pi^2*(8 + 48*x + 12*x^2 + 151*x^3))/(72*x) + 
        (116 + 1617*x - 3678*x^2 + 1837*x^3)/(432*x) + 
        (Pi^2*(-1 - 2*x - 2*x^2)*HPL[{-1}, x])/4 + 
        ((Pi^2*x*(-3 + 2*x))/12 + (354 - 48*x - 457*x^2)/72)*HPL[{0}, x] + 
        ((Pi^2*(1 - 2*x + 2*x^2))/8 - (88 - 210*x + 75*x^2 + 74*x^3)/(72*x))*
         HPL[{1}, x] + ep*Lv^3*(-((-1 + x)*(4 + 7*x + 31*x^2))/(24*x) + 
          ((1 + 4*x)*HPL[{0}, x])/4 + ((-1 + 2*x - 2*x^2)*HPL[{1}, x])/4) + 
        ((-5 + 8*x - 65*x^2)*HPL[{2}, x])/4 + 
        ((3 + 14*x - 2*x^2)*HPL[{3}, x])/2 + 
        ((1 + 2*x + 2*x^2)*HPL[{-2, 0}, x])/2 + 
        ((1 + x)*(1 + 4*x)*HPL[{-1, 0}, x])/2 + (1 + 2*x + 2*x^2)*
         HPL[{-1, 2}, x] + ((-15 - 12*x - 346*x^2)*HPL[{0, 0}, x])/24 - 
        ((-16 - 18*x - 120*x^2 + 151*x^3)*HPL[{1, 0}, x])/(12*x) - 
        ((-1 + x)*(8 + 14*x + 77*x^2)*HPL[{1, 1}, x])/(6*x) + 
        Lv^2*(-((-1 + x)*(4 + 7*x + 31*x^2))/(24*x) + ((1 + 4*x)*HPL[{0}, x])/
           4 + ((-1 + 2*x - 2*x^2)*HPL[{1}, x])/4 + 
          ep*((Pi^2*(-3 - 26*x + 6*x^2))/24 - (116 - 207*x - 342*x^2 + 
              469*x^3)/(144*x) + ((-15 + 9*x - 104*x^2)*HPL[{0}, x])/12 - 
            ((-8 - 11*x - 54*x^2 + 71*x^3)*HPL[{1}, x])/(8*x) - 
            (3*(-1 - 10*x + 2*x^2)*HPL[{2}, x])/4 + 
            ((1 + 2*x + 2*x^2)*HPL[{-1, 0}, x])/2 + 
            ((5 + 16*x)*HPL[{0, 0}, x])/4 - (3*(1 - 2*x + 2*x^2)*
              HPL[{1, 0}, x])/2 - (5*(1 - 2*x + 2*x^2)*HPL[{1, 1}, x])/2)) + 
        ((-6 + 12*x - 11*x^2)*HPL[{1, 2}, x])/2 + 
        ((3 + 21*x - 4*x^2)*HPL[{2, 0}, x])/2 + 
        ((1 + 22*x - 6*x^2)*HPL[{2, 1}, x])/2 + (-1 - 2*x - 2*x^2)*
         HPL[{-1, -1, 0}, x] + (3*(1 + 2*x + 2*x^2)*HPL[{-1, 0, 0}, x])/2 + 
        ((9 + 20*x)*HPL[{0, 0, 0}, x])/4 + (-3 + 6*x - 5*x^2)*
         HPL[{1, 0, 0}, x] + ((-15 + 30*x - 28*x^2)*HPL[{1, 1, 0}, x])/4 - 
        (13*(1 - 2*x + 2*x^2)*HPL[{1, 1, 1}, x])/4 + 
        ((9 + 8*x + 8*x^2)*Zeta[3])/4 + Lv*((Pi^2*(-1 - 8*x + 2*x^2))/12 - 
          (44 - 87*x - 12*x^2 + 73*x^3)/(72*x) + 
          ((-3 + 2*x - 28*x^2)*HPL[{0}, x])/4 - ((-1 + x)*(8 + 17*x + 71*x^2)*
            HPL[{1}, x])/(12*x) + ((1 + 10*x - 2*x^2)*HPL[{2}, x])/2 + 
          ((1 + 2*x + 2*x^2)*HPL[{-1, 0}, x])/2 + (1 + 3*x)*HPL[{0, 0}, x] + 
          (-1 + 2*x - 2*x^2)*HPL[{1, 0}, x] - (3*(1 - 2*x + 2*x^2)*
            HPL[{1, 1}, x])/2 + ep*((Pi^2*(36 + 225*x + 12*x^2 + 751*x^3))/
             (144*x) - (-164 - 2697*x + 24*x^2 + 2945*x^3)/(432*x) + 
            (Pi^2*(-1 - 2*x - 2*x^2)*HPL[{-1}, x])/2 + 
            ((678 + 204*x - 1063*x^2)/72 + (Pi^2*(-1 - 28*x + 12*x^2))/24)*
             HPL[{0}, x] + ((17*Pi^2*(1 - 2*x + 2*x^2))/24 - (116 - 288*x - 
                297*x^2 + 496*x^3)/(36*x))*HPL[{1}, x] + 
            ((-45 + 66*x - 479*x^2)*HPL[{2}, x])/12 + 
            ((7 + 38*x - 6*x^2)*HPL[{3}, x])/2 + (1 + 2*x + 2*x^2)*
             HPL[{-2, 0}, x] + ((1 + 8*x + 6*x^2)*HPL[{-1, 0}, x])/2 + 
            2*(1 + 2*x + 2*x^2)*HPL[{-1, 2}, x] + ((-33 - 18*x - 386*x^2)*
              HPL[{0, 0}, x])/12 - ((-20 - 30*x - 150*x^2 + 191*x^3)*
              HPL[{1, 0}, x])/(6*x) - ((-8 - 11*x - 64*x^2 + 80*x^3)*
              HPL[{1, 1}, x])/(2*x) + (-9 + 18*x - 17*x^2)*HPL[{1, 2}, x] - 
            3*(-1 - 9*x + 2*x^2)*HPL[{2, 0}, x] + (1 + 34*x - 10*x^2)*
             HPL[{2, 1}, x] - 2*(1 + 2*x + 2*x^2)*HPL[{-1, -1, 0}, x] + 
            3*(1 + 2*x + 2*x^2)*HPL[{-1, 0, 0}, x] + (5 + 12*x)*
             HPL[{0, 0, 0}, x] + (-7 + 14*x - 12*x^2)*HPL[{1, 0, 0}, x] + 
            ((-21 + 42*x - 40*x^2)*HPL[{1, 1, 0}, x])/2 - 
            (25*(1 - 2*x + 2*x^2)*HPL[{1, 1, 1}, x])/2 + (4 - x + 5*x^2)*
             Zeta[3])) + ep*((Pi^4*(-40 - 77*x - 18*x^2))/360 - 
          (-2236 - 31611*x + 30324*x^2 + 4171*x^3)/(2592*x) + 
          (Pi^2*(-332 - 1293*x - 1668*x^2 + 5691*x^3))/(864*x) + 
          (Pi^2*(-1 - 2*x - 2*x^2)*HPL[{-2}, x])/2 + 
          ((759 + 312*x - 1369*x^2)/36 + (Pi^2*(-13 - 142*x + 34*x^2))/24)*
           HPL[{2}, x] + ((-51 - 18*x - 799*x^2)*HPL[{3}, x])/12 + 
          ((17 + 54*x - 6*x^2)*HPL[{4}, x])/2 + (1 + 2*x + 2*x^2)*
           HPL[{-3, 0}, x] + ((1 + 8*x + 6*x^2)*HPL[{-2, 0}, x])/2 + 
          2*(1 + 2*x + 2*x^2)*HPL[{-2, 2}, x] + Pi^2*(1 + 2*x + 2*x^2)*
           HPL[{-1, -1}, x] + ((1 + x)*(2 + 7*x) - (9*Pi^2*(1 + 2*x + 2*x^2))/
             8)*HPL[{-1, 0}, x] + (1 + 8*x + 6*x^2)*HPL[{-1, 2}, x] + 
          6*(1 + 2*x + 2*x^2)*HPL[{-1, 3}, x] + 
          ((933 + 102*x - 2719*x^2)/72 + (Pi^2*(-11 - 30*x + 6*x^2))/12)*
           HPL[{0, 0}, x] + ((Pi^2*(17 - 34*x + 32*x^2))/24 - 
            (188 - 297*x - 1008*x^2 + 1153*x^3)/(36*x))*HPL[{1, 0}, x] + 
          ((Pi^2*(43 - 86*x + 82*x^2))/24 - (464 - 1026*x - 1827*x^2 + 
              2470*x^3)/(72*x))*HPL[{1, 1}, x] - 
          ((-40 - 60*x - 333*x^2 + 409*x^3)*HPL[{1, 2}, x])/(6*x) + 
          ((-29 + 58*x - 50*x^2)*HPL[{1, 3}, x])/2 + 
          ((-87 + 174*x - 728*x^2)*HPL[{2, 0}, x])/12 + 
          ((-45 + 99*x - 515*x^2)*HPL[{2, 1}, x])/6 + (3 + 60*x - 17*x^2)*
           HPL[{2, 2}, x] + ((9 + 61*x - 12*x^2)*HPL[{3, 0}, x])/2 + 
          (5 + 42*x - 10*x^2)*HPL[{3, 1}, x] - 2*(1 + 2*x + 2*x^2)*
           HPL[{-2, -1, 0}, x] + 3*(1 + 2*x + 2*x^2)*HPL[{-2, 0, 0}, x] - 
          2*(1 + 2*x + 2*x^2)*HPL[{-1, -2, 0}, x] + (-1 - 8*x - 6*x^2)*
           HPL[{-1, -1, 0}, x] - 4*(1 + 2*x + 2*x^2)*HPL[{-1, -1, 2}, x] + 
          (3*(1 + 8*x + 6*x^2)*HPL[{-1, 0, 0}, x])/2 + 2*(1 + 2*x + 2*x^2)*
           HPL[{-1, 2, 0}, x] + 4*(1 + 2*x + 2*x^2)*HPL[{-1, 2, 1}, x] + 
          ((-35 - 36*x - 346*x^2)*HPL[{0, 0, 0}, x])/8 - 
          ((-56 - 96*x - 468*x^2 + 569*x^3)*HPL[{1, 0, 0}, x])/(12*x) - 
          ((-24 - 29*x - 234*x^2 + 270*x^3)*HPL[{1, 1, 0}, x])/(4*x) - 
          ((-32 - 35*x - 296*x^2 + 350*x^3)*HPL[{1, 1, 1}, x])/(4*x) + 
          ((-47 + 94*x - 88*x^2)*HPL[{1, 1, 2}, x])/2 - 
          (5*(13 - 26*x + 24*x^2)*HPL[{1, 2, 0}, x])/4 + 
          ((-41 + 82*x - 77*x^2)*HPL[{1, 2, 1}, x])/2 + 
          ((5 + 89*x - 24*x^2)*HPL[{2, 0, 0}, x])/2 + 
          ((-1 + 116*x - 40*x^2)*HPL[{2, 1, 0}, x])/2 + 
          ((-1 + 146*x - 50*x^2)*HPL[{2, 1, 1}, x])/2 + 4*(1 + 2*x + 2*x^2)*
           HPL[{-1, -1, -1, 0}, x] - 6*(1 + 2*x + 2*x^2)*HPL[{-1, -1, 0, 0}, 
            x] + (9*(1 + 2*x + 2*x^2)*HPL[{-1, 0, 0, 0}, x])/2 + 
          ((29 + 62*x)*HPL[{0, 0, 0, 0}, x])/4 + 
          ((-37 + 74*x - 60*x^2)*HPL[{1, 0, 0, 0}, x])/4 + 
          ((-77 + 154*x - 136*x^2)*HPL[{1, 1, 0, 0}, x])/4 + 
          (-25 + 50*x - 48*x^2)*HPL[{1, 1, 1, 0}, x] - 
          (111*(1 - 2*x + 2*x^2)*HPL[{1, 1, 1, 1}, x])/4 + 
          ((-3 + 69*x + 218*x^2)*Zeta[3])/6 + HPL[{-1}, x]*
           ((Pi^2*(-1 - 8*x - 6*x^2))/4 - (15*(1 + 2*x + 2*x^2)*Zeta[3])/2) + 
          HPL[{0}, x]*((4650 - 2040*x - 8273*x^2)/432 + 
            (Pi^2*(87 + 150*x + 1034*x^2))/144 + ((-3 - 19*x + 10*x^2)*
              Zeta[3])/2) + HPL[{1}, x]*((Pi^2*(-104 - 93*x - 882*x^2 + 1049*
                x^3))/(144*x) - (-164 - 2886*x + 294*x^2 + 2837*x^3)/
             (216*x) + ((21 - 42*x + 34*x^2)*Zeta[3])/4))) + 
      CF*TF*((-157 + 442*x - 305*x^2)/32 + (Pi^2*(1 - 20*x + 16*x^2))/48 + 
        ((-31 + 201*x - 174*x^2)/16 + (Pi^2*(7 - 14*x + 24*x^2))/24)*
         HPL[{0}, x] + ((-26 + 135*x - 88*x^2)/8 + (Pi^2*(1 - 2*x + 2*x^2))/
           3)*HPL[{1}, x] + ep*Lv^3*((3*(-1 + 4*x))/16 - 
          (3*(1 - 2*x + 4*x^2)*HPL[{0}, x])/8 - (3*(1 - 2*x + 2*x^2)*
            HPL[{1}, x])/4) + ((-7 + 48*x - 48*x^2)*HPL[{2}, x])/4 + 
        (-3 + 6*x - 10*x^2)*HPL[{3}, x] + (1 - 2*x + 2*x^2)*HPL[{-2, 0}, x] - 
        (1 + x)*(-1 + 3*x)*HPL[{-1, 0}, x] + 
        ((-11 + 60*x - 4*x^2)*HPL[{0, 0}, x])/16 + 
        ((-23 + 72*x - 56*x^2)*HPL[{1, 0}, x])/8 + 
        ((-23 + 80*x - 63*x^2)*HPL[{1, 1}, x])/4 + 
        Lv^2*((3*(-1 + 4*x))/16 - (3*(1 - 2*x + 4*x^2)*HPL[{0}, x])/8 - 
          (3*(1 - 2*x + 2*x^2)*HPL[{1}, x])/4 + 
          ep*((-19 + 61*x - 13*x^2)/16 + (Pi^2*(1 - 2*x + 11*x^2))/12 + 
            ((-3 + 52*x - 53*x^2)*HPL[{0}, x])/8 + ((-18 + 86*x - 53*x^2)*
              HPL[{1}, x])/8 + ((-7 + 14*x - 23*x^2)*HPL[{2}, x])/2 - 
            (11*(1 - 2*x + 4*x^2)*HPL[{0, 0}, x])/8 - 
            (7*(1 - 2*x + 2*x^2)*HPL[{1, 0}, x])/2 - 
            (13*(1 - 2*x + 2*x^2)*HPL[{1, 1}, x])/2)) + 
        ((-9 + 18*x - 20*x^2)*HPL[{1, 2}, x])/2 + 
        ((-11 + 22*x - 40*x^2)*HPL[{2, 0}, x])/4 - 
        (3*(7 - 14*x + 22*x^2)*HPL[{2, 1}, x])/4 + 
        ((-17 + 34*x - 52*x^2)*HPL[{0, 0, 0}, x])/8 + 
        ((1 - 2*x - 6*x^2)*HPL[{1, 0, 0}, x])/4 + 
        ((-7 + 14*x - 16*x^2)*HPL[{1, 1, 0}, x])/2 - 
        (35*(1 - 2*x + 2*x^2)*HPL[{1, 1, 1}, x])/4 + 
        ((13 - 26*x + 34*x^2)*Zeta[3])/4 + Lv*((-12 + 34*x - 11*x^2)/8 + 
          (Pi^2*(1 - 2*x + 8*x^2))/12 + ((-5 + 40*x - 46*x^2)*HPL[{0}, x])/
           8 + ((-8 + 34*x - 23*x^2)*HPL[{1}, x])/4 + 
          ((-5 + 10*x - 16*x^2)*HPL[{2}, x])/2 + (-1 + 2*x - 4*x^2)*
           HPL[{0, 0}, x] - (5*(1 - 2*x + 2*x^2)*HPL[{1, 0}, x])/2 - 
          (9*(1 - 2*x + 2*x^2)*HPL[{1, 1}, x])/2 + 
          ep*((-245 + 704*x - 415*x^2)/32 + (Pi^2*(-51 - 124*x + 84*x^2))/
             96 + ((-73 + 376*x - 394*x^2)/16 + (Pi^2*(41 - 82*x + 140*x^2))/
               48)*HPL[{0}, x] + ((-13 + 73*x - 51*x^2)/2 + 
              (11*Pi^2*(1 - 2*x + 2*x^2))/8)*HPL[{1}, x] + 
            ((-5 + 118*x - 109*x^2)*HPL[{2}, x])/4 + (-8 + 16*x - 27*x^2)*
             HPL[{3}, x] + 2*(1 - 2*x + 2*x^2)*HPL[{-2, 0}, x] - 
            2*(1 + x)*(-1 + 3*x)*HPL[{-1, 0}, x] - (5*(1 - 16*x + 2*x^2)*
              HPL[{0, 0}, x])/8 + ((-10 + 44*x - 31*x^2)*HPL[{1, 0}, x])/2 + 
            ((-21 + 98*x - 71*x^2)*HPL[{1, 1}, x])/2 - 2*(7 - 14*x + 15*x^2)*
             HPL[{1, 2}, x] - 4*(2 - 4*x + 7*x^2)*HPL[{2, 0}, x] + 
            ((-31 + 62*x - 98*x^2)*HPL[{2, 1}, x])/2 + (-5 + 10*x - 16*x^2)*
             HPL[{0, 0, 0}, x] + ((-3 + 6*x - 14*x^2)*HPL[{1, 0, 0}, x])/2 - 
            2*(6 - 12*x + 13*x^2)*HPL[{1, 1, 0}, x] - 
            (55*(1 - 2*x + 2*x^2)*HPL[{1, 1, 1}, x])/2 + 
            ((15 - 30*x + 44*x^2)*Zeta[3])/2)) + 
        ep*((-541 + 2624*x - 2123*x^2)/64 + (Pi^4*(-22 + 44*x - 47*x^2))/
           360 + (Pi^2*(24 - 152*x + 57*x^2))/96 + Pi^2*(-1 + 2*x - 2*x^2)*
           HPL[{-2}, x] + Pi^2*(1 + x)*(-1 + 3*x)*HPL[{-1}, x] + 
          ((-71 + 430*x - 394*x^2)/8 + (11*Pi^2*(5 - 10*x + 16*x^2))/24)*
           HPL[{2}, x] + ((-4 + 98*x - 31*x^2)*HPL[{3}, x])/4 + 
          ((-25 + 50*x - 74*x^2)*HPL[{4}, x])/2 + 2*(1 - 2*x + 2*x^2)*
           HPL[{-3, 0}, x] + (1 - 6*x - 6*x^2)*HPL[{-2, 0}, x] + 
          4*(1 - 2*x + 2*x^2)*HPL[{-2, 2}, x] - 2*(1 + x)*(-2 + 7*x)*
           HPL[{-1, 0}, x] - 4*(1 + x)*(-1 + 3*x)*HPL[{-1, 2}, x] + 
          ((Pi^2*(13 - 26*x + 36*x^2))/8 + (-44 + 273*x + 62*x^2)/16)*
           HPL[{0, 0}, x] + ((-7 + 215*x - 190*x^2)/8 + 
            (Pi^2*(41 - 82*x + 86*x^2))/24)*HPL[{1, 0}, x] + 
          ((-92 + 625*x - 504*x^2)/8 + (Pi^2*(83 - 166*x + 174*x^2))/24)*
           HPL[{1, 1}, x] + ((-19 + 86*x - 68*x^2)*HPL[{1, 2}, x])/2 + 
          (-7 + 14*x - 22*x^2)*HPL[{1, 3}, x] + ((13 + 208*x - 212*x^2)*
            HPL[{2, 0}, x])/8 + ((-6 + 228*x - 215*x^2)*HPL[{2, 1}, x])/4 + 
          ((-31 + 62*x - 114*x^2)*HPL[{2, 2}, x])/2 + 
          ((-35 + 70*x - 128*x^2)*HPL[{3, 0}, x])/4 + 
          ((-35 + 70*x - 114*x^2)*HPL[{3, 1}, x])/2 - 4*(1 - 2*x + 2*x^2)*
           HPL[{-2, -1, 0}, x] + 6*(1 - 2*x + 2*x^2)*HPL[{-2, 0, 0}, x] + 
          4*(1 + x)*(-1 + 3*x)*HPL[{-1, -1, 0}, x] - 6*(1 + x)*(-1 + 3*x)*
           HPL[{-1, 0, 0}, x] + ((11 + 152*x + 124*x^2)*HPL[{0, 0, 0}, x])/
           16 + ((-15 + 36*x - 40*x^2)*HPL[{1, 0, 0}, x])/8 + 
          (-10 + 39*x - 31*x^2)*HPL[{1, 1, 0}, x] + 
          ((-84 + 370*x - 287*x^2)*HPL[{1, 1, 1}, x])/4 + 
          ((-47 + 94*x - 106*x^2)*HPL[{1, 1, 2}, x])/2 + 
          ((-23 + 46*x - 56*x^2)*HPL[{1, 2, 0}, x])/2 - 
          (5*(11 - 22*x + 24*x^2)*HPL[{1, 2, 1}, x])/2 - 
          (7*(1 - 2*x + 12*x^2)*HPL[{2, 0, 0}, x])/4 + 
          ((-27 + 54*x - 106*x^2)*HPL[{2, 1, 0}, x])/2 - 
          (3*(41 - 82*x + 130*x^2)*HPL[{2, 1, 1}, x])/4 + 
          ((-53 + 106*x - 156*x^2)*HPL[{0, 0, 0, 0}, x])/8 + 
          ((9 - 18*x - 10*x^2)*HPL[{1, 0, 0, 0}, x])/4 - 
          9*x^2*HPL[{1, 1, 0, 0}, x] + ((-43 + 86*x - 94*x^2)*
            HPL[{1, 1, 1, 0}, x])/2 - (209*(1 - 2*x + 2*x^2)*
            HPL[{1, 1, 1, 1}, x])/4 + ((51 - 248*x + 4*x^2)*Zeta[3])/8 + 
          HPL[{1}, x]*((-187 + 1088*x - 859*x^2)/16 + 
            (Pi^2*(20 - 166*x + 117*x^2))/48 + 4*(2 - 4*x + 5*x^2)*Zeta[3]) + 
          HPL[{0}, x]*((-188 + 1336*x - 1553*x^2)/32 + 
            (Pi^2*(7 - 80*x + 10*x^2))/32 + ((43 - 86*x + 120*x^2)*Zeta[3])/
             4))))
 
rsigc["NNLO", "qp"] = 
    CF*TF*x*coeff[
       "Sum[k,QB]*(Cif[qi,qk]*(vi^2+ai^2)+Cif[qj,qk]*(vj^2+aj^2))*alphas^2/Pi\
^2"]*(((-1 + x)*(-116 - 1895*x + 703*x^2))/(432*x) + 
       (Pi^2*(8 + 33*x + 21*x^2 + 16*x^3))/(72*x) + 
       ((Pi^2*(-1 - x))/12 + (345 - 48*x + 20*x^2)/72)*HPL[{0}, x] + 
       ep*Lv^3*(-((-1 + x)*(4 + 7*x + 4*x^2))/(24*x) + ((1 + x)*HPL[{0}, x])/
          4) + ((-1 + x)*(22 - 17*x + 22*x^2)*HPL[{1}, x])/(18*x) + 
       ((-3 - 6*x - 4*x^2)*HPL[{2}, x])/2 + 2*(1 + x)*HPL[{3}, x] - 
       (5*(3 + 15*x + 8*x^2)*HPL[{0, 0}, x])/24 + 
       Lv^2*(-((-1 + x)*(4 + 7*x + 4*x^2))/(24*x) + ((1 + x)*HPL[{0}, x])/4 + 
         ep*((Pi^2*(-1 - x))/4 + ((-1 + x)*(29 - 16*x + 11*x^2))/(36*x) + 
           ((-30 - 45*x - 28*x^2)*HPL[{0}, x])/24 - 
           ((-1 + x)*(4 + 7*x + 4*x^2)*HPL[{1}, x])/(4*x) + 
           (3*(1 + x)*HPL[{2}, x])/2 + (5*(1 + x)*HPL[{0, 0}, x])/4)) - 
       ((-1 + x)*(16 + 37*x + 16*x^2)*HPL[{1, 0}, x])/(12*x) - 
       ((-1 + x)*(4 + 7*x + 4*x^2)*HPL[{1, 1}, x])/(3*x) + 
       (5*(1 + x)*HPL[{2, 0}, x])/2 + 2*(1 + x)*HPL[{2, 1}, x] + 
       (9*(1 + x)*HPL[{0, 0, 0}, x])/4 + (1 + x)*Zeta[3] + 
       ep*((-11*Pi^4*(1 + x))/120 + ((-1 + x)*(-2236 - 43567*x + 4901*x^2))/
          (2592*x) - (Pi^2*(166 + 1101*x + 105*x^2 + 30*x^3))/(432*x) + 
         ((Pi^2*(-1 + x)*(52 + 109*x + 52*x^2))/(72*x) + 
           ((-1 + x)*(-164 - 3185*x + 187*x^2))/(216*x))*HPL[{1}, x] + 
         ((-5*Pi^2*(1 + x))/4 + (651 + 6*x + 44*x^2)/36)*HPL[{2}, x] + 
         ((-33 - 90*x - 44*x^2)*HPL[{3}, x])/6 + 10*(1 + x)*HPL[{4}, x] + 
         ((-7*Pi^2*(1 + x))/6 + (906 + 21*x + 80*x^2)/72)*HPL[{0, 0}, x] + 
         ((-1 + x)*(188 - 73*x + 62*x^2)*HPL[{1, 0}, x])/(36*x) + 
         (2*(-1 + x)*(29 - 16*x + 11*x^2)*HPL[{1, 1}, x])/(9*x) - 
         (4*(-1 + x)*(5 + 11*x + 5*x^2)*HPL[{1, 2}, x])/(3*x) + 
         ((-99 - 153*x - 80*x^2)*HPL[{2, 0}, x])/12 + 
         ((-30 - 45*x - 28*x^2)*HPL[{2, 1}, x])/3 + 
         12*(1 + x)*HPL[{2, 2}, x] + (15*(1 + x)*HPL[{3, 0}, x])/2 + 
         10*(1 + x)*HPL[{3, 1}, x] + ((-35 - 99*x - 40*x^2)*
           HPL[{0, 0, 0}, x])/8 - ((-1 + x)*(56 + 143*x + 56*x^2)*
           HPL[{1, 0, 0}, x])/(12*x) - (6*(-1 + x)*(1 + x)^2*
           HPL[{1, 1, 0}, x])/x - (2*(-1 + x)*(4 + 7*x + 4*x^2)*
           HPL[{1, 1, 1}, x])/x + (19*(1 + x)*HPL[{2, 0, 0}, x])/2 + 
         10*(1 + x)*HPL[{2, 1, 0}, x] + 12*(1 + x)*HPL[{2, 1, 1}, x] + 
         (29*(1 + x)*HPL[{0, 0, 0, 0}, x])/4 + ((-3 + 45*x + 40*x^2)*Zeta[3])/
          12 + HPL[{0}, x]*((4974 - 4362*x - 173*x^2)/432 + 
           (Pi^2*(45 + 216*x + 116*x^2))/144 - (9*(1 + x)*Zeta[3])/2)) + 
       Lv*((Pi^2*(-1 - x))/6 + ((-1 + x)*(22 - 17*x + 22*x^2))/(36*x) + 
         ((-3 - 6*x - 4*x^2)*HPL[{0}, x])/4 - ((-1 + x)*(4 + 7*x + 4*x^2)*
           HPL[{1}, x])/(6*x) + (1 + x)*HPL[{2}, x] + 
         (1 + x)*HPL[{0, 0}, x] + ep*(((-1 + x)*(-164 - 3185*x + 187*x^2))/
            (432*x) + (Pi^2*(36 + 183*x + 117*x^2 + 76*x^3))/(144*x) + 
           ((-7*Pi^2*(1 + x))/24 + (651 + 6*x + 44*x^2)/72)*HPL[{0}, x] + 
           ((-1 + x)*(29 - 16*x + 11*x^2)*HPL[{1}, x])/(9*x) + 
           ((-30 - 45*x - 28*x^2)*HPL[{2}, x])/6 + 5*(1 + x)*HPL[{3}, x] + 
           ((-33 - 90*x - 44*x^2)*HPL[{0, 0}, x])/12 - 
           (2*(-1 + x)*(5 + 11*x + 5*x^2)*HPL[{1, 0}, x])/(3*x) - 
           ((-1 + x)*(4 + 7*x + 4*x^2)*HPL[{1, 1}, x])/x + 
           6*(1 + x)*HPL[{2, 0}, x] + 6*(1 + x)*HPL[{2, 1}, x] + 
           5*(1 + x)*HPL[{0, 0, 0}, x] + (1 + x)*Zeta[3]))) + 
     CF*TF*x*coeff["Cif[qi,qi]*vi*vj*(delta[i,q]*delta[j,b]+delta[i,b]*delta[\
j,q]-delta[i,q]*delta[j,q]-delta[i,b]*delta[j,b])*alphas^2/Pi^2"]*
      (10*(-1 + x) + (Pi^2*(10 + x))/12 - (Pi^2*(2 + 2*x + x^2)*HPL[{-1}, x])/
        (2*x) + (-5 + (5*Pi^2)/3 - 4*x)*HPL[{0}, x] - 
       (Pi^2*(4 - 6*x + 3*x^2)*HPL[{1}, x])/(12*x) + 
       (10 + x)*HPL[{-2, 0}, x] + 5*(1 + x)*HPL[{-1, 0}, x] - 
       (13*x*HPL[{0, 0}, x])/2 + ((5 - 4*x)*HPL[{1, 0}, x])/2 + 
       ((10 - x)*HPL[{2, 0}, x])/2 - (6*(2 + 2*x + x^2)*HPL[{-1, -1, 0}, x])/
        x + (5*(2 + 2*x + x^2)*HPL[{-1, 0, 0}, x])/x - 
       2*x*HPL[{0, 0, 0}, x] - ((4 - 6*x + 3*x^2)*HPL[{1, 0, 0}, x])/x - 
       ((4 - 6*x + 3*x^2)*HPL[{1, 1, 0}, x])/(2*x) + 2*(5 + x)*Zeta[3] + 
       ep*Lv*(20*(-1 + x) + (Pi^2*(10 + x))/6 - 
         (Pi^2*(2 + 2*x + x^2)*HPL[{-1}, x])/x + ((10*Pi^2)/3 - 2*(5 + 4*x))*
          HPL[{0}, x] - (Pi^2*(4 - 6*x + 3*x^2)*HPL[{1}, x])/(6*x) + 
         2*(10 + x)*HPL[{-2, 0}, x] + 10*(1 + x)*HPL[{-1, 0}, x] - 
         13*x*HPL[{0, 0}, x] + (5 - 4*x)*HPL[{1, 0}, x] + 
         (10 - x)*HPL[{2, 0}, x] - (12*(2 + 2*x + x^2)*HPL[{-1, -1, 0}, x])/
          x + (10*(2 + 2*x + x^2)*HPL[{-1, 0, 0}, x])/x - 
         4*x*HPL[{0, 0, 0}, x] - (2*(4 - 6*x + 3*x^2)*HPL[{1, 0, 0}, x])/x - 
         ((4 - 6*x + 3*x^2)*HPL[{1, 1, 0}, x])/x + 4*(5 + x)*Zeta[3]) + 
       ep*((4*Pi^4*(-11 + x))/45 + (93*(-1 + x))/2 + (Pi^2*(80 + 33*x))/12 + 
         (Pi^2*(-42 - 5*x)*HPL[{-2}, x])/4 + ((Pi^2*(-10 + x))/6 - 
           4*(5 + 4*x))*HPL[{2}, x] - 26*x*HPL[{3}, x] - 8*x*HPL[{4}, x] + 
         (22 + 3*x)*HPL[{-3, 0}, x] + ((22 + 21*x)*HPL[{-2, 0}, x])/2 + 
         4*(10 + x)*HPL[{-2, 2}, x] + (13*Pi^2*(2 + 2*x + x^2)*
           HPL[{-1, -1}, x])/(2*x) + ((40 + 41*x)/2 - 
           (31*Pi^2*(2 + 2*x + x^2))/(6*x))*HPL[{-1, 0}, x] + 
         20*(1 + x)*HPL[{-1, 2}, x] + (20*(2 + 2*x + x^2)*HPL[{-1, 3}, x])/
          x + (-11 - 42*x + (Pi^2*(11 + 5*x))/3)*HPL[{0, 0}, x] + 
         ((40 - 39*x)/4 + (5*Pi^2*(4 - 6*x + 3*x^2))/(12*x))*HPL[{1, 0}, x] - 
         2*(-5 + 4*x)*HPL[{1, 2}, x] - (4*(4 - 6*x + 3*x^2)*HPL[{1, 3}, x])/
          x + ((22 + 3*x)*HPL[{2, 0}, x])/4 - 2*(-10 + x)*HPL[{2, 2}, x] + 
         ((22 + x)*HPL[{3, 0}, x])/2 + (-46 - 7*x)*HPL[{-2, -1, 0}, x] + 
         ((130 + 17*x)*HPL[{-2, 0, 0}, x])/2 - 
         (14*(2 + 2*x + x^2)*HPL[{-1, -2, 0}, x])/x + 
         (-23 - 17*x)*HPL[{-1, -1, 0}, x] - (24*(2 + 2*x + x^2)*
           HPL[{-1, -1, 2}, x])/x + (5*(13 + 11*x)*HPL[{-1, 0, 0}, x])/2 - 
         (4*(2 + 2*x + x^2)*HPL[{-1, 2, 0}, x])/x - 24*x*HPL[{0, 0, 0}, x] + 
         ((25 - 17*x)*HPL[{1, 0, 0}, x])/2 + ((20 - 13*x)*HPL[{1, 1, 0}, x])/
          4 - (2*(4 - 6*x + 3*x^2)*HPL[{1, 1, 2}, x])/x - 
         (3*(4 - 6*x + 3*x^2)*HPL[{1, 2, 0}, x])/(2*x) - 
         (5*(-10 + x)*HPL[{2, 0, 0}, x])/2 + (10 - x)*HPL[{2, 1, 0}, x] + 
         (30*(2 + 2*x + x^2)*HPL[{-1, -1, -1, 0}, x])/x - 
         (41*(2 + 2*x + x^2)*HPL[{-1, -1, 0, 0}, x])/x + 
         (18*(2 + 2*x + x^2)*HPL[{-1, 0, 0, 0}, x])/x - 
         7*x*HPL[{0, 0, 0, 0}, x] - (11*(4 - 6*x + 3*x^2)*HPL[{1, 0, 0, 0}, 
            x])/(2*x) - (9*(4 - 6*x + 3*x^2)*HPL[{1, 1, 0, 0}, x])/(2*x) - 
         (2*(4 - 6*x + 3*x^2)*HPL[{1, 1, 1, 0}, x])/x + 
         3*(7 + 11*x)*Zeta[3] + HPL[{0}, x]*(-27 + 23*x + (Pi^2*(11 + 32*x))/
            6 + 2*(21 + 5*x)*Zeta[3]) + HPL[{-1}, x]*((Pi^2*(-21 - 19*x))/4 - 
           (37*(2 + 2*x + x^2)*Zeta[3])/x) + HPL[{1}, x]*
          (40*(-1 + x) + (Pi^2*(-20 + 19*x))/24 + 
           (7*(4 - 6*x + 3*x^2)*Zeta[3])/(2*x))))
 
rsigc["NNLO", "qq"] = 
    x*coeff["delta[i,j]*Sum[k,QB](Cif[qi,qk]*(vi^2+ai^2))*alphas^2/Pi^2"]*
      (CA*CF*((Pi^2*(-1 + x)^2)/8 + ((-1 + x)*(-15 + 13*x))/16 + 
         ((7 - 6*x)/8 + (Pi^2*(-1 + x)^2)/12)*HPL[{0}, x] + 
         (Pi^2*(-1 + x)^2*HPL[{1}, x])/12 + (3*(-1 + x)^2*HPL[{0, 0}, x])/4 + 
         (3*(-1 + x)^2*HPL[{1, 0}, x])/4 + ((-1 + x)^2*HPL[{2, 0}, x])/2 + 
         ((-1 + x)^2*HPL[{0, 0, 0}, x])/2 + ((-1 + x)^2*HPL[{1, 0, 0}, x])/
          2 + ((-1 + x)^2*HPL[{1, 1, 0}, x])/2 + 
         ep*Lv*((Pi^2*(-1 + x)^2)/4 + ((-1 + x)*(-15 + 13*x))/8 + 
           ((7 - 6*x)/4 + (Pi^2*(-1 + x)^2)/6)*HPL[{0}, x] + 
           (Pi^2*(-1 + x)^2*HPL[{1}, x])/6 + (3*(-1 + x)^2*HPL[{0, 0}, x])/
            2 + (3*(-1 + x)^2*HPL[{1, 0}, x])/2 + (-1 + x)^2*HPL[{2, 0}, x] + 
           (-1 + x)^2*HPL[{0, 0, 0}, x] + (-1 + x)^2*HPL[{1, 0, 0}, x] + 
           (-1 + x)^2*HPL[{1, 1, 0}, x]) + ep*((17*Pi^2*(-1 + x))/24 - 
           (53*Pi^4*(-1 + x)^2)/360 + ((-1 + x)*(53 + 93*x))/32 + 
           ((7 - 6*x)*HPL[{2}, x])/2 + 3*(-1 + x)^2*HPL[{3}, x] + 
           2*(-1 + x)^2*HPL[{4}, x] + ((15 - 8*x)/8 - (Pi^2*(-1 + x)^2)/4)*
            HPL[{0, 0}, x] + (-(Pi^2*(-1 + x)^2)/4 + (-3 + 5*x)/4)*
            HPL[{1, 0}, x] + 3*(-1 + x)^2*HPL[{1, 2}, x] + 
           2*(-1 + x)^2*HPL[{1, 3}, x] + ((-1 - 5*x + 2*x^2)*HPL[{2, 0}, x])/
            2 + 2*(-1 + x)^2*HPL[{2, 2}, x] + ((-1 + x)^2*HPL[{3, 0}, x])/2 + 
           ((7 - 28*x + 13*x^2)*HPL[{0, 0, 0}, x])/4 + 
           ((13 - 28*x + 13*x^2)*HPL[{1, 0, 0}, x])/4 + 
           ((-2 + x)*(-1 + 2*x)*HPL[{1, 1, 0}, x])/2 + 2*(-1 + x)^2*
            HPL[{1, 1, 2}, x] + ((-1 + x)^2*HPL[{1, 2, 0}, x])/2 + 
           (7*(-1 + x)^2*HPL[{2, 0, 0}, x])/2 + 2*(-1 + x)^2*
            HPL[{2, 1, 0}, x] + 2*(-1 + x)^2*HPL[{0, 0, 0, 0}, x] + 
           2*(-1 + x)^2*HPL[{1, 0, 0, 0}, x] + 
           (7*(-1 + x)^2*HPL[{1, 1, 0, 0}, x])/2 + 2*(-1 + x)^2*
            HPL[{1, 1, 1, 0}, x] + (3*(-3 + x)*(1 + x)*Zeta[3])/4 + 
           HPL[{1}, x]*(((-1 + x)*(-15 + 13*x))/4 + (Pi^2*(-4 + 7*x - 4*x^2))/
              12 - (5*(-1 + x)^2*Zeta[3])/2) + HPL[{0}, x]*
            ((Pi^2*(-7 + 7*x - 4*x^2))/12 + (19 - 136*x + 52*x^2)/16 - 
             (5*(-1 + x)^2*Zeta[3])/2))) + CF^2*(-(Pi^2*(-1 + x)^2)/4 - 
         ((-1 + x)*(-15 + 13*x))/8 + (-(Pi^2*(-1 + x)^2)/6 + (-7 + 6*x)/4)*
          HPL[{0}, x] - (Pi^2*(-1 + x)^2*HPL[{1}, x])/6 - 
         (3*(-1 + x)^2*HPL[{0, 0}, x])/2 - (3*(-1 + x)^2*HPL[{1, 0}, x])/2 - 
         (-1 + x)^2*HPL[{2, 0}, x] - (-1 + x)^2*HPL[{0, 0, 0}, x] - 
         (-1 + x)^2*HPL[{1, 0, 0}, x] - (-1 + x)^2*HPL[{1, 1, 0}, x] + 
         ep*Lv*(-(Pi^2*(-1 + x)^2)/2 - ((-1 + x)*(-15 + 13*x))/4 + 
           (-(Pi^2*(-1 + x)^2)/3 + (-7 + 6*x)/2)*HPL[{0}, x] - 
           (Pi^2*(-1 + x)^2*HPL[{1}, x])/3 - 3*(-1 + x)^2*HPL[{0, 0}, x] - 
           3*(-1 + x)^2*HPL[{1, 0}, x] - 2*(-1 + x)^2*HPL[{2, 0}, x] - 
           2*(-1 + x)^2*HPL[{0, 0, 0}, x] - 2*(-1 + x)^2*HPL[{1, 0, 0}, x] - 
           2*(-1 + x)^2*HPL[{1, 1, 0}, x]) + ep*((-17*Pi^2*(-1 + x))/12 + 
           (53*Pi^4*(-1 + x)^2)/180 - ((-1 + x)*(53 + 93*x))/16 + 
           (-7 + 6*x)*HPL[{2}, x] - 6*(-1 + x)^2*HPL[{3}, x] - 
           4*(-1 + x)^2*HPL[{4}, x] + ((Pi^2*(-1 + x)^2)/2 + (-15 + 8*x)/4)*
            HPL[{0, 0}, x] + ((3 - 5*x)/2 + (Pi^2*(-1 + x)^2)/2)*
            HPL[{1, 0}, x] - 6*(-1 + x)^2*HPL[{1, 2}, x] - 
           4*(-1 + x)^2*HPL[{1, 3}, x] + (1 + 5*x - 2*x^2)*HPL[{2, 0}, x] - 
           4*(-1 + x)^2*HPL[{2, 2}, x] - (-1 + x)^2*HPL[{3, 0}, x] + 
           ((-7 + 28*x - 13*x^2)*HPL[{0, 0, 0}, x])/2 + 
           ((-13 + 28*x - 13*x^2)*HPL[{1, 0, 0}, x])/2 - (-2 + x)*(-1 + 2*x)*
            HPL[{1, 1, 0}, x] - 4*(-1 + x)^2*HPL[{1, 1, 2}, x] - 
           (-1 + x)^2*HPL[{1, 2, 0}, x] - 7*(-1 + x)^2*HPL[{2, 0, 0}, x] - 
           4*(-1 + x)^2*HPL[{2, 1, 0}, x] - 4*(-1 + x)^2*HPL[{0, 0, 0, 0}, 
             x] - 4*(-1 + x)^2*HPL[{1, 0, 0, 0}, x] - 7*(-1 + x)^2*
            HPL[{1, 1, 0, 0}, x] - 4*(-1 + x)^2*HPL[{1, 1, 1, 0}, x] - 
           (3*(-3 + x)*(1 + x)*Zeta[3])/2 + HPL[{1}, x]*
            (-((-1 + x)*(-15 + 13*x))/2 + (Pi^2*(4 - 7*x + 4*x^2))/6 + 
             5*(-1 + x)^2*Zeta[3]) + HPL[{0}, x]*((-19 + 136*x - 52*x^2)/8 + 
             (Pi^2*(7 - 7*x + 4*x^2))/6 + 5*(-1 + x)^2*Zeta[3])))) + 
     x*coeff["(Cif[qi,qj]*(vi^2+ai^2)+Cif[qj,qi]*(vj^2+aj^2))*alphas^2/Pi^2"]*
      (CF^2*((17*(-1 + x))/8 + (Pi^2*(7 + 3*x))/24 - 
         (Pi^2*(4 + 3*x^2)*HPL[{-1}, x])/(6*(1 + x)) + 
         ((-9 + 7*x)/8 + (Pi^2*(5 + 4*x^2))/(12*(1 + x)))*HPL[{0}, x] + 
         2*(-1 + x)*HPL[{1}, x] + (-1 - x)*HPL[{2}, x] - 
         ((1 + x^2)*HPL[{3}, x])/(1 + x) + ((2 + x^2)*HPL[{-2, 0}, x])/
          (1 + x) + ((1 + x)*HPL[{-1, 0}, x])/2 + 
         (2*(1 + x^2)*HPL[{-1, 2}, x])/(1 + x) + ((-1 - 3*x)*HPL[{0, 0}, x])/
          2 + ep*Lv^2*(-1 + x + (Pi^2*(1 + x^2))/(12*(1 + x)) + 
           ((-1 - x)*HPL[{0}, x])/2 + ((1 + x^2)*HPL[{-1, 0}, x])/(1 + x) - 
           ((1 + x^2)*HPL[{0, 0}, x])/(2*(1 + x))) + ((1 - x)*HPL[{1, 0}, x])/
          2 + ((1 + x^2)*HPL[{2, 0}, x])/(2*(1 + x)) - 
         (2*(2 + x^2)*HPL[{-1, -1, 0}, x])/(1 + x) + 
         ((4 + 3*x^2)*HPL[{-1, 0, 0}, x])/(1 + x) - 
         ((5 + 3*x^2)*HPL[{0, 0, 0}, x])/(4*(1 + x)) + 
         ((15 + 11*x^2)*Zeta[3])/(4*(1 + x)) + 
         Lv*(-1 + x + (Pi^2*(1 + x^2))/(12*(1 + x)) + ((-1 - x)*HPL[{0}, x])/
            2 + ((1 + x^2)*HPL[{-1, 0}, x])/(1 + x) - 
           ((1 + x^2)*HPL[{0, 0}, x])/(2*(1 + x)) + 
           ep*((13*(-1 + x))/4 + (Pi^2*(3 + 5*x + x^2))/(6*(1 + x)) - 
             (Pi^2*(4 + 3*x^2)*HPL[{-1}, x])/(3*(1 + x)) + 
             ((-7 + 9*x)/4 + (Pi^2*(5 + 4*x^2))/(6*(1 + x)))*HPL[{0}, x] + 
             4*(-1 + x)*HPL[{1}, x] - 2*(1 + x)*HPL[{2}, x] - 
             (2*(1 + x^2)*HPL[{3}, x])/(1 + x) + (2*(2 + x^2)*HPL[{-2, 0}, 
                x])/(1 + x) + (2*x*HPL[{-1, 0}, x])/(1 + x) + 
             (4*(1 + x^2)*HPL[{-1, 2}, x])/(1 + x) - ((1 + 8*x + 5*x^2)*HPL[
                {0, 0}, x])/(2*(1 + x)) + (1 - x)*HPL[{1, 0}, x] + 
             ((1 + x^2)*HPL[{2, 0}, x])/(1 + x) - (4*(2 + x^2)*HPL[
                {-1, -1, 0}, x])/(1 + x) + (2*(4 + 3*x^2)*HPL[{-1, 0, 0}, x])/
              (1 + x) - ((5 + 3*x^2)*HPL[{0, 0, 0}, x])/(2*(1 + x)) + 
             ((15 + 11*x^2)*Zeta[3])/(2*(1 + x)))) + 
         ep*((Pi^2*(12 - 11*x))/24 + (19*(-1 + x))/4 - (Pi^4*(24 + 37*x^2))/
            (360*(1 + x)) - (Pi^2*(25 + 12*x^2)*HPL[{-2}, x])/(12*(1 + x)) + 
           ((13*(-1 + x))/2 + (Pi^2*(-1 + x))/6)*HPL[{1}, x] + 
           ((-7 + 9*x)/2 - (Pi^2*(1 + x^2))/(6*(1 + x)))*HPL[{2}, x] - 
           ((1 + 8*x + 5*x^2)*HPL[{3}, x])/(1 + x) - 
           ((5 + 3*x^2)*HPL[{4}, x])/(1 + x) + ((5 + 2*x^2)*HPL[{-3, 0}, x])/
            (1 + x) - ((1 - x + x^2)*HPL[{-2, 0}, x])/(1 + x) + 
           (4*(2 + x^2)*HPL[{-2, 2}, x])/(1 + x) + 
           (Pi^2*(13 + 6*x^2)*HPL[{-1, -1}, x])/(3*(1 + x)) + 
           ((2 + 3*x)/2 - (3*Pi^2*(5 + 3*x^2))/(4*(1 + x)))*HPL[{-1, 0}, x] + 
           (4*x*HPL[{-1, 2}, x])/(1 + x) + (4*(4 + 3*x^2)*HPL[{-1, 3}, x])/
            (1 + x) + ((-13 + x)/8 + (Pi^2*(35 + 17*x^2))/(24*(1 + x)))*
            HPL[{0, 0}, x] + (9*(-1 + x)*HPL[{1, 0}, x])/2 + 
           8*(-1 + x)*HPL[{1, 1}, x] - 2*(-1 + x)*HPL[{1, 2}, x] - 
           ((5 + 10*x + 7*x^2)*HPL[{2, 0}, x])/(2*(1 + x)) - 
           4*(1 + x)*HPL[{2, 1}, x] + (2*(1 + x^2)*HPL[{2, 2}, x])/(1 + x) - 
           ((1 + 3*x^2)*HPL[{3, 0}, x])/(2*(1 + x)) - 
           (4*(1 + x^2)*HPL[{3, 1}, x])/(1 + x) - 
           ((9 + 4*x^2)*HPL[{-2, -1, 0}, x])/(1 + x) + 
           ((23 + 12*x^2)*HPL[{-2, 0, 0}, x])/(2*(1 + x)) - 
           (2*(5 + 2*x^2)*HPL[{-1, -2, 0}, x])/(1 + x) + 
           (2*(1 - x + x^2)*HPL[{-1, -1, 0}, x])/(1 + x) - 
           (8*(2 + x^2)*HPL[{-1, -1, 2}, x])/(1 + x) - 
           ((1 - 5*x + x^2)*HPL[{-1, 0, 0}, x])/(1 + x) + 
           (2*(1 + 2*x^2)*HPL[{-1, 2, 0}, x])/(1 + x) + 
           (8*(1 + x^2)*HPL[{-1, 2, 1}, x])/(1 + x) - 
           ((-1 + 22*x + 15*x^2)*HPL[{0, 0, 0}, x])/(4*(1 + x)) - 
           (5*(-1 + x)*HPL[{1, 0, 0}, x])/2 + (1 - x)*HPL[{1, 1, 0}, x] + 
           (5*(1 + x^2)*HPL[{2, 0, 0}, x])/(2*(1 + x)) + 
           ((1 + x^2)*HPL[{2, 1, 0}, x])/(1 + x) + 
           (4*(5 + 2*x^2)*HPL[{-1, -1, -1, 0}, x])/(1 + x) - 
           (2*(13 + 6*x^2)*HPL[{-1, -1, 0, 0}, x])/(1 + x) + 
           ((13 + 9*x^2)*HPL[{-1, 0, 0, 0}, x])/(1 + x) - 
           ((15 + 7*x^2)*HPL[{0, 0, 0, 0}, x])/(4*(1 + x)) + 
           ((-1 + 13*x + x^2)*Zeta[3])/(2*(1 + x)) + HPL[{0}, x]*
            ((-23 + 37*x)/8 + (Pi^2*(-1 + 28*x + 11*x^2))/(24*(1 + x)) + 
             (3*(9 + 4*x^2)*Zeta[3])/(2*(1 + x))) + HPL[{-1}, x]*
            ((Pi^2*(1 - 5*x + x^2))/(6*(1 + x)) - ((28 + 15*x^2)*Zeta[3])/
              (1 + x)))) + CA*CF*((Pi^2*(-7 - 3*x))/48 - (17*(-1 + x))/16 + 
         (Pi^2*(4 + 3*x^2)*HPL[{-1}, x])/(12*(1 + x)) + 
         ((9 - 7*x)/16 - (Pi^2*(5 + 4*x^2))/(24*(1 + x)))*HPL[{0}, x] + 
         (1 - x)*HPL[{1}, x] + ((1 + x)*HPL[{2}, x])/2 + 
         ((1 + x^2)*HPL[{3}, x])/(2*(1 + x)) - ((2 + x^2)*HPL[{-2, 0}, x])/
          (2*(1 + x)) + ((-1 - x)*HPL[{-1, 0}, x])/4 - 
         ((1 + x^2)*HPL[{-1, 2}, x])/(1 + x) + ((1 + 3*x)*HPL[{0, 0}, x])/4 + 
         ep*Lv^2*((1 - x)/2 - (Pi^2*(1 + x^2))/(24*(1 + x)) + 
           ((1 + x)*HPL[{0}, x])/4 - ((1 + x^2)*HPL[{-1, 0}, x])/
            (2*(1 + x)) + ((1 + x^2)*HPL[{0, 0}, x])/(4*(1 + x))) + 
         ((-1 + x)*HPL[{1, 0}, x])/4 - ((1 + x^2)*HPL[{2, 0}, x])/
          (4*(1 + x)) + ((2 + x^2)*HPL[{-1, -1, 0}, x])/(1 + x) - 
         ((4 + 3*x^2)*HPL[{-1, 0, 0}, x])/(2*(1 + x)) + 
         ((5 + 3*x^2)*HPL[{0, 0, 0}, x])/(8*(1 + x)) - 
         ((15 + 11*x^2)*Zeta[3])/(8*(1 + x)) + 
         Lv*((1 - x)/2 - (Pi^2*(1 + x^2))/(24*(1 + x)) + 
           ((1 + x)*HPL[{0}, x])/4 - ((1 + x^2)*HPL[{-1, 0}, x])/
            (2*(1 + x)) + ((1 + x^2)*HPL[{0, 0}, x])/(4*(1 + x)) + 
           ep*((-13*(-1 + x))/8 - (Pi^2*(3 + 5*x + x^2))/(12*(1 + x)) + 
             (Pi^2*(4 + 3*x^2)*HPL[{-1}, x])/(6*(1 + x)) + 
             ((7 - 9*x)/8 - (Pi^2*(5 + 4*x^2))/(12*(1 + x)))*HPL[{0}, x] - 
             2*(-1 + x)*HPL[{1}, x] + (1 + x)*HPL[{2}, x] + 
             ((1 + x^2)*HPL[{3}, x])/(1 + x) - ((2 + x^2)*HPL[{-2, 0}, x])/
              (1 + x) - (x*HPL[{-1, 0}, x])/(1 + x) - 
             (2*(1 + x^2)*HPL[{-1, 2}, x])/(1 + x) + ((1 + 8*x + 5*x^2)*HPL[
                {0, 0}, x])/(4*(1 + x)) + ((-1 + x)*HPL[{1, 0}, x])/2 - 
             ((1 + x^2)*HPL[{2, 0}, x])/(2*(1 + x)) + 
             (2*(2 + x^2)*HPL[{-1, -1, 0}, x])/(1 + x) - 
             ((4 + 3*x^2)*HPL[{-1, 0, 0}, x])/(1 + x) + 
             ((5 + 3*x^2)*HPL[{0, 0, 0}, x])/(4*(1 + x)) - 
             ((15 + 11*x^2)*Zeta[3])/(4*(1 + x)))) + 
         ep*((-19*(-1 + x))/8 + (Pi^2*(-12 + 11*x))/48 + (Pi^4*(24 + 37*x^2))/
            (720*(1 + x)) + (Pi^2*(25 + 12*x^2)*HPL[{-2}, x])/(24*(1 + x)) + 
           ((Pi^2*(1 - x))/12 - (13*(-1 + x))/4)*HPL[{1}, x] + 
           ((7 - 9*x)/4 + (Pi^2*(1 + x^2))/(12*(1 + x)))*HPL[{2}, x] + 
           ((1 + 8*x + 5*x^2)*HPL[{3}, x])/(2*(1 + x)) + 
           ((5 + 3*x^2)*HPL[{4}, x])/(2*(1 + x)) - 
           ((5 + 2*x^2)*HPL[{-3, 0}, x])/(2*(1 + x)) + 
           ((1 - x + x^2)*HPL[{-2, 0}, x])/(2*(1 + x)) - 
           (2*(2 + x^2)*HPL[{-2, 2}, x])/(1 + x) - 
           (Pi^2*(13 + 6*x^2)*HPL[{-1, -1}, x])/(6*(1 + x)) + 
           ((-2 - 3*x)/4 + (3*Pi^2*(5 + 3*x^2))/(8*(1 + x)))*
            HPL[{-1, 0}, x] - (2*x*HPL[{-1, 2}, x])/(1 + x) - 
           (2*(4 + 3*x^2)*HPL[{-1, 3}, x])/(1 + x) + 
           ((13 - x)/16 - (Pi^2*(35 + 17*x^2))/(48*(1 + x)))*HPL[{0, 0}, x] - 
           (9*(-1 + x)*HPL[{1, 0}, x])/4 - 4*(-1 + x)*HPL[{1, 1}, x] + 
           (-1 + x)*HPL[{1, 2}, x] + ((5 + 10*x + 7*x^2)*HPL[{2, 0}, x])/
            (4*(1 + x)) + 2*(1 + x)*HPL[{2, 1}, x] - 
           ((1 + x^2)*HPL[{2, 2}, x])/(1 + x) + ((1 + 3*x^2)*HPL[{3, 0}, x])/
            (4*(1 + x)) + (2*(1 + x^2)*HPL[{3, 1}, x])/(1 + x) + 
           ((9 + 4*x^2)*HPL[{-2, -1, 0}, x])/(2*(1 + x)) - 
           ((23 + 12*x^2)*HPL[{-2, 0, 0}, x])/(4*(1 + x)) + 
           ((5 + 2*x^2)*HPL[{-1, -2, 0}, x])/(1 + x) - 
           ((1 - x + x^2)*HPL[{-1, -1, 0}, x])/(1 + x) + 
           (4*(2 + x^2)*HPL[{-1, -1, 2}, x])/(1 + x) + 
           ((1 - 5*x + x^2)*HPL[{-1, 0, 0}, x])/(2*(1 + x)) - 
           ((1 + 2*x^2)*HPL[{-1, 2, 0}, x])/(1 + x) - 
           (4*(1 + x^2)*HPL[{-1, 2, 1}, x])/(1 + x) + 
           ((-1 + 22*x + 15*x^2)*HPL[{0, 0, 0}, x])/(8*(1 + x)) + 
           (5*(-1 + x)*HPL[{1, 0, 0}, x])/4 + ((-1 + x)*HPL[{1, 1, 0}, x])/
            2 - (5*(1 + x^2)*HPL[{2, 0, 0}, x])/(4*(1 + x)) - 
           ((1 + x^2)*HPL[{2, 1, 0}, x])/(2*(1 + x)) - 
           (2*(5 + 2*x^2)*HPL[{-1, -1, -1, 0}, x])/(1 + x) + 
           ((13 + 6*x^2)*HPL[{-1, -1, 0, 0}, x])/(1 + x) - 
           ((13 + 9*x^2)*HPL[{-1, 0, 0, 0}, x])/(2*(1 + x)) + 
           ((15 + 7*x^2)*HPL[{0, 0, 0, 0}, x])/(8*(1 + x)) - 
           ((-1 + 13*x + x^2)*Zeta[3])/(4*(1 + x)) + HPL[{0}, x]*
            ((23 - 37*x)/16 - (Pi^2*(-1 + 28*x + 11*x^2))/(48*(1 + x)) - 
             (3*(9 + 4*x^2)*Zeta[3])/(4*(1 + x))) + HPL[{-1}, x]*
            (-(Pi^2*(1 - 5*x + x^2))/(12*(1 + x)) + ((28 + 15*x^2)*Zeta[3])/
              (2*(1 + x))))))
 
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)."
