]>
Commit | Line | Data |
---|---|---|
1 | #include "isajet/pilot.h" | |
2 | SUBROUTINE SSTPBF | |
3 | C----------------------------------------------------------------------- | |
4 | C | |
5 | C Calculate the top branching ratios. | |
6 | C Source: H. Baer (modified by F. Paige) | |
7 | C | |
8 | C----------------------------------------------------------------------- | |
9 | #if defined(CERNLIB_IMPNONE) | |
10 | IMPLICIT NONE | |
11 | #endif | |
12 | #include "isajet/sslun.inc" | |
13 | #include "isajet/ssmode.inc" | |
14 | #include "isajet/sssm.inc" | |
15 | #include "isajet/sspar.inc" | |
16 | #include "isajet/sstype.inc" | |
17 | C | |
18 | COMPLEX ZI,ZONE,ZA,ZB,ZPP,ZPM,ZAUIZ,ZBUIZ | |
19 | REAL SSXLAM,G,AL2,BE2,TANB,COTB,GTBW,GTBH,BWLL,GF,BWQQ,PI,SR2 | |
20 | REAL WID,AS,BS,MZIZ,CS2THW,GP,FT,FB,SNZI,THIZ | |
21 | REAL SINT,COST,SINB,COSB,AWI,BWI,AMW1,AMW2,SNWI | |
22 | REAL THX,THY,XM,YM,BETA,ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS | |
23 | DOUBLE PRECISION SSMQCD | |
24 | INTEGER IZ,ISZIZ(4) | |
25 | DATA ZONE/(1.,0.)/,ZI/(0.,1.)/ | |
26 | C | |
27 | PI=4.*ATAN(1.) | |
28 | SR2=SQRT(2.) | |
29 | G=SQRT(4*PI*ALFAEM/SN2THW) | |
30 | GP=G*SQRT(SN2THW/(1.-SN2THW)) | |
31 | TANB=1./RV2V1 | |
32 | COTB=1./TANB | |
33 | BETA=ATAN(TANB) | |
34 | CS2THW=1.-SN2THW | |
35 | ASMB=SUALFS(AMBT**2,.36,AMTP,3) | |
36 | MBMB=AMBT*(1.-4*ASMB/3./PI) | |
37 | MBQ=SSMQCD(DBLE(MBMB),DBLE(AMTP)) | |
38 | ASMT=SUALFS(AMTP**2,.36,AMTP,3) | |
39 | MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))* | |
40 | $(ASMT/PI)**2) | |
41 | MTQ=SSMQCD(DBLE(MTMT),DBLE(AMTP)) | |
42 | FB=G*MBQ/SR2/AMW/COS(BETA) | |
43 | FT=G*MTQ/SR2/AMW/SIN(BETA) | |
44 | SINT=SIN(THETAT) | |
45 | COST=COS(THETAT) | |
46 | SINB=SIN(THETAB) | |
47 | COSB=COS(THETAB) | |
48 | ISZIZ(1)=ISZ1 | |
49 | ISZIZ(2)=ISZ2 | |
50 | ISZIZ(3)=ISZ3 | |
51 | ISZIZ(4)=ISZ4 | |
52 | XM=1./TAN(GAMMAL) | |
53 | YM=1./TAN(GAMMAR) | |
54 | THX=SIGN(1.,XM) | |
55 | THY=SIGN(1.,YM) | |
56 | C | |
57 | C W decays | |
58 | C | |
59 | GF=1.16E-5 | |
60 | GTBW=GF*AMTP**3*SSXLAM(1.,AMW**2/AMTP**2,AMBT**2/AMTP**2)* | |
61 | $((1.-AMBT**2/AMTP**2)**2+AMW**2/AMTP**2*(1.+AMBT**2/AMTP**2) | |
62 | $-2*AMW**4/AMTP**4)/(8.*PI*SR2) | |
63 | BWQQ=3./9. | |
64 | BWLL=1./9. | |
65 | CALL SSSAVE(IDTP,BWQQ*GTBW,IDUP,-IDDN,IDBT,0,0) | |
66 | CALL SSSAVE(IDTP,BWQQ*GTBW,IDCH,-IDST,IDBT,0,0) | |
67 | CALL SSSAVE(IDTP,BWLL*GTBW,-IDE,IDNE,IDBT,0,0) | |
68 | CALL SSSAVE(IDTP,BWLL*GTBW,-IDMU,IDNM,IDBT,0,0) | |
69 | CALL SSSAVE(IDTP,BWLL*GTBW,-IDTAU,IDNT,IDBT,0,0) | |
70 | C | |
71 | C H+ decays | |
72 | C | |
73 | AL2=(G/2/SR2/AMW*(AMBT*TANB+AMTP*COTB))**2 | |
74 | BE2=(G/2/SR2/AMW*(AMBT*TANB-AMTP*COTB))**2 | |
75 | IF (AMTP.GT.(AMBT+AMHC)) THEN | |
76 | GTBH=AMTP/16./PI*((AL2+BE2) | |
77 | $ *(1.+AMBT**2/AMTP**2-AMHC**2/AMTP**2) | |
78 | $ +2*(AL2-BE2)*AMBT/AMTP) | |
79 | $ *SQRT(SSXLAM(1.,AMHC**2/AMTP**2,AMBT**2/AMTP**2)) | |
80 | CALL SSSAVE(IDTP,GTBH,ISHC,IDBT,0,0,0) | |
81 | END IF | |
82 | C | |
83 | C t->t_1 + z_i decays | |
84 | DO 100 IZ=1,4 | |
85 | MZIZ=ABS(AMZISS(IZ)) | |
86 | SNZI=SIGN(1.,AMZISS(IZ)) | |
87 | IF (SNZI.EQ.1.) THEN | |
88 | THIZ=0. | |
89 | ELSE | |
90 | THIZ=1. | |
91 | END IF | |
92 | ZAUIZ=ZI**(THIZ-1.)*SNZI* | |
93 | $(-G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ)) | |
94 | ZBUIZ=ZI**(THIZ-1.)*4*GP*ZMIXSS(4,IZ)/3./SR2 | |
95 | ZPP=ZI**THIZ | |
96 | ZPM=(-ZI)**THIZ | |
97 | ZA=((ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*COST- | |
98 | $(ZI*ZBUIZ-ZPM*FT*ZMIXSS(1,IZ))*SINT)/2. | |
99 | ZB=((-ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*COST- | |
100 | $(ZI*ZBUIZ+ZPM*FT*ZMIXSS(1,IZ))*SINT)/2. | |
101 | AS=ZA*CONJG(ZA) | |
102 | BS=ZB*CONJG(ZB) | |
103 | IF (AMTP.GT.(AMT1SS+MZIZ)) THEN | |
104 | WID=(AS*((AMTP+MZIZ)**2-AMT1SS**2)+BS* | |
105 | $((AMTP-MZIZ)**2-AMT1SS**2))/16./PI/AMTP* | |
106 | $SQRT(SSXLAM(1.,AMT1SS**2/AMTP**2,MZIZ**2/AMTP**2)) | |
107 | CALL SSSAVE(IDTP,WID,ISZIZ(IZ),ISTP1,0,0,0) | |
108 | END IF | |
109 | 100 CONTINUE | |
110 | C | |
111 | C t -> sb_1 + sW_i | |
112 | C | |
113 | AMW1=ABS(AMW1SS) | |
114 | AMW2=ABS(AMW2SS) | |
115 | IF (AMTP.GT.(AMB1SS+AMW1)) THEN | |
116 | SNWI=SIGN(1.,AMW1SS) | |
117 | AWI=-G*SIN(GAMMAL)*COSB+FB*COS(GAMMAL)*SINB | |
118 | BWI=-FT*(-SNWI)*COS(GAMMAR) | |
119 | WID=AMTP*((AWI**2+BWI**2*COSB**2)*(1.+AMW1**2/AMTP**2 | |
120 | $-AMB1SS**2/AMTP**2)+4*AMW1/AMTP*AWI*BWI*COST)/32./PI* | |
121 | $SQRT(SSXLAM(1.,AMW1**2/AMTP**2,AMB1SS**2/AMTP**2)) | |
122 | CALL SSSAVE(IDTP,WID,ISW1,ISBT1,0,0,0) | |
123 | END IF | |
124 | c | |
125 | IF (AMTP.GT.(AMB1SS+AMW2)) THEN | |
126 | SNWI=SIGN(1.,AMW2SS) | |
127 | AWI=-G*THX*COS(GAMMAL)*COSB-FB*THX*SIN(GAMMAL)*SINB | |
128 | BWI=FT*(-SNWI)*THY*SIN(GAMMAR) | |
129 | WID=AMTP*((AWI**2+BWI**2*COSB**2)*(1.+AMW2**2/AMTP**2 | |
130 | $-AMB1SS**2/AMTP**2)+4*AMW2/AMTP*AWI*BWI*COST)/32./PI* | |
131 | $SQRT(SSXLAM(1.,AMW2**2/AMTP**2,AMB1SS**2/AMTP**2)) | |
132 | CALL SSSAVE(IDTP,WID,ISW2,ISBT1,0,0,0) | |
133 | END IF | |
134 | C | |
135 | C t -> sb_2 + sW_i | |
136 | C | |
137 | IF (AMTP.GT.(AMB2SS+AMW1)) THEN | |
138 | SNWI=SIGN(1.,AMW1SS) | |
139 | AWI=-G*SIN(GAMMAL)*SINB-FB*COS(GAMMAL)*COSB | |
140 | BWI=-FT*(-SNWI)*COS(GAMMAR) | |
141 | WID=AMTP*((AWI**2+BWI**2*SINB**2)*(1.+AMW1**2/AMTP**2 | |
142 | $-AMB2SS**2/AMTP**2)+4*AMW1/AMTP*AWI*BWI*COST)/32./PI* | |
143 | $SQRT(SSXLAM(1.,AMW1**2/AMTP**2,AMB2SS**2/AMTP**2)) | |
144 | CALL SSSAVE(IDTP,WID,ISW1,ISBT2,0,0,0) | |
145 | END IF | |
146 | c | |
147 | IF (AMTP.GT.(AMB2SS+AMW2)) THEN | |
148 | SNWI=SIGN(1.,AMW2SS) | |
149 | AWI=-G*THX*COS(GAMMAL)*SINB+FB*THX*SIN(GAMMAL)*COSB | |
150 | BWI=FT*(-SNWI)*THY*SIN(GAMMAR) | |
151 | WID=AMTP*((AWI**2+BWI**2*SINB**2)*(1.+AMW2**2/AMTP**2 | |
152 | $-AMB2SS**2/AMTP**2)+4*AMW2/AMTP*AWI*BWI*COST)/32./PI* | |
153 | $SQRT(SSXLAM(1.,AMW2**2/AMTP**2,AMB2SS**2/AMTP**2)) | |
154 | CALL SSSAVE(IDTP,WID,ISW2,ISBT2,0,0,0) | |
155 | END IF | |
156 | C | |
157 | C | |
158 | C Normalize branching ratios | |
159 | C | |
160 | CALL SSNORM(IDTP) | |
161 | C | |
162 | RETURN | |
163 | END |