1 #include "isajet/pilot.h"
4 C Calculate d(sigma)/d(pt**2)d(y1)d(y2) for
5 C Wh, WH, Zh, ZH, hA, HA and H+H- production in SUSY
7 C SIGMA = cross section summed over types allowed by
9 C SIGS(I) = partial cross section for I1 + I2 --> I3 + I4
10 C INOUT(I) = IOPAK**3*I4 + IOPAK**2*I3 + IOPAK*I2 +I1
12 C Extra factor of 1/2 needed for nonidentical final jets.
13 C Y=-log(tan(theta/2)) gives jacobean P1*P2/E1*E2
16 #if defined(CERNLIB_IMPNONE)
19 #include "isajet/itapes.inc"
20 #include "isajet/const.inc"
21 #include "isajet/jetpar.inc"
22 #include "isajet/jetsig.inc"
23 #include "isajet/primar.inc"
24 #include "isajet/q1q2.inc"
25 #include "isajet/qcdpar.inc"
26 #include "isajet/wcon.inc"
27 #include "isajet/sspar.inc"
31 EQUIVALENCE (S,SHAT),(T,THAT),(U,UHAT)
32 REAL SIG,S,T,U,FAC,AMW,AMZ,AMW2,AMZ2,E1,E2,EQ1
33 REAL QFCN,STRUC,SIGHW,SCFAC,BETA,SINW,COS2W
34 REAL PROPZ,PROPW,GV(2),GA(2),AMH,GAMW,GAMZ
35 INTEGER IS2UD(25),IQ,IH,I,IQ1,IQ2,IFLQ
38 C IS2UD: Susy jettype -> u/d code
39 DATA IS2UD/0,1,1,2,2,2,2,1,1,2,2,1,1,1,1,2,2,2,2,1,1,2,2,1,1/
42 QFCN(IQ,IH)=STRUC(X(IH),QSQ,IQ,IDIN(IH))/X(IH)
66 SCFAC=SIN(ALFAH+BETA)**2
69 SCFAC=COS(ALFAH+BETA)**2
73 C Wh, WH production via W-*
75 IF (GOQ(79,1).AND.GOQ(IH,2)) THEN
76 CALL TWOKIN(0.,0.,AMW,AMH)
77 IF(X1.GE.1..OR.X2.GE.1.) GO TO 100
78 E1=SQRT(P(1)**2+AMW**2)
79 E2=SQRT(P(2)**2+AMH**2)
81 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
82 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
83 SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
84 $ PROPW*TBRWW(3,1)*SCFAC
85 SIG=.5*SIGHW*FAC*QFCN(3,1)*QFCN(4,2)
86 CALL SIGFIL(SIG,3,4,79,IH)
87 SIG=.5*SIGHW*FAC*QFCN(4,1)*QFCN(3,2)
88 CALL SIGFIL(SIG,4,3,79,IH)
89 SIG=.5*SIGHW*FAC*QFCN(9,1)*QFCN(6,2)
90 CALL SIGFIL(SIG,9,6,79,IH)
91 SIG=.5*SIGHW*FAC*QFCN(6,1)*QFCN(9,2)
92 CALL SIGFIL(SIG,6,9,79,IH)
96 IF (GOQ(IH,1).AND.GOQ(79,2)) THEN
97 CALL TWOKIN(0.,0.,AMH,AMW)
98 IF(X1.GE.1..OR.X2.GE.1.) GO TO 110
99 E1=SQRT(P(1)**2+AMH**2)
100 E2=SQRT(P(2)**2+AMW**2)
102 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
103 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
104 SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
105 $ PROPW*TBRWW(3,2)*SCFAC
106 SIG=.5*SIGHW*FAC*QFCN(3,1)*QFCN(4,2)
107 CALL SIGFIL(SIG,3,4,IH,79)
108 SIG=.5*SIGHW*FAC*QFCN(4,1)*QFCN(3,2)
109 CALL SIGFIL(SIG,4,3,IH,79)
110 SIG=.5*SIGHW*FAC*QFCN(9,1)*QFCN(6,2)
111 CALL SIGFIL(SIG,9,6,IH,79)
112 SIG=.5*SIGHW*FAC*QFCN(6,1)*QFCN(9,2)
113 CALL SIGFIL(SIG,6,9,IH,79)
118 C Wh, WH production via W+*
120 IF (GOQ(78,1).AND.GOQ(IH,2)) THEN
121 CALL TWOKIN(0.,0.,AMW,AMH)
122 IF(X1.GE.1..OR.X2.GE.1.) GO TO 120
123 E1=SQRT(P(1)**2+AMW**2)
124 E2=SQRT(P(2)**2+AMH**2)
126 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
127 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
128 SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
129 $ PROPW*TBRWW(2,1)*SCFAC
130 SIG=.5*SIGHW*FAC*QFCN(2,1)*QFCN(5,2)
131 CALL SIGFIL(SIG,2,5,78,IH)
132 SIG=.5*SIGHW*FAC*QFCN(5,1)*QFCN(2,2)
133 CALL SIGFIL(SIG,5,2,78,IH)
134 SIG=.5*SIGHW*FAC*QFCN(8,1)*QFCN(7,2)
135 CALL SIGFIL(SIG,8,7,78,IH)
136 SIG=.5*SIGHW*FAC*QFCN(7,1)*QFCN(8,2)
137 CALL SIGFIL(SIG,7,8,78,IH)
141 IF (GOQ(IH,1).AND.GOQ(78,2)) THEN
142 CALL TWOKIN(0.,0.,AMH,AMW)
143 IF(X1.GE.1..OR.X2.GE.1.) GO TO 130
144 E1=SQRT(P(1)**2+AMH**2)
145 E2=SQRT(P(2)**2+AMW**2)
147 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
148 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
149 SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
150 $ PROPW*TBRWW(2,2)*SCFAC
151 SIG=.5*SIGHW*FAC*QFCN(2,1)*QFCN(5,2)
152 CALL SIGFIL(SIG,2,5,IH,78)
153 SIG=.5*SIGHW*FAC*QFCN(5,1)*QFCN(2,2)
154 CALL SIGFIL(SIG,5,2,IH,78)
155 SIG=.5*SIGHW*FAC*QFCN(8,1)*QFCN(7,2)
156 CALL SIGFIL(SIG,8,7,IH,78)
157 SIG=.5*SIGHW*FAC*QFCN(7,1)*QFCN(8,2)
158 CALL SIGFIL(SIG,7,8,IH,78)
162 C Zh, ZH production via Z*
164 IF (GOQ(80,1).AND.GOQ(IH,2)) THEN
165 CALL TWOKIN(0.,0.,AMZ,AMH)
166 IF(X1.GE.1..OR.X2.GE.1.) GO TO 200
167 E1=SQRT(P(1)**2+AMZ2)
168 E2=SQRT(P(2)**2+AMH**2)
170 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
171 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
175 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 210
176 SIG=GF**2*AMZ**8*(GV(IFLQ)**2+GA(IFLQ)**2)*
177 $ (S/AMZ2+(1.-T/AMZ2)*(1.-U/AMZ2))/PROPZ*TBRWW(4,1)*SCFAC
178 SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
179 CALL SIGFIL(SIG,IQ1,IQ2,80,IH)
183 C hZ, HZ production via Z*
185 IF (GOQ(IH,1).AND.GOQ(80,2)) THEN
186 CALL TWOKIN(0.,0.,AMH,AMZ)
187 IF(X1.GE.1..OR.X2.GE.1.) GO TO 220
188 E1=SQRT(P(1)**2+AMH**2)
189 E2=SQRT(P(2)**2+AMZ2)
191 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
192 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
196 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 230
197 SIG=GF**2*AMZ**8*(GV(IFLQ)**2+GA(IFLQ)**2)*
198 $ (S/AMZ2+(1.-T/AMZ2)*(1.-U/AMZ2))/PROPZ*TBRWW(4,2)*SCFAC
199 SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
200 CALL SIGFIL(SIG,IQ1,IQ2,IH,80)
205 C Next, do Ah and AH production
207 IF (GOQ(83,1).AND.GOQ(IH,2)) THEN
208 CALL TWOKIN(0.,0.,AMHA,AMH)
209 IF(X1.GE.1..OR.X2.GE.1.) GO TO 240
210 E1=SQRT(P(1)**2+AMHA**2)
211 E2=SQRT(P(2)**2+AMH**2)
213 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
214 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
218 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 250
219 SIG=GF**2*AMZ**4*(GV(IFLQ)**2+GA(IFLQ)**2)*
220 $ ((AMHA**2+U-T-AMH**2)*(AMHA**2+T-U-AMH**2)-
221 $ S*(2*AMHA**2+2*AMH**2-S))/PROPZ*SCFAC
222 SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
223 CALL SIGFIL(SIG,IQ1,IQ2,83,IH)
227 IF (GOQ(IH,1).AND.GOQ(83,2)) THEN
228 CALL TWOKIN(0.,0.,AMH,AMHA)
229 IF(X1.GE.1..OR.X2.GE.1.) GO TO 260
230 E1=SQRT(P(1)**2+AMH**2)
231 E2=SQRT(P(2)**2+AMHA**2)
233 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
234 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
238 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 270
239 SIG=GF**2*AMZ**4*(GV(IFLQ)**2+GA(IFLQ)**2)*
240 $ ((AMHA**2+U-T-AMH**2)*(AMHA**2+T-U-AMH**2)-
241 $ S*(2*AMHA**2+2*AMH**2-S))/PROPZ*SCFAC
242 SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
243 CALL SIGFIL(SIG,IQ1,IQ2,IH,83)
249 C Next, do H+H- production
251 IF (GOQ(84,1).AND.GOQ(85,2)) THEN
252 CALL TWOKIN(0.,0.,AMHC,AMHC)
253 IF(X1.GE.1..OR.X2.GE.1.) GO TO 300
254 E1=SQRT(P(1)**2+AMHC**2)
255 E2=SQRT(P(2)**2+AMHC**2)
257 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
258 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
267 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 310
268 SIG=((4*PI*ALFA)**2*EQ1**2/S/S+32*PI*ALFA*EQ1*GF*AMZ**2*
269 $ COS2W*GV(IFLQ)*(S-AMZ**2)/S/PROPZ/SQRT2+8*GF**2*
270 $ AMZ**4*COS2W**2*(GV(IFLQ)**2+GA(IFLQ)**2)/PROPZ)*
271 $ ((U-T)*(T-U)-S*(4*AMHC**2-S))
272 SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
273 CALL SIGFIL(SIG,IQ1,IQ2,84,85)
277 IF (GOQ(85,1).AND.GOQ(84,2)) THEN
278 CALL TWOKIN(0.,0.,AMHC,AMHC)
279 IF(X1.GE.1..OR.X2.GE.1.) GO TO 320
280 E1=SQRT(P(1)**2+AMHC**2)
281 E2=SQRT(P(2)**2+AMHC**2)
283 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
284 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
293 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 330
294 SIG=((4*PI*ALFA)**2*EQ1**2/S/S+32*PI*ALFA*EQ1*GF*AMZ**2*
295 $ COS2W*GV(IFLQ)*(S-AMZ**2)/S/PROPZ/SQRT2+8*GF**2*
296 $ AMZ**4*COS2W**2*(GV(IFLQ)**2+GA(IFLQ)**2)/PROPZ)*
297 $ ((U-T)*(T-U)-S*(4*AMHC**2-S))
298 SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
299 CALL SIGFIL(SIG,IQ1,IQ2,85,84)