ISAJET/code/nogood.F

   1 #include "isajet/pilot.h"
   2       LOGICAL FUNCTION NOGOOD(KK)
   3 C
   4 C          Insure proper distribution and check kinematics.
   5 C          Select jet types.
   6 C
   7 #if defined(CERNLIB_IMPNONE)
   8       IMPLICIT NONE
   9 #endif
  10 #include "isajet/itapes.inc"
  11 #include "isajet/keys.inc"
  12 #include "isajet/wcon.inc"
  13 #include "isajet/const.inc"
  14 #include "isajet/wsig.inc"
  15 #include "isajet/wgen.inc"
  16 #include "isajet/dylim.inc"
  17 #include "isajet/jetlim.inc"
  18 #include "isajet/jetpar.inc"
  19 #include "isajet/jetsig.inc"
  20 #include "isajet/ptpar.inc"
  21 #include "isajet/hcon.inc"
  22 #include "isajet/xmssm.inc"
  23 C
  24       REAL RANF,SIGINV,SUM,TRY,BRANCH
  25       INTEGER KK,I,II,K,IFL
  26 C
  27       NOGOOD=.TRUE.
  28       GO TO (1,2,3,4,5,6),KK
  29 C
  30 C          TWOJET, SUPERSYM, WPAIR or PHOTON events
  31 C
  32 1     CONTINUE
  33       IF(KEYS(1)) THEN
  34         CALL SIGQCD
  35       ELSEIF(KEYS(5)) THEN
  36         CALL SIGSSY
  37       ELSEIF(KEYS(6)) THEN
  38         CALL SIGWW
  39       ELSEIF(KEYS(8)) THEN
  40         CALL SIGGAM
  41       ELSEIF(KEYS(10)) THEN
  42         CALL SIGWH
  43       ENDIF
  44       IF(SIGMA.LE.0) RETURN
  45       IF(SIGMAX*RANF().GT.SIGMA) RETURN
  46       NOGOOD=.FALSE.
  47       SIGINV=1./SIGMA
  48       SUM=0.
  49       TRY=RANF()
  50       DO 100 I=1,NSIGS
  51         SUM=SUM+SIGS(I)*SIGINV
  52         IF(SUM.LT.TRY) GO TO 100
  53 C          Find reaction
  54         ISIGS=I
  55         SIGEVT=SIGS(ISIGS)
  56         II=INOUT(I)
  57         DO 110 K=1,2
  58         INITYP(K)=MOD(II,IOPAK)
  59 110     II=II/IOPAK
  60         DO 120 K=1,2
  61         JETTYP(K)=MOD(II,IOPAK)
  62 120     II=II/IOPAK
  63         RETURN
  64 100   CONTINUE
  65       RETURN
  66 C
  67 C          DRELLYAN events--test of SIGDY
  68 C
  69 2     CONTINUE
  70       IF(KEYS(3)) THEN
  71         CALL SIGDY
  72       ELSEIF(KEYS(7).AND..NOT.GOMSSM) THEN
  73         CALL SIGH
  74       ELSEIF(KEYS(7).AND.GOMSSM) THEN
  75         CALL SIGHSS
  76       ELSEIF(KEYS(9)) THEN
  77         CALL SIGTC
  78       ELSEIF(KEYS(11)) THEN
  79         CALL SIGKKG
  80       ENDIF
  81       IF(SIGMA.LE.0.) RETURN
  82       IF(SIGSL(KSEL)*RANF().GT.SIGMA) RETURN
  83       NOGOOD=.FALSE.
  84       SIGINV=1./SIGMA
  85       SUM=0.
  86       TRY=RANF()
  87 C          Find reaction.
  88       DO 200 I=1,NSIGS
  89         SUM=SUM+SIGS(I)*SIGINV
  90         IF(SUM.LT.TRY) GO TO 200
  91         ISIGS=I
  92         SIGEVT=SIGS(ISIGS)
  93         GO TO 210
  94 200   CONTINUE
  95 C          Unpack INOUT to find JETTYP and INITYP
  96 210   IF(KEYS(3).OR.KEYS(11)) THEN
  97         II=INOUT(I)
  98         DO 220 K=1,2
  99         INITYP(K)=MOD(II,IOPAK)
 100 220     II=II/IOPAK
 101         JWTYP=MOD(II,IOPAK)
 102         II=II/IOPAK
 103         JETTYP(3)=MOD(II,IOPAK)
 104       ELSEIF(KEYS(7).OR.KEYS(9)) THEN
 105         II=INOUT(ISIGS)
 106         DO 230 I=1,2
 107         INITYP(I)=MOD(II,IOPAK)
 108 230     II=II/IOPAK
 109         DO 240 I=1,2
 110         JETTYP(I)=MOD(II,IOPAK)
 111 240     II=II/IOPAK
 112       ENDIF
 113       RETURN
 114 C
 115 C          DRELLYAN events--test of SIGDY2
 116 C
 117 3     CONTINUE
 118       IF(KEYS(3)) THEN
 119         CALL SIGDY2
 120         IFL=JETTYP(1)/2
 121         BRANCH=(AQ(IFL,JWTYP)**2+BQ(IFL,JWTYP)**2)/COUT(JWTYP)
 122       ELSEIF(KEYS(7).AND..NOT.GOMSSM) THEN
 123         CALL SIGH2
 124         BRANCH=1.
 125       ELSEIF(KEYS(7).AND.GOMSSM) THEN
 126         SIGLLQ=SIGMA/(4*PI)
 127         NOGOOD=.FALSE.
 128         RETURN
 129       ELSEIF(KEYS(9)) THEN
 130         CALL SIGTC2
 131         BRANCH=1.
 132       ENDIF
 133       IF(SIGLLQ.GT.SIGS(ISIGS)*BRANCH*3.*RANF()/(4.*PI))
 134      1NOGOOD=.FALSE.
 135       RETURN
 136 C
 137 C          DRELLYAN events--test of kinematics
 138 C
 139 4     CONTINUE
 140       DO 400 I=1,2
 141         IF(P(I).LT.PMIN(I).OR.P(I).GT.PMAX(I)) GO TO 410
 142         IF(PT(I).LT.PTMIN(I).OR.PT(I).GT.PTMAX(I)) GO TO 410
 143         IF(YJ(I).LT.YJMIN(I).OR.YJ(I).GT.YJMAX(I)) GO TO 410
 144         IF(PHI(I).LT.PHIMIN(I).OR.PHI(I).GT.PHIMAX(I)) GO TO 410
 145 400   CONTINUE
 146       NOGOOD=.FALSE.
 147 410   RETURN
 148 C
 149 5     CONTINUE
 150 6     CONTINUE
 151       RETURN
 152 C
 153       END