Radiation Transport Calculations by a Monte Carlo Method. Hideo Hirayama and Yoshihito Namito KEK, High Energy Accelerator Research Organization
|
|
- Mary Gaines
- 5 years ago
- Views:
Transcription
1 Radiation Transport Calculations by a Monte Carlo Method Hideo Hirayama and Yoshihito Namito KEK, High Energy Accelerator Research Organization
2 Monte Carlo Method A method used to solve a problem with random numbers is called a Monte Carlo Method. A name of Monte Carlo was introduced by J. von Neumann and S. M. Ulam. Random numbers are key tool for the Monte Carlo method
3 Random numbers Use a dice, a roulette etc. very slow. Use a table of random numbers. A table of random numbers has been well examined concerning its statistical characteristics. It is required to store a whole table in computer data storage. It is currently is not very fast to produce random numbers. Use physical random numbers like the decay of a radioisotope. It is not easy to digitalize, and has a weakness concerning stability and reproducibility.
4 Pseudo random Numbers Produce random numbers successively from a seed random number, R 0, using a recurrence formula in the form of R n+1 =f(r n ) Pseudo random number residuals by a divider, m. There are m different integers at most and, therefore, pseudo random numbers have a limited period. God pseudo random numbers Fast to create A long sequence Reproducibility Good statistical characteristics It is possible to create pseudo random number between 0 and 1 by dividing by m.
5 Another Type of Random Number Generator RANECU random number generator F. James, A review of pseudorandom number generators, Comp. Phys Comm 60 (1990) Periodicity Can be easily controlled by 2 seeds (iseed1 and iseed2) Use as the random number generator of egs5 Marasaglia-Zaman random number generator G. Masaglia and A. Zaman, A New Class of Random Number Generator, Annals of Applied Probability 1(1991) Long periodicity ~ Portable to all 32-bit machines
6 Pseudo Random Number Linear congruence methods proposed by D. H. Lehmer is most widely used. R n+1 =mod(ar n +b,m) a, b and m are positive integer and m is the length of the integer value allowed in the compiler. Name a b m RANDU SLAC RAN SLAC RAN
7 Production of pseudo random numbers using pocket calculator Produce 10 random numbers for R 0 =3, a=5, b=0 and m=16. Confirm that the same random sequence appears from some point. What is a sequence in this case? Check for a different R 0.
8 n R n R n *5 R n+1 =mod(r n *a,m) 0 R 0 =3 15 R 1 =15-INT(15/16)*16= mod(r n *a,m)= R n *a-int(rn*a/m) INT(15/16)=0
9 Calculation π using random numbers Select 2 random numbers from an arbitrary place in Table 2 (created using SLAC RAN6) in order Count the number of pair which satisfy R 2 2 = ξ + η 1.0 AREA= π /4 π= 4 x (probability within this area)
10 Discrete probability process If a probability variable (x) takes on discrete value (x i ) with probabilities (p i ) such that x=x i if Cumulative probability function (PDF) 1, ) ( 1 = = = n i x n p i F + = + = = < = ) ( ) ( i j j i i j j i p x F p x F η
11 Discrete type sampling(1) Example) Sample reaction when photoelectric is 30%, Compton scattering 50% and pair production 20% Photoelectric Compton scattering Pair creation
12 Discrete type sampling (2) Random number η Probability of Compton scattering Photoelectric Compton Scattering Pair creation Cumulative probability function (pdf)
13 Example determination of photon reaction The probability of the photoelectric effect, Compton scattering and pair creation at a photon interaction are P photo, P Compt and P pair, respectively. P photo +P Compt + P pair =1 η P photo, photoelectric P photo η < Pphoto + PCompt, Compton scattering P photo P η, pair creation + Compt
14 Sampling method (direct sampling) A probability distribution function (PDF:f(x)) is defined over the range [a,b]. b ( f ( x) dx = 1) Its cumulative probability function (CDF:F(x)) a a x F ( x) = f ( xi ) dx By its definition, we can map F(x) onto a range of random variable, η, where 0 η 1. Having mapped the random number onto F(x), we may invert the equation to give x = F 1 ( η) The way to determine x by solving the above equation is called a direct method. i
15 F(x) 1 η 0 a x b x
16 Example-determination of flight distance If the interaction probability of a particle per unit distance is Σ t, the probability that a first interaction occurs between l and l+dl is p( l) dl = e Σ l t Σ t dl η = P( l) = 1 l = Σ t l 0 p( l 1 ) dl = 1 e ln(1 η) = λ ln(1 η) l : flight distance λ: mean free path 1 η is equivalent to η 1 Σ l t l = λ ln( η)
17
18 Red:in User Code Black:egs5 Infinite geometry Photon Initial Condition:Energy Position Direction E 0, X 0, Y 0,Z 0, U 0, V 0, W 0
19 Determine Interaction Point l=-ln(δ)/µ x=x 0 +u 0 l, y=y 0 +v 0 l. z=z 0 +w 0 l l Initial Condition:Energy Position Direction e 0, x 0, y 0, z 0, u 0, v 0, w 0
20 Determine Interaction Type Photoelectric:a, Compton scattering:b, Pair creation:c δ<=a/(a+b+c):photoelectric a/(a+b+c)< δ<=(a+c)/(a+b+c):compton scattering δ>(a+c)/(a+b+c):pair creation x=x 0 +u 0 l, y=y 0 +v 0 l. z=z 0 +w 0 l l Initial Condition:Energy Position Direction e 0, x 0, y 0, z 0, u 0, v 0, w 0
21 Electron Follows produced particle Determination of energy and direction of each particle Photon l x=x 0 +u 0 l, y=y 0 +v 0 l. z=z 0 +w 0 l Initial Condition:Energy Position Direction e 0, x 0, y 0, z 0, u 0, v 0, w 0
22 d>l:move to Interaction point d<=l:move a distance d Region boundary Same material:distance to Interaction point=l-d Different material:determine Interaction point Calculate a straight-line distance to the boundary (HOWFAR) l Interaction point: x=x 0 +u 0 l, y=y 0 +v 0 l. z=z 0 +w 0 l d Initial Condition:Energy Position Direction e 0, x 0, y 0, z 0, u 0, v 0, w 0
23 Region boundary Calculate a straight-line distance to the boundary (HOWFAR) l Store Information (AUSGAB) Particle moves:energy deposition Track length Boundary crossing Photoelectric:Energy deposition Below cut-off IDISC=1 IAUSFL(IARG)=1:optional d Initial Condition:Energy Position Direction e 0, x 0, y 0, z 0, u 0, v 0, w 0
24 AUSGAB Region boundary Store Information HOWFAR Calculate a straight-line distance to the boundary Set IDISCT=1 l d Main program Material assignment, Geometry Initial condition:energy Position Direction E 0, X 0, Y 0,Z 0, U 0, V 0, W 0 Analysis of results and output them
25 Determine Interaction Point l=-ln(δ)/σ It is difficult to treat an elastic scattering of electrons (or positrons) like photons due to its large number. l 5 x l 6 Charged particle lose a part of its l 4 x kinetic energy when it moves some x x l distance via ionization or excitation 3 l Energy deposition at each step Electron and positron True track length x Stopping power (de/dx) Initial Condition:Energy Position Direction e 0, x 0, y 0, z 0, u 0, v 0, w 0 l x x l 7 l 8 x Condensed History Technique Divide flight path to many small steps and apply Multiple Scattering Model to obtain true track length and change of direction and displacement due to many elastic scattering
26 Photon transport using a pocket calculator (Fig.1) Suppose that material A having thickness of 50 cm exist as shown Fig MeV photons enter this system from the left end, the mean free path is 20 cm, the ratio of the photoelectric effect and Compton scattering is 1:1 and, a scattered photon does not change its energy or direction. Exp.1 First random number λ=-20.0ln(0.2336)= (cm)<50.0(cm) Next random number is (<0.5). photoelectric effect Exp. 2 Next random number is λ=-20.0ln( )= (cm)<50.0(cm) Next random number is (>0.5). Compton scattering Next random number is λ=-20.0ln( )= (cm)< (cm)
27 Single layer No. d(cm) Random l(cm) d>l d l Random Photo. Compt. number number Exp * * Exp * * * * * * *
28 Fig. 1 Trajectories for a single layer 0cm 10cm 20cm 30cm 40cm 50cm Exp. 1 Exp. 2 C C P C
29 Photon transport using a pocket calculator (Fig.2) Suppose that material B (10cm) exits behind material A (40cm) as shown in Fig MeV photons enter this system from the left end, the mean free path and the ratio of the photoelectric effect and Compton scattering in medium A are same as the previous case, the mean free path of medium B is 3cm, ratio of the photoelectric effect and Compton scattering of medium B is 3:1 and, a scattered photon does not change its energy or direction. Exp.1 First random number λ=-20.0ln( )= (cm)<40.0(cm) Next random number is (>0.5). Compton scattering Next random number is λ=-20.0ln(0.2336)= (cm)> (cm) Move to boundary (40.0cm) Next random number is λ=-3.0ln( )= (cm)<10.0(cm)
30 Medium A No. d(cm) Random number l(cm) d>l d l Random number Photo. Compt Exp * * * Medium B d(cm) Random number l(cm) d>l d l Random number Photo. Compt * * * * * *
31 Fig. 2 Trajectories in double layers. 0cm 10cm 20cm 30cm 40cm 50cm Ex. 1 C C C P Medium A Medium B
32 Complex, but more realistic, example of photon transport Consider 10cm aluminum plate shown in Fig.3. Suppose that 0.5 MeV photon enters this system from left end. Photon is scattered with equal probability for each 45 at Compton scattering for all photon energies. The photon energy after scattering is calculated by E = E0 E (1 cosθ )
33 Complex, but more realistic, example of photon transport Suppose that the azimuthal angle after Compton scattering is 0º or 180º with an equal probability. 0º is 90º left from the particle direction and 180º is 90º right. Use the mean free path (mfp) and branching ratio for each photon energy in Fig. 4 and 5. Set the cutoff energy of photons to 0.05 MeV.
34
35
36
Monte Carlo radiation transport codes
Monte Carlo radiation transport codes How do they work? Michel Maire (Lapp/Annecy) 16/09/2011 introduction to Monte Carlo radiation transport codes 1 Outline From simplest case to complete process : Decay
More informationMonte Carlo radiation transport codes
Monte Carlo radiation transport codes How do they work? Michel Maire (Lapp/Annecy) 23/05/2007 introduction to Monte Carlo radiation transport codes 1 Decay in flight (1) An unstable particle have a time
More informationPassage of particles through matter
Passage of particles through matter Alexander Khanov PHYS6260: Experimental Methods is HEP Oklahoma State University September 11, 2017 Delta rays During ionization, the energy is transferred to electrons
More informationHow does the computer generate observations from various distributions specified after input analysis?
1 How does the computer generate observations from various distributions specified after input analysis? There are two main components to the generation of observations from probability distributions.
More informationThe Monte Carlo Method in Medical Radiation Physics
Monte Carlo in Medical Physics The Monte Carlo Method in Medical Radiation Physics P Andreo, Professor of Medical Radiation Physics Stockholm University @ Karolinska Univ Hospital, Stockholm, Sweden ICRM
More information2. Passage of Radiation Through Matter
2. Passage of Radiation Through Matter Passage of Radiation Through Matter: Contents Energy Loss of Heavy Charged Particles by Atomic Collision (addendum) Cherenkov Radiation Energy loss of Electrons and
More informationEnrico Fermi and the FERMIAC. Mechanical device that plots 2D random walks of slow and fast neutrons in fissile material
Monte Carlo History Statistical sampling Buffon s needles and estimates of π 1940s: neutron transport in fissile material Origin of name People: Ulam, von Neuman, Metropolis, Teller Other areas of use:
More informationMonte Carlo Techniques
Physics 75.502 Part III: Monte Carlo Methods 40 Monte Carlo Techniques Monte Carlo refers to any procedure that makes use of random numbers. Monte Carlo methods are used in: Simulation of natural phenomena
More informationOutline. Chapter 6 The Basic Interactions between Photons and Charged Particles with Matter. Photon interactions. Photoelectric effect
Chapter 6 The Basic Interactions between Photons and Charged Particles with Matter Radiation Dosimetry I Text: H.E Johns and J.R. Cunningham, The physics of radiology, 4 th ed. http://www.utoledo.edu/med/depts/radther
More informationAt the conclusion of this lesson the trainee will be able to: a) Write a typical equation for the production of each type of radiation.
RADIOACTIVITY - SPONTANEOUS NUCLEAR PROCESSES OBJECTIVES At the conclusion of this lesson the trainee will be able to: 1. For~, p and 7 decays a) Write a typical equation for the production of each type
More informationBethe-Block. Stopping power of positive muons in copper vs βγ = p/mc. The slight dependence on M at highest energies through T max
Bethe-Block Stopping power of positive muons in copper vs βγ = p/mc. The slight dependence on M at highest energies through T max can be used for PID but typically de/dx depend only on β (given a particle
More informationMean Intensity. Same units as I ν : J/m 2 /s/hz/sr (ergs/cm 2 /s/hz/sr) Function of position (and time), but not direction
MCRT: L5 MCRT estimators for mean intensity, absorbed radiation, radiation pressure, etc Path length sampling: mean intensity, fluence rate, number of absorptions Random number generators J Mean Intensity
More informationChapter 4: Monte Carlo Methods. Paisan Nakmahachalasint
Chapter 4: Monte Carlo Methods Paisan Nakmahachalasint Introduction Monte Carlo Methods are a class of computational algorithms that rely on repeated random sampling to compute their results. Monte Carlo
More informationDetecting high energy photons. Interactions of photons with matter Properties of detectors (with examples)
Detecting high energy photons Interactions of photons with matter Properties of detectors (with examples) Interactions of high energy photons with matter Cross section/attenution length/optical depth Photoelectric
More informationDetectors in Nuclear Physics: Monte Carlo Methods. Dr. Andrea Mairani. Lectures I-II
Detectors in Nuclear Physics: Monte Carlo Methods Dr. Andrea Mairani Lectures I-II INTRODUCTION Sampling from a probability distribution Sampling from a probability distribution X λ Sampling from a probability
More informationIII. Energy Deposition in the Detector and Spectrum Formation
1 III. Energy Deposition in the Detector and Spectrum Formation a) charged particles Bethe-Bloch formula de 4πq 4 z2 e 2m v = NZ ( ) dx m v ln ln 1 0 2 β β I 0 2 2 2 z, v: atomic number and velocity of
More informationEEE4101F / EEE4103F Radiation Interactions & Detection
EEE4101F / EEE4103F Radiation Interactions & Detection 1. Interaction of Radiation with Matter Dr. Steve Peterson 5.14 RW James Department of Physics University of Cape Town steve.peterson@uct.ac.za March
More informationOutline. Radiation Interactions. Spurs, Blobs and Short Tracks. Introduction. Radiation Interactions 1
Outline Radiation Interactions Introduction Interaction of Heavy Charged Particles Interaction of Fast Electrons Interaction of Gamma Rays Interactions of Neutrons Radiation Exposure & Dose Sources of
More informationEmphasis on what happens to emitted particle (if no nuclear reaction and MEDIUM (i.e., atomic effects)
LECTURE 5: INTERACTION OF RADIATION WITH MATTER All radiation is detected through its interaction with matter! INTRODUCTION: What happens when radiation passes through matter? Emphasis on what happens
More informationRandom numbers and generators
Chapter 2 Random numbers and generators Random numbers can be generated experimentally, like throwing dice or from radioactive decay measurements. In numerical calculations one needs, however, huge set
More information2008 Winton. Review of Statistical Terminology
1 Review of Statistical Terminology 2 Formal Terminology An experiment is a process whose outcome is not known with certainty The experiment s sample space S is the set of all possible outcomes. A random
More informationQuantitative Assessment of Scattering Contributions in MeV-Industrial X-ray Computed Tomography
11th European Conference on Non-Destructive Testing (ECNDT 2014), October 6-10, 2014, Prague, Czech Republic More Info at Open Access Database www.ndt.net/?id=16530 Quantitative Assessment of Scattering
More informationMonte Carlo Radiation Transfer I
Monte Carlo Radiation Transfer I Monte Carlo Photons and interactions Sampling from probability distributions Optical depths, isotropic emission, scattering Monte Carlo Basics Emit energy packet, hereafter
More informationReview of Statistical Terminology
Review of Statistical Terminology An experiment is a process whose outcome is not known with certainty. The experiment s sample space S is the set of all possible outcomes. A random variable is a function
More informationEEE4106Z Radiation Interactions & Detection
EEE4106Z Radiation Interactions & Detection 2. Radiation Detection Dr. Steve Peterson 5.14 RW James Department of Physics University of Cape Town steve.peterson@uct.ac.za May 06, 2015 EEE4106Z :: Radiation
More informationCHARGED PARTICLE INTERACTIONS
CHARGED PARTICLE INTERACTIONS Background Charged Particles Heavy charged particles Charged particles with Mass > m e α, proton, deuteron, heavy ion (e.g., C +, Fe + ), fission fragment, muon, etc. α is
More informationInteraction of Particles and Matter
MORE CHAPTER 11, #7 Interaction of Particles and Matter In this More section we will discuss briefly the main interactions of charged particles, neutrons, and photons with matter. Understanding these interactions
More informationBasics of Monte Carlo Simulation
Basics of Monte Carlo Simulation condensed from a presentation by M. Asai Dennis Wright Geant4 Tutorial at Lund University 3 September 2018 Contents Historical review of Monte Carlo methods Basics of Monte
More informationChapter Four (Interaction of Radiation with Matter)
Al-Mustansiriyah University College of Science Physics Department Fourth Grade Nuclear Physics Dr. Ali A. Ridha Chapter Four (Interaction of Radiation with Matter) Different types of radiation interact
More informationChapter NP-4. Nuclear Physics. Particle Behavior/ Gamma Interactions TABLE OF CONTENTS INTRODUCTION OBJECTIVES 1.0 IONIZATION
Chapter NP-4 Nuclear Physics Particle Behavior/ Gamma Interactions TABLE OF CONTENTS INTRODUCTION OBJECTIVES 1.0 IONIZATION 2.0 ALPHA PARTICLE INTERACTIONS 3.0 BETA INTERACTIONS 4.0 GAMMA INTERACTIONS
More informationPhysics 736. Experimental Methods in Nuclear-, Particle-, and Astrophysics. Lecture 4
Physics 736 Experimental Methods in Nuclear-, Particle-, and Astrophysics Lecture 4 Karsten Heeger heeger@wisc.edu Homework Homework is posted on course website http://neutrino.physics.wisc.edu/teaching/phys736/
More informationStatistische Methoden der Datenanalyse. Kapitel 3: Die Monte-Carlo-Methode
1 Statistische Methoden der Datenanalyse Kapitel 3: Die Monte-Carlo-Methode Professor Markus Schumacher Freiburg / Sommersemester 2009 Basiert auf Vorlesungen und Folien von Glen Cowan und Abbildungen
More informationThe interaction of radiation with matter
Basic Detection Techniques 2009-2010 http://www.astro.rug.nl/~peletier/detectiontechniques.html Detection of energetic particles and gamma rays The interaction of radiation with matter Peter Dendooven
More informationElectromagnetic and hadronic showers development. G. Gaudio, M. Livan The Art of Calorimetry Lecture II
Electromagnetic and hadronic showers development 1 G. Gaudio, M. Livan The Art of Calorimetry Lecture II Summary (Z dependence) Z Z 4 5 Z(Z + 1) Z Z(Z + 1) 2 A simple shower 3 Electromagnetic Showers Differences
More informationParticle Interactions in Detectors
Particle Interactions in Detectors Dr Peter R Hobson C.Phys M.Inst.P. Department of Electronic and Computer Engineering Brunel University, Uxbridge Peter.Hobson@brunel.ac.uk http://www.brunel.ac.uk/~eestprh/
More information3 Radioactivity - Spontaneous Nuclear Processes
3 Radioactivity - Spontaneous Nuclear Processes Becquerel was the first to detect radioactivity. In 1896 he was carrying out experiments with fluorescent salts (which contained uranium) and found that
More informationParticle-Matter Interactions
Particle-Matter Interactions to best detect radiations and particles we must know how they behave inside the materials 8/30/2010 PHYS6314 Prof. Lou 1 Stable Particles Visible to a Detector Hadrons (Baryon/Meson)
More informationDETECTORS. I. Charged Particle Detectors
DETECTORS I. Charged Particle Detectors A. Scintillators B. Gas Detectors 1. Ionization Chambers 2. Proportional Counters 3. Avalanche detectors 4. Geiger-Muller counters 5. Spark detectors C. Solid State
More informationShell Atomic Model and Energy Levels
Shell Atomic Model and Energy Levels (higher energy, deeper excitation) - Radio waves: Not absorbed and pass through tissue un-attenuated - Microwaves : Energies of Photos enough to cause molecular rotation
More informationS.A. Teukolsky, Computers in Physics, Vol. 6, No. 5,
Physics 75.502 Part III: Monte Carlo Methods 139 Part III: Monte Carlo Methods Topics: Introduction Random Number generators Special distributions General Techniques Multidimensional simulation References:
More informationScintillation Detector
Scintillation Detector Introduction The detection of ionizing radiation by the scintillation light produced in certain materials is one of the oldest techniques on record. In Geiger and Marsden s famous
More informationCopyright 2008, University of Chicago, Department of Physics. Experiment VI. Gamma Ray Spectroscopy
Experiment VI Gamma Ray Spectroscopy 1. GAMMA RAY INTERACTIONS WITH MATTER In order for gammas to be detected, they must lose energy in the detector. Since gammas are electromagnetic radiation, we must
More informationDrift plane. substrate (20ÉIm polyimide) 200ÉIm. Back strip (180ÉIm width) Base (Ceramic) Anode strip (10ÉIm width) Cathode strip (100ÉIm width)
Proceedings of the Second International Workshop on EGS, 8.-1. August, Tsukuba, Japan KEK Proceedings -, pp.11-17 Development of Gamma-Ray Direction Detector Based on MSGC T. Nagayoshi 1, H. Kubo 1, A.
More informationInteraction of Ionizing Radiation with Matter
Type of radiation charged particles photonen neutronen Uncharged particles Charged particles electrons (β - ) He 2+ (α), H + (p) D + (d) Recoil nuclides Fission fragments Interaction of ionizing radiation
More informationOffice of Nonproliferation and Verification Research and Development University and Industry Technical Interchange (UITI2011) Review Meeting
Office of Nonproliferation and Verification Research and Development niversity and Industry Technical Interchange (ITI2011) Review Meeting Modeling of SNM Fission Signatures and December 7, 2011 Gennady
More informationNeutrino detection. Kate Scholberg, Duke University International Neutrino Summer School Sao Paulo, Brazil, August 2015
Neutrino detection Kate Scholberg, Duke University International Neutrino Summer School Sao Paulo, Brazil, August 2015 Sources of wild neutrinos The Big Bang The Atmosphere (cosmic rays) Super novae AGN's,
More informationMonte Carlo sampling. FLUKA Beginner s Course
Monte Carlo sampling FLUKA Beginner s Course Overview: General concepts: Phase space Monte Carlo foundations Simulation vs. integration Sampling techniques discrete by inversion by rejection Results and
More informationMotivation. g-spectroscopy deals with g-ray detection and is one of the most relevant methods to investigate excited states in nuclei.
Motivation Spins and excited states of double-magic nucleus 16 O Decay spectra are caused by electro-magnetic transitions. g-spectroscopy deals with g-ray detection and is one of the most relevant methods
More informationModern Methods of Data Analysis - WS 07/08
Modern Methods of Data Analysis Lecture III (29.10.07) Contents: Overview & Test of random number generators Random number distributions Monte Carlo The terminology Monte Carlo-methods originated around
More informationA Short Course on Geant4 Simulation Toolkit. Introduction.
A Short Course on Geant4 Simulation Toolkit Introduction http://cern.ch/geant4 The full set of lecture notes of this Geant4 Course is available at http://www.ge.infn.it/geant4/events/nss2003/geant4course.html
More informationPolarised Geant4 Applications at the ILC
Polarised Geant4 Applications at the ILC Andreas Schälicke, Karim Laihem 2 and Pavel Starovoitov - DESY Platanenallee 6, 578 Zeuthen - Germany 2- RWTH Aachen - Phys. Inst. IIIB Physikzentrum, 5256 Aachen-
More informationParticle Detectors. How to See the Invisible
Particle Detectors How to See the Invisible Which Subatomic Particles are Seen? Which particles live long enough to be visible in a detector? 2 Which Subatomic Particles are Seen? Protons Which particles
More informationWeek 2: Chap. 2 Interaction of Radiation
Week 2: Chap. 2 Interaction of Radiation Introduction -- Goals, roll back the fog -- General Nomenclature -- Decay Equations -- Laboratory Sources Interaction of Radiation with Matter -- Charged Particles
More informationPhysics of Particle Beams. Hsiao-Ming Lu, Ph.D., Jay Flanz, Ph.D., Harald Paganetti, Ph.D. Massachusetts General Hospital Harvard Medical School
Physics of Particle Beams Hsiao-Ming Lu, Ph.D., Jay Flanz, Ph.D., Harald Paganetti, Ph.D. Massachusetts General Hospital Harvard Medical School PTCOG 53 Education Session, Shanghai, 2014 Dose External
More informationBlackbody Radiation. Rayleigh-Jeans law was an attempt to explain blackbody radiation based on classical ideas:
Blackbody Radiation A Blackbody is an ideal system that absorbs all radiation incident on it. Emission of radiation by a blackbody is independent of the properties of its wall, but depends only on its
More informationInteraction of charged particles and photons with matter
Interaction of charged particles and photons with matter Robert Miyaoka, Ph.D. Old Fisheries Center, Room 200 rmiyaoka@u.washington.edu Passage of radiation through matter depends on Type of radiation
More informationComparison of available measurements of the absolute air-fluorescence yield*
Comparison of available measurements of the absolute air-fluorescence yield* J. Rosado, F. Blanco and F. Arqueros Universidad Complutense de Madrid * J. Rosado et al., Astropart. Phys. 34 (2010) 164 Outline
More informationInteraction of Electron and Photons with Matter
Interaction of Electron and Photons with Matter In addition to the references listed in the first lecture (of this part of the course) see also Calorimetry in High Energy Physics by Richard Wigmans. (Oxford
More informationPHYS 571 Radiation Physics
PHYS 571 Radiation Physics Prof. Gocha Khelashvili http://blackboard.iit.edu login Interaction of Electrons with Matter The Plan Interactions of Electrons with Matter Energy-Loss Mechanism Collisional
More informationAiro International Research Journal October, 2015 Volume VI, ISSN:
1 INTERACTION BETWEEN CHARGED PARTICLE AND MATTER Kamaljeet Singh NET Qualified Declaration of Author: I hereby declare that the content of this research paper has been truly made by me including the title
More informationMCRT: L4 A Monte Carlo Scattering Code
MCRT: L4 A Monte Carlo Scattering Code Plane parallel scattering slab Optical depths & physical distances Emergent flux & intensity Internal intensity moments Constant density slab, vertical optical depth
More informationPhysics of particles. H. Paganetti PhD Massachusetts General Hospital & Harvard Medical School
Physics of particles H. Paganetti PhD Massachusetts General Hospital & Harvard Medical School Introduction Dose The ideal dose distribution ideal Dose: Energy deposited Energy/Mass Depth [J/kg] [Gy] Introduction
More informationNeutron Transport Calculations Using Monte-Carlo Methods. Sean Lourette Fairport High School Advisor: Christian Stoeckl
Neutron Transport Calculations Using Monte-Carlo Methods Sean Lourette Fairport High School Advisor: Christian Stoeckl Laboratory for Laser Energetics University of Rochester Summer High School Research
More informationNeutron pulse height analysis (R405n)
Neutron pulse height analysis (R405n) Y. Satou April 6, 2011 Abstract A pulse height analysis was made for the neutron counter hodoscope used in R405n. By normalizing the pulse height distributions measured
More informationHow does the computer generate observations from various distributions specified after input analysis?
1 How does the computer generate observations from various distributions specified after input analysis? There are two main components to the generation of observations from probability distributions.
More informationPhysics of Radiotherapy. Lecture II: Interaction of Ionizing Radiation With Matter
Physics of Radiotherapy Lecture II: Interaction of Ionizing Radiation With Matter Charge Particle Interaction Energetic charged particles interact with matter by electrical forces and lose kinetic energy
More informationFrom Random Numbers to Monte Carlo. Random Numbers, Random Walks, Diffusion, Monte Carlo integration, and all that
From Random Numbers to Monte Carlo Random Numbers, Random Walks, Diffusion, Monte Carlo integration, and all that Random Walk Through Life Random Walk Through Life If you flip the coin 5 times you will
More informationOverview: In this experiment we study the decay of a radioactive nucleus, Cesium 137. Figure 1: The Decay Modes of Cesium 137
Radioactivity (Part I and Part II) 7-MAC Objectives: To measure the absorption of beta and gamma rays To understand the concept of half life and to measure the half life of Ba 137* Apparatus: Radioactive
More informationCHAPTER 3 Prelude to Quantum Theory. Observation of X Rays. Thomson s Cathode-Ray Experiment. Röntgen s X-Ray Tube
CHAPTER Prelude to Quantum Theory.1 Discovery of the X Ray and the Electron. Determination of Electron Charge. Line Spectra.4 Quantization.5 Blackbody Radiation.6 Photoelectric Effect.7 X-Ray Production.8
More informationCalculation of the Stopping Power for Intermediate Energy Positrons. Önder Kabaday, and M. Çaǧatay Tufan
CHINESE JOURNAL OF PHYSICS VOL. 44, NO. 4 AUGUST 006 Calculation of the Stopping Power for Intermediate Energy Positrons Hasan Gümüş, Önder Kabaday, and M. Çaǧatay Tufan Department of Physics, Faculty
More informationDetector Simulation. Mihaly Novak CERN PH/SFT
Detector Simulation Mihaly Novak CERN PH/SFT CERN Summer Student Program, 1 August 2017 Foreword This lecture is aimed to offer a simple and general introduction to detector simulation. Geant4 will be
More informationLecture 20. Randomness and Monte Carlo. J. Chaudhry. Department of Mathematics and Statistics University of New Mexico
Lecture 20 Randomness and Monte Carlo J. Chaudhry Department of Mathematics and Statistics University of New Mexico J. Chaudhry (UNM) CS 357 1 / 40 What we ll do: Random number generators Monte-Carlo integration
More informationGeneral Overview of Gas Filled Detectors
GAS-FILLED DETECTOR General Overview of Gas Filled Detectors Gas-Filled Detectors Ion chamber Proportional counter G-M (Geiger-Miller) counter Diagram of a Generic Gas-Filled Detector A Anode High-voltage
More informationAppendix A2. Particle Accelerators and Detectors The Large Hadron Collider (LHC) in Geneva, Switzerland on the Border of France.
Appendix A. Particle Accelerators and Detectors The Large Hadron Collider (LHC) in Geneva, Switzerland on the Border of France. Prepared by: Arash Akbari-Sharbaf Why Build Accelerators? Probe deeper From
More informationChapiter VII: Ionization chamber
Chapiter VII: Ionization chamber 1 Types of ionization chambers Sensitive volume: gas (most often air direct measurement of exposure) ionization chamber Sensitive volume: semiconductor (silicon, germanium,
More informationChapter V: Cavity theories
Chapter V: Cavity theories 1 Introduction Goal of radiation dosimetry: measure of the dose absorbed inside a medium (often assimilated to water in calculations) A detector (dosimeter) never measures directly
More informationOverview: In this experiment we will study the decay of a radioactive nucleus, Cesium. Figure 1: The Decay Modes of Cesium 137
Radioactivity (Part I and Part II) Objectives: To measure the absorption of beta and gamma rays To understand the concept of half life and to measure the half life of Ba 137* Apparatus: Radioactive source,
More informationPhysics of Radiography
Physics of Radiography Yao Wang Polytechnic Institute of NYU Brooklyn, NY 11201 Based on J L Prince and J M Links Medical Imaging Signals and Based on J. L. Prince and J. M. Links, Medical Imaging Signals
More informationPhysics 111 Homework Solutions Week #9 - Friday
Physics 111 Homework Solutions Week #9 - Friday Tuesday, March 1, 2011 Chapter 24 Questions 246 The Compton shift in wavelength for the proton and the electron are given by Δλ p = h ( 1 cosφ) and Δλ e
More information1.E Neutron Energy (MeV)
Proceedings of the Second International Workshop on EGS, 8.-12. August 2000, Tsukuba, Japan KEK Proceedings 200-20, pp.130-134 Measurements of Photoneutron Spectra from Thick Pb Target Bombarded by 1.2
More informationAlgorithms and Networking for Computer Games
Algorithms and Networking for Computer Games Chapter 2: Random Numbers http://www.wiley.com/go/smed What are random numbers good for (according to D.E. Knuth) simulation sampling numerical analysis computer
More informationDavid B. Cassidy. Department of Physics and Astronomy, University of California, Riverside, USA. Varenna, July 09
Experimental production of many- positron systems: L2, techniques David B. Cassidy Department of Physics and Astronomy, University of California, Riverside, USA cassidy@physics.ucr.edu Varenna, July 09
More information11. Point Detector Tally
and suppose that we want half the source to be in the range 005 < E < 1 and the other half to be in the range 1 < E < 20 Then the input is SPn -5 a SIn 005 1 20 SBn C 0 5 1 MCNP breaks up the function
More informationCHAPTER 2 INTERACTION OF RADIATION WITH MATTER
CHAPTER 2 INTERACTION OF RADIATION WITH MATTER 2.1 Introduction When gamma radiation interacts with material, some of the radiation will be absorbed by the material. There are five mechanisms involve in
More informationPhysics 736. Experimental Methods in Nuclear-, Particle-, and Astrophysics. Lecture 3
Physics 736 Experimental Methods in Nuclear-, Particle-, and Astrophysics Lecture 3 Karsten Heeger heeger@wisc.edu Review of Last Lecture a colleague shows you this data... what type of reaction is this?
More informationRandom processes and probability distributions. Phys 420/580 Lecture 20
Random processes and probability distributions Phys 420/580 Lecture 20 Random processes Many physical processes are random in character: e.g., nuclear decay (Poisson distributed event count) P (k, τ) =
More informationShielded Scintillator for Neutron Characterization
Shielded Scintillator for Neutron Characterization A Thesis Submitted in Partial Fulfillment of the Requirements for Graduation with Research Distinction in Engineering Physics By Patrick X. Belancourt
More informationUniversity of Oslo. Department of Physics. Interaction Between Ionizing Radiation And Matter, Part 2 Charged-Particles.
Interaction Between Ionizing Radiation And Matter, Part Charged-Particles Audun Sanderud Excitation / ionization Incoming charged particle interact with atom/molecule: Ionization Excitation Ion pair created
More informationChapter V: Interactions of neutrons with matter
Chapter V: Interactions of neutrons with matter 1 Content of the chapter Introduction Interaction processes Interaction cross sections Moderation and neutrons path For more details see «Physique des Réacteurs
More informationActinide Target Station & 238
Actinide Target Station & 238 U Photofission Yield Pierre Bricault SEEC, 25-26 March 2008 TRIUMF CANADA S NATIONAL LABORATORY FOR PARTICLE AND NUCLEAR PHYSICS LABORATOIRE NATIONAL CANADIEN POUR LA RECHERCHE
More informationelectrons out of, or ionize, material in their paths as they pass. Such radiation is known as
Detecting radiation It is always possible to detect charged particles moving through matter because they rip electrons out of, or ionize, material in their paths as they pass. Such radiation is known as
More informationRESPONSE CALCULATIONS OF THE CdZnTe DETECTOR USING EGS4. James C. Liu, W. R. Nelson and R. Seefred
SLAC-PUB-8606 (September 5, 2000) RESPONSE CALCULATIONS OF THE CdZnTe DETECTOR USING EGS4 James C. Liu, W. R. Nelson and R. Seefred Stanford Linear Accelerator Center, MS 48, P. O. Box 4349, Stanford,
More informationPHYS 3650L - Modern Physics Laboratory
PHYS 3650L - Modern Physics Laboratory Laboratory Advanced Sheet Photon Attenuation 1. Objectives. The objectives of this laboratory exercise are: a. To measure the mass attenuation coefficient at a gamma
More information05/11/2013. Nuclear Fuel Cycle Ionizing radiation. Typical decay energies. Radiation with energy > 100 ev. Ionize an atom < 15eV
Nuclear Fuel Cycle 2013 Lecture 4: Interaction of Ionizing Radiation with Matter Ionizing radiation Radiation with energy > 100 ev Ionize an atom < 15eV Break a bond 1-5 ev Typical decay energies α: 4-9
More informationComparative Analysis of Nuclear Cross Sections in Monte Carlo Methods for Medical Physics Applications
Comparative Analysis of Nuclear Cross Sections in Monte Carlo Methods for Medical Physics Applications Christopher T. Myers 1 Georgia Institute of Technology Bernadette L. Kirk 2 Luiz C. Leal 2 Oak Ridge
More information1 The Cathode Rays experiment is associated. with: Millikan A B. Thomson. Townsend. Plank Compton
1 The Cathode Rays experiment is associated with: A B C D E Millikan Thomson Townsend Plank Compton 1 2 The electron charge was measured the first time in: A B C D E Cathode ray experiment Photoelectric
More informationParticle Detectors. Summer Student Lectures 2010 Werner Riegler, CERN, History of Instrumentation History of Particle Physics
Particle Detectors Summer Student Lectures 2010 Werner Riegler, CERN, werner.riegler@cern.ch History of Instrumentation History of Particle Physics The Real World of Particles Interaction of Particles
More informationNuclear Physics and Astrophysics
Nuclear Physics and Astrophysics PHY-30 Dr. E. Rizvi Lecture 4 - Detectors Binding Energy Nuclear mass MN less than sum of nucleon masses Shows nucleus is a bound (lower energy) state for this configuration
More informationFall Quarter 2010 UCSB Physics 225A & UCSD Physics 214 Homework 1
Fall Quarter 2010 UCSB Physics 225A & UCSD Physics 214 Homework 1 Problem 2 has nothing to do with what we have done in class. It introduces somewhat strange coordinates called rapidity and pseudorapidity
More informationPhotonuclear Reactions and Nuclear Transmutation. T. Tajima 1 and H. Ejiri 2
Draft Photonuclear Reactions and Nuclear Transmutation T. Tajima 1 and H. Ejiri 2 1) Kansai JAERI 2) JASRI/SPring-8, Mikazuki-cho, Sayou-gun, Hyougo, 679-5198 JAPAN Abstract Photonuclear reactions are
More information