ISAJET/code/evol11.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE EVOL11
   3 C----------------------------------------------------------------------
   4 C-
   5 C-   Purpose and Methods :
   6 C-        Setup for process 11 (EXTRADIM)
   7 C-        Lorentz frames and perform initial and final QCD jet
   8 C-        evolution in leading-log approximation.
   9 C-
  10 C----------------------------------------------------------------------
  11 #if defined(CERNLIB_IMPNONE)
  12       IMPLICIT NONE
  13 #endif
  14 #include "isajet/primar.inc"
  15 #include "isajet/jetpar.inc"
  16 #include "isajet/pjets.inc"
  17 #include "isajet/jetset.inc"
  18 #include "isajet/jwork.inc"
  19 #include "isajet/jwork2.inc"
  20 #include "isajet/q1q2.inc"
  21 #include "isajet/frame.inc"
  22 #include "isajet/wcon.inc"
  23 C
  24       INTEGER K,NJFINL,J
  25       REAL EVOLMS
  26 C----------------------------------------------------------------------
  27 C
  28 C          Add recoil jet (jet 3)
  29       NJSET=NJSET+1
  30       N0JETS=NJSET
  31       JORIG(NJSET)=JPACK*3
  32       JTYPE(NJSET)=IDJETS(3)
  33       JDCAY(NJSET)=0
  34       DO 105 K=1,5
  35 105   PJSET(K,NJSET)=PJETS(K,3)
  36       IFRAME(3)=1
  37
  38 C          Add W (=KKG)
  39       NJSET=NJSET+1
  40       N0W=NJSET
  41       JORIG(NJSET)=0
  42       JTYPE(NJSET)=IDENTW
  43       JDCAY(NJSET)=0
  44       DO 120 K=1,5
  45 120   PJSET(K,NJSET)=QWJET(K)
  46 C
  47 C          Set flags and maximum off-shell masses and generate
  48 C          initial QCD parton shower.
  49 C
  50       CALL ISTRAD(1.0)
  51       IF(NJSET.LT.0) RETURN
  52 C
  53 C          Final state evolution.
  54 C          Define Lorentz frames and JMATCH pointers for jet evolution
  55 C          and fragmentation.
  56 C
  57       CALL IFRAMS(N0JETS,N0W,1,.FALSE.)
  58 C
  59 C          Set maximum off-shell masses and JDCAY flags.
  60 C
  61       NJFINL=N0JETS
  62       DO 321 J=N0JETS,N0W
  63         IF(IABS(JTYPE(J)).LT.10) THEN
  64           PJSET(5,J)=EVOLMS(J,WFUDGE)
  65           JDCAY(J)=-1
  66         ENDIF
  67 321   CONTINUE
  68 C
  69 C          Produce final-state QCD parton cascade
  70 C
  71       CALL QCDJET(NJFINL)
  72 C
  73       RETURN
  74       END