]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | SUBROUTINE SSHCC | |
3 | C----------------------------------------------------------------------- | |
4 | C Calculates the decay widths of all neutral Higgses into all | |
5 | C possible pairs of charginos. | |
6 | C | |
7 | C Bisset's CHGINO | |
8 | C----------------------------------------------------------------------- | |
9 | #if defined(CERNLIB_IMPNONE) | |
10 | IMPLICIT NONE | |
11 | #endif | |
12 | #include "isajet/sspar.inc" | |
13 | #include "isajet/sssm.inc" | |
14 | #include "isajet/sstype.inc" | |
15 | C | |
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 | |
18 | $,DWID,TEMP,TEMPXY | |
19 | DOUBLE PRECISION MHI(3),IDHI(3),SHP(3),SHM(3),SH(3),PH(3) | |
20 | DOUBLE PRECISION SSDLAM | |
21 | REAL WID | |
22 | INTEGER NUMH,IDHHA | |
23 | C | |
24 | PI=4.*ATAN(1.D0) | |
25 | SR2=SQRT(2.D0) | |
26 | XM=1./TAN(GAMMAL) | |
27 | THETX=SIGN(1.D0,XM) | |
28 | YM=1./TAN(GAMMAR) | |
29 | THETY=SIGN(1.D0,YM) | |
30 | SGL=1/(DSQRT(1+XM**2)) | |
31 | CGL=SGL*XM | |
32 | SGR=1/(DSQRT(1+YM**2)) | |
33 | CGR=SGR*YM | |
34 | MW1=DBLE(ABS(AMW1SS)) | |
35 | MW2=DBLE(ABS(AMW2SS)) | |
36 | THETM=SIGN(1.,AMW1SS) | |
37 | THETP=SIGN(1.,AMW2SS) | |
38 | G2=4*PI*ALFAEM/SN2THW | |
39 | GP2=4*PI*ALFAEM/(1-SN2THW) | |
40 | BETA=ATAN(1.0/RV2V1) | |
41 | ALPHA=ALFAH | |
42 | C | |
43 | C Arrays for loops | |
44 | C | |
45 | MHI(1)=AMHL | |
46 | MHI(2)=AMHH | |
47 | MHI(3)=AMHA | |
48 | IDHI(1)=ISHL | |
49 | IDHI(2)=ISHH | |
50 | IDHI(3)=ISHA | |
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 | |
55 | SHP(1)=SHP(1)*TEMPXY | |
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 | |
59 | PH(1)=-SH(1) | |
60 | T1=THETX*CGL*CGR*COS(ALPHA)*THETM | |
61 | SH(1)=SH(1)+T1 | |
62 | PH(1)=PH(1)-T1 | |
63 | T1=THETY*SGL*SGR*COS(ALPHA)*THETP | |
64 | SH(1)=SH(1)-T1 | |
65 | PH(1)=PH(1)-T1 | |
66 | T1=THETY*CGL*CGR*SIN(ALPHA)*THETP | |
67 | SH(1)=SH(1)+T1 | |
68 | PH(1)=PH(1)+T1 | |
69 | SH(1)=0.5*SH(1) | |
70 | PH(1)=0.5*PH(1) | |
71 | SHP(2)=COS(ALPHA)*CGR*SGL-SIN(ALPHA)*CGL*SGR | |
72 | SHP(2)=SHP(2)*TEMPXY | |
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 | |
76 | PH(2)=-SH(2) | |
77 | T1=THETX*CGL*CGR*SIN(ALPHA)*THETM | |
78 | SH(2)=SH(2)-T1 | |
79 | PH(2)=PH(2)+T1 | |
80 | T1=THETY*SGL*SGR*SIN(ALPHA)*THETP | |
81 | SH(2)=SH(2)+T1 | |
82 | PH(2)=PH(2)+T1 | |
83 | T1=THETY*CGL*CGR*COS(ALPHA)*THETP | |
84 | SH(2)=SH(2)+T1 | |
85 | PH(2)=PH(2)+T1 | |
86 | SH(2)=0.5*SH(2) | |
87 | PH(2)=0.5*PH(2) | |
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 | |
93 | PH(3)=-SH(3) | |
94 | T1=THETX*CGL*CGR*COS(BETA)*THETM | |
95 | SH(3)=SH(3)+T1 | |
96 | PH(3)=PH(3)-T1 | |
97 | T1=THETY*SGL*SGR*COS(BETA)*THETP | |
98 | SH(3)=SH(3)+T1 | |
99 | PH(3)=PH(3)+T1 | |
100 | T1=THETY*CGL*CGR*SIN(BETA)*THETP | |
101 | SH(3)=SH(3)-T1 | |
102 | PH(3)=PH(3)-T1 | |
103 | SH(3)=0.5*SH(3) | |
104 | PH(3)=0.5*PH(3) | |
105 | C | |
106 | C Loop over neutral Higgs | |
107 | C | |
108 | DO 100 NUMH=1,3 | |
109 | MH=MHI(NUMH) | |
110 | IDHHA=IDHI(NUMH) | |
111 | C w1 + w1 | |
112 | M1=ABS(AMW1SS) | |
113 | M2=M1 | |
114 | IF(MH.GT.M1+M2) THEN | |
115 | LAMB=SSDLAM(MH**2,M1**2,M2**2) | |
116 | TEMP=1-4*M1**2/MH**2 | |
117 | DWID=G2*MH*SHM(NUMH)**2/(4.0*PI) | |
118 | DWID=DWID*SQRT(TEMP**3) | |
119 | WID=DWID | |
120 | CALL SSSAVE(IDHHA,WID,ISW1,-ISW1,0,0,0) | |
121 | ENDIF | |
122 | C w2 + w2 | |
123 | M1=ABS(AMW2SS) | |
124 | M2=M1 | |
125 | IF(MH.GT.M1+M2) THEN | |
126 | TEMP=1-4*M1**2/MH**2 | |
127 | DWID=G2*MH*SHP(NUMH)**2/(4*PI) | |
128 | DWID=DWID*SQRT(TEMP**3) | |
129 | WID=DWID | |
130 | CALL SSSAVE(IDHHA,WID,ISW2,-ISW2,0,0,0) | |
131 | ENDIF | |
132 | C w1 + w2 | |
133 | M1=ABS(AMW1SS) | |
134 | M2=ABS(AMW2SS) | |
135 | IF(MH.GT.M1+M2) THEN | |
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)) | |
140 | WID=DWID | |
141 | CALL SSSAVE(IDHHA,WID,ISW1,-ISW2,0,0,0) | |
142 | CALL SSSAVE(IDHHA,WID,-ISW1,ISW2,0,0,0) | |
143 | ENDIF | |
144 | 100 CONTINUE | |
145 | C | |
146 | RETURN | |
147 | END |