]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ISAJET/code/evol02.F
- Reset TProcessID count after each event
[u/mrichter/AliRoot.git] / ISAJET / code / evol02.F
CommitLineData
0795afa3 1#include "isajet/pilot.h"
2 SUBROUTINE EVOL02
3C----------------------------------------------------------------------
4C-
5C- Purpose and Methods :
6C- Setup for process 2 (E+E-)
7C- Lorentz frames and perform initial and final QCD jet
8C- evolution in leading-log approximation.
9C-
10C- Created 13-AUG-1991 Frank E. Paige,Serban D. Protopopescu
11C-
12C----------------------------------------------------------------------
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
26C----------------------------------------------------------------------
27C
28C Copy momenta from /PJETS/ to /JETSET/
29 N0JETS=NJSET+1
30 CALL IPJSET
31 NJSAVE=NJSET
32C
33C Final state evolution.
34C Define Lorentz frames and JMATCH pointers for jet evolution
35C and fragmentation.
36C
37 CALL IFRAMS(N0JETS,NJSAVE,1,.FALSE.)
38C
39C Set maximum off-shell masses and JDCAY flags.
40C
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
49310 CONTINUE
50C
51C Produce final-state QCD parton cascade
52C
53 CALL QCDJET(NJFINL)
54C
55 RETURN
56 END