Basic MonteCarlo concepts
|
|
- Judith Armstrong
- 5 years ago
- Views:
Transcription
1 Departamento de Física Atómica, Molecular y Nuclear Universidad de Sevilla Introductory tutorial to Geant4 ITA/IEav, São José dos Campos, SP, Brazil July August 1 st, 2014
2 Index 1 Probability distributions Binomial Poisson Normal
3 Index 1 Probability distributions Binomial Poisson Normal 2 Generation of a continuous random variable Inverse transform method Acceptance-rejection method
4 Index 1 Probability distributions Binomial Poisson Normal 2 Generation of a continuous random variable Inverse transform method Acceptance-rejection method 3 The MonteCarlo method Central Limit Theorem Concept
5 Index The MonteCarlo method 1 Probability distributions Binomial Poisson Normal 2 Generation of a continuous random variable Inverse transform method Acceptance-rejection method 3 The MonteCarlo method Central Limit Theorem Concept Normal distribution
6 Normal distribution It is a discrete distribution which was introduced by the swiss mathematician Jacob Bernoulli ( ). Let s consider a discrete variable, which has two possible results, with probabilities p and q = 1 p. The probability to obtain x results of the first type after N measurements of the variable is P(x) = N! x!(n x)! px q (N x) with the following mean value and standard deviation x = Np σ = Npq = Np(1 p) Therefore it s a discrete distribution defined by two parameters: N and p.
7 Index The MonteCarlo method 1 Probability distributions Binomial Poisson Normal 2 Generation of a continuous random variable Inverse transform method Acceptance-rejection method 3 The MonteCarlo method Central Limit Theorem Concept Normal distribution
8 Normal distribution It was introduced by Siméon Denis Poisson ( ) during its studies on criminal matter. It describes a discrete random variable x which tallies the number of occurrences of a given type which take place during a time interval. It s the limiting case of a binomial distribution when N, p 0, being x = Np its only parameter. The Poisson distribution is discrete and is fully characterized by it s mean value x. p x (x) = x x x! e x
9 Normal distribution Let s consider the disintegration of a radiactive atomic nucleus with probability per unit time λ. Ultimately, we will resort to differential calculus, therefore being dt the differential time interval and t the total elapsed time. t = N dt N dt 0 each sampling of the random variable (disintegration: does it happen? or doesn t it happen?) corresponds to a time interval dt with probability of success (i.e. it happens) which evidently fulfills p = λ dt p 0
10 Normal distribution In these conditions ( N and p 0 ) the binomial distribution tends to the Poisson distribution (also known as Law of rare events). The condition to satisfy this statistics is that the probability of success in each sampling be small, which is obviously accomplished since each samplig corresponts to a differential time interval dt. Let p x (t) be the probability of x successes (disintegrations) in a given interval [0, t]. In this case from which it follows dp x dt which solution is p x (t + dt) = p x (t)(1 λdt) + p x 1 (t)λdt = p x(t + dt) p x (t) dt p x (t) = (λt)x e λt x! = λ [p x 1 (t) p x (t)]
11 Evidently the average number of successes is Normal distribution x = λt which leads to the the most usual representation p x (x) = x x It s straightforward to demonstrate x! e x x 2 = x 2 + x i.e. σ = x. The Poisson distribution is discrete and fully characterized by its mean value x. The rare events condition has been implicitly applied when the probability of simulataneous occurrence of two disintegrations in the same time interval dt has been neglected.
12 Normal distribution In our differential treatment this condition is met, since p = λdt 0 For this reason the Poisson distribution plays a key rôle in the simulations by the MonteCarlo method of processes in which the probability of occurrence per unit time or unit lenght is known (for instance, the radiactive disintegration and processes taking place during radiation transport through matter).
13 Index The MonteCarlo method 1 Probability distributions Binomial Poisson Normal 2 Generation of a continuous random variable Inverse transform method Acceptance-rejection method 3 The MonteCarlo method Central Limit Theorem Concept Normal distribution
14 Normal distribution The normal distribution (also called Gauss distribution) p x,σ (x) = 1 e (x x)2 2σ 2 2πσ is characterized by two parameters: the mean value x and the standard deviation σ. The Poisson distribution tends to the normal for large values of the mean value x, therefore becoming : p x (x) = x x x! e x 1 e (x x)2 2x 2πx The standard normal distribution corresponds to the situation where x = 0 and σ = 1 p estndar (x) = 1 2π e x
15 Calculations with Mathematica Normal distribution Figura: x = 5 Figura: x = 100
16 Index The MonteCarlo method 1 Probability distributions Binomial Poisson Normal 2 Generation of a continuous random variable Inverse transform method Acceptance-rejection method 3 The MonteCarlo method Central Limit Theorem Concept
17 The goal is to generate values of ξ with probability density p(x), defined in the interval [x 1, x 2 ]. It can be shown that ξ x 1 p(x)dx = γ (1) where γ is uniformly in [0, 1]. Therefore, provided a random number generator is available, the problem reduces to solve the above integral equation.
18 PROOF let s define x y(x) = p(x)dx x 1 which is monotonically increasing (y (x) = p(x) > 0 ) with y(x 1 ) = 0 and y(x 2 ) = 1. It follows that y(x) = valor has a unique solution If x (a, b) y(x) [y(a), y(b)] P [a < ξ < b] = P [y(a) < γ < y(b)] Since γ is uniformly distributed: y(b) y(a) P [y(a) < γ < y(b)] = y(x 2 ) y(x 1 ) = b y(b) y(a) = p(x)dx a Therefore P [a < ξ < b] = b a p(x)dx
19 Index The MonteCarlo method 1 Probability distributions Binomial Poisson Normal 2 Generation of a continuous random variable Inverse transform method Acceptance-rejection method 3 The MonteCarlo method Central Limit Theorem Concept
20 The inverse transform method may be unpractical if: The integral of p(x) cannot be solved analitically. The probability density p(x) is known grafically or is tabulated. Method: let s consider ξ defined in [x 1, x 2 ] and its probability density p(x) M 0, then the values of ξ can be sampled as follows: 1 A point Γ of coordinates (η, η ) is built from two values γ and γ of the random variable uniformly distributed in [0, 1] : η = x 1 + γ (x 2 x 1 ) η = γ M 0 2 If the point Γ is below the curve y = p(x), then ξ = η. Otherwise a new pair of values (γ, γ ) is sampled and the procedure is repeated until the above condition is fulfilled.
21 Index The MonteCarlo method 1 Probability distributions Binomial Poisson Normal 2 Generation of a continuous random variable Inverse transform method Acceptance-rejection method 3 The MonteCarlo method Central Limit Theorem Concept
22 Generation of a gaussian distribution (inverse transform method) Figure: Inverse transform method. CPU time Mathematica: 103 seconds
23 Generation of a gaussian distribution (acceptance-reject method) Figure: Acceptance-rejection method. CPU time Mathematica: 10 seconds
24 2D uniform distribution 1 = R 2π P(r, φ)drdφ = R 2π Uniform distribution: p(r, φ) r drdφ p(r, φ) = C C = 1 πr 2 1 = C R 2π r drdφ 0 0 P(r, φ)drdφ = r drdφ π R 2 Figura: Polar coordinates The distribution in r: [ 2π ] P r (r)dr = P(r, φ)dφ dr = 2rdr 0 R 2 The distribution in φ: [ R ] P φ(φ)dφ = P(r, φ)dr dφ = dφ 0 2π
25 2D uniform distribution (cont.) Therefore for sampling the r and φ values, being γ a continuous uniformly distributed random variable in [0, 1] (random number generator in mathematical libraries available in all major scientific languajes, as FORTRAN and C++): r 0 φ 0 P r (r )dr = γ = 2 r R 2 P φ (φ )dφ = γ = 1 2π 0 φ 0 r dr = r 2 R 2 r = R γ dφ = φ 2π φ = 2πγ
26 Mathematica calculations Figura: Uniform 2D distribution Figura: Pseudouniform 2D distribution
27 3D uniform distribution Figura: Spherical coordinates The distribution in r: [ π ] 2π P r (r)dr = P(r, θ, φ)dθdφ dr 0 0 = 3r 2 dr R 3 The distribution in θ: [ R ] 2π P θ(θ)dr = P(r, θ, φ)drdφ dθ 0 0 = sen(θ)dθ 2 1 = R π 2π 0 0 R π 2π P(r, θ, φ)drdθdφ = p(r, θ, φ) r 2 sen(θ) drdθdφ Uniform distribution: p(r, θ, φ) = C 1 = C R π 2π r 2 sen(θ) drdθdφ C = 3 4πR 3 P(r, θ, φ)drdθdφ = 3r 2 sen(θ) drdθdφ 4π R 3 The distribution in φ: [ R ] π P φ(φ)dφ = P(r, θ, φ)drdθ dφ 0 0 = dφ 2π
28 3D uniform distribution (cont.) Therefore, for sampling the r, θ and φ values, being γ a continuous uniformly distributed random variable in [0, 1] (random number generator generador in mathematical libraries available in all major scientific languajes, as FORTRAN and C++): r 0 P r (r )dr = γ = 3 r R 3 r 2 dr = r 3 0 R 3 r = Rγ1/3 θ 0 P θ (θ )dθ = γ = 1 2 φ 0 θ 0 sen(θ )dθ = P φ (φ )dφ = γ = 1 2π 1 2(1 cos(θ)) θ = acos(2γ 1) φ 0 dφ = φ 2π φ = 2πγ
29 Mathematica calculations Figura: Uniform distribution Figura: Pseudouniform distribution
30 Index The MonteCarlo method 1 Probability distributions Binomial Poisson Normal 2 Generation of a continuous random variable Inverse transform method Acceptance-rejection method 3 The MonteCarlo method Central Limit Theorem Concept
31 Theorem The MonteCarlo method Given N independent, equal, random variables: ξ 1, ξ 2,...ξ N ξ 1 = ξ 2 =... = ξ N = m ξ 1 = ξ 2 =... = ξ N = σ and being ρ their sum N N ρ N = ξ i ρ N = ξ i = Nm i=1 i=1 ρ N = Nσ for large N, ρ is a random normal variable (i.e. distributed acoording to Gauss law) with the above parameters P[a ρ N b] b a 1 p ρn (x) dx con p ρn (x) = e (x ρ N ) 2 ρ 2 N 2π ρn 2
32 Index The MonteCarlo method 1 Probability distributions Binomial Poisson Normal 2 Generation of a continuous random variable Inverse transform method Acceptance-rejection method 3 The MonteCarlo method Central Limit Theorem Concept
33 It is a direct application of the Central Limit Theorem. The goal is to calculate a magnitude m (for instance, the number of neutrons which get through a given material thickness) The method consists in: 1 Take N equal and independent random variables (which ultimately are the N Monte Carlo samplings or histories). ξ 1, ξ 2,...ξ N ξ 1 = ξ 2 =... = ξ N = m 2 Calculated their mean value: ξ 1 = ξ 2 =... = ξ N = σ ξ = ρ N N N = i=1 ξi N According to the Central Limit Theorem, it follows: P[Nm 3σ N ρ N Nm + 3σ N] P[m 3 σ N ρ N N m + 3 σ N ] P[ 1 i=n ξ i m < 3 σ ] N i=1 N
34 Basically, a MonteCarlo calculation is a numerical experiment which uses random number generators to simulate the physical quantities, which are sampled according to probability distributions (Physics models)
35 Evidently, the larger the number N of samplings, the closer will be the experimental average (by MonteCarlo) to the theoretical mean value, which strictly would ask for N. In a typical application a huge number of events takes place, for instance a beam of particles with intensity p/s impinges on a a macroscopic target, which contains of the order of the Avogadro number ( ) of scattering centers. As a direct consequence of the Central Limit Theorem the quality of the Montecarlo simulation is proportional to N, since the relative error goes as σ/ N, where σ is the standad deviation; from this follows the limitation coming from the available computational power.
36 Index The MonteCarlo method 1 Probability distributions Binomial Poisson Normal 2 Generation of a continuous random variable Inverse transform method Acceptance-rejection method 3 The MonteCarlo method Central Limit Theorem Concept
37 Calculation of π The MonteCarlo method The numerical integration by geometrical (acceptance-rejection) method was at the origin of the MonteCarlo calculations. We will start with a simple geometrical example for calculating the value of π. From the figure, it s obvious that the area of the cicle within one of the quadrants is r=1 es π/4.
38 Let s imagine that we throw darts against the figure. After a high enough number of throws, the number of hits (impacts within the circle) will be approx. proportional to the area of the circle. In other words, the ratio of hits over total number of throws will be π/4. If we actually carry ouit this experiment, we woud verify that a very large number of throws (well above 1000) is needed in order to reach an acceptable estimation of π. For this very reason it s much better to carry the experiment virtually with a computer by using use a random number generator withing this very simple algorithm. The quality of the estimation will depend on the quality of the random number generator and the number of iterations (virtual throws).
39 Distribution law from Poisson statistics If Σ is the probability per unit lenght for a certain process to occur, that leads to the Poisson distribution of the number of occurrences which take place within the length interval [0, l]: p x (l) = (Σl)x e Σl x! Nevertheless, in a numerical simulation by the MonteCarlo method what we do really need is the occurrence probability of the first success within the interval [l, l + dl]: Distribution law of the lenght where I 1 (l)dl = p 0 (l) Σdl = e Σl Σdl p 0 (l) is the non-occurrence probability within the inverval [0, l]. Σdl is the occurrence probability of a single success within the interval [l, l + dl].
40 Distribution law from Poisson statistics (cont.) By applying the inverse transform method: Sampling algorithm l = 1 Σ lnγ where γ is a random variable uniformly distributed in [0, 1]. The same applies for distribution of occurrence times (for instance, in radiactive disintegration) provided that the disintegration constant (probability per unit time) is known.
41 Neutron transport The MonteCarlo method In this simple example only two processes are considered: elastic scattering absorption Therefore there are 3 possibilities after transportation through the slab: reflection (probability: p ) absorption (probability: p 0 ) transmission (probability: p + ) The macroscopic cross section is defined as the probability per unit path length of a neutron experiencing a collision. Σ = σρ, where σ is the cross section of the process and ρ is the density of targets. Figura: Neutrons through a slab Σ = 1 λ = Σ abs + Σ disp (2) being λ the mean free path of a neutron in the slab material.
42 sampling of λ λ: random variable with distribution function sampling of θ θ: random variable with distribution function: p(λ) = 1 λ e 1 λ λ = Σe Σλ By applying the inverse transform method: λ = 1 Σ lnγ P θ (θ)dr = sen(θ)dθ 2 By applying the inverse transform method: θ = acos(1 2γ) Method (after the k-th collision the neutron is in position z k ) 1 Sampling of λ: λ k = 1 Σ lnγ 2 Sampling of θ: θ k = acos(2γ 1) 3 The new position is calculated: z k+1 = z k + λ k cos(θ k ) 4 If z k+1 > h transmitted neutron N + = N If z k+1 < 0 reflected neutron N = N If 0 < z k+1 < h the neutron is still inside: If γ < Σ abs Σ N0 = N (captured neutron) Otherwise return to point 1 and repeat the process Finally: p + = N + /N ; p 0 = N 0 /N ; p = N /N.
43 Potential problem A problem may arise if after N repetitions few particles (if any) survive (bad statistics). Solution: weighting method 1 Instead of considering individual neutrons, an initial bunch made of ω 0 neutrons is created In a given interaction, (k), in average 2 ω k 1 Σ abs Σ Σabs ω k 1 Σ Σdisp ω k 1 Σ are absorbed. is are scattered. is added to N0 (absorbed particles) and the process is Σ repeated with the remaining bunch, with weight ω disp k 1 Σ, moving in the same direction. Σ All former formulae apply, susbtituting ω k+1 = ω disp k Σ as the part of the bunch which survives after each collision (ω k is the weight of the neutron. the initial weight is ω 0 = 1. In this way a trajectory never ends in an absorption: there is always a remnant, which is a fractionary quantity, obviously proportional to the transmission probability. In these circumstances, this method noticeably improves the efficiency of the simulation.
Monte 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 informationPHYS 352. Charged Particle Interactions with Matter. Intro: Cross Section. dn s. = F dω
PHYS 352 Charged Particle Interactions with Matter Intro: Cross Section cross section σ describes the probability for an interaction as an area flux F number of particles per unit area per unit time dσ
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 informationElastic scattering. Elastic scattering
Elastic scattering Now we have worked out how much energy is lost when a neutron is scattered through an angle, θ We would like to know how much energy, on average, is lost per collision In order to do
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 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 informationInteraction theory Photons. Eirik Malinen
Interaction theory Photons Eirik Malinen Introduction Interaction theory Dosimetry Radiation source Ionizing radiation Atoms Ionizing radiation Matter - Photons - Charged particles - Neutrons Ionizing
More informationME 595M: Monte Carlo Simulation
ME 595M: Monte Carlo Simulation T.S. Fisher Purdue University 1 Monte Carlo Simulation BTE solutions are often difficult to obtain Closure issues force the use of BTE moments Inherent approximations Ballistic
More informationLectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 1. Title: Neutron Life Cycle
Lectures on Nuclear Power Safety Lecture No 1 Title: Neutron Life Cycle Department of Energy Technology KTH Spring 2005 Slide No 1 Outline of the Lecture Infinite Multiplication Factor, k Four Factor Formula
More informationRandom variables. DS GA 1002 Probability and Statistics for Data Science.
Random variables DS GA 1002 Probability and Statistics for Data Science http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall17 Carlos Fernandez-Granda Motivation Random variables model numerical quantities
More informationMonte Carlo Sampling
Monte Carlo Sampling Sampling from PDFs: given F(x) in analytic or tabulated form, generate a random number ξ in the range (0,1) and solve the equation to get the randomly sampled value X X ξ = F(x)dx
More informationScattering of an α Particle by a Radioactive Nucleus
EJTP 3, No. 1 (6) 93 33 Electronic Journal of Theoretical Physics Scattering of an α Particle by a Radioactive Nucleus E. Majorana Written 198 published 6 Abstract: In the following we reproduce, translated
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 1 (2/3/04) Overview -- Interactions, Distributions, Cross Sections, Applications
.54 Neutron Interactions and Applications (Spring 004) Chapter 1 (/3/04) Overview -- Interactions, Distributions, Cross Sections, Applications There are many references in the vast literature on nuclear
More informationOlber s Paradox Solved
Olber s Paradox Solved by Ronald Koster (http://www.ronaldkoster.net) ersion.2, 22-7-3 Introduction Why is the night sky dark? According to Olber s paradox the night (and day) sky should be infinitely
More informationB. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 2015 Sept.-Dec.
2: Fission and Other Neutron Reactions B. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 2015 Sept.-Dec. 2015 September 1 Contents Concepts: Fission and other
More informationUNIVERSITY OF CAPE TOWN. EEE4101F / EEE4103F: Basic Nuclear Physics Problem Set 04. Due 12:00 (!) Wednesday 8 April 2015
UNIVERSITY OF CAPE TOWN EEE4101F / EEE4103F: Basic Nuclear Physics Problem Set 04 Due 12:00 (!) Wednesday 8 April 2015 1) Explain the concept of a cross-section, using language that a high-school teacher
More informationAtmospheric Refraction. And Related Phenomena
Atmospheric Refraction And Related Phenomena Overview 1. Refractive invariant 2. The atmosphere and its index of refraction 3. Ray tracing through atmosphere 4. Sunset distortions 5. Green flash Spherical
More informationStatistics, Data Analysis, and Simulation SS 2013
Statistics, Data Analysis, and Simulation SS 213 8.128.73 Statistik, Datenanalyse und Simulation Dr. Michael O. Distler Mainz, 23. April 213 What we ve learned so far Fundamental
More informationToday: Fundamentals of Monte Carlo
Today: Fundamentals of Monte Carlo What is Monte Carlo? Named at Los Alamos in 940 s after the casino. Any method which uses (pseudo)random numbers as an essential part of the algorithm. Stochastic - not
More informationProbability and Measure
Probability and Measure Robert L. Wolpert Institute of Statistics and Decision Sciences Duke University, Durham, NC, USA Convergence of Random Variables 1. Convergence Concepts 1.1. Convergence of Real
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 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 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 informationApplied Statistics I
Applied Statistics I (IMT224β/AMT224β) Department of Mathematics University of Ruhuna A.W.L. Pubudu Thilan Department of Mathematics University of Ruhuna Applied Statistics I(IMT224β/AMT224β) 1/158 Chapter
More informationDS-GA 1002 Lecture notes 2 Fall Random variables
DS-GA 12 Lecture notes 2 Fall 216 1 Introduction Random variables Random variables are a fundamental tool in probabilistic modeling. They allow us to model numerical quantities that are uncertain: the
More informationCross-Sections for Neutron Reactions
22.05 Reactor Physics Part Four Cross-Sections for Neutron Reactions 1. Interactions: Cross-sections deal with the measurement of interactions between moving particles and the material through which they
More informationProbability Distributions - Lecture 5
Probability Distributions - Lecture 5 1 Introduction There are a number of mathematical models of probability density functions that represent the behavior of physical systems. In this lecture we explore
More informationToday: Fundamentals of Monte Carlo
Today: Fundamentals of Monte Carlo What is Monte Carlo? Named at Los Alamos in 1940 s after the casino. Any method which uses (pseudo)random numbers as an essential part of the algorithm. Stochastic -
More information7 Convergence in R d and in Metric Spaces
STA 711: Probability & Measure Theory Robert L. Wolpert 7 Convergence in R d and in Metric Spaces A sequence of elements a n of R d converges to a limit a if and only if, for each ǫ > 0, the sequence a
More informationMATH 280 Multivariate Calculus Fall Integration over a curve
dr dr y MATH 28 Multivariate Calculus Fall 211 Integration over a curve Given a curve C in the plane or in space, we can (conceptually) break it into small pieces each of which has a length ds. In some
More information2. The Steady State and the Diffusion Equation
2. The Steady State and the Diffusion Equation The Neutron Field Basic field quantity in reactor physics is the neutron angular flux density distribution: Φ( r r, E, r Ω,t) = v(e)n( r r, E, r Ω,t) -- distribution
More informationBMIR Lecture Series on Probability and Statistics Fall, 2015 Uniform Distribution
Lecture #5 BMIR Lecture Series on Probability and Statistics Fall, 2015 Department of Biomedical Engineering and Environmental Sciences National Tsing Hua University s 5.1 Definition ( ) A continuous random
More informationDecays and Scattering. Decay Rates Cross Sections Calculating Decays Scattering Lifetime of Particles
Decays and Scattering Decay Rates Cross Sections Calculating Decays Scattering Lifetime of Particles 1 Decay Rates There are THREE experimental probes of Elementary Particle Interactions - bound states
More informationStatistical Methods in Particle Physics
Statistical Methods in Particle Physics Lecture 3 October 29, 2012 Silvia Masciocchi, GSI Darmstadt s.masciocchi@gsi.de Winter Semester 2012 / 13 Outline Reminder: Probability density function Cumulative
More informationToday: Fundamentals of Monte Carlo
Today: Fundamentals of Monte Carlo What is Monte Carlo? Named at Los Alamos in 1940 s after the casino. Any method which uses (pseudo)random numbers as an essential part of the algorithm. Stochastic -
More informationIntroduction to Elementary Particle Physics I
Physics 56400 Introduction to Elementary Particle Physics I Lecture 2 Fall 2018 Semester Prof. Matthew Jones Cross Sections Reaction rate: R = L σ The cross section is proportional to the probability of
More informationSolutions: Homework 7
Solutions: Homework 7 Ex. 7.1: Frustrated Total Internal Reflection a) Consider light propagating from a prism, with refraction index n, into air, with refraction index 1. We fix the angle of incidence
More informationLecture 15 Random variables
Lecture 15 Random variables Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn No.1
More informationProbability and Distributions
Probability and Distributions What is a statistical model? A statistical model is a set of assumptions by which the hypothetical population distribution of data is inferred. It is typically postulated
More informationfrom which follow by application of chain rule relations y = y (4) ˆL z = i h by constructing θ , find also ˆL x ˆL y and
9 Scattering Theory II 9.1 Partial wave analysis Expand ψ in spherical harmonics Y lm (θ, φ), derive 1D differential equations for expansion coefficients. Spherical coordinates: x = r sin θ cos φ (1) y
More informationAssignment 4. u n+1 n(n + 1) i(i + 1) = n n (n + 1)(n + 2) n(n + 2) + 1 = (n + 1)(n + 2) 2 n + 1. u n (n + 1)(n + 2) n(n + 1) = n
Assignment 4 Arfken 5..2 We have the sum Note that the first 4 partial sums are n n(n + ) s 2, s 2 2 3, s 3 3 4, s 4 4 5 so we guess that s n n/(n + ). Proving this by induction, we see it is true for
More informationExploring Monte Carlo Methods
Exploring Monte Carlo Methods William L Dunn J. Kenneth Shultis AMSTERDAM BOSTON HEIDELBERG LONDON NEW YORK OXFORD PARIS SAN DIEGO SAN FRANCISCO SINGAPORE SYDNEY TOKYO ELSEVIER Academic Press Is an imprint
More informationProbability Density Functions
Probability Density Functions Probability Density Functions Definition Let X be a continuous rv. Then a probability distribution or probability density function (pdf) of X is a function f (x) such that
More informationA possible wavefunction for two particles (1 and 2) moving in one dimension is. πx 2 x 0 2x
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.044 Statistical Physics I Spring Term 03 Problem Set # Due in hand-in box by :40 PM, Wednesday, February 0 Problem : Two Quantum Particles A possible
More informationReactor Physics: Basic Definitions and Perspectives. Table of Contents
Reactor Physics - Basic Definitions and Perspectives Reactor Physics: Basic Definitions and Perspectives prepared by Wm. J. Garland, Professor, Department of Engineering Physics, McMaster University, Hamilton,
More informationChing-Han Hsu, BMES, National Tsing Hua University c 2015 by Ching-Han Hsu, Ph.D., BMIR Lab. = a + b 2. b a. x a b a = 12
Lecture 5 Continuous Random Variables BMIR Lecture Series in Probability and Statistics Ching-Han Hsu, BMES, National Tsing Hua University c 215 by Ching-Han Hsu, Ph.D., BMIR Lab 5.1 1 Uniform Distribution
More informationElastic Scattering. R = m 1r 1 + m 2 r 2 m 1 + m 2. is the center of mass which is known to move with a constant velocity (see previous lectures):
Elastic Scattering In this section we will consider a problem of scattering of two particles obeying Newtonian mechanics. The problem of scattering can be viewed as a truncated version of dynamic problem
More informationEFFECTS OF PERMEABILITY ON SOUND ABSORPTION AND SOUND INSULATION PERFORMANCE OF ACOUSTIC CEILING PANELS
EFFECTS OF PERMEABILITY ON SOUND ABSORPTION AND SOUND INSULATION PERFORMANCE OF ACOUSTIC CEILING PANELS Kento Hashitsume and Daiji Takahashi Graduate School of Engineering, Kyoto University email: kento.hashitsume.ku@gmail.com
More informationMath 221 Examination 2 Several Variable Calculus
Math Examination Spring Instructions These problems should be viewed as essa questions. Before making a calculation, ou should explain in words what our strateg is. Please write our solutions on our own
More informationSTOCHASTIC & DETERMINISTIC SOLVERS
STOCHASTIC & DETERMINISTIC SOLVERS Outline Spatial Scales of Optical Technologies & Mathematical Models Prototype RTE Problems 1-D Transport in Slab Geometry: Exact Solution Stochastic Models: Monte Carlo
More informationLecture 11: Probability, Order Statistics and Sampling
5-75: Graduate Algorithms February, 7 Lecture : Probability, Order tatistics and ampling Lecturer: David Whitmer cribes: Ilai Deutel, C.J. Argue Exponential Distributions Definition.. Given sample space
More informationChapter 4. Repeated Trials. 4.1 Introduction. 4.2 Bernoulli Trials
Chapter 4 Repeated Trials 4.1 Introduction Repeated indepentent trials in which there can be only two outcomes are called Bernoulli trials in honor of James Bernoulli (1654-1705). As we shall see, Bernoulli
More informationApplied Nuclear Physics (Fall 2006) Lecture 19 (11/22/06) Gamma Interactions: Compton Scattering
.101 Applied Nuclear Physics (Fall 006) Lecture 19 (11//06) Gamma Interactions: Compton Scattering References: R. D. Evans, Atomic Nucleus (McGraw-Hill New York, 1955), Chaps 3 5.. W. E. Meyerhof, Elements
More informationChap 2.1 : Random Variables
Chap 2.1 : Random Variables Let Ω be sample space of a probability model, and X a function that maps every ξ Ω, toa unique point x R, the set of real numbers. Since the outcome ξ is not certain, so is
More informationThe next three lectures will address interactions of charged particles with matter. In today s lecture, we will talk about energy transfer through
The next three lectures will address interactions of charged particles with matter. In today s lecture, we will talk about energy transfer through the property known as stopping power. In the second lecture,
More informationHow Electronics Started! And JLab Hits the Wall!
How Electronics Started! And JLab Hits the Wall! In electronics, a vacuum diode or tube is a device used to amplify, switch, otherwise modify, or create an electrical signal by controlling the movement
More informationProbability Density Functions
Statistical Methods in Particle Physics / WS 13 Lecture II Probability Density Functions Niklaus Berger Physics Institute, University of Heidelberg Recap of Lecture I: Kolmogorov Axioms Ingredients: Set
More informationNeutron interactions and dosimetry. Eirik Malinen Einar Waldeland
Neutron interactions and dosimetry Eirik Malinen Einar Waldeland Topics 1. Neutron interactions 1. Scattering 2. Absorption 2. Neutron dosimetry 3. Applications The neutron Uncharged particle, mass close
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 informationChapter 2: Random Variables
ECE54: Stochastic Signals and Systems Fall 28 Lecture 2 - September 3, 28 Dr. Salim El Rouayheb Scribe: Peiwen Tian, Lu Liu, Ghadir Ayache Chapter 2: Random Variables Example. Tossing a fair coin twice:
More informationProbability Distributions
The ASTRO509-3 dark energy puzzle Probability Distributions I have had my results for a long time: but I do not yet know how I am to arrive at them. Johann Carl Friedrich Gauss 1777-1855 ASTR509 Jasper
More informationKnown probability distributions
Known probability distributions Engineers frequently wor with data that can be modeled as one of several nown probability distributions. Being able to model the data allows us to: model real systems design
More informationMonte 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 informationin Computer Simulations for Bioinformatics
Bi04a_ Unit 04a: Stochastic Processes and their Applications in Computer Simulations for Bioinformatics Stochastic Processes and their Applications in Computer Simulations for Bioinformatics Basic Probability
More informationNotes 19 Gradient and Laplacian
ECE 3318 Applied Electricity and Magnetism Spring 218 Prof. David R. Jackson Dept. of ECE Notes 19 Gradient and Laplacian 1 Gradient Φ ( x, y, z) =scalar function Φ Φ Φ grad Φ xˆ + yˆ + zˆ x y z We can
More informationFigure 10.1: Recording when the event E occurs
10 Poisson Processes Let T R be an interval. A family of random variables {X(t) ; t T} is called a continuous time stochastic process. We often consider T = [0, 1] and T = [0, ). As X(t) is a random variable
More informationM378K In-Class Assignment #1
The following problems are a review of M6K. M7K In-Class Assignment # Problem.. Complete the definition of mutual exclusivity of events below: Events A, B Ω are said to be mutually exclusive if A B =.
More informationMcGill University April Calculus 3. Tuesday April 29, 2014 Solutions
McGill University April 4 Faculty of Science Final Examination Calculus 3 Math Tuesday April 9, 4 Solutions Problem (6 points) Let r(t) = (t, cos t, sin t). i. Find the velocity r (t) and the acceleration
More informationNumerical Methods for Data Analysis
Michael O. Distler distler@uni-mainz.de Bosen (Saar), August 29 - September 3, 2010 Fundamentals Probability distributions Expectation values, error propagation Parameter estimation Regression analysis
More informationThe arithmetic geometric mean (Agm)
The arithmetic geometric mean () Pictures by David Lehavi This is part of exercise 5 question 4 in myc> calculus class Fall 200: a = c a n+ =5a n 2 lim n an = a = a an + bn a n+ = 2 b = b b n+ = a nb n
More informationMCRT L10: Scattering and clarification of astronomy/medical terminology
MCRT L10: Scattering and clarification of astronomy/medical terminology What does the scattering? Shape of scattering Sampling from scattering phase functions Co-ordinate frames Refractive index changes
More informationProbability Distributions Columns (a) through (d)
Discrete Probability Distributions Columns (a) through (d) Probability Mass Distribution Description Notes Notation or Density Function --------------------(PMF or PDF)-------------------- (a) (b) (c)
More informationLECTURE NOTES FYS 4550/FYS EXPERIMENTAL HIGH ENERGY PHYSICS AUTUMN 2013 PART I A. STRANDLIE GJØVIK UNIVERSITY COLLEGE AND UNIVERSITY OF OSLO
LECTURE NOTES FYS 4550/FYS9550 - EXPERIMENTAL HIGH ENERGY PHYSICS AUTUMN 2013 PART I PROBABILITY AND STATISTICS A. STRANDLIE GJØVIK UNIVERSITY COLLEGE AND UNIVERSITY OF OSLO Before embarking on the concept
More informationChapter 1: Useful definitions
Chapter 1: Useful definitions 1.1. Cross-sections (review) The Nuclear and Radiochemistry class listed as a prerequisite is a good place to start. The understanding of a cross-section being fundamentai
More informationThe effect of the spectator charge on the charged pion spectra in peripheral ultrarelativistic heavy-ion collisions
The effect of the spectator charge on the charged pion spectra in peripheral ultrarelativistic heavy-ion collisions Antoni Szczurek and Andrzej Rybicki INSTITUTE OF NUCLEAR PHYSICS POLISH ACADEMY OF SCIENCES
More informationSolve for an unknown rate of change using related rates of change.
Objectives: Solve for an unknown rate of change using related rates of change. 1. Draw a diagram. 2. Label your diagram, including units. If a quantity in the diagram is not changing, label it with a number.
More informationModern Discrete Probability Branching processes
Modern Discrete Probability IV - Branching processes Review Sébastien Roch UW Madison Mathematics November 15, 2014 1 Basic definitions 2 3 4 Galton-Watson branching processes I Definition A Galton-Watson
More informationSTAT/MATH 395 A - PROBABILITY II UW Winter Quarter Moment functions. x r p X (x) (1) E[X r ] = x r f X (x) dx (2) (x E[X]) r p X (x) (3)
STAT/MATH 395 A - PROBABILITY II UW Winter Quarter 07 Néhémy Lim Moment functions Moments of a random variable Definition.. Let X be a rrv on probability space (Ω, A, P). For a given r N, E[X r ], if it
More information10 BIVARIATE DISTRIBUTIONS
BIVARIATE DISTRIBUTIONS After some discussion of the Normal distribution, consideration is given to handling two continuous random variables. The Normal Distribution The probability density function f(x)
More informationFigure 21:The polar and Cartesian coordinate systems.
Figure 21:The polar and Cartesian coordinate systems. Coordinate systems in R There are three standard coordinate systems which are used to describe points in -dimensional space. These coordinate systems
More informationLecture 2 Binomial and Poisson Probability Distributions
Binomial Probability Distribution Lecture 2 Binomial and Poisson Probability Distributions Consider a situation where there are only two possible outcomes (a Bernoulli trial) Example: flipping a coin James
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 informationPhysics 115A: Statistical Physics
Physics 115A: Statistical Physics Prof. Clare Yu email: cyu@uci.edu phone: 949-824-6216 Office: 210E RH Fall 2013 LECTURE 1 Introduction So far your physics courses have concentrated on what happens to
More informationElectrostatics. Chapter Maxwell s Equations
Chapter 1 Electrostatics 1.1 Maxwell s Equations Electromagnetic behavior can be described using a set of four fundamental relations known as Maxwell s Equations. Note that these equations are observed,
More informationProbability. Carlo Tomasi Duke University
Probability Carlo Tomasi Due University Introductory concepts about probability are first explained for outcomes that tae values in discrete sets, and then extended to outcomes on the real line 1 Discrete
More informationNumerical Methods I Monte Carlo Methods
Numerical Methods I Monte Carlo Methods Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Dec. 9th, 2010 A. Donev (Courant Institute) Lecture
More informationOptical Imaging Chapter 5 Light Scattering
Optical Imaging Chapter 5 Light Scattering Gabriel Popescu University of Illinois at Urbana-Champaign Beckman Institute Quantitative Light Imaging Laboratory http://light.ece.uiuc.edu Principles of Optical
More information1 CHAPTER 5 ABSORPTION, SCATTERING, EXTINCTION AND THE EQUATION OF TRANSFER
1 CHAPTER 5 ABSORPTION, SCATTERING, EXTINCTION AND THE EQUATION OF TRANSFER 5.1 Introduction As radiation struggles to make its way upwards through a stellar atmosphere, it may be weakened by absorption
More informationStochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno
Stochastic Processes M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno 1 Outline Stochastic (random) processes. Autocorrelation. Crosscorrelation. Spectral density function.
More informationElectrodynamics Exam Solutions
Electrodynamics Exam Solutions Name: FS 215 Prof. C. Anastasiou Student number: Exercise 1 2 3 4 Total Max. points 15 15 15 15 6 Points Visum 1 Visum 2 The exam lasts 18 minutes. Start every new exercise
More information14 Supernovae (short overview) introduc)on to Astrophysics, C. Bertulani, Texas A&M-Commerce 1
14 Supernovae (short overview) introduc)on to Astrophysics, C. Bertulani, Texas A&M-Commerce 1 The core-collapse of a supernova The core of a pre-supernova is made of nuclei in the iron-mass range A ~
More informationStatistics, Data Analysis, and Simulation SS 2015
Statistics, Data Analysis, and Simulation SS 2015 08.128.730 Statistik, Datenanalyse und Simulation Dr. Michael O. Distler Mainz, 27. April 2015 Dr. Michael O. Distler
More informationGeneral Random Variables
1/65 Chia-Ping Chen Professor Department of Computer Science and Engineering National Sun Yat-sen University Probability A general random variable is discrete, continuous, or mixed. A discrete random variable
More informationContinuous Distributions
Continuous Distributions 1.8-1.9: Continuous Random Variables 1.10.1: Uniform Distribution (Continuous) 1.10.4-5 Exponential and Gamma Distributions: Distance between crossovers Prof. Tesler Math 283 Fall
More informationNumerical integration - I. M. Peressi - UniTS - Laurea Magistrale in Physics Laboratory of Computational Physics - Unit V
Numerical integration - I M. Peressi - UniTS - Laurea Magistrale in Physics Laboratory of Computational Physics - Unit V deterministic methods in 1D equispaced points (trapezoidal, Simpson...), others...
More informationSystem Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models
System Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models Fatih Cavdur fatihcavdur@uludag.edu.tr March 20, 2012 Introduction Introduction The world of the model-builder
More informationAccelerated Observers
Accelerated Observers In the last few lectures, we ve been discussing the implications that the postulates of special relativity have on the physics of our universe. We ve seen how to compute proper times
More informationMTH739U/P: Topics in Scientific Computing Autumn 2016 Week 6
MTH739U/P: Topics in Scientific Computing Autumn 16 Week 6 4.5 Generic algorithms for non-uniform variates We have seen that sampling from a uniform distribution in [, 1] is a relatively straightforward
More informationAdaptive Rejection Sampling with fixed number of nodes
Adaptive Rejection Sampling with fixed number of nodes L. Martino, F. Louzada Institute of Mathematical Sciences and Computing, Universidade de São Paulo, Brazil. Abstract The adaptive rejection sampling
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More information