]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.2 1996/09/30 13:30:10 ravndal | |
6 | * Bremsstrahlung extension for higher energies | |
7 | * | |
8 | * Revision 1.1.1.1 1995/10/24 10:21:22 cernlib | |
9 | * Geant | |
10 | * | |
11 | * | |
12 | #include "geant321/pilot.h" | |
13 | *CMZ : 3.21/04 24/02/95 14.50.31 by S.Giani | |
14 | *-- Author : | |
15 | FUNCTION GBRSGE(ZZ,T,BCUT) | |
16 | C. | |
17 | C. ****************************************************************** | |
18 | C. * * | |
19 | C. * Calculates cross-section in current material * | |
20 | C. * for discrete(hard) electron BREMSSTRAHLUNG. * | |
21 | C. * (SIG in barn/atom) * | |
22 | C. * * | |
23 | C. * ==>Called by : GBRSGA * | |
24 | C. * Author L.Urban ********* * | |
25 | C. * * | |
26 | C. * correction for T> 100 GeV ! (by L.Urban on 23/09/96) * | |
27 | C. ****************************************************************** | |
28 | C. | |
29 | #include "geant321/gconsp.inc" | |
30 | #include "geant321/gcmate.inc" | |
31 | DIMENSION C(100),C1(60),C2(40) | |
32 | EQUIVALENCE (C(1),C1(1)),(C(61),C2(1)) | |
33 | SAVE C | |
34 | #if !defined(CERNLIB_BETHE) | |
35 | PARAMETER (AKSI=1.80,ALFA=0.98,VS=0.0001) | |
36 | DATA C1/ 0.430748E-02, 0.576058E-02,-0.122564E-02, 0.114843E-03 | |
37 | + ,-0.489452E-05, 0.795991E-07, 0.326746E-02,-0.132872E-02 | |
38 | + , 0.217197E-03,-0.179769E-04, 0.766114E-06,-0.125603E-07 | |
39 | + , 0.326452E-02,-0.175331E-02, 0.415488E-03,-0.507652E-04 | |
40 | + , 0.297569E-05,-0.651741E-07, 0.847189E-03,-0.433923E-03 | |
41 | + , 0.116672E-03,-0.166799E-04, 0.110237E-05,-0.263383E-07 | |
42 | + , 0.846052E-04,-0.415764E-04, 0.129610E-04,-0.212844E-05 | |
43 | + , 0.152871E-06,-0.384393E-08, 0.300838E-05,-0.136833E-05 | |
44 | + , 0.507296E-06,-0.943623E-07, 0.720305E-08,-0.187210E-09 | |
45 | + , 0.448230E-01,-0.210048E-01, 0.379434E-02,-0.328431E-03 | |
46 | + , 0.136710E-04,-0.220593E-06,-0.539248E-02, 0.330244E-02 | |
47 | + ,-0.733726E-03, 0.732312E-04,-0.336810E-05, 0.583913E-07 | |
48 | + ,-0.106983E-02, 0.378021E-03,-0.384854E-04, 0.978156E-06 | |
49 | + , 0.410622E-07,-0.174250E-08,-0.117501E-04,-0.983887E-05 | |
50 | + , 0.239644E-05,-0.190104E-06, 0.619226E-08,-0.680932E-10/ | |
51 | DATA C2/ 0.168074E-03,-0.934609E-04, 0.141293E-04,-0.854216E-06 | |
52 | + , 0.183287E-07, 0.932144E-04,-0.234926E-04, 0.136656E-05 | |
53 | + , 0.351109E-07,-0.330189E-08, 0.174523E-04, 0.253854E-05 | |
54 | + ,-0.171643E-05, 0.183074E-06,-0.566331E-08, 0.111970E-05 | |
55 | + , 0.112776E-05,-0.386924E-06, 0.367597E-07,-0.108504E-08 | |
56 | + , 0.171604E-07, 0.738801E-07,-0.218761E-07, 0.199032E-08 | |
57 | + ,-0.576173E-10,-0.105531E-03, 0.362995E-04,-0.433334E-05 | |
58 | + , 0.207664E-06,-0.330250E-08,-0.168293E-05,-0.773204E-06 | |
59 | + , 0.227974E-06,-0.159385E-07, 0.321958E-09, 0.167046E-05 | |
60 | + ,-0.440761E-06, 0.396377E-07,-0.151053E-08, 0.215624E-10/ | |
61 | #endif | |
62 | #if defined(CERNLIB_BETHE) | |
63 | PARAMETER (AKSI=1.80,ALFA=1.00,VS=0.0001) | |
64 | DATA C1/ 0.111394E-01, 0.138592E-02,-0.274910E-03, 0.198389E-04 | |
65 | + ,-0.472291E-06,-0.668760E-09, 0.839876E-02,-0.449619E-02 | |
66 | + , 0.917722E-03,-0.874992E-04, 0.390015E-05,-0.653499E-07 | |
67 | + , 0.464718E-02,-0.260321E-02, 0.582582E-03,-0.631769E-04 | |
68 | + , 0.325801E-05,-0.638109E-07, 0.101243E-02,-0.532089E-03 | |
69 | + , 0.129100E-03,-0.161260E-04, 0.958116E-06,-0.212606E-07 | |
70 | + , 0.926006E-04,-0.455288E-04, 0.124755E-04,-0.182510E-05 | |
71 | + , 0.122078E-06,-0.294182E-08, 0.308111E-05,-0.134712E-05 | |
72 | + , 0.436230E-06,-0.751770E-07, 0.552327E-08,-0.140541E-09 | |
73 | + , 0.292552E-02,-0.830719E-03, 0.210705E-04, 0.103750E-04 | |
74 | + ,-0.953318E-06, 0.236453E-07, 0.110907E-02,-0.219463E-03 | |
75 | + , 0.128517E-04,-0.554575E-06, 0.507378E-07,-0.182214E-08 | |
76 | + ,-0.639866E-03, 0.209918E-03,-0.250183E-04, 0.138030E-05 | |
77 | + ,-0.358845E-07, 0.366305E-09, 0.300095E-04,-0.133668E-04 | |
78 | + , 0.198606E-05,-0.133100E-06, 0.414062E-08,-0.485929E-10/ | |
79 | DATA C2/ 0.938677E-04,-0.613470E-04, 0.981984E-05,-0.646289E-06 | |
80 | + , 0.150731E-07, 0.331764E-04, 0.447248E-05,-0.313021E-05 | |
81 | + , 0.329777E-06,-0.997210E-08, 0.876484E-08, 0.110042E-04 | |
82 | + ,-0.313406E-05, 0.280040E-06,-0.794116E-08,-0.103097E-05 | |
83 | + , 0.220550E-05,-0.569158E-06, 0.491359E-07,-0.136987E-08 | |
84 | + ,-0.781963E-07, 0.123162E-06,-0.302825E-07, 0.255485E-08 | |
85 | + ,-0.702300E-10,-0.193213E-03, 0.640349E-04,-0.720586E-05 | |
86 | + , 0.329673E-06,-0.514682E-08, 0.362138E-04,-0.142783E-04 | |
87 | + , 0.184063E-05,-0.967983E-07, 0.179386E-08, 0.496130E-06 | |
88 | + , 0.210105E-06,-0.538512E-07, 0.372388E-08,-0.812734E-10/ | |
89 | #endif | |
90 | C. | |
91 | DATA THIGH,CHIGH/100.,50./ | |
92 | C. ------------------------------------------------------------------ | |
93 | C. | |
94 | GBRSGE=0. | |
95 | * | |
96 | IF(BCUT.LE.0.) GOTO 99 | |
97 | IF(BCUT.GE.T) GOTO 99 | |
98 | * | |
99 | IF(T.GT.THIGH) THEN | |
100 | TT=THIGH | |
101 | IF(BCUT.GE.THIGH) THEN | |
102 | EC=CHIGH | |
103 | ELSE | |
104 | EC=BCUT | |
105 | ENDIF | |
106 | ELSE | |
107 | TT=T | |
108 | EC=BCUT | |
109 | ENDIF | |
110 | * | |
111 | E=TT+EMASS | |
112 | * | |
113 | C | |
114 | X=LOG(E/EMASS) | |
115 | Y=LOG(E*VS/EC) | |
116 | C | |
117 | S=0. | |
118 | YY=1. | |
119 | DO 30 I=1,2 | |
120 | XX=1. | |
121 | DO 20 J=1,6 | |
122 | K=6*I+J-6 | |
123 | S=S+C(K)*XX*YY | |
124 | XX=XX*X | |
125 | 20 CONTINUE | |
126 | YY=YY*Y | |
127 | 30 CONTINUE | |
128 | DO 50 I=3,6 | |
129 | XX=1. | |
130 | DO 40 J=1,6 | |
131 | K=6*I+J-6 | |
132 | IF(Y.LE.0.) THEN | |
133 | S=S+C(K)*XX*YY | |
134 | ELSE | |
135 | S=S+C(K+24)*XX*YY | |
136 | ENDIF | |
137 | XX=XX*X | |
138 | 40 CONTINUE | |
139 | YY=YY*Y | |
140 | 50 CONTINUE | |
141 | SS=0. | |
142 | YY=1. | |
143 | DO 70 I=1,2 | |
144 | XX=1. | |
145 | DO 60 J=1,5 | |
146 | K=5*I+J+55 | |
147 | SS=SS+C(K)*XX*YY | |
148 | XX=XX*X | |
149 | 60 CONTINUE | |
150 | YY=YY*Y | |
151 | 70 CONTINUE | |
152 | DO 90 I=3,5 | |
153 | XX=1. | |
154 | DO 80 J=1,5 | |
155 | K=5*I+J+55 | |
156 | IF(Y.LE.0.) THEN | |
157 | SS=SS+C(K)*XX*YY | |
158 | ELSE | |
159 | SS=SS+C(K+15)*XX*YY | |
160 | ENDIF | |
161 | XX=XX*X | |
162 | 80 CONTINUE | |
163 | YY=YY*Y | |
164 | 90 CONTINUE | |
165 | C | |
166 | S=S+Z*SS | |
167 | IF(S.LE.0.) GOTO 99 | |
168 | C | |
169 | FAC=ZZ*(ZZ+AKSI)*E*E*(LOG(TT/EC))**ALFA/(TT*(E+EMASS)) | |
170 | IF(FAC.LE.0.) GOTO 99 | |
171 | C | |
172 | GBRSGE=FAC*S | |
173 | * | |
174 | IF(T.GT.THIGH) THEN | |
175 | ||
176 | IF(BCUT.LT.THIGH) THEN | |
177 | RAT=BCUT/T | |
178 | S=(-LOG(RAT)-2./3.+RAT-RAT*RAT/3.) | |
179 | RAT=BCUT/TT | |
180 | S=S/(-LOG(RAT)-2./3.+RAT-RAT*RAT/3.) | |
181 | ELSE | |
182 | RAT=BCUT/T | |
183 | S=(-LOG(RAT)-2./3.+RAT-RAT*RAT/3.) | |
184 | RAT=EC/TT | |
185 | S=S/(-LOG(RAT)-2./3.+RAT-RAT*RAT/3.) | |
186 | ENDIF | |
187 | ||
188 | GBRSGE=GBRSGE*S | |
189 | ENDIF | |
190 | * | |
191 | 99 RETURN | |
192 | * | |
193 | END |