1 #include "isajet/pilot.h"
4 C Compute the integrated MSSM Higgs cross section
5 C d(sigma)/d(QMW**2)d(YW)
6 C Since SUSY Higgs are always narrow, can use the widths to
7 C determine couplings and ignore interference with continuum.
9 C SIGMA = cross section summed over quark types allowed by
10 C JETTYPE and WTYPE cards.
11 C SIGS(I) = partial cross section for I1 + I2 --> I3 + I4.
12 C INOUT(I) = IOPAK**3*I4 + IOPAK**2*I3 + IOPAK*I2 + I1
13 C using JETTYPE code from LISTSS.
15 C Ver 7.18: Correct GOQ's and include TBRWW for W/Z modes.
17 #if defined(CERNLIB_IMPNONE)
20 #include "isajet/itapes.inc"
21 #include "isajet/qcdpar.inc"
22 #include "isajet/jetpar.inc"
23 #include "isajet/primar.inc"
24 #include "isajet/q1q2.inc"
25 #include "isajet/jetsig.inc"
26 #include "isajet/qsave.inc"
27 #include "isajet/wcon.inc"
28 #include "isajet/const.inc"
29 #include "isajet/jetlim.inc"
30 #include "isajet/hcon.inc"
34 REAL AM1,AM2,S,T,U,Q2SAVE,YHAT,EY,ANEFF,QMW2,QZW,EHAT,SIG0,SIG,
36 INTEGER JT1,JT2,I,J,IH,IQ,I1,I2,JTGL,JTOFF
37 EQUIVALENCE (S,SHAT),(T,THAT),(U,UHAT),(X(1),X1)
39 C Kinematics (identical to Drell-Yan)
42 QTMW=SQRT(QMW2+QTW**2)
45 QW=SQRT(QZW**2+QTW**2)
49 IF(ABS(CTHW).LT.1.) THEN
64 ANEFF=4.+QSQ/(QSQ+AMASS(5)**2)+QSQ/(QSQ+AMASS(6)**2)
65 ALFQSQ=12.*PI/((33.-ANEFF)*ALOG(QSQ/ALAM2))
78 IF(X1.GE.1..OR.X2.GE.1.) RETURN
80 C Compute structure functions
84 120 QSAVE(IQ,IH)=STRUC(X(IH),QSQ,IQ,IDIN(IH))/X(IH)
92 SIG0=PI*HMASS**2/(8*S**2)*HGAMSS(JTGL,JTGL)*X1*X2*UNITS
93 $/((S-HMASS**2)**2+(HMASS*HGAM)**2)
94 SIG0=SIG0*QSAVE(1,1)*QSAVE(1,2)
97 IF(HGAMSS(I,J).EQ.0) GO TO 210
98 IF(.NOT.(GOQ(I,1).AND.GOQ(J,2))) GO TO 210
100 C Include W/Z branching ratios
101 IF((I.GE.78.AND.I.LE.80).AND.(J.GE.78.AND.J.LE.80)) THEN
102 SIG=SIG*TBRWW(I-76,1)*TBRWW(J-76,2)
104 CALL SIGFIL(SIG,JTGL,JTGL,I,J)
111 C Note I1,I2 run over quarks; JT1,JT2,I,J over LISTSS
118 IF(HGAMSS(JT1,JT2).LE.0) GO TO 310
119 SIG0=4*PI*HMASS**2/(9*S**2)*HGAMSS(JT1,JT2)*X1*X2*UNITS
120 $ /((S-HMASS**2)**2+(HMASS*HGAM)**2)
121 SIG0=SIG0*QSAVE(I1,1)*QSAVE(I2,2)
122 C Decay partial cross sections
125 IF(HGAMSS(I,J).EQ.0) GO TO 330
126 IF(.NOT.(GOQ(I,1).AND.GOQ(J,2))) GO TO 330
128 C Include W/Z branching ratios
129 IF((I.GE.78.AND.I.LE.80).AND.(J.GE.78.AND.J.LE.80)) THEN
130 SIG=SIG*TBRWW(I-76,1)*TBRWW(J-76,2)
132 CALL SIGFIL(SIG,JT1,JT2,I,J)