- Page 1: THE EGS5 CODE SYSTEM 1Hideo Hirayam
- Page 7 and 8: B.4.9 Output of Results (Step 9) .
- Page 9 and 10: List of Figures2.1 Program flow and
- Page 11 and 12: C.4 Subprogram relationships in PEG
- Page 13 and 14: B.5 Variable descriptions for COMMO
- Page 15 and 16: PREFACEIn the nineteen years since
- Page 18 and 19: of-print and this history might hav
- Page 20 and 21: constructed SHOWER3 in modular form
- Page 22 and 23: It soon became clear that, in the t
- Page 24 and 25: subroutine ELECTR.13. The control l
- Page 26 and 27: - Geometry subprograms and correspo
- Page 28 and 29: indicates that a much slower simula
- Page 30 and 31: Although the transport mechanics al
- Page 32 and 33: Photoelectron angular distributionI
- Page 34 and 35: scattering is included (as an optio
- Page 36 and 37: Chapter 2RADIATION TRANSPORT INEGS5
- Page 38 and 39: that the reader is already acquaint
- Page 40 and 41: so thatP r{ ˆx < x } = G(h(x)) = F
- Page 42 and 43: 2.3 Simulating the Physical Process
- Page 44 and 45: 1. Compute λ at the current locati
- Page 46 and 47: 2.6 General Implementation NotesWe
- Page 48 and 49: UserControlDataUSER CODEInformation
- Page 50 and 51: Table 2.1: Symbols used in EGS5 and
- Page 52 and 53:
Table 2.1 (cont.)Math FORTRAN Progr
- Page 54 and 55:
then we see that for k ≫ k c , F
- Page 56 and 57:
where δ = 272 Z −1/3 ∆ as befo
- Page 58 and 59:
For φ rad in the completely screen
- Page 60 and 61:
Now suppose that the f(x) in Equati
- Page 62 and 63:
Also, recall from Table 2.1 the def
- Page 64 and 65:
We proceed now by using Equation 2.
- Page 66 and 67:
Then to sample from f 1j (E), first
- Page 68 and 69:
d˘Σ P airdE=(ZA + Z B − (Z F if
- Page 70 and 71:
The various parameters needed to sa
- Page 72 and 73:
used for medical linacs (6-50 MeV,
- Page 74 and 75:
where Θ ± is the e ± scattering
- Page 76 and 77:
and N g is a normalization factor t
- Page 78 and 79:
are then given byP 1 = (E 1 , 0, 0,
- Page 80 and 81:
Thus˘Σ Compt, T otal (˘k 0 ) =(E
- Page 82 and 83:
2.10 Møller ScatteringThe form of
- Page 84 and 85:
2.11 Bhabha ScatteringThe different
- Page 86 and 87:
integrating over azimuth, and chang
- Page 88 and 89:
Using the value of C 2 and rearrang
- Page 90 and 91:
over the range of energy transfers
- Page 92 and 93:
Table 2.2: Default atomic numbers,
- Page 94 and 95:
Table 2.3: Default Sternheimer dens
- Page 96 and 97:
and x 0 , x 1 , a, m s are paramete
- Page 98 and 99:
As a final comment on continuous en
- Page 100 and 101:
The general formula for the f (i) i
- Page 102 and 103:
Scott’s formula (7.25) for χ α
- Page 104 and 105:
Figure 2.6: Plots of Molière funct
- Page 106 and 107:
Figure 2.7: Plot of Equation 2.290
- Page 108 and 109:
at large angles for low energies an
- Page 110 and 111:
and〈cos 2 Θ〉 GS = 1 3 [1 + 2 e
- Page 112 and 113:
2.15 Transport Mechanics in EGS5As
- Page 114 and 115:
xφFinalDirectionΘ(1−ζ) t∆xt
- Page 116 and 117:
FinalDirectionxφΘt( (1 − ζ) K1
- Page 118 and 119:
FinalDirectionxEnergy HingesφΘt(
- Page 120 and 121:
2.15.5 EGS5 Transport Mechanics Alg
- Page 122 and 123:
tally, as in EGS4. Depositing energ
- Page 124 and 125:
2.15.6 Electron Step-Size Selection
- Page 126 and 127:
where P Eh (t) is the probability t
- Page 128 and 129:
despite being able to take much lon
- Page 130 and 131:
Different scattering anglesL (= D)D
- Page 132 and 133:
Fractional Energy Deposited10.80.60
- Page 134 and 135:
10.10.01Electrons on Ti20MeV10MeV7M
- Page 136 and 137:
particle on its very first track, i
- Page 138 and 139:
50PWLF-LEM40Kβ1Kβ2ExactPWLFGMFP (
- Page 140 and 141:
Equation 2.397. Additionally, the v
- Page 142 and 143:
Table 2.7: L x-ray energies and rep
- Page 144 and 145:
Table 2.8: Total cross section (10
- Page 146 and 147:
are reasonable, but at lower initia
- Page 148 and 149:
Implementation in EGS5 If the user
- Page 150 and 151:
XSke=Ze0Ok0YFigure 2.24: Direction
- Page 152 and 153:
Transformation to the laboratory sy
- Page 154 and 155:
Table 2.9: Formulas used in various
- Page 156 and 157:
Chapter 3A SERIES OF SHORT EGS5TUTO
- Page 158 and 159:
! usersc contains emaxecommon/geom/
- Page 160 and 161:
uin=0.0vin=0.0win=1.0! Moving along
- Page 162 and 163:
integer iargreal*8 angle,ekine! Arg
- Page 164 and 165:
! In regon with source! This must b
- Page 166 and 167:
In addition, this user code include
- Page 168 and 169:
HATCH-call comes nextEGS SUCCESSFUL
- Page 170 and 171:
subroutine ausgab(iarg)implicit non
- Page 172 and 173:
accumulated over the given distance
- Page 174 and 175:
! ====================call counters
- Page 176 and 177:
! Define initial variables for 2 Me
- Page 178 and 179:
total=0.0do i=1,3total=total+escore
- Page 180 and 181:
Following is the output produced by
- Page 182 and 183:
-----------Total fraction of energy
- Page 184 and 185:
include ’include/randomm.f’! bo
- Page 186 and 187:
!----------------------------------
- Page 188 and 189:
end ifend dowrite(6,160) ein*1000.,
- Page 190 and 191:
HATCH-call comes nextRAYLEIGH OPTIO
- Page 192 and 193:
! backwards in region 1!! | |! Regi
- Page 194 and 195:
! Step 1: Initialization!----------
- Page 196 and 197:
iraylr(2)=1! Turn on rayleigh scatt
- Page 198 and 199:
pnorm(1,2) =0.0pnorm(2,2) =0.0pnorm
- Page 200 and 201:
eturnend!--------------------------
- Page 202 and 203:
EGS SUCCESSFULLY ’HATCHED’ FOR
- Page 204 and 205:
Chapter 4ADVANCED EGS5 USER CODESIn
- Page 206 and 207:
wt(np)=wt(np)/lspltdo isplt=1,nsplt
- Page 208 and 209:
ekenp=e(np)-RMcall randomset(rnnolp
- Page 210 and 211:
where ẋ and ẍ (for example) are
- Page 212 and 213:
Figure 4.2: UCBEND simulation at 8.
- Page 214 and 215:
Figure 4.4: UCBEND simulation at 8.
- Page 216 and 217:
Figure 4.5: Geometry and particle t
- Page 218 and 219:
KEY for EGS5 Flow Diagramssubroutin
- Page 220 and 221:
subroutineaphi(br)Version051219-143
- Page 222 and 223:
subroutinebhabhaVersion051219-1435i
- Page 224 and 225:
subroutinebremsVersion051219-1435ib
- Page 226 and 227:
4 5"Rejection sampling loop"callran
- Page 228 and 229:
¥ne¥ne¥ne¥ne¥ne1 2 3 45yesiaus
- Page 230 and 231:
§ne§ne9 10111213callausgab(iarg)c
- Page 232 and 233:
12rejf3 = 1. - br*sinthe/(1. + br*b
- Page 234 and 235:
¨7ese = eig - esg + RMiprofr(irloc
- Page 236 and 237:
subroutineedgbinVersion051219-1435i
- Page 238 and 239:
subroutineeiiVersion051219-1435ieii
- Page 240 and 241:
4 5jnp=1jnp=jnp+1 jnp>neispl yesiq(
- Page 242 and 243:
1 2 3 4 5 6dedx0 = pdedx1(lelke,med
- Page 244 and 245:
1415 16 17 18 19 2021 22k1s0 = k1s0
- Page 246 and 247:
3132 33 34 35 36 37 3839 40dnear(np
- Page 248 and 249:
5051 52 53 54 55 56 57tinel.eq.0.no
- Page 250 and 251:
65 66 67 6869 70 71 72iausfl(iarg+1
- Page 252 and 253:
77 78 79return(to shower)"Ecut disc
- Page 254 and 255:
subroutinehatchVersion060318-1555ih
- Page 256 and 257:
6 7 8 9 1011 12theta = wid*float(is
- Page 258 and 259:
18j=1j=j+1i=1i=i+1j >nmednoi >nregn
- Page 260 and 261:
20 21nolok(im) = 1nm = nm + 1write(
- Page 262 and 263:
24 25irayl .eq.1.or.irayl .eq.2.or.
- Page 264 and 265:
29stopwrite(6,5007) emaxeyesemaxe.e
- Page 266 and 267:
35ecutmn = 1.d20vacdst = vacdst*dun
- Page 268 and 269:
subroutinekaugerVersion051219-1435i
- Page 270 and 271:
subroutinekshellVersion051219-1435i
- Page 272 and 273:
subroutinelaugerVersion051219-1435i
- Page 274 and 275:
subroutinelshell(ll)Version051219-1
- Page 276 and 277:
subroutinelxray(ll)Version051219-14
- Page 278 and 279:
subroutinemollerVersion051219-1435i
- Page 280 and 281:
subroutinemscatVersion060313-1005im
- Page 282 and 283:
34costhe = 1.d0 - 2.d0 * xmusinthe
- Page 284 and 285:
subroutinepairVersion051219-1435ipa
- Page 286 and 287:
4567ichrg.eq.1yescalluphi(2,1)nonp
- Page 288 and 289:
13 14theta=RM/eig"Polar angle ism/E
- Page 290 and 291:
123eigk.le.embind(iz)noyestcros(i)
- Page 292 and 293:
9 10 11nxray.ge.1noyesiphot=1iphot=
- Page 294 and 295:
1516171819costhe = 2.D0*rnnow - 1.D
- Page 296 and 297:
12345eig.lt.0.15noiij=1iij=iij+1yes
- Page 298 and 299:
111213141516yesmedium.ne.medoldnoia
- Page 300 and 301:
2223242526iausfl(iarg+1)ne0noyescal
- Page 302 and 303:
subroutinerayleiVersion051219-1435i
- Page 304 and 305:
1 2 3ie=1ie=ie+1ie>nne(im)noyesz2w
- Page 306 and 307:
7 8 9 10 11 12 13 14 15 16 17k1ez(m
- Page 308 and 309:
subroutinermsfitVersion060314-0855o
- Page 310 and 311:
3.not.useGStrueclose(17)returnfalse
- Page 312 and 313:
6k=1k=k+1k>nmednoyeslok(k) = 05read
- Page 314 and 315:
9 10 11 12 13 14 15 16 17 19read(17
- Page 316 and 317:
1 2calluphi(3,2)"e+,e- or photon"ir
- Page 318 and 319:
1iarg=22iausfl(iarg+1)ne0noyescalla
- Page 320 and 321:
1kount = kount + 1kount.ge.igiganoy
- Page 322 and 323:
subroutinerluxgoVersion051219-1435l
- Page 324 and 325:
2in24 = mod(kount, nskip+24)in24 =
- Page 326 and 327:
subroutinerluxoutVersion051219-1435
- Page 328 and 329:
B.1 IntroductionVersion 5 of the EG
- Page 330 and 331:
+---------+ +-------------+| User |
- Page 332 and 333:
Table B.3: Variable descriptions fo
- Page 334 and 335:
Table B.7: Variable descriptions fo
- Page 336 and 337:
Table B.10: Variable descriptions f
- Page 338 and 339:
Table B.15: Variable descriptions f
- Page 340 and 341:
Second, because of the way PEGS and
- Page 342 and 343:
ECUT and PCUT The ECUT and PCUT arr
- Page 344 and 345:
the scattering strength at ECUT for
- Page 346 and 347:
IBRSPL and NBRSPL The parameters IB
- Page 348 and 349:
vi=0.0wi=1.0iri=2wti=1.0ncases=10id
- Page 350 and 351:
in AUSGAB would keep a running tota
- Page 352 and 353:
inside the shower call loop.B.4.9 O
- Page 354 and 355:
| |Region | Region | Region| |1 | 2
- Page 356 and 357:
end ifreturnendNote that a number o
- Page 358 and 359:
Table B.19: IARG values, IAUSFL ind
- Page 360 and 361:
B.7 UCSAMPL5 — An Example of a
- Page 362 and 363:
! ====================call counters
- Page 364 and 365:
emaxe = ei + RMend if! photon! ----
- Page 366 and 367:
240 FORMAT(//,’ Total energy frac
- Page 368 and 369:
if (w(np).gt.0.0) thendeltaz=(zthic
- Page 370 and 371:
C.1 IntroductionPEGS (Preprocessor
- Page 372 and 373:
+------+| PEGS |+------+| | +------
- Page 374 and 375:
+------------+ +------+| BLOCK DATA
- Page 376 and 377:
+------+ +------+|BREMTM|----------
- Page 378 and 379:
+------+|PAIRTU|+------+| |+ +| |V
- Page 380 and 381:
NameDIFFEREFUNSGFUNSHPLT1LAYMIXPLOT
- Page 382 and 383:
convention) the word is prefixed wi
- Page 384 and 385:
NameBREMDRBREMFRBREMDZBREMFZBREMRMB
- Page 386 and 387:
NamePBR1PBR2PDEDXPHOTTEPHOTTZPSIGQD
- Page 388 and 389:
Table C.7: ELEM option input data l
- Page 390 and 391:
Table C.9: COMP option input data l
- Page 392 and 393:
Table C.11: ENER option input data
- Page 394 and 395:
Table C.13: DECK option input data
- Page 396 and 397:
IRAYL flag is set to 1 or 2 in NAME
- Page 398 and 399:
parameters for the pair production
- Page 400 and 401:
The following are examples of sets
- Page 402 and 403:
Note that the coding of EGS is desi
- Page 404 and 405:
IDFNAMALINDLOGALKEADFMOLTable C.19:
- Page 406 and 407:
C.4 Concluding RemarksIn the previo
- Page 408 and 409:
D.1 Installation of EGS5The current
- Page 410 and 411:
thenCOMPILER="user_defined_compiler
- Page 412 and 413:
egs5run script has started=========
- Page 414 and 415:
working directory is /home/user/egs
- Page 416 and 417:
The EGS5 distribution archive conta
- Page 418 and 419:
tions (in a subdirectory called peg
- Page 420 and 421:
eeldx015.tab eeldx053.tab eeldx091.
- Page 422 and 423:
iff_b.dat iff_k.dat iff_m.dat iff_w
- Page 424 and 425:
that are useful in performing input
- Page 426 and 427:
[13] J. E. Augustin, A. M. Boyarski
- Page 428 and 429:
[44] E. Casnati, A. Tartari, and C.
- Page 430 and 431:
[75] H. Hirayama, W. R. Nelson, A.
- Page 432 and 433:
[108] G. Z. Molière. Theorie der S
- Page 434 and 435:
[140] F. Rohrlich and B. C. Carlson
- Page 436 and 437:
[172] Y. S. Tsai. Pair production a
- Page 438 and 439:
cutoff energy, 26, 28, 30, 31, 38,
- Page 440:
RLUXDAT, 320, 330, 331run5again, 39