]>
Commit | Line | Data |
---|---|---|
e74335a4 | 1 | * $Id$ |
2 | ||
3 | C********************************************************************* | |
4 | ||
5 | SUBROUTINE PYINRE_HIJING | |
6 | ||
7 | C...Calculates full and effective widths of guage bosons, stores masses | |
8 | C...and widths, rescales coefficients to be used for resonance | |
9 | C...production generation. | |
10 | #include "ludat1_hijing.inc" | |
11 | #include "ludat2_hijing.inc" | |
12 | #include "ludat3_hijing.inc" | |
13 | #include "pysubs_hijing.inc" | |
14 | #include "pypars_hijing.inc" | |
15 | #include "pyint1_hijing.inc" | |
16 | #include "pyint2_hijing.inc" | |
17 | #include "pyint4_hijing.inc" | |
18 | #include "pyint6_hijing.inc" | |
19 | DIMENSION WDTP(0:40),WDTE(0:40,0:5) | |
20 | ||
21 | C...Calculate full and effective widths of gauge bosons. | |
22 | AEM=PARU(101) | |
23 | XW=PARU(102) | |
24 | DO 100 I=21,40 | |
25 | DO 100 J=0,40 | |
26 | WIDP(I,J)=0. | |
27 | 100 WIDE(I,J)=0. | |
28 | ||
29 | C...W+/-: | |
30 | WMAS=PMAS(24,1) | |
31 | WFAC=AEM/(24.*XW)*WMAS | |
32 | CALL PYWIDT_HIJING(24,WMAS,WDTP,WDTE) | |
33 | WIDS(24,1)=((WDTE(0,1)+WDTE(0,2))*(WDTE(0,1)+WDTE(0,3))+ | |
34 | &(WDTE(0,1)+WDTE(0,2)+WDTE(0,1)+WDTE(0,3))*(WDTE(0,4)+WDTE(0,5))+ | |
35 | &2.*WDTE(0,4)*WDTE(0,5))/WDTP(0)**2 | |
36 | WIDS(24,2)=(WDTE(0,1)+WDTE(0,2)+WDTE(0,4))/WDTP(0) | |
37 | WIDS(24,3)=(WDTE(0,1)+WDTE(0,3)+WDTE(0,4))/WDTP(0) | |
38 | DO 110 I=0,40 | |
39 | WIDP(24,I)=WFAC*WDTP(I) | |
40 | 110 WIDE(24,I)=WFAC*WDTE(I,0) | |
41 | ||
42 | C...H+/-: | |
43 | HCMAS=PMAS(37,1) | |
44 | HCFAC=AEM/(8.*XW)*(HCMAS/WMAS)**2*HCMAS | |
45 | CALL PYWIDT_HIJING(37,HCMAS,WDTP,WDTE) | |
46 | WIDS(37,1)=((WDTE(0,1)+WDTE(0,2))*(WDTE(0,1)+WDTE(0,3))+ | |
47 | &(WDTE(0,1)+WDTE(0,2)+WDTE(0,1)+WDTE(0,3))*(WDTE(0,4)+WDTE(0,5))+ | |
48 | &2.*WDTE(0,4)*WDTE(0,5))/WDTP(0)**2 | |
49 | WIDS(37,2)=(WDTE(0,1)+WDTE(0,2)+WDTE(0,4))/WDTP(0) | |
50 | WIDS(37,3)=(WDTE(0,1)+WDTE(0,3)+WDTE(0,4))/WDTP(0) | |
51 | DO 120 I=0,40 | |
52 | WIDP(37,I)=HCFAC*WDTP(I) | |
53 | 120 WIDE(37,I)=HCFAC*WDTE(I,0) | |
54 | ||
55 | C...Z0: | |
56 | ZMAS=PMAS(23,1) | |
57 | ZFAC=AEM/(48.*XW*(1.-XW))*ZMAS | |
58 | CALL PYWIDT_HIJING(23,ZMAS,WDTP,WDTE) | |
59 | WIDS(23,1)=((WDTE(0,1)+WDTE(0,2))**2+ | |
60 | &2.*(WDTE(0,1)+WDTE(0,2))*(WDTE(0,4)+WDTE(0,5))+ | |
61 | &2.*WDTE(0,4)*WDTE(0,5))/WDTP(0)**2 | |
62 | WIDS(23,2)=(WDTE(0,1)+WDTE(0,2)+WDTE(0,4))/WDTP(0) | |
63 | WIDS(23,3)=0. | |
64 | DO 130 I=0,40 | |
65 | WIDP(23,I)=ZFAC*WDTP(I) | |
66 | 130 WIDE(23,I)=ZFAC*WDTE(I,0) | |
67 | ||
68 | C...H0: | |
69 | HMAS=PMAS(25,1) | |
70 | HFAC=AEM/(8.*XW)*(HMAS/WMAS)**2*HMAS | |
71 | CALL PYWIDT_HIJING(25,HMAS,WDTP,WDTE) | |
72 | WIDS(25,1)=((WDTE(0,1)+WDTE(0,2))**2+ | |
73 | &2.*(WDTE(0,1)+WDTE(0,2))*(WDTE(0,4)+WDTE(0,5))+ | |
74 | &2.*WDTE(0,4)*WDTE(0,5))/WDTP(0)**2 | |
75 | WIDS(25,2)=(WDTE(0,1)+WDTE(0,2)+WDTE(0,4))/WDTP(0) | |
76 | WIDS(25,3)=0. | |
77 | DO 140 I=0,40 | |
78 | WIDP(25,I)=HFAC*WDTP(I) | |
79 | 140 WIDE(25,I)=HFAC*WDTE(I,0) | |
80 | ||
81 | C...Z'0: | |
82 | ZPMAS=PMAS(32,1) | |
83 | ZPFAC=AEM/(48.*XW*(1.-XW))*ZPMAS | |
84 | CALL PYWIDT_HIJING(32,ZPMAS,WDTP,WDTE) | |
85 | WIDS(32,1)=((WDTE(0,1)+WDTE(0,2)+WDTE(0,3))**2+ | |
86 | &2.*(WDTE(0,1)+WDTE(0,2))*(WDTE(0,4)+WDTE(0,5))+ | |
87 | &2.*WDTE(0,4)*WDTE(0,5))/WDTP(0)**2 | |
88 | WIDS(32,2)=(WDTE(0,1)+WDTE(0,2)+WDTE(0,4))/WDTP(0) | |
89 | WIDS(32,3)=0. | |
90 | DO 150 I=0,40 | |
91 | WIDP(32,I)=ZPFAC*WDTP(I) | |
92 | 150 WIDE(32,I)=ZPFAC*WDTE(I,0) | |
93 | ||
94 | C...R: | |
95 | RMAS=PMAS(40,1) | |
96 | RFAC=0.08*RMAS/((MSTP(1)-1)*(1.+6.*(1.+ULALPS_HIJING(RMAS**2) | |
97 | $ /PARU(1)))) | |
98 | CALL PYWIDT_HIJING(40,RMAS,WDTP,WDTE) | |
99 | WIDS(40,1)=((WDTE(0,1)+WDTE(0,2))*(WDTE(0,1)+WDTE(0,3))+ | |
100 | &(WDTE(0,1)+WDTE(0,2)+WDTE(0,1)+WDTE(0,3))*(WDTE(0,4)+WDTE(0,5))+ | |
101 | &2.*WDTE(0,4)*WDTE(0,5))/WDTP(0)**2 | |
102 | WIDS(40,2)=(WDTE(0,1)+WDTE(0,2)+WDTE(0,4))/WDTP(0) | |
103 | WIDS(40,3)=(WDTE(0,1)+WDTE(0,3)+WDTE(0,4))/WDTP(0) | |
104 | DO 160 I=0,40 | |
105 | WIDP(40,I)=WFAC*WDTP(I) | |
106 | 160 WIDE(40,I)=WFAC*WDTE(I,0) | |
107 | ||
108 | C...Q: | |
109 | KFLQM=1 | |
110 | DO 170 I=1,MIN(8,MDCY(21,3)) | |
111 | IDC=I+MDCY(21,2)-1 | |
112 | IF(MDME(IDC,1).LE.0) GOTO 170 | |
113 | KFLQM=I | |
114 | 170 CONTINUE | |
115 | MINT(46)=KFLQM | |
116 | KFPR(81,1)=KFLQM | |
117 | KFPR(81,2)=KFLQM | |
118 | KFPR(82,1)=KFLQM | |
119 | KFPR(82,2)=KFLQM | |
120 | ||
121 | C...Set resonance widths and branching ratios in JETSET. | |
122 | DO 180 I=1,6 | |
123 | IF(I.LE.3) KC=I+22 | |
124 | IF(I.EQ.4) KC=32 | |
125 | IF(I.EQ.5) KC=37 | |
126 | IF(I.EQ.6) KC=40 | |
127 | PMAS(KC,2)=WIDP(KC,0) | |
128 | PMAS(KC,3)=MIN(0.9*PMAS(KC,1),10.*PMAS(KC,2)) | |
129 | DO 180 J=1,MDCY(KC,3) | |
130 | IDC=J+MDCY(KC,2)-1 | |
131 | BRAT(IDC)=WIDE(KC,J)/WIDE(KC,0) | |
132 | 180 CONTINUE | |
133 | ||
134 | C...Special cases in treatment of gamma*/Z0: redefine process name. | |
135 | IF(MSTP(43).EQ.1) THEN | |
136 | PROC(1)='f + fb -> gamma*' | |
137 | ELSEIF(MSTP(43).EQ.2) THEN | |
138 | PROC(1)='f + fb -> Z0' | |
139 | ELSEIF(MSTP(43).EQ.3) THEN | |
140 | PROC(1)='f + fb -> gamma*/Z0' | |
141 | ENDIF | |
142 | ||
143 | C...Special cases in treatment of gamma*/Z0/Z'0: redefine process name. | |
144 | IF(MSTP(44).EQ.1) THEN | |
145 | PROC(141)='f + fb -> gamma*' | |
146 | ELSEIF(MSTP(44).EQ.2) THEN | |
147 | PROC(141)='f + fb -> Z0' | |
148 | ELSEIF(MSTP(44).EQ.3) THEN | |
149 | PROC(141)='f + fb -> Z''0' | |
150 | ELSEIF(MSTP(44).EQ.4) THEN | |
151 | PROC(141)='f + fb -> gamma*/Z0' | |
152 | ELSEIF(MSTP(44).EQ.5) THEN | |
153 | PROC(141)='f + fb -> gamma*/Z''0' | |
154 | ELSEIF(MSTP(44).EQ.6) THEN | |
155 | PROC(141)='f + fb -> Z0/Z''0' | |
156 | ELSEIF(MSTP(44).EQ.7) THEN | |
157 | PROC(141)='f + fb -> gamma*/Z0/Z''0' | |
158 | ENDIF | |
159 | ||
160 | RETURN | |
161 | END |