Models/ggHFullLoop: functions.f

File functions.f, 1.2 KB (added by dbfranzosi, 5 years ago)

For newer versions of MG.

Line 
1
2
3      double complex function myGGH(ss)
4      implicit none     
5      include 'input.inc' ! include all model parameter
6      INCLUDE 'coupl.inc'
7      complex*16 c0
8      external c0
9      real*8 mb2, mt2,ss
10     
11      mt2=mdl_mt**2
12      mb2=mdl_mb**2
13           
14      myGGH = 3d0/ss*(mdl_mt*mdl_ymt*(2d0 + (4d0*mt2 - ss)*c0(ss,mt2))
15     &   +mdl_mb*mdl_ymb*(2d0 + (4d0*mb2 - ss)*c0(ss,mb2)))
16      return
17      end
18     
19      double complex function myGGA(ss)
20      implicit none     
21      include 'input.inc' ! include all model parameter
22      INCLUDE 'coupl.inc'
23      complex*16 c0
24      external c0
25      real*8 mb2, mt2,ss
26     
27      mt2=mdl_mt**2
28      mb2=mdl_mb**2
29           
30      myGGA = -2d0*(mdl_mb*mdl_ymb*c0(ss,mb2)+mdl_mt*mdl_ymt*c0(ss,mt2))
31      return
32      end
33     
34      double complex function c0(ss,mass2)
35      implicit none
36     
37      COMPLEX*16 IMAG1
38      PARAMETER (IMAG1=(0D0,1D0))     
39      real*8 pi,xx,ss,mass2
40      parameter (pi=3.141592653589793d0)
41     
42      xx=4d0*mass2/ss
43     
44      if (xx.le.1d0) then
45        c0=1d0/(2d0*ss)*
46     &   (log((sqrt(1d0-xx)+1d0)/(1d0-sqrt(1d0-xx)))-imag1*pi)**2
47      else
48        c0=-2d0/ss*dasin(1d0/sqrt(xx))**2
49      endif 
50     
51      end
52
53     
54     
55