1 #include "isajet/pilot.h"
4 C Calculate d(sigma)/d(pt**2)d(y1)d(y2) for WH and ZH
5 C associated production.
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/hcon.inc"
28 #include "isajet/xmssm.inc"
32 EQUIVALENCE (S,SHAT),(T,THAT),(U,UHAT)
33 REAL SIG,S,T,U,FAC,AMW,AMZ,AMW2,AMZ2,E1,E2
35 REAL PROPZ,PROPW,GV(2),GA(2),AMH,AMH2,GAMW,GAMZ
36 INTEGER IS2UD(25),IQ,IH,I,IQ1,IQ2,IFLQ
39 C IS2UD: Susy jettype -> u/d code
40 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/
43 QFCN(IQ,IH)=STRUC(X(IH),QSQ,IQ,IDIN(IH))/X(IH)
68 C WH production via W-*
70 IF (GOQ(28,1).AND.GOQ(30,2)) THEN
71 CALL TWOKIN(0.,0.,AMW,AMH)
72 IF(X1.GE.1..OR.X2.GE.1.) GO TO 100
73 E1=SQRT(P(1)**2+AMW**2)
74 E2=SQRT(P(2)**2+AMH**2)
76 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
77 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
78 SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
80 SIG=.5*SIGHW*FAC*QFCN(3,1)*QFCN(4,2)
81 CALL SIGFIL(SIG,3,4,28,30)
82 SIG=.5*SIGHW*FAC*QFCN(4,1)*QFCN(3,2)
83 CALL SIGFIL(SIG,4,3,28,30)
84 SIG=.5*SIGHW*FAC*QFCN(9,1)*QFCN(6,2)
85 CALL SIGFIL(SIG,9,6,28,30)
86 SIG=.5*SIGHW*FAC*QFCN(6,1)*QFCN(9,2)
87 CALL SIGFIL(SIG,6,9,28,30)
91 IF (GOQ(30,1).AND.GOQ(28,2)) THEN
92 CALL TWOKIN(0.,0.,AMH,AMW)
93 IF(X1.GE.1..OR.X2.GE.1.) GO TO 110
94 E1=SQRT(P(1)**2+AMH**2)
95 E2=SQRT(P(2)**2+AMW**2)
97 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
98 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
99 SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
101 SIG=.5*SIGHW*FAC*QFCN(3,1)*QFCN(4,2)
102 CALL SIGFIL(SIG,3,4,30,28)
103 SIG=.5*SIGHW*FAC*QFCN(4,1)*QFCN(3,2)
104 CALL SIGFIL(SIG,4,3,30,28)
105 SIG=.5*SIGHW*FAC*QFCN(9,1)*QFCN(6,2)
106 CALL SIGFIL(SIG,9,6,30,28)
107 SIG=.5*SIGHW*FAC*QFCN(6,1)*QFCN(9,2)
108 CALL SIGFIL(SIG,6,9,30,28)
113 C WH production via W+*
115 IF (GOQ(27,1).AND.GOQ(30,2)) THEN
116 CALL TWOKIN(0.,0.,AMW,AMH)
117 IF(X1.GE.1..OR.X2.GE.1.) GO TO 120
118 E1=SQRT(P(1)**2+AMW**2)
119 E2=SQRT(P(2)**2+AMH**2)
121 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
122 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
123 SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
125 SIG=.5*SIGHW*FAC*QFCN(2,1)*QFCN(5,2)
126 CALL SIGFIL(SIG,2,5,27,30)
127 SIG=.5*SIGHW*FAC*QFCN(5,1)*QFCN(2,2)
128 CALL SIGFIL(SIG,5,2,27,30)
129 SIG=.5*SIGHW*FAC*QFCN(8,1)*QFCN(7,2)
130 CALL SIGFIL(SIG,8,7,27,30)
131 SIG=.5*SIGHW*FAC*QFCN(7,1)*QFCN(8,2)
132 CALL SIGFIL(SIG,7,8,27,30)
136 IF (GOQ(30,1).AND.GOQ(27,2)) THEN
137 CALL TWOKIN(0.,0.,AMH,AMW)
138 IF(X1.GE.1..OR.X2.GE.1.) GO TO 130
139 E1=SQRT(P(1)**2+AMH**2)
140 E2=SQRT(P(2)**2+AMW**2)
142 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
143 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
144 SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
146 SIG=.5*SIGHW*FAC*QFCN(2,1)*QFCN(5,2)
147 CALL SIGFIL(SIG,2,5,30,27)
148 SIG=.5*SIGHW*FAC*QFCN(5,1)*QFCN(2,2)
149 CALL SIGFIL(SIG,5,2,30,27)
150 SIG=.5*SIGHW*FAC*QFCN(8,1)*QFCN(7,2)
151 CALL SIGFIL(SIG,8,7,30,27)
152 SIG=.5*SIGHW*FAC*QFCN(7,1)*QFCN(8,2)
153 CALL SIGFIL(SIG,7,8,30,27)
157 C ZH production via Z*
159 IF (GOQ(29,1).AND.GOQ(30,2)) THEN
160 CALL TWOKIN(0.,0.,AMZ,AMH)
161 IF(X1.GE.1..OR.X2.GE.1.) GO TO 200
162 E1=SQRT(P(1)**2+AMZ2)
163 E2=SQRT(P(2)**2+AMH2)
165 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
166 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
170 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 210
171 SIG=GF**2*AMZ**8*(GV(IFLQ)**2+GA(IFLQ)**2)*
172 $ (S/AMZ2+(1.-T/AMZ2)*(1.-U/AMZ2))/PROPZ*TBRWW(4,1)
173 SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
174 CALL SIGFIL(SIG,IQ1,IQ2,29,30)
178 C HZ production via Z*
180 IF (GOQ(30,1).AND.GOQ(29,2)) THEN
181 CALL TWOKIN(0.,0.,AMH,AMZ)
182 IF(X1.GE.1..OR.X2.GE.1.) GO TO 220
183 E1=SQRT(P(1)**2+AMH2)
184 E2=SQRT(P(2)**2+AMZ2)
186 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
187 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
191 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 230
192 SIG=GF**2*AMZ**8*(GV(IFLQ)**2+GA(IFLQ)**2)*
193 $ (S/AMZ2+(1.-T/AMZ2)*(1.-U/AMZ2))/PROPZ*TBRWW(4,2)
194 SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
195 CALL SIGFIL(SIG,IQ1,IQ2,30,29)