ISAJET/isasusy/ssalfs.F

   1 #include "isajet/pilot.h"
   2       DOUBLE PRECISION FUNCTION SSALFS(Q2)
   3 C-----------------------------------------------------------------------
   4 C     Strong coupling formula from page 201 of Barger and Phillips:
   5 C     (using ALQCD4 for 4 flavor Lambda)
   6 C
   7 C     Bisset's STRCPLH
   8 C-----------------------------------------------------------------------
   9 #include "isajet/sssm.inc"
  10 C
  11       DOUBLE PRECISION Q2,AS,TH5,TH6,PI
  12       DATA PI/3.14159265D0/
  13 C
  14       TH5=4*AMBT**2
  15       TH6=4*AMTP**2
  16       IF (Q2.LE.TH5)THEN
  17         AS=12*PI/(25*LOG(Q2/ALQCD4**2))
  18       ELSE IF(Q2.GT.TH5.AND.Q2.LE.TH6) THEN
  19         AS=25*LOG(Q2/ALQCD4**2)-2*LOG(Q2/TH5)
  20         AS=12*PI/AS
  21       ELSEIF(Q2.GT.TH6)THEN
  22         AS=25*LOG(Q2/ALQCD4**2)
  23         AS=AS-2*(LOG(Q2/TH5)+LOG(Q2/TH6))
  24         AS=12*PI/AS
  25       ENDIF
  26       SSALFS=AS
  27       RETURN
  28       END