1 | REAL FUNCTION GBRSGM(Z,E)
|
---|
2 |
|
---|
3 | C-----------------------------------------------------------------------
|
---|
4 | C CALCULATES MUON BREMSSTRAHLUNG CROSS SECTIONS
|
---|
5 | C
|
---|
6 | C THIS SUBROUTINE IS TAKEN FROM GEANT321 PACKAGE (WITH MODIFICATIONS)
|
---|
7 | C CALCULATES CROSS-SECTION IN CURRENT MATERIAL FOR DISCRETE(HARD) MUON
|
---|
8 | C BREMSSTRAHLUNG. (SIG IN BARN/ATOM)
|
---|
9 | C FOR A DESCRIPTION SEE: CERN PROGRAM LIBRARY LONG WRITEUP W5013 (1993)
|
---|
10 | C THIS FUNCTION IS CALLED FROM BOX2
|
---|
11 | C ARGUMENTS:
|
---|
12 | C Z (R4) = ATOMIC NUMBER OF PENETRATET MATERIAL
|
---|
13 | C E (R4) = TOTAL ENERGY OF MUON
|
---|
14 | C
|
---|
15 | C AUTHOR : L.URBAN
|
---|
16 | C MODIFIED: D. HECK IK3 FZK KARLSRUHE
|
---|
17 | C-----------------------------------------------------------------------
|
---|
18 |
|
---|
19 | IMPLICIT NONE
|
---|
20 | *KEEP,MUPART.
|
---|
21 | COMMON /MUPART/ AMUPAR,BCUT,CMUON,FMUBRM,FMUORG
|
---|
22 | DOUBLE PRECISION AMUPAR(14),BCUT,CMUON(11)
|
---|
23 | LOGICAL FMUBRM,FMUORG
|
---|
24 | *KEEP,PAM.
|
---|
25 | COMMON /PAM/ PAMA,SIGNUM
|
---|
26 | DOUBLE PRECISION PAMA(6000),SIGNUM(6000)
|
---|
27 | *KEEP,RUNPAR.
|
---|
28 | COMMON /RUNPAR/ FIXHEI,THICK0,HILOECM,HILOELB,
|
---|
29 | * STEPFC,NRRUN,NSHOW,PATAPE,MONIIN,
|
---|
30 | * MONIOU,MDEBUG,NUCNUC,
|
---|
31 | * CETAPE,
|
---|
32 | * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
|
---|
33 | * N1STTR,MDBASE,
|
---|
34 | * DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
|
---|
35 | * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
|
---|
36 | * ,GHEISH,GHESIG
|
---|
37 | COMMON /RUNPAC/ DSN,HOST,USER
|
---|
38 | DOUBLE PRECISION FIXHEI,THICK0,HILOECM,HILOELB
|
---|
39 | REAL STEPFC
|
---|
40 | INTEGER NRRUN,NSHOW,PATAPE,MONIIN,MONIOU,MDEBUG,NUCNUC,
|
---|
41 | * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
|
---|
42 | * N1STTR,MDBASE
|
---|
43 | INTEGER CETAPE
|
---|
44 | CHARACTER*79 DSN
|
---|
45 | CHARACTER*20 HOST,USER
|
---|
46 |
|
---|
47 | LOGICAL DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
|
---|
48 | * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
|
---|
49 | * ,GHEISH,GHESIG
|
---|
50 | *KEND.
|
---|
51 |
|
---|
52 | REAL C(52),AKSI,ALFA,E,ECMAX,FAC,GAM,
|
---|
53 | * S,SS,X,XX,Y,YY,Z
|
---|
54 | INTEGER I,J,K
|
---|
55 | SAVE C
|
---|
56 | DATA AKSI/2.30/, ALFA/1.06/, GAM/0.63/
|
---|
57 | DATA C/ 0.949313E-07,-0.819600E-07, 0.529075E-07,-0.832023E-08
|
---|
58 | + , 0.539299E-09,-0.127042E-10,-0.165784E-08,-0.307788E-07
|
---|
59 | + , 0.977905E-08,-0.113658E-08, 0.574481E-10,-0.106221E-11
|
---|
60 | + , 0.968339E-09,-0.108640E-08,-0.177634E-09, 0.889497E-10
|
---|
61 | + ,-0.876878E-11, 0.264303E-12, 0.216263E-08,-0.152680E-08
|
---|
62 | + , 0.380989E-09,-0.455274E-10, 0.264172E-11,-0.596016E-13
|
---|
63 | + , 0.444927E-09,-0.272978E-09, 0.645634E-10,-0.748783E-11
|
---|
64 | + , 0.424890E-12,-0.940837E-14, 0.162289E-10,-0.362486E-11
|
---|
65 | + ,-0.576652E-12, 0.211269E-12,-0.185482E-13, 0.522065E-15
|
---|
66 | + ,-0.215590E-09, 0.112204E-09,-0.819133E-11, 0.145128E-12
|
---|
67 | + ,-0.206029E-09, 0.559940E-10,-0.483350E-11, 0.134252E-12
|
---|
68 | + ,-0.368469E-10, 0.999457E-11,-0.904967E-12, 0.272717E-13
|
---|
69 | + ,-0.303446E-11, 0.853429E-12,-0.785466E-13, 0.236435E-14/
|
---|
70 | C-----------------------------------------------------------------------
|
---|
71 |
|
---|
72 | GBRSGM = 0.
|
---|
73 | IF ( E-PAMA(5) .LE. BCUT ) RETURN
|
---|
74 | ECMAX = E - CMUON(10) * Z**0.333333
|
---|
75 | IF ( ECMAX .LE. BCUT ) RETURN
|
---|
76 | X = LOG(E/PAMA(5))
|
---|
77 |
|
---|
78 | S = 0.
|
---|
79 | YY = 1.
|
---|
80 | DO 30 I = 1,6
|
---|
81 | XX = 1.
|
---|
82 | DO 20 J = 1,6
|
---|
83 | K = 6*I + J - 6
|
---|
84 | S = S + C(K) * XX * YY
|
---|
85 | XX = XX * X
|
---|
86 | 20 CONTINUE
|
---|
87 | YY = YY * CMUON(11)
|
---|
88 | 30 CONTINUE
|
---|
89 | SS = 0.
|
---|
90 | YY = 1.
|
---|
91 | DO 50 I = 1,4
|
---|
92 | XX = 1.
|
---|
93 | DO 40 J = 1,4
|
---|
94 | K = 4*I + J + 32
|
---|
95 | SS = SS + C(K) * XX * YY
|
---|
96 | XX = XX * X
|
---|
97 | 40 CONTINUE
|
---|
98 | YY = YY * CMUON(11)
|
---|
99 | 50 CONTINUE
|
---|
100 | S = S + Z * SS
|
---|
101 | IF ( S .LE. 0. ) RETURN
|
---|
102 |
|
---|
103 | FAC = LOG(ECMAX/BCUT)
|
---|
104 | IF ( FAC .LE. 0. ) RETURN
|
---|
105 | FAC = Z * ( Z + AKSI * (1.+GAM*LOG(Z)) ) * FAC**ALFA
|
---|
106 | GBRSGM = FAC * S
|
---|
107 |
|
---|
108 | * IF ( DEBUG ) WRITE(MDEBUG,444) Z,E,GBRSGM
|
---|
109 | * 444 FORMAT(' GBRSGM: Z=',F3.0,' E=',1P,E10.4,' GBRSGM=',E10.4)
|
---|
110 |
|
---|
111 | 99 RETURN
|
---|
112 | END
|
---|