ISAJET/code/evol05.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE EVOL05
   3 C----------------------------------------------------------------------
   4 C-
   5 C-   Purpose and Methods :
   6 C-        Setup for process 5 (SUPERSYM)
   7 C-        Lorentz frames and perform initial and final QCD jet
   8 C-        evolution in leading-log approximation.
   9 C-
  10 C-   Created  13-AUG-1991   Frank E. Paige,Serban D. Protopopescu
  11 C-
  12 C----------------------------------------------------------------------
  13 #include "isajet/primar.inc"
  14 #include "isajet/jetpar.inc"
  15 #include "isajet/pjets.inc"
  16 #include "isajet/jetset.inc"
  17 #include "isajet/jwork.inc"
  18 #include "isajet/jwork2.inc"
  19 #include "isajet/frame.inc"
  20       REAL    EVOLMS
  21       INTEGER I,K,J,NJSAVE,NJFINL,JTABS
  22 C----------------------------------------------------------------------
  23 C
  24 C          Copy momenta from /PJETS/ to /JETSET/
  25       N0JETS=NJSET+1
  26       CALL IPJSET
  27       NJSAVE=NJSET
  28 C
  29 C          Set flags and maximum off-shell masses and generate
  30 C          initial QCD parton shower.
  31 C
  32       CALL ISTRAD(1.0)
  33 C
  34       IF(NJSET.LT.0) RETURN
  35 C
  36 C
  37 C          Final state evolution.
  38 C          Define Lorentz frames and JMATCH pointers for jet evolution
  39 C          and fragmentation.
  40 C
  41       CALL IFRAMS(N0JETS,NJSAVE,1,.FALSE.)
  42 C
  43 C          Set maximum off-shell masses and JDCAY flags.
  44 C
  45         NJFINL=N0JETS
  46         DO 325 J=N0JETS,NJSAVE
  47           JTABS=IABS(JTYPE(J))
  48           IF(JTABS.GT.20.AND.JTABS.LT.30) THEN
  49           PJSET(5,J)=EVOLMS(J,1.0)
  50           JDCAY(J)=-1
  51         ENDIF
  52 325   CONTINUE
  53 C
  54 C          Produce final-state QCD parton cascade
  55 C
  56       CALL QCDJET(NJFINL)
  57 C
  58       RETURN
  59       END