& -TWO*(U*CI(K)+T*CI(J)+U1*CI(K1)+T1*CI(J1)-U*T*DI(I))/S**2 )
RETURN
+ END
C-----------------------------------------------------------------------
- ENTRY HWHIG2(S,T,U,EH2,EQ2,I,J,K,I1,J1,K1)
+ FUNCTION HWHIG2(S,T,U,EH2,EQ2,I,J,K,I1,J1,K1)
C-----------------------------------------------------------------------
C ++- helicity amplitude for: g+g --> g+H
+C-----------------------------------------------------------------------
+C-----------------------------------------------------------------------
+
+ IMPLICIT NONE
+
+ DOUBLE COMPLEX HWHIG2,BI(4),CI(7),DI(3)
+
+ DOUBLE PRECISION S,T,U,EH2,EQ2,S1,T1,U1,ONE,TWO,FOUR,HALF
+
+ INTEGER I,J,K,I1,J1,K1
+
+ COMMON/CINTS/BI,CI,DI
+
+ PARAMETER (ONE =1.D0, TWO =2.D0, FOUR =4.D0, HALF =0.5D0)
+
C-----------------------------------------------------------------------
S1=S-EH2
& -HALF*(EH2-FOUR*EQ2)*(S*T*DI(3)+U*S*DI(2)+U*T*DI(1)) )/(S*T*U)
RETURN
+ END
C-----------------------------------------------------------------------
- ENTRY HWHIG5(S,T,U,EH2,EQ2,I,J,K,I1,J1,K1)
+ FUNCTION HWHIG5(S,T,U,EH2,EQ2,I,J,K,I1,J1,K1)
C-----------------------------------------------------------------------
C Amplitude for: q+qbar --> g+H
+C-----------------------------------------------------------------------
+C-----------------------------------------------------------------------
+
+ IMPLICIT NONE
+
+ DOUBLE COMPLEX HWHIG5,BI(4),CI(7),DI(3)
+
+ DOUBLE PRECISION S,T,U,EH2,EQ2,S1,T1,U1,ONE,TWO,FOUR,HALF
+
+ INTEGER I,J,K,I1,J1,K1
+
+ COMMON/CINTS/BI,CI,DI
+
+ PARAMETER (ONE =1.D0, TWO =2.D0, FOUR =4.D0, HALF =0.5D0)
+
C-----------------------------------------------------------------------
HWHIG5=DCMPLX(TWO+TWO*S/(S-EH2))*BI(I)+DCMPLX(FOUR*EQ2-U-T)*CI(K)