ISAJET/isasusy/sssave.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE SSSAVE(IIN,GAM,IOUT1,IOUT2,IOUT3,IOUT4,IOUT5)
   3 C-----------------------------------------------------------------------
   4 C     Store a SUSY decay mode in /SSMODE/
   5 C     Ver 7.14: Increment NSSMOD only after test
   6 C-----------------------------------------------------------------------
   7 #if defined(CERNLIB_IMPNONE)
   8       IMPLICIT NONE
   9 #endif
  10 C
  11 #include "isajet/sslun.inc"
  12 #include "isajet/ssmode.inc"
  13 C
  14       INTEGER IIN,IOUT1,IOUT2,IOUT3,IOUT4,IOUT5,I
  15       REAL GAM
  16 C
  17       IF (GAM.LE.0.) THEN
  18         IF(GAM.LT.0.) THEN
  19           WRITE(LOUT,1000) IIN,IOUT1,IOUT2,IOUT3,IOUT4,IOUT5,GAM
  20 1000      FORMAT(' WARNING: SSSAVE: ',I5,' --> ',5I5,E14.5)
  21         ENDIF
  22         LSSMOD=.FALSE.
  23         GO TO 999
  24       ENDIF
  25       NSSMOD=NSSMOD+1
  26       LSSMOD=.TRUE.
  27       IF(NSSMOD.GT.MXSS) THEN
  28         WRITE(LOUT,*) 'SSSAVE: TOO MANY MODES, IIN = ',IIN
  29         STOP99
  30       ENDIF
  31       ISSMOD(NSSMOD)=IIN
  32       JSSMOD(1,NSSMOD)=IOUT1
  33       JSSMOD(2,NSSMOD)=IOUT2
  34       JSSMOD(3,NSSMOD)=IOUT3
  35       JSSMOD(4,NSSMOD)=IOUT4
  36       JSSMOD(5,NSSMOD)=IOUT5
  37       GSSMOD(NSSMOD)=GAM
  38       BSSMOD(NSSMOD)=0.
  39       MSSMOD(NSSMOD)=0
  40 C          Check that quarks and gluons appear at end of list.
  41       DO 100 I=1,4
  42         IF(IABS(JSSMOD(I,NSSMOD)).LE.9.AND.
  43      $  IABS(JSSMOD(I+1,NSSMOD)).GT.9) THEN
  44           WRITE(LOUT,1100) IIN,IOUT1,IOUT2,IOUT3,IOUT4,IOUT5
  45 1100      FORMAT(' WARNING: SSSAVE: BAD ORDER: ',I5,' --> ',5I5)
  46           STOP99
  47         ENDIF
  48 100   CONTINUE
  49 C
  50 999   RETURN
  51       END