* dse_HHHH_gl.F
* generated 25-Sep-2020 15:35
* this file is part of FeynHiggs
* please do not edit directly

#include "externals.h"
#include "types.h"
#include "debug.h"


	subroutine dse_HHHH_gl(se,k2)
	implicit none
	ComplexType se
	RealType k2

#include "FH.h"
#include "looptools.h"


#ifdef DETAILED_DEBUG
	DHIGGS "dse_HHHH_gl k2 =", k2 ENDL
#endif

	se = -(3/(64.D0*Pi)*
     &      (Alfa1L*S2B**2*Mf2(tM1,3)*
     &         (4*(B1(k2,Mf2(tM1,3),Mf2(tM1,3)) + 
     &              k2*DB1(k2,Mf2(tM1,3),Mf2(tM1,3))) + 
     &           8*DB0(k2,Mf2(tM1,3),Mf2(tM1,3))*Mf2(tM1,3) - 
     &           DB0(k2,MStgl2(2),MStgl2(2))*
     &            (2*Mf(tM1,3) - YtglC*UCStgl(1,3) - 
     &               Ytgl*UCStglC(1,3))**2 - 
     &           DB0(k2,MStgl2(1),MStgl2(1))*
     &            (2*Mf(tM1,3) + YtglC*UCStgl(1,3) + 
     &               Ytgl*UCStglC(1,3))**2 - 
     &           2*DB0(k2,MStgl2(1),MStgl2(2))*
     &            (YtglC*UCStgl(3,3) + Ytgl*UCStgl(3,4))*
     &            (Ytgl*UCStglC(3,3) + YtglC*UCStglC(3,4))))/
     &       (MW2*SB2**2*SW2))

#ifdef DETAILED_DEBUG
	DHIGGS "dse_HHHH_gl  t/st =", se ENDL
#endif

	if( mssmpart .le. 1 ) return

	se = se - 3/(16.D0*Pi)*
     &     (Alfa1L*SB2*Mf2(bM1,3)*
     &        (4*(B1(k2,Mf2(bM1,3),Mf2(bM1,3)) + 
     &             k2*DB1(k2,Mf2(bM1,3),Mf2(bM1,3))) + 
     &          8*DB0(k2,Mf2(bM1,3),Mf2(bM1,3))*Mf2(bM1,3) - 
     &          DB0(k2,MSbgl2(2),MSbgl2(2))*
     &           (2*Mf(bM1,3) - YbglC*UCSbgl(1,3) - 
     &              Ybgl*UCSbglC(1,3))**2 - 
     &          DB0(k2,MSbgl2(1),MSbgl2(1))*
     &           (2*Mf(bM1,3) + YbglC*UCSbgl(1,3) + 
     &              Ybgl*UCSbglC(1,3))**2 - 
     &          2*DB0(k2,MSbgl2(1),MSbgl2(2))*
     &           (YbglC*UCSbgl(3,3) + Ybgl*UCSbgl(3,4))*
     &           (Ybgl*UCSbglC(3,3) + YbglC*UCSbglC(3,4))))/
     &      (CB2*MW2*SW2)

#ifdef DETAILED_DEBUG
	DHIGGS "dse_HHHH_gl +b/sb =", se ENDL
#endif

	end
