ISAJET/isasusy/sspole.F

   1 #include "isajet/pilot.h"
   2       REAL FUNCTION SSPOLE(MGMS,MUSQ,AS)
   3 C*********************************************************************
   4 C* Computes the on-shell (pole) gluino mass for given running (MSbar)*
   5 C* gluino mass, defined at scale MUSQ, and given alpha_s (AS). The   *
   6 C* squark masses are stored in the SQUARK COMMON block.              *
   7 C* This function needs the complex functions B0 and B1.              *
   8 C* Contributed by M. Drees; modified by H. Baer                      *
   9 C                                                                    *
  10 C  Version 7.30: Cast COMPLEX*16 to REAL*8 in standard way. :-(      *
  11 C*********************************************************************
  12 #if defined(CERNLIB_IMPNONE)
  13       IMPLICIT NONE
  14 #endif
  15 #include "isajet/sssm.inc"
  16 #include "isajet/sspar.inc"
  17 #include "isajet/sstmp.inc"
  18 #include "isajet/ssinf.inc"
  19       REAL MGMS,MUSQ,AS,MGSQ,FAC
  20       DOUBLE PRECISIONDMUSQ,DFAC
  21       COMPLEX*16 SSB1
  22       DMUSQ=MUSQ
  23       XLAM = LOG(DMUSQ)
  24       MGSQ = MGMS*MGMS
  25 C
  26 C          Cast COMPLEX*16 to REAL*8:
  27 C
  28       DFAC = -8*SSB1(MGSQ,0.,AMULSS) - (SSB1(MGSQ,AMTP,AMT1SS)+
  29      $SSB1(MGSQ,AMTP,AMT2SS)+SSB1(MGSQ,4.0,AMB1SS)+
  30      $SSB1(MGSQ,4.0,AMB2SS))
  31       DFAC = DFAC + 12.D0 + 9.D0*LOG(DMUSQ/MGSQ)
  32       FAC=DFAC
  33       SSPOLE = MGMS*(1.0 + .0796*AS*FAC )
  34       RETURN
  35       END