2/2/2016 1 A Basic Monte Carlo Course (Electron Gamma Shower) 1 Speaker: Dr. Joel Y.C Cheung Date: 24 th Jan 2015 Venue: Queen Elizabeth Hospital, HK Monte Carlo – A Statistical Calculation, problem solving 2 Probability of drawing 6 ? 1) Perform an actual experiment ! 2) Using the Monte Carlo technique! 3 Monte Carlo - π Calculation 4 Radius = 1 unit Calculation of π π π π using the Monte Carlo method x 2 + y 2 = 1 History run = 19 ù = A r ea of the square â n o: of no: of = 4 â 19 11 = 2:32 π 5 Radius = 1 unit Calculation of π π π π using the Monte Carlo method – cont. History run = ~200 ù = A r ea of the square â n o: of no: of = 4 â 200 150 = 3:0 π 6
15
Embed
Monte Carlo - ππππ Calculation - HKRPS Basic Course in Monte Carlo... · 2/2/2016 2 Monte Carlo – Sales Forecast 7 Sales Forecast Monte Carlo Simulation 8 Monte Carlo in Radiation
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
2/2/2016
1
A Basic Monte Carlo Course
(Electron Gamma Shower)
1
Speaker: Dr. Joel Y.C Cheung
Date: 24th Jan 2015
Venue: Queen Elizabeth Hospital, HK
Monte Carlo –
A Statistical Calculation,
problem solving
2
Probability of drawing 6 ?
1) Perform an actual experiment !
2) Using the Monte Carlo technique!
3
Monte Carlo - ππππ Calculation
4
Radius = 1 unit
Calculation of π π π π using the Monte Carlo method
x2 + y2 = 1
History run = 19
ù = Area of the square âno: ofno: of
= 4 â1911
= 2:32
ππππx
x
5
Radius = 1 unit
Calculation of π π π π using the Monte Carlo method – cont.
History run = ~200
ù = Area of the square âno: ofno: of
= 4 â200150
= 3:0
ππππx
x
6
2/2/2016
2
Monte Carlo –
Sales Forecast
7
Sales Forecast Monte Carlo Simulation
8
Monte Carlo in Radiation Physics
9
Monte Carlo in Radiation Physics
Monte Carlo modeling of particle transport problems in medical and radiation physics gives more advantages than other techniques. Experiments can be done without setting up the physical situation, and results of some “impossible” experiments can be obtained.
e.g. scoring the numbers of created particles or calculating the relative OPFs of narrow beams.
10
Photon or Charge Particle
Before InteractionPosition ? (X,Y,Z)Moving Direction ? (U,V,W)Energy ? (E)
After Interaction / Next StepPosition ? (X,Y,Z)Moving Direction ? (U,V,W)Energy ? (E) Energy Deposition (EDEP)
Determine by Random Numbers ?
11
Radiation Interactions
GNU Free Documentation & public domain
12
2/2/2016
3
Monte Carlo Packages
13
The availability of standardized Monte Carlo packages such as:
•EGSNRC•BEAM•EGS•ETRAN/ITS•PENELOPE•MCNP•GEANT
along with the development of more powerful and inexpensive computers has allowed more widespread use of the technique.
Monte Carlo packages
14
EGS4 – Electron Gamma
Shower
15
EGS4 (Electron Gamma Shower) code
Standford Linear Accelerator Center- by Nelson, Harayama and Rogers.
16
http://rcwww.kek.jp/research/egs/
17
Energy < 100keV
http://rcwww.kek.jp/research/egs/epub.html
18
2/2/2016
4
http://rcwww.kek.jp/research/egs/egs5.html
19
EGSnrc Code
20
http://irs.inms.nrc.ca/
21 22
23
Fortran Compilers
24
2/2/2016
5
Fortran Compilers
25
EGS4 Code
26
EGS4 Structure
27
EGS4 Distribution / Fortran Compiler
28
EGS4 folders on PC
Creating material file(s) for radiation interaction