]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | SUBROUTINE SSHWW | |
3 | C----------------------------------------------------------------------- | |
4 | C Calculate HL, HH -> WW, ZZ, using either the on-shell matrix | |
5 | C element if kinematically allowed or the WW* or ZZ* matrix | |
6 | C element from Eqn.(6) for Keung and Marciano (PRD. 84: 248). | |
7 | C For the latter, save the mode as W(Z) f fbar, and require that | |
8 | C MH > MW + 2 * MB. | |
9 | C | |
10 | C Bisset's GBDCY | |
11 | C----------------------------------------------------------------------- | |
12 | #if defined(CERNLIB_IMPNONE) | |
13 | IMPLICIT NONE | |
14 | #endif | |
15 | #include "isajet/sspar.inc" | |
16 | #include "isajet/sssm.inc" | |
17 | #include "isajet/sstmp.inc" | |
18 | #include "isajet/sstype.inc" | |
19 | C | |
20 | EXTERNAL SSHWW1,SSHWW2 | |
21 | DOUBLE PRECISION SSHWW1,SSHWW2 | |
22 | DOUBLE PRECISION PI,SR2,G2,BETA,ALPHA,SW2,CW2,CAB2,SAB2,MW,MZ | |
23 | $,MH,COUPL,LOWER,UPPER,FWW1,FWW2,FWW3,FWW,DWID,FZZ | |
24 | DOUBLE PRECISION SSDINT,SSDLAM | |
25 | REAL WID | |
26 | REAL BRZN,BRZL,BRZU,BRZD,BRWL,BRWQ | |
27 | INTEGER IDHHA,IH | |
28 | C Hard wired Z branching ratios | |
29 | DATA BRZN,BRZL,BRZU,BRZD/.06839,.03442,.11792,.15191/ | |
30 | DATA BRWL,BRWQ/.11111,.33333/ | |
31 | C | |
32 | C Mass matrix parameters | |
33 | C | |
34 | PI=4*ATAN(1.D0) | |
35 | SR2=SQRT(2.D0) | |
36 | G2=4*PI*ALFAEM/SN2THW | |
37 | BETA=ATAN(1.0/RV2V1) | |
38 | ALPHA=ALFAH | |
39 | SW2=SN2THW | |
40 | CW2=1.-SN2THW | |
41 | CAB2=(DCOS(ALPHA+BETA))**2 | |
42 | SAB2=1.0-CAB2 | |
43 | MW=AMW | |
44 | MZ=AMZ | |
45 | C | |
46 | C WW* and ZZ* decays | |
47 | C | |
48 | DO 100 IH=1,2 | |
49 | IF(IH.EQ.1) THEN | |
50 | MH=AMHL | |
51 | IDHHA=ISHL | |
52 | COUPL=SAB2 | |
53 | ELSE | |
54 | MH=AMHH | |
55 | IDHHA=ISHH | |
56 | COUPL=CAB2 | |
57 | ENDIF | |
58 | C H -> W + W* -> W + f + fbar | |
59 | TMP(1)=MH | |
60 | IF(MH.GT.MW+2*AMBT.AND.MH.LE.2*MW) THEN | |
61 | LOWER=2*MW/MH | |
62 | UPPER=1+MW**2/MH**2 | |
63 | IF (LOWER.LT.0.998D0) THEN | |
64 | IF (UPPER.LE.1.001D0) THEN | |
65 | FWW1=SSDINT(LOWER,SSHWW1,0.998D0) | |
66 | FWW2=SSDINT(0.998D0,SSHWW1,UPPER) | |
67 | FWW=FWW1+FWW2 | |
68 | ELSEIF(UPPER.GT.1.001D0) THEN | |
69 | FWW1=SSDINT(LOWER,SSHWW1,0.998D0) | |
70 | FWW2=SSDINT(0.998D0,SSHWW1,1.001D0) | |
71 | FWW3=SSDINT(1.001D0,SSHWW1,UPPER) | |
72 | FWW=FWW1+FWW2+FWW3 | |
73 | ENDIF | |
74 | ELSE IF (0.998D0.LT.LOWER.AND.LOWER.LT.1.001D0) THEN | |
75 | IF (UPPER.LE.1.001D0) THEN | |
76 | FWW=SSDINT(LOWER,SSHWW1,UPPER) | |
77 | ELSEIF(UPPER.GT.1.001D0) THEN | |
78 | FWW1=SSDINT(LOWER,SSHWW1,1.001D0) | |
79 | FWW2=SSDINT(1.001D0,SSHWW1,UPPER) | |
80 | FWW=FWW1+FWW2 | |
81 | ENDIF | |
82 | ELSE IF (LOWER.GT.1.001D0) THEN | |
83 | FWW=SSDINT(LOWER,SSHWW1,UPPER) | |
84 | END IF | |
85 | DWID=3*(G2**2)*MH*FWW/(512.0*PI**3) | |
86 | WID=DWID*COUPL | |
87 | CALL SSSAVE(IDHHA,0.5*BRWL*WID,IDW,IDE,-IDNE,0,0) | |
88 | CALL SSSAVE(IDHHA,0.5*BRWL*WID,IDW,IDMU,-IDNM,0,0) | |
89 | CALL SSSAVE(IDHHA,0.5*BRWL*WID,IDW,IDTAU,-IDNT,0,0) | |
90 | CALL SSSAVE(IDHHA,0.5*BRWQ*WID,IDW,-IDUP,IDDN,0,0) | |
91 | CALL SSSAVE(IDHHA,0.5*BRWQ*WID,IDW,-IDCH,IDST,0,0) | |
92 | CALL SSSAVE(IDHHA,0.5*BRWL*WID,-IDW,-IDE,IDNE,0,0) | |
93 | CALL SSSAVE(IDHHA,0.5*BRWL*WID,-IDW,-IDMU,IDNM,0,0) | |
94 | CALL SSSAVE(IDHHA,0.5*BRWL*WID,-IDW,-IDTAU,IDNT,0,0) | |
95 | CALL SSSAVE(IDHHA,0.5*BRWQ*WID,-IDW,IDUP,-IDDN,0,0) | |
96 | CALL SSSAVE(IDHHA,0.5*BRWQ*WID,-IDW,IDCH,-IDST,0,0) | |
97 | ENDIF | |
98 | C H -> Z + Z* -> Z + f + fbar | |
99 | IF(MH.GT.MZ+2*AMBT.AND.MH.LE.2*MZ) THEN | |
100 | LOWER=2*MZ/MH | |
101 | UPPER=1+MZ**2/MH**2 | |
102 | FZZ=SSDINT(LOWER,SSHWW2,UPPER) | |
103 | DWID=7.0-40*SW2/3+160*SW2**2/9 | |
104 | DWID=DWID/CW2**2 | |
105 | DWID=DWID*G2**2*MH*FZZ/(2048*PI**3) | |
106 | WID=DWID*COUPL | |
107 | CALL SSSAVE(IDHHA,BRZN*WID,IDZ,IDNE,-IDNE,0,0) | |
108 | CALL SSSAVE(IDHHA,BRZN*WID,IDZ,IDNM,-IDNM,0,0) | |
109 | CALL SSSAVE(IDHHA,BRZN*WID,IDZ,IDNT,-IDNT,0,0) | |
110 | CALL SSSAVE(IDHHA,BRZL*WID,IDZ,IDE,-IDE,0,0) | |
111 | CALL SSSAVE(IDHHA,BRZL*WID,IDZ,IDMU,-IDMU,0,0) | |
112 | CALL SSSAVE(IDHHA,BRZL*WID,IDZ,IDTAU,-IDTAU,0,0) | |
113 | CALL SSSAVE(IDHHA,BRZU*WID,IDZ,IDUP,-IDUP,0,0) | |
114 | CALL SSSAVE(IDHHA,BRZU*WID,IDZ,IDCH,-IDCH,0,0) | |
115 | CALL SSSAVE(IDHHA,BRZD*WID,IDZ,IDDN,-IDDN,0,0) | |
116 | CALL SSSAVE(IDHHA,BRZD*WID,IDZ,IDST,-IDST,0,0) | |
117 | CALL SSSAVE(IDHHA,BRZD*WID,IDZ,IDBT,-IDBT,0,0) | |
118 | ENDIF | |
119 | 100 CONTINUE | |
120 | C | |
121 | C HH -> WW, ZZ | |
122 | C If these are allowed, the WW* and ZZ* are not. | |
123 | C | |
124 | MH=AMHH | |
125 | IF(MH.GT.2*MW) THEN | |
126 | DWID=3+(MH/MW)**4/4-(MH/MW)**2 | |
127 | DWID=DWID*G2*CAB2*MW**2/(16.0*PI*MH**3) | |
128 | WID=DWID*SQRT(SSDLAM(MH**2,MW**2,MW**2)) | |
129 | CALL SSSAVE(ISHH,WID,IDW,-IDW,0,0,0) | |
130 | ENDIF | |
131 | IF(MH.GT.2*MZ) THEN | |
132 | DWID=3+(MH/MZ)**4/4-(MH/MZ)**2 | |
133 | DWID=DWID*G2*CAB2*MW**2/(16.0*PI*MH**3)/(2.0*CW2**2) | |
134 | WID=DWID*SQRT(SSDLAM(MH**2,MZ**2,MZ**2)) | |
135 | CALL SSSAVE(ISHH,WID,IDZ,IDZ,0,0,0) | |
136 | ENDIF | |
137 | C | |
138 | RETURN | |
139 | END |