ISAJET/code/evol01.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE EVOL01
   3 C----------------------------------------------------------------------
   4 C-
   5 C-   Purpose and Methods :
   6 C-        Setup for process 1 (TWOJET)
   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 #if defined(CERNLIB_IMPNONE)
  14       IMPLICIT NONE
  15 #endif
  16 #include "isajet/primar.inc"
  17 #include "isajet/jetpar.inc"
  18 #include "isajet/pjets.inc"
  19 #include "isajet/jetset.inc"
  20 #include "isajet/jwork.inc"
  21 #include "isajet/jwork2.inc"
  22 #include "isajet/frame.inc"
  23       REAL    EVOLMS
  24       INTEGER I,K,J,NJSAVE,NJFINL
  25 C----------------------------------------------------------------------
  26 C
  27 C          Copy momenta from /PJETS/ to /JETSET/
  28       N0JETS=NJSET+1
  29       CALL IPJSET
  30       NJSAVE=NJSET
  31 C
  32 C          Set flags and maximum off-shell masses and generate
  33 C          initial QCD parton shower.
  34 C
  35       CALL ISTRAD(1.0)
  36 C
  37       IF(NJSET.LT.0) RETURN
  38 C
  39 C          Final state evolution.
  40 C          Define Lorentz frames and JMATCH pointers for jet evolution
  41 C          and fragmentation.
  42 C
  43       CALL IFRAMS(N0JETS,NJSAVE,1,.FALSE.)
  44 C
  45 C          Set maximum off-shell masses and JDCAY flags.
  46 C
  47       NJFINL=N0JETS
  48       DO 310 J=N0JETS,NJSAVE
  49         IF(IABS(JTYPE(J)).LT.10) THEN
  50           PJSET(5,J)=EVOLMS(J,1.0)
  51           JDCAY(J)=-1
  52         ENDIF
  53 310   CONTINUE
  54 C
  55 C          Produce final-state QCD parton cascade
  56 C
  57       CALL QCDJET(NJFINL)
  58 C
  59       RETURN
  60       END