1 #include "isajet/pilot.h"
3 C----------------------------------------------------------------------
5 C- Purpose and Methods :
6 C- Call for each process a subroutine to set up
7 C- Lorentz frames and perform initial and final QCD jet
8 C- evolution in leading-log approximation.
10 C- Created 13-AUG-1991 Frank E. Paige,Serban D. Protopopescu
12 C----------------------------------------------------------------------
13 #if defined(CERNLIB_IMPNONE)
16 #include "isajet/primar.inc"
17 #include "isajet/jetpar.inc"
18 #include "isajet/pjets.inc"
19 #include "isajet/pinits.inc"
20 #include "isajet/jetset.inc"
21 #include "isajet/jwork.inc"
22 #include "isajet/jwork2.inc"
23 #include "isajet/keys.inc"
24 #include "isajet/frame.inc"
26 INTEGER I,K,J,JJET,IFR
27 C----------------------------------------------------------------------
34 C Copy momenta from /PINITS/ to /JETSET/
38 JORIG(NJSET)=JPACK*(10+I)
39 JTYPE(NJSET)=IDINIT(I)
40 JDCAY(NJSET)=JPACK*I+I
42 105 PJSET(K,NJSET)=PINITS(K,I)
46 C Handle each process separately
48 IF(KEYS(1).OR.KEYS(8)) THEN
56 ELSEIF(KEYS(6).OR.KEYS(10)) THEN
58 ELSEIF(KEYS(7).OR.KEYS(9)) THEN
68 C Boost /JETSET/ partons back to PP COM
75 IF(JJET.GT.10) GO TO 500
76 IF(IDJETS(JJET).EQ.10.AND.KEYS(6)) GO TO 500
81 505 BP=BP+FRAME(K,IFR)*PJSET(K,J)
84 510 PJSET(K,J)=PJSET(K,J)+FRAME(K,IFR)*PJSET(4,J)/FRAME(5,IFR)
85 1 +FRAME(K,IFR)*BP/(FRAME(4,IFR)+FRAME(5,IFR))
86 PJSET(4,J)=FRAME(4,IFR)*PJSET(4,J)/FRAME(5,IFR)+BP
91 IF(JDCAY(J).EQ.JPACK*J+J) THEN
92 JJET=JORIG(J)/JPACK-10
93 PINCOM=.5*(PJSET(4,J)+ABS(PJSET(3,J)))
94 PBEAM(JJET)=HALFE-PINCOM
98 C Check for zero energy partons