ISAJET/code/evol02.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE EVOL02
   3 C----------------------------------------------------------------------
   4 C-
   5 C-   Purpose and Methods :
   6 C-        Setup for process 2 (E+E-)
   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/keys.inc"
  23 #include "isajet/frame.inc"
  24       REAL    EVOLMS
  25       INTEGER I,K,J,NJSAVE,NJFINL
  26 C----------------------------------------------------------------------
  27 C
  28 C          Copy momenta from /PJETS/ to /JETSET/
  29       N0JETS=NJSET+1
  30       CALL IPJSET
  31       NJSAVE=NJSET
  32 C
  33 C          Final state evolution.
  34 C          Define Lorentz frames and JMATCH pointers for jet evolution
  35 C          and fragmentation.
  36 C
  37       CALL IFRAMS(N0JETS,NJSAVE,1,.FALSE.)
  38 C
  39 C          Set maximum off-shell masses and JDCAY flags.
  40 C
  41       NJFINL=N0JETS
  42       DO 310 J=N0JETS,NJSAVE
  43         IF((IABS(JTYPE(J)).LT.10).OR.
  44      $     (IABS(JTYPE(J)).GE.21.AND.IABS(JTYPE(J)).LE.29).OR.
  45      $     (IABS(JTYPE(J)).GE.41.AND.IABS(JTYPE(J)).LE.46))THEN
  46           PJSET(5,J)=EVOLMS(J,1.0)
  47           JDCAY(J)=-1
  48         ENDIF
  49 310   CONTINUE
  50 C
  51 C          Produce final-state QCD parton cascade
  52 C
  53       CALL QCDJET(NJFINL)
  54 C
  55       RETURN
  56       END