1 #include "isajet/pilot.h"
3 C-----------------------------------------------------------------------
4 C Calculates the decay widths of all neutral Higgses into all
5 C possible pairs of charginos.
8 C-----------------------------------------------------------------------
9 #if defined(CERNLIB_IMPNONE)
12 #include "isajet/sspar.inc"
13 #include "isajet/sssm.inc"
14 #include "isajet/sstype.inc"
16 DOUBLE PRECISION PI,SR2,XM,THETX,YM,THETY,SGL,CGL,SGR,CGR
17 $,MW1,MW2,THETM,THETP,G2,GP2,BETA,ALPHA,T1,MH,M1,M2,LAMB
19 DOUBLE PRECISION MHI(3),IDHI(3),SHP(3),SHM(3),SH(3),PH(3)
20 DOUBLE PRECISION SSDLAM
30 SGL=1/(DSQRT(1+XM**2))
32 SGR=1/(DSQRT(1+YM**2))
39 GP2=4*PI*ALFAEM/(1-SN2THW)
51 C The following came from Bisset's MASZIN, but with L,H,P
52 C replaced by a generic H and a subscript.
53 TEMPXY=0.5*THETX*THETY*(-THETP)
54 SHP(1)=SIN(ALPHA)*CGR*SGL+COS(ALPHA)*CGL*SGR
56 SHM(1)=SIN(ALPHA)*SGR*CGL+COS(ALPHA)*SGL*CGR
57 SHM(1)=SHM(1)*0.5*THETM
58 SH(1)=-THETX*SGR*SGL*SIN(ALPHA)*THETM
60 T1=THETX*CGL*CGR*COS(ALPHA)*THETM
63 T1=THETY*SGL*SGR*COS(ALPHA)*THETP
66 T1=THETY*CGL*CGR*SIN(ALPHA)*THETP
71 SHP(2)=COS(ALPHA)*CGR*SGL-SIN(ALPHA)*CGL*SGR
73 SHM(2)=COS(ALPHA)*SGR*CGL-SIN(ALPHA)*SGL*CGR
74 SHM(2)=SHM(2)*0.5*THETM
75 SH(2)=-THETX*SGR*SGL*COS(ALPHA)*THETM
77 T1=THETX*CGL*CGR*SIN(ALPHA)*THETM
80 T1=THETY*SGL*SGR*SIN(ALPHA)*THETP
83 T1=THETY*CGL*CGR*COS(ALPHA)*THETP
88 SHP(3)=SIN(BETA)*CGR*SGL+COS(BETA)*CGL*SGR
89 SHP(3)=SHP(3)*0.5*THETX*THETY*(-THETP)
90 SHM(3)=SIN(BETA)*SGR*CGL+COS(BETA)*SGL*CGR
91 SHM(3)=SHM(3)*0.5*THETM
92 SH(3)=-THETX*SGR*SGL*SIN(BETA)*THETM
94 T1=THETX*CGL*CGR*COS(BETA)*THETM
97 T1=THETY*SGL*SGR*COS(BETA)*THETP
100 T1=THETY*CGL*CGR*SIN(BETA)*THETP
106 C Loop over neutral Higgs
115 LAMB=SSDLAM(MH**2,M1**2,M2**2)
117 DWID=G2*MH*SHM(NUMH)**2/(4.0*PI)
118 DWID=DWID*SQRT(TEMP**3)
120 CALL SSSAVE(IDHHA,WID,ISW1,-ISW1,0,0,0)
127 DWID=G2*MH*SHP(NUMH)**2/(4*PI)
128 DWID=DWID*SQRT(TEMP**3)
130 CALL SSSAVE(IDHHA,WID,ISW2,-ISW2,0,0,0)
136 LAMB=SSDLAM(MH**2,M1**2,M2**2)
137 DWID=PH(NUMH)**2*(MH**2-(M1-M2)**2)
138 DWID=DWID+SH(NUMH)**2*(MH**2-(M1+M2)**2)
139 DWID=DWID*G2*SQRT(LAMB)/(16.0*PI*(MH**3))
141 CALL SSSAVE(IDHHA,WID,ISW1,-ISW2,0,0,0)
142 CALL SSSAVE(IDHHA,WID,-ISW1,ISW2,0,0,0)