]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ISAJET/code/evol01.F
EffC++ warnings corrected.
[u/mrichter/AliRoot.git] / ISAJET / code / evol01.F
CommitLineData
0795afa3 1#include "isajet/pilot.h"
2 SUBROUTINE EVOL01
3C----------------------------------------------------------------------
4C-
5C- Purpose and Methods :
6C- Setup for process 1 (TWOJET)
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/frame.inc"
23 REAL EVOLMS
24 INTEGER I,K,J,NJSAVE,NJFINL
25C----------------------------------------------------------------------
26C
27C Copy momenta from /PJETS/ to /JETSET/
28 N0JETS=NJSET+1
29 CALL IPJSET
30 NJSAVE=NJSET
31C
32C Set flags and maximum off-shell masses and generate
33C initial QCD parton shower.
34C
35 CALL ISTRAD(1.0)
36C
37 IF(NJSET.LT.0) RETURN
38C
39C Final state evolution.
40C Define Lorentz frames and JMATCH pointers for jet evolution
41C and fragmentation.
42C
43 CALL IFRAMS(N0JETS,NJSAVE,1,.FALSE.)
44C
45C Set maximum off-shell masses and JDCAY flags.
46C
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
53310 CONTINUE
54C
55C Produce final-state QCD parton cascade
56C
57 CALL QCDJET(NJFINL)
58C
59 RETURN
60 END