ISAJET/isasusy/sugfrz.F

   1 #include "isajet/pilot.h"
   2 C------------------------------------------------------------------
   3       SUBROUTINE SUGFRZ(Q,G,G0,IG)
   4 C------------------------------------------------------------------
   5 C
   6 C     Freeze out final soft breaking parameters
   7 C
   8 #if defined(CERNLIB_IMPNONE)
   9       IMPLICIT NONE
  10 #endif
  11 #include "isajet/sslun.inc"
  12 #include "isajet/sugpas.inc"
  13       DIMENSION G(29),G0(29)
  14       INTEGER IG(29)
  15       REAL Q,MT
  16       REAL G,G0,TANB
  17       INTEGER I
  18 C
  19       TANB=XTANB
  20       MT=AMT
  21       DO 200 I=1,5
  22         G0(I)=G(I)
  23 200   CONTINUE
  24       IF (Q.LT.MT.AND.IG(6).EQ.0) THEN
  25         G3MT=G(3)
  26         G0(6)=G(6)
  27         IG(6)=1
  28       END IF
  29 C          Freeze out running gluino mass at MGL
  30       DO 210 I=7,12
  31         IF (Q.LT.ABS(G(I)).AND.IG(I).EQ.0) THEN
  32           G0(I)=G(I)
  33           IG(I)=1
  34         ELSE IF (IG(I).EQ.0) THEN
  35           G0(I)=G(I)
  36         END IF
  37 210   CONTINUE
  38 C          Freeze out Higgs paremeters at HIGFRZ
  39       DO 211 I=13,14
  40         IF (Q.LT.HIGFRZ.AND.IG(I).EQ.0) THEN
  41           G0(I)=G(I)
  42           IG(I)=1
  43           G0(I+12)=G(I+12)
  44           IG(I+12)=1
  45         ELSE IF (IG(I).EQ.0) THEN
  46           G0(I)=G(I)
  47           G0(I+12)=G(I+12)
  48         END IF
  49 211   CONTINUE
  50 C          Freeze out rest at own masses
  51       DO 220 I=15,24
  52 C        IF (G(I).LT.0.) THEN
  53 C          G(I)=0.
  54 C          NOGOOD=1
  55 C          GO TO 100
  56 C        END IF
  57         IF (Q.LT.SQRT(ABS(G(I))).AND.IG(I).EQ.0) THEN
  58           G0(I)=G(I)
  59           IG(I)=1
  60         ELSE IF (IG(I).EQ.0) THEN
  61           G0(I)=G(I)
  62         END IF
  63 220   CONTINUE
  64 C          Freeze our N_R parameters at Majorana mass scale
  65       DO 230 I=27,29
  66         IF (G(I).NE.0.) G0(I)=G(I)
  67         IF (Q.LT.AMNRMJ.AND.IG(I).EQ.0.) THEN
  68           IG(I)=1
  69         END IF
  70 230   CONTINUE
  71 100   RETURN
  72       END