1 #include "isajet/pilot.h"
3 C-----------------------------------------------------------------------
5 C Calculate the top branching ratios.
6 C Source: H. Baer (modified by F. Paige)
8 C-----------------------------------------------------------------------
9 #if defined(CERNLIB_IMPNONE)
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"
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
25 DATA ZONE/(1.,0.)/,ZI/(0.,1.)/
29 G=SQRT(4*PI*ALFAEM/SN2THW)
30 GP=G*SQRT(SN2THW/(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))*
41 MTQ=SSMQCD(DBLE(MTMT),DBLE(AMTP))
42 FB=G*MBQ/SR2/AMW/COS(BETA)
43 FT=G*MTQ/SR2/AMW/SIN(BETA)
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)
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)
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)
86 SNZI=SIGN(1.,AMZISS(IZ))
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
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.
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)
115 IF (AMTP.GT.(AMB1SS+AMW1)) THEN
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)
125 IF (AMTP.GT.(AMB1SS+AMW2)) THEN
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)
137 IF (AMTP.GT.(AMB2SS+AMW1)) THEN
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)
147 IF (AMTP.GT.(AMB2SS+AMW2)) THEN
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)
158 C Normalize branching ratios