1 #include "isajet/pilot.h"
3 C----------------------------------------------------------------------
5 C- Purpose and Methods :
6 C- Setup for process 7 (HIGGS)
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/frame.inc"
25 INTEGER I,K,J,NJSAVE,NJFINL,JTRUE
26 DOUBLE PRECISION DPASS(5),DSUM(5)
28 C----------------------------------------------------------------------
30 C Copy momenta from /PJETS/ to /JETSET/
34 C Add extra momenta for WPAIR
35 IDABS1=IABS(IDJETS(1))
36 IDABS2=IABS(IDJETS(2))
37 IF(IDABS1.EQ.80.OR.IDABS1.EQ.90.OR.IDABS2.EQ.80.OR.
42 JORIG(NJSET)=JPACK*JPAIR(J)
43 JTYPE(NJSET)=IDPAIR(J)
46 135 PJSET(K,NJSET)=PPAIR(K,J)
51 JDCAY(N0JETS+JET-1)=JTRUE*JPACK+JTRUE+1
56 C Set flags and maximum off-shell masses and generate
57 C initial QCD parton shower.
59 IF(IABS(IDINIT(1)).LT.80) THEN
64 C Special initial state evolution for W-W fusion.
72 PJSET(5,J)=-PJSET(5,JET)
74 CALL QCDINI(NJSET-3,NJSET-1)
79 C Final state evolution.
80 C Define Lorentz frames and JMATCH pointers for jet evolution
88 CALL DBLVEC(PJSET(1,N0JETS),DSUM)
89 CALL DBLVEC(PJSET(1,N0JETS+1),DPASS)
91 231 DSUM(K)=DSUM(K)+DPASS(K)
92 DSUM(5)=DSQRT(DSUM(4)**2-DSUM(1)**2-DSUM(2)**2-DSUM(3)**2)
100 FRAME(K,I)=PJSET(K,N0JETS+I-1)
106 C Set up and generate final state QCD parton shower.
107 C Boost PJSET with -FRAME.
112 IF(JET.GT.10) GO TO 240
113 C Do this boost in double precision for 32-bit machines
114 CALL DBOOST(-1,FRAME(1,JET),PJSET(1,J))
117 C Set maximum off-shell masses and JDCAY flags.
121 DO 340 J=N0JETS,NJSAVE
122 IF(IABS(JTYPE(J)).LT.10) THEN
123 PJSET(5,J)=EVOLMS(J,1.0)
130 IF(IABS(JTYPE(N0PAIR+J-1)).LT.10) THEN
131 PJSET(5,N0PAIR+J-1)=PJETS(5,JPAIR(J))
137 C Produce final-state QCD parton cascade