1 #include "isajet/pilot.h"
4 C Compute the integrated technirho cross section
5 C d(sigma)/d(qmw**2)d(yw) = d(sigma)/d(qmw**2)*f(x1)*f(x2)/scm
6 C including W-technirho mixing from EHLQ 6.22 and 6.23 and
7 C elastic resonance in longitudinal WW fusion.
9 C Use WTYPE for control with
13 C SIGMA = cross section summed over allowed types.
14 C SIGS(I) = partial cross section for I1 + I2 --> I3 + I4.
15 C INOUT(I) = IOPAK**3*I4 + IOPAK**2*I3 + IOPAK*I2 + I1
18 #if defined(CERNLIB_IMPNONE)
21 #include "isajet/itapes.inc"
22 #include "isajet/qcdpar.inc"
23 #include "isajet/jetpar.inc"
24 #include "isajet/primar.inc"
25 #include "isajet/q1q2.inc"
26 #include "isajet/jetsig.inc"
27 #include "isajet/qsave.inc"
28 #include "isajet/wcon.inc"
29 #include "isajet/const.inc"
30 #include "isajet/jetlim.inc"
31 #include "isajet/hcon.inc"
32 #include "isajet/tcpar.inc"
34 REAL AMQCUR(6),WTHELI(4),FINT(9),X(2)
35 EQUIVALENCE (S,SHAT),(T,THAT),(U,UHAT),(X(1),X1)
37 REAL ACOSH,Z,ATANH,AMASS,QMW2,QMZ,EHAT,ANEFF,Q2SAVE,YHAT,EY,AMW,
38 $AMZ,STRUC,STRUCW,WM,ZM,PWWCM,SIG0,S,T,U,FACINV,RATZ,Q1L,Q1R,SIG1,
40 INTEGER I,IH,IQ,IW,IQ1,IQ2,IQ3,IQ4,IRHO,LISTW(4)
42 DATA AMQCUR/.005,.009,.175,1.25,4.50,30./
43 DATA LISTW/10,80,-80,90/
44 DATA MATCHT/0,0,0,0, 0,29,0,27, 0,0,29,28, 0,28,27,0/
47 ACOSH(Z)=ALOG(Z+SQRT(Z**2-1.))
48 ATANH(Z)=.5*ALOG((1.+Z)/(1.-Z))
50 C Kinematics (identical to Drell-Yan)
54 QTMW=SQRT(QMW2+QTW**2)
57 QW=SQRT(QZW**2+QTW**2)
61 IF(ABS(CTHW).LT.1.) THEN
76 ANEFF=4.+QSQ/(QSQ+AMASS(5)**2)+QSQ/(QSQ+AMASS(6)**2)
77 ALFQSQ=12.*PI/((33.-ANEFF)*ALOG(QSQ/ALAM2))
91 IF(X1.GE.1..OR.X2.GE.1.) RETURN
95 C Compute structure functions
99 120 QSAVE(IQ,IH)=STRUC(X(IH),QSQ,IQ,IDIN(IH))/X(IH)
104 IF(QMW.GT.2.*AMW) THEN
105 QSAVE(25+IW,IH)=STRUCW(X(IH),IW,IDIN(IH))/X(IH)
112 C qk + qb --> technirho0
114 IF(.NOT.((GOQ(27,1).AND.GOQ(28,2)).OR.(GOQ(28,1).AND.GOQ(27,2))))
118 IF(QMW.LE.2.*AMW) GO TO 300
119 PWWCM=.5*SQRT(QMW**2-4.*WM**2)
120 SIG0=PI*ALFA**2/(72.*SIN2W*S)*(2.*PWWCM/QMW)**3*X1*X2*UNITS
121 SIG0=SIG0*TCMRHO**2/((S-TCMRHO**2)**2+TCMRHO**2*TCGRHO**2)
125 IF(IQ2.EQ.0) GO TO 210
126 FACINV=2.*SQRT(SIN2W*(1.-SIN2W))
128 Q1L=AQ(IQ1/2,4)*FACINV
129 Q1R=BQ(IQ1/2,4)*FACINV
130 SIG1=.25*SIG0*(1.-RATZ*Q1L/(Q1R*(1.-SIN2W))
131 $ +RATZ**2*(Q1L**2+Q1R**2)/(4.*(1-SIN2W)**2))
132 $ *QSAVE(IQ1,1)*QSAVE(IQ2,2)
136 IF(GOQ(IQ3,1).AND.GOQ(IQ4,2)) THEN
137 SIG=SIG1*TBRWW(IQ3-25,1)*TBRWW(IQ4-25,2)
138 CALL SIGFIL(SIG,IQ1,IQ2,IQ3,IQ4)
143 C W+ + W- -> technirho0 -> W+ + W-
145 SIG0=12*PI/PWWCM**2*TCGRHO**2*X1*X2*UNITS
146 $/((S-TCMRHO**2)**2+TCMRHO**2*TCGRHO**2)
150 SIG1=.25*SIG0*QSAVE(IQ1,1)*QSAVE(IQ2,2)
154 IF(GOQ(IQ3,1).AND.GOQ(IQ4,2)) THEN
155 SIG=SIG1*TBRWW(IQ3-25,1)*TBRWW(IQ4-25,2)
156 CALL SIGFIL(SIG,IQ1,IQ2,IQ3,IQ4)
161 C q + qbar -> technirho+-
163 300 IF(.NOT.((GOQ(27,1).AND.GOQ(29,2)).OR.(GOQ(28,1).AND.GOQ(29,2))
164 $.OR.(GOQ(29,1).AND.GOQ(27,2)).OR.(GOQ(29,1).AND.GOQ(28,2))))
168 IF(QMW.LE.WM+ZM) GO TO 400
169 PWWCM=SQRT((S-WM**2-ZM**2)**2-4.*WM**2*ZM**2)/(2.*QMW)
170 SIG0=PI*ALFA**2/(144.*SIN2W)*S/(S-WM**2)**2*(2.*PWWCM/QMW)**3
172 SIG0=SIG0*TCMRHO**2/((S-TCMRHO**2)**2+TCMRHO**2*TCGRHO**2)
177 IF(IQ2.EQ.0) GO TO 320
178 SIG1=.25*SIG0*QSAVE(IQ1,1)*QSAVE(IQ2,2)
181 IQ4=MATCHT(IQ3-25,IRHO)
182 IF(IQ4.EQ.0) GO TO 330
183 IF(GOQ(IQ3,1).AND.GOQ(IQ4,2)) THEN
184 SIG=SIG1*TBRWW(IQ3-25,1)*TBRWW(IQ4-25,2)
185 CALL SIGFIL(SIG,IQ1,IQ2,IQ3,IQ4)
191 C W+- + Z0 -> technirho+- -> W+- + Z0
193 SIG0=12*PI/PWWCM**2*TCGRHO**2*X1*X2*UNITS
194 $/((S-TCMRHO**2)**2+TCMRHO**2*TCGRHO**2)
198 IQ2=MATCHT(IQ1-25,IRHO)
199 IF(IQ2.EQ.0) GO TO 350
200 SIG1=.25*SIG0*QSAVE(IQ1,1)*QSAVE(IQ2,2)
203 IQ4=MATCHT(IQ3-25,IRHO)
204 IF(IQ4.EQ.0) GO TO 360
205 IF(GOQ(IQ3,1).AND.GOQ(IQ4,2)) THEN
206 SIG=SIG1*TBRWW(IQ3-25,1)*TBRWW(IQ4-25,2)
207 CALL SIGFIL(SIG,IQ1,IQ2,IQ3,IQ4)