update vel series

feynamps/ggH_++.m

-(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
-(Alfas*EL*MT^2*(4*MH^2 - (MH^2 - 4*MT^2)*
-    Log[(-MH^2 + 2*MT^2 + MH*Sqrt[MH^2 - 4*MT^2])/(2*MT^2)]^2)*
-  Mat[SUNT[Glu1, Glu2, 0, 0]])/(4*MH^2*MW*Pi*SW)

feynamps/ggH_++_vel.m

-(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
-(Alfas*EL*MT^2*(4*MH^2 - (MH^2 - 4*MT^2)*
-    Log[(-MH^2 + 2*MT^2 + MH*Sqrt[MH^2 - 4*MT^2])/(2*MT^2)]^2)*
-  Mat[SUNT[Glu1, Glu2, 0, 0]])/(4*MH^2*MW*Pi*SW)

feynamps/ggH_+-.m

-(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
-0

feynamps/ggH_+-_vel.m

-(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
-0

feynamps/ggH_-+.m

-(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
-0

feynamps/ggH_-+_vel.m

-(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
-0

feynamps/ggH_--.m

-(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
-(Alfas*EL*MT^2*(4*MH^2 - (MH^2 - 4*MT^2)*
-    Log[(-MH^2 + 2*MT^2 + MH*Sqrt[MH^2 - 4*MT^2])/(2*MT^2)]^2)*
-  Mat[SUNT[Glu1, Glu2, 0, 0]])/(4*MH^2*MW*Pi*SW)

feynamps/ggH_--_vel.m

-(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
-(Alfas*EL*MT^2*(4*MH^2 - (MH^2 - 4*MT^2)*
-    Log[(-MH^2 + 2*MT^2 + MH*Sqrt[MH^2 - 4*MT^2])/(2*MT^2)]^2)*
-  Mat[SUNT[Glu1, Glu2, 0, 0]])/(4*MH^2*MW*Pi*SW)

feynamps/ggH_LR.m

 (* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
-(Alfas*EL*MT^2*Mat[SUNT[Glu1, Glu2, 0, 0]]*
-  ((2 + Eps^(-1) + DiscB[MH^2, Sqrt[MT^2], Sqrt[MT^2]] + Log[Mu^2/MT^2] - 
-     4*(3/4 + DiscB[MH^2, Sqrt[MT^2], Sqrt[MT^2]]/4 + 
-       (MT^2*Log[(-MH^2 + 2*MT^2 + Sqrt[-(MH^2*(-MH^2 + 4*MT^2))])/(2*MT^2)]^
-          2)/(4*MH^2) + (Eps^(-1) + Log[Mu^2/MT^2])/4))*Pair[e[1], e[2]] + 
-   4*(1/(2*MH^2) + 
-     (MT^2*Log[(-MH^2 + 2*MT^2 + Sqrt[-(MH^2*(-MH^2 + 4*MT^2))])/(2*MT^2)]^2)/
-      (2*MH^4))*Pair[e[1], k[2]]*Pair[e[2], k[1]] - 
-   (Log[(-MH^2 + 2*MT^2 + Sqrt[-(MH^2*(-MH^2 + 4*MT^2))])/(2*MT^2)]^2*
-     (-(MH^2*Pair[e[1], e[2]]) + 2*Pair[e[1], k[2]]*Pair[e[2], k[1]]))/
+(Alfas*c*EL*MT^2*
+  ((-1 + 2*Log[MT] - (MT^2*Log[(-MH^2 + 2*MT^2 + MH*Sqrt[MH^2 - 4*MT^2])/
+          (2*MT^2)]^2)/MH^2 - Log[Mu^2] + Log[Mu^2/MT^2])*Pair[e[1], e[2]] + 
+   (2*(MH^2 + MT^2*Log[(-MH^2 + 2*MT^2 + MH*Sqrt[MH^2 - 4*MT^2])/(2*MT^2)]^2)*
+     Pair[e[1], k[2]]*Pair[e[2], k[1]])/MH^4 + 
+   (Log[(-MH^2 + 2*MT^2 + MH*Sqrt[MH^2 - 4*MT^2])/(2*MT^2)]^2*
+     (MH^2*Pair[e[1], e[2]] - 2*Pair[e[1], k[2]]*Pair[e[2], k[1]]))/
     (4*MH^2)))/(MW*Pi*SW)

feynamps/ggH_LRS.m

 (* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
-(Alfas*EL*MT^2*Mat[SUNT[Glu1, Glu2, 0, 0]]*
-  ((2 + Eps^(-1) + DiscB[MH^2, Sqrt[MT^2], Sqrt[MT^2]] + Log[Mu^2/MT^2] - 
-     4*(3/4 + DiscB[MH^2, Sqrt[MT^2], Sqrt[MT^2]]/4 + 
-       (MT^2*Log[(-MH^2 + 2*MT^2 + Sqrt[-(MH^2*(-MH^2 + 4*MT^2))])/(2*MT^2)]^
-          2)/(4*MH^2) + (Eps^(-1) + Log[Mu^2/MT^2])/4))*Pair[e[1], e[2]] + 
-   4*(1/(2*MH^2) + 
-     (MT^2*Log[(-MH^2 + 2*MT^2 + Sqrt[-(MH^2*(-MH^2 + 4*MT^2))])/(2*MT^2)]^2)/
-      (2*MH^4))*Pair[e[1], k[2]]*Pair[e[2], k[1]] - 
-   (Log[(-MH^2 + 2*MT^2 + Sqrt[-(MH^2*(-MH^2 + 4*MT^2))])/(2*MT^2)]^2*
-     (-(MH^2*Pair[e[1], e[2]]) + 2*Pair[e[1], k[2]]*Pair[e[2], k[1]]))/
+(Alfas*c*EL*MT^2*
+  ((-1 + 2*Log[MT] - (MT^2*Log[(-MH^2 + 2*MT^2 + MH*Sqrt[MH^2 - 4*MT^2])/
+          (2*MT^2)]^2)/MH^2 - Log[Mu^2] + Log[Mu^2/MT^2])*Pair[e[1], e[2]] + 
+   (2*(MH^2 + MT^2*Log[(-MH^2 + 2*MT^2 + MH*Sqrt[MH^2 - 4*MT^2])/(2*MT^2)]^2)*
+     Pair[e[1], k[2]]*Pair[e[2], k[1]])/MH^4 + 
+   (Log[(-MH^2 + 2*MT^2 + MH*Sqrt[MH^2 - 4*MT^2])/(2*MT^2)]^2*
+     (MH^2*Pair[e[1], e[2]] - 2*Pair[e[1], k[2]]*Pair[e[2], k[1]]))/
     (4*MH^2)))/(MW*Pi*SW)

feynamps/ggH_PV.m

 (* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
-(Alfas*EL*MT^2*Mat[SUNT[Glu1, Glu2, 0, 0]]*
-  (-((-(MH^2*Pair[e[1], e[2]]) + 2*Pair[e[1], k[2]]*Pair[e[2], k[1]])*
+(Alfas*c*EL*MT^2*
+  (((MH^2*Pair[e[1], e[2]] - 2*Pair[e[1], k[2]]*Pair[e[2], k[1]])*
      PVC[0, 0, 0, 0, MH^2, 0, Sqrt[MT^2], Sqrt[MT^2], Sqrt[MT^2]])/2 + 
    4*Pair[e[1], k[2]]*Pair[e[2], k[1]]*PVC[0, 1, 1, 0, MH^2, 0, Sqrt[MT^2], 
      Sqrt[MT^2], Sqrt[MT^2]] + Pair[e[1], e[2]]*