]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | SUBROUTINE EVOL03 | |
3 | C---------------------------------------------------------------------- | |
4 | C- | |
5 | C- Purpose and Methods : | |
6 | C- Setup for process 3 (DRELLYAN) | |
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/q1q2.inc" | |
23 | #include "isajet/frame.inc" | |
24 | #include "isajet/wcon.inc" | |
25 | REAL EVOLMS,BP | |
26 | INTEGER I,K,J,NJFINL | |
27 | C---------------------------------------------------------------------- | |
28 | C | |
29 | C Add W momentum and recoil jets | |
30 | N0JETS=NJSET+1 | |
31 | IF(.NOT.STDDY) THEN | |
32 | DO 101 I=3,NJET | |
33 | NJSET=NJSET+1 | |
34 | JORIG(NJSET)=JPACK*I | |
35 | JTYPE(NJSET)=IDJETS(I) | |
36 | JDCAY(NJSET)=0 | |
37 | DO 105 K=1,5 | |
38 | 105 PJSET(K,NJSET)=PJETS(K,I) | |
39 | IFRAME(I)=1 | |
40 | 101 CONTINUE | |
41 | NJSET=NJSET+1 | |
42 | N0W=NJSET | |
43 | JORIG(NJSET)=0 | |
44 | JTYPE(NJSET)=IDENTW | |
45 | JDCAY(NJSET)=(N0W+1)*JPACK+N0W+2 | |
46 | DO 120 K=1,5 | |
47 | 120 PJSET(K,NJSET)=QWJET(K) | |
48 | ENDIF | |
49 | C | |
50 | C Add W decays | |
51 | DO 110 I=1,2 | |
52 | NJSET=NJSET+1 | |
53 | JORIG(NJSET)=JPACK*I | |
54 | JTYPE(NJSET)=IDJETS(I) | |
55 | JDCAY(NJSET)=0 | |
56 | DO 115 K=1,5 | |
57 | 115 PJSET(K,NJSET)=PJETS(K,I) | |
58 | IFRAME(I)=2 | |
59 | IF(STDDY) IFRAME(I)=1 | |
60 | 110 CONTINUE | |
61 | C | |
62 | C Set flags and maximum off-shell masses and generate | |
63 | C initial QCD parton shower. | |
64 | C | |
65 | CALL ISTRAD(WFUDGE) | |
66 | C | |
67 | IF(NJSET.LT.0) RETURN | |
68 | C | |
69 | C Final state evolution. | |
70 | C Define Lorentz frames and JMATCH pointers for jet evolution | |
71 | C and fragmentation. | |
72 | C | |
73 | IF(STDDY) THEN | |
74 | CALL IFRAMS(3,4,1,.FALSE.) | |
75 | ELSE | |
76 | CALL IFRAMS(N0W+1,N0W+2,2,.FALSE.) | |
77 | CALL IFRAMS(N0JETS,N0W,1,.FALSE.) | |
78 | ENDIF | |
79 | C | |
80 | C Set maximum off-shell masses and JDCAY flags. | |
81 | C | |
82 | IF(STDDY) THEN | |
83 | NJFINL=3 | |
84 | DO 310 J=3,4 | |
85 | IF(IABS(JTYPE(J)).LT.10) THEN | |
86 | PJSET(5,J)=QMW | |
87 | JDCAY(J)=-1 | |
88 | ENDIF | |
89 | 310 CONTINUE | |
90 | ELSE | |
91 | NJFINL=N0JETS | |
92 | DO 320 J=N0W+1,N0W+2 | |
93 | IF(IABS(JTYPE(J)).LT.10) THEN | |
94 | PJSET(5,J)=QMW | |
95 | JDCAY(J)=-1 | |
96 | ENDIF | |
97 | 320 CONTINUE | |
98 | C Need fudge factor for DRELLYAN | |
99 | DO 321 J=N0JETS,N0W | |
100 | IF(IABS(JTYPE(J)).LT.10) THEN | |
101 | PJSET(5,J)=EVOLMS(J,WFUDGE) | |
102 | JDCAY(J)=-1 | |
103 | ENDIF | |
104 | 321 CONTINUE | |
105 | ENDIF | |
106 | C | |
107 | C Produce final-state QCD parton cascade | |
108 | C | |
109 | CALL QCDJET(NJFINL) | |
110 | C | |
111 | C Reset FRAME using W momentum modified by evolution | |
112 | IF(.NOT.STDDY) THEN | |
113 | BP=0. | |
114 | DO 400 K=1,3 | |
115 | 400 BP=BP+FRAME(K,1)*PJSET(K,N0W) | |
116 | BP=BP/FRAME(5,1) | |
117 | DO 410 K=1,3 | |
118 | FRAME(K,2)=PJSET(K,N0W)+FRAME(K,1)*PJSET(4,N0W)/FRAME(5,1) | |
119 | $ +FRAME(K,1)*BP/(FRAME(4,1)+FRAME(5,1)) | |
120 | 410 CONTINUE | |
121 | FRAME(4,2)=FRAME(4,1)*PJSET(4,N0W)/FRAME(5,1)+BP | |
122 | ENDIF | |
123 | C | |
124 | RETURN | |
125 | END |