1 #include "isajet/pilot.h"
4 C COMPUTE THE INTEGRATED WEINBERG-SALAM HIGGS CROSS SECTION
5 C D(SIGMA)/D(QMW**2)D(YW)
7 C SIGMA = CROSS SECTION SUMMED OVER QUARK TYPES ALLOWED BY
8 C JETTYPE3 AND WTYPE CARDS.
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
13 C VER. 7.14: CHECK INITIAL QUARK MASS IS ALLOWED
15 #include "isajet/itapes.inc"
16 #include "isajet/qcdpar.inc"
17 #include "isajet/jetpar.inc"
18 #include "isajet/primar.inc"
19 #include "isajet/q1q2.inc"
20 #include "isajet/jetsig.inc"
21 #include "isajet/qsave.inc"
22 #include "isajet/wcon.inc"
23 #include "isajet/const.inc"
24 #include "isajet/jetlim.inc"
25 #include "isajet/hcon.inc"
27 DIMENSION AMQCUR(6),LISTW(4),WTHELI(4),FINT(9)
29 EQUIVALENCE (S,SHAT),(T,THAT),(U,UHAT),(X(1),X1)
30 #if defined(CERNLIB_DOUBLE)
31 DOUBLE PRECISION C,TERM,SUM,FINT,ZLIM
33 DATA AMQCUR/.005,.009,.175,1.25,4.50,30./
34 DATA LISTW/10,80,-80,90/
35 C WTHELI ARE WEIGHTS OF HELICITY AMPLITUDES IN SIGMA.
36 DATA WTHELI/1.,2.,2.,4./
39 ACOSH(Z)=ALOG(Z+SQRT(Z**2-1.))
40 ATANH(Z)=.5*ALOG((1.+Z)/(1.-Z))
42 C KINEMATICS (IDENTICAL TO DRELL-YAN)
46 QTMW=SQRT(QMW2+QTW**2)
49 QW=SQRT(QZW**2+QTW**2)
53 IF(ABS(CTHW).LT.1.) THEN
68 ANEFF=4.+QSQ/(QSQ+AMASS(5)**2)+QSQ/(QSQ+AMASS(6)**2)
69 ALFQSQ=12.*PI/((33.-ANEFF)*ALOG(QSQ/ALAM2))
83 IF(X1.GE.1..OR.X2.GE.1.) RETURN
85 C COMPUTE STRUCTURE FUNCTIONS
88 120 QSAVE(IQ,IH)=STRUC(X(IH),QSQ,IQ,IDIN(IH))/X(IH)
93 IF(QMW.GT.2.*AMW) THEN
94 QSAVE(25+IW,IH)=STRUCW(X(IH),IW,IDIN(IH))/X(IH)
101 C CALCULATE HIGGS-GLUON-GLUON COUPLING FOR GIVEN Q**2
106 IF(AMQ.LE.0.) GO TO 150
109 ETAR=ETAR+.5*RQ*(1.+(1.-RQ)*ASIN(1./SQRT(RQ))**2)
111 RQLOG=ALOG((1.+SQRT(1.-RQ))/(1.-SQRT(1.-RQ)))
112 PHIR=.25*(RQLOG**2-PI**2)
113 ETAR=ETAR+.5*RQ*(1.+(RQ-1.)*PHIR)
115 ETAI=ETAI+.5*RQ*(RQ-1.)*PHII
118 ETAHGG=ETAR**2+ETAI**2
122 SIG0=GF*ALFQSQ**2/(32.*PI*SQRT2)*ETAHGG*X1*X2*UNITS
123 SIG0=SIG0*S/(PI*HMASS*((S-HMASS**2)**2+(HMASS*HGAM)**2))
124 SIG0=SIG0*QSAVE(1,1)*QSAVE(1,2)
127 IF(GOQ(IQ1,1).AND.GOQ(IQ2,2)) THEN
129 IF(IQ1.GT.25) SIG=SIG*TBRWW(IQ1-25,1)*TBRWW(IQ2-25,2)
130 CALL SIGFIL(SIG,1,1,IQ1,IQ2)
136 SIG0=PI*GF/(3.*SQRT2*HMASS**2)*X1*X2*UNITS
137 SIG0=SIG0*S/(PI*HMASS*((S-HMASS**2)**2+(HMASS*HGAM)**2))
141 IF(QMW.LE.2*AMQ) GO TO 210
142 SIG1=SIG0*AMQ**2*QSAVE(IQ1,1)*QSAVE(IQ2,2)
145 IF(GOQ(IQ3,1).AND.GOQ(IQ4,2)) THEN
147 IF(IQ3.GT.25) SIG=SIG*TBRWW(IQ3-25,1)*TBRWW(IQ4-25,2)
148 CALL SIGFIL(SIG,IQ1,IQ2,IQ3,IQ4)
153 C W+W FUSION AND W+W->W+W IN EFFECTIVE W APPROXIMATION WITH
154 C ANGULAR DISTRIBUTION CUT OFF BY PTMIN.
155 C Z0 Z0 FINAL STATE HAS SYMMETRY FACTOR OF .5
157 IF(QMW.LE.2.*AMASS(80)) GO TO 500
161 IF(.NOT.((GOQ(27,1).AND.GOQ(28,2)).OR.(GOQ(28,1).AND.GOQ(27,2))))
164 PWWCM=.5*SQRT(QMW**2-4.*WM**2)
165 STHLIM=PTMIN(1)/PWWCM
167 ZLIM=SQRT(1.-STHLIM**2)
173 C SUM CROSS SECTION TERMS. I,J RUN OVER AMPLITUDE TERMS.
174 C L RUNS OVER HELICITY STATES. N RUNS OVER POWERS.
175 C REMEMBER THAT L=4 IS MISSING SIN(THETA)/SQRT(2)
179 CALL SIGINT(FINT,ZLIM,ADWWWW(1,I),ADWWWW(2,I),ADWWWW(1,J),
188 314 C=C+ANWWWW(K+1,I,L)*ANWWWW(N-K+1,J,L)
192 TERM=TERM+.5*C*FINT(N+1)-.5*C*FINT(N+3)
194 TERM=TERM+C*FINT(N+1)
200 C ADD INTEGRAL OF IMAGINARY PART SQUARED.
201 SUM=SUM+2.*ZLIM*(WTHELI(1)*AIWWWW(1)**2+WTHELI(2)*AIWWWW(2)**2
202 $+WTHELI(3)*AIWWWW(3)**2+WTHELI(4)*AIWWWW(4)**2)
204 SIG0=SUM/(32.*PI*S*SCM)*UNITS
205 SIG1=.5*SIG0*QSAVE(27,1)*QSAVE(28,2)
206 IF(GOQ(27,1).AND.GOQ(28,2)) THEN
207 SIG=SIG1*TBRWW(2,1)*TBRWW(3,2)
208 CALL SIGFIL(SIG,27,28,27,28)
210 IF(GOQ(28,1).AND.GOQ(27,2)) THEN
211 SIG=SIG1*TBRWW(3,1)*TBRWW(2,2)
212 CALL SIGFIL(SIG,27,28,28,27)
214 SIG1=.5*SIG0*QSAVE(28,1)*QSAVE(27,2)
215 IF(GOQ(27,1).AND.GOQ(28,2)) THEN
216 SIG=SIG1*TBRWW(2,1)*TBRWW(3,2)
217 CALL SIGFIL(SIG,28,27,27,28)
219 IF(GOQ(28,1).AND.GOQ(27,2)) THEN
220 SIG=SIG1*TBRWW(3,1)*TBRWW(2,2)
221 CALL SIGFIL(SIG,28,27,28,27)
227 IF(QMW.LE.2.*AMASS(90)) GO TO 500
229 C SUM CROSS SECTION TERMS. I,J RUN OVER AMPLITUDE TERMS.
230 C L RUNS OVER HELICITY STATES. N RUNS OVER POWERS.
231 C REMEMBER THAT L=4 IS MISSING SIN(THETA)/SQRT(2)
235 CALL SIGINT(FINT,ZLIM,ADWWWW(1,I),ADWWWW(2,I),ADWWWW(1,J),
244 324 C=C+ANWWWW(K+1,I,L)*ANWWWW(N-K+1,J,L)
248 TERM=TERM+.5*C*FINT(N+1)-.5*C*FINT(N+3)
250 TERM=TERM+C*FINT(N+1)
256 C ADD INTEGRAL OF IMAGINARY PART SQUARED.
257 SUM=SUM+2.*ZLIM*(WTHELI(1)*AIWWWW(1)**2+WTHELI(2)*AIWWWW(2)**2
258 $+WTHELI(3)*AIWWWW(3)**2+WTHELI(4)*AIWWWW(4)**2)
260 SIG0=SUM/(32.*PI*S*SCM)*UNITS
261 SIG1=.5*SIG0*QSAVE(29,1)*QSAVE(29,2)
262 IF(GOQ(27,1).AND.GOQ(28,2)) THEN
263 SIG=SIG1*TBRWW(2,1)*TBRWW(3,2)
264 CALL SIGFIL(SIG,29,29,27,28)
266 IF(GOQ(28,1).AND.GOQ(27,2)) THEN
267 SIG=SIG1*TBRWW(3,1)*TBRWW(2,2)
268 CALL SIGFIL(SIG,29,29,28,27)
273 400 IF(QMW.LE.2.*AMASS(90)) GO TO 500
274 IF(.NOT.(GOQ(29,1).AND.GOQ(29,2))) GO TO 500
276 PWWCM=.5*SQRT(QMW**2-4.*WM**2)
277 STHLIM=PTMIN(1)/PWWCM
279 ZLIM=SQRT(1.-STHLIM**2)
285 C SUM CROSS SECTION TERMS. I,J RUN OVER AMPLITUDE TERMS.
286 C L RUNS OVER HELICITY STATES. N RUNS OVER POWERS.
287 C REMEMBER THAT L=4 IS MISSING SIN(THETA)/SQRT(2)
291 CALL SIGINT(FINT,ZLIM,ADWWWW(1,I),ADWWWW(2,I),ADWWWW(1,J),
300 414 C=C+ANWWWW(K+1,I,L)*ANWWWW(N-K+1,J,L)
304 TERM=TERM+.5*C*FINT(N+1)-.5*C*FINT(N+3)
306 TERM=TERM+C*FINT(N+1)
312 C ADD INTEGRAL OF IMAGINARY PART SQUARED.
313 SUM=SUM+2.*ZLIM*(WTHELI(1)*AIWWWW(1)**2+WTHELI(2)*AIWWWW(2)**2
314 $+WTHELI(3)*AIWWWW(3)**2+WTHELI(4)*AIWWWW(4)**2)
316 SIG0=SUM/(32.*PI*S*SCM)*UNITS
318 SIG0=SIG0*TBRWW(4,1)*TBRWW(4,2)
319 SIG=SIG0*QSAVE(27,1)*QSAVE(28,2)
320 CALL SIGFIL(SIG,27,28,29,29)
321 SIG=SIG0*QSAVE(28,1)*QSAVE(27,2)
322 CALL SIGFIL(SIG,28,27,29,29)
328 C SUM CROSS SECTION TERMS. I,J RUN OVER AMPLITUDE TERMS.
329 C L RUNS OVER HELICITY STATES. N RUNS OVER POWERS.
330 C REMEMBER THAT L=4 IS MISSING SIN(THETA)/SQRT(2)
334 CALL SIGINT(FINT,ZLIM,ADWWWW(1,I),ADWWWW(2,I),ADWWWW(1,J),
343 424 C=C+ANWWWW(K+1,I,L)*ANWWWW(N-K+1,J,L)
347 TERM=TERM+.5*C*FINT(N+1)-.5*C*FINT(N+3)
349 TERM=TERM+C*FINT(N+1)
355 C ADD INTEGRAL OF IMAGINARY PART SQUARED.
356 SUM=SUM+2.*ZLIM*(WTHELI(1)*AIWWWW(1)**2+WTHELI(2)*AIWWWW(2)**2
357 $+WTHELI(3)*AIWWWW(3)**2+WTHELI(4)*AIWWWW(4)**2)
359 SIG0=SUM/(32.*PI*S*SCM)*UNITS
361 SIG0=SIG0*TBRWW(4,1)*TBRWW(4,2)
362 SIG=SIG0*QSAVE(29,1)*QSAVE(29,2)
363 CALL SIGFIL(SIG,29,29,29,29)