source: trunk/MagicSoft/Simulation/Corsika/Mmcs/gbrsgm.f

Last change on this file was 286, checked in by harald, 25 years ago
This is the start point for further developments of the Magic Monte Carlo Simulation written by Jose Carlos Gonzales. Now it is under control of one CVS repository for the whole collaboration. Everyone should use this CVS repository for further developments.
File size: 4.0 KB
Line 
1 REAL FUNCTION GBRSGM(Z,E)
2
3C-----------------------------------------------------------------------
4C CALCULATES MUON BREMSSTRAHLUNG CROSS SECTIONS
5C
6C THIS SUBROUTINE IS TAKEN FROM GEANT321 PACKAGE (WITH MODIFICATIONS)
7C CALCULATES CROSS-SECTION IN CURRENT MATERIAL FOR DISCRETE(HARD) MUON
8C BREMSSTRAHLUNG. (SIG IN BARN/ATOM)
9C FOR A DESCRIPTION SEE: CERN PROGRAM LIBRARY LONG WRITEUP W5013 (1993)
10C THIS FUNCTION IS CALLED FROM BOX2
11C ARGUMENTS:
12C Z (R4) = ATOMIC NUMBER OF PENETRATET MATERIAL
13C E (R4) = TOTAL ENERGY OF MUON
14C
15C AUTHOR : L.URBAN
16C MODIFIED: D. HECK IK3 FZK KARLSRUHE
17C-----------------------------------------------------------------------
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/
70C-----------------------------------------------------------------------
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
Note: See TracBrowser for help on using the repository browser.