| 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 |
|
|---|
