Decadimento b e bb: codici e tools di calcolo
|
|
- Hilary Campbell
- 6 years ago
- Views:
Transcription
1 Decadimento b e bb: codici e tools di calcolo Il sito nucleonica.net: Decay engine Dosimetry and shielding++ Beta dose rate Gamma spectrum generator++ Wespa++ (peak analysis) Virtual Cloud Chamber Databases: Nucleosynthesis simulator NACRE EXFOR EXPACS Codici di calcolo analitici: Talys Neutrino Oscillation Simulator Codici MC: Metodo MC Generatori decadimento Fluka Geant4
2 Nucleonica.net: decay engine
3 Nucleonica.net: decay engine
4 Nucleonica.net: dosimetry and shielding++
5 Nucleonica.net: dosimetry and shielding++
6 Nucleonica.net: dosimetry and shielding++
7 Nucleonica.net: beta dose rate
8 Nucleonica.net: beta dose rate
9 Nucleonica.net: gamma spectrum generator++
10 Nucleonica.net: gamma spectrum generator++
11 Nucleonica.net: WESPA++ (Web Spectrum Analyzer)
12 Nucleonica.net: WESPA++ (Web Spectrum Analyzer)
13 Nucleonica.net: Virtual Cloud Chamber
14 Nucleonica.net: Virtual Cloud Chamber electron-positron pairs are created through the interaction of 10 MeV gamma photons incident on lead. 3 MeV positrons (red) from the radioactive source are blocked by a lead shield (green).
15 Nucleosynthesis simulator
16 Nucleosynthesis simulator
17 NACRE-II
18 NACRE-II
19 EXFOR
20 EXPACS
21 TALYS
22 NOS: Neutrino oscillation Simulator
23 NOS: Neutrino oscillation Simulator
24 Il metodo montecarlo (MC)
25 Il metodo montecarlo (MC)
26 Il metodo montecarlo (MC)
27 Il metodo montecarlo (MC)
28 Particle transport Particle transport is a typical physical process described by probabilities (cross sections = interaction probabilities per unit distance) Therefore it lends itself naturally to be simulated by Monte Carlo Many applications, especially in high energy physics and medicine, are based on simulations where the history of each particle (trajectory, interactions) is reproduced in detail However in other types of application, typically shielding design, the user is interested only in the expectation values of some quantities (fluence and dose) at some space point or region, which are calculated as solutions of a mathematical equation This equation (the Boltzmann equation), describes the statistical distribution of particles in phase space and therefore does indeed represent a physical stochastic process But in order to estimate the desired expectation values it is not necessary that the Monte Carlo process be identical to it
29 Particle transport Monte Carlo Assumptions made by most MC codes: Static, homogeneous, isotropic, amorphous media and geometry Problems: e.g. moving targets, atmosphere must be represented by discrete layers of uniform density, radioactive decay may take place in a geometry different from that in which the radionuclides were produced Markovian process: the fate of a particle depends only on its actual present properties, not on previous events or histories Particles do not interact with each other Problem: e.g. the Chudakov effect (charges cancelling in e + e pairs) Particles interact with individual electrons / atoms / nuclei / molecules Problem: invalid at low energies Material properties are not affected by particle reactions Problem: e.g. burnup
30 Particle transport Application of Monte Carlo to particle transport and interaction: Each particle is followed on its path through matter At each step the occurrence and outcome of interactions are decided by random selection from the appropriate probability distributions All the secondaries issued from the same primary are stored in a stack or bank and are transported before a new history is started The accuracy and reliability of a Monte Carlo depend on the models or data on which the probability distribution functions are based Statistical accuracy of results depends on the number of histories" Statistical convergence can be accelerated by biasing" techniques.
31 Analog vs non-analog simulation
32 Analogue vs non-analogue What is analogue simulation? Sample using natural probability distribution, N(x) Predicts mean with correct fluctuations Can be inefficient for certain applications What is non-analogue/event biased simulation? Cheat - apply artificial biasing probability distribution, B(x) in place of natural one, N(x) B(x) enhances production of whatever it is that is interesting To get meaningful results, must apply a weight correction Predicts same analogue mean with smaller variance Increases efficiency of the Monte Carlo Doesn t predict correct fluctuations Should be used with care
33 Monte Carlo Flavors: microscopic Microscopic Analog Monte Carlo Uses theoretical models to describe physical processes whenever possible Samples from actual physical phase space distributions Predicts average quantities and all statistical moments of any order Preserves correlations (provided the physics is correct, of course!) Reproduces fluctuations (provided... see above) Is (almost) safe and (sometimes) can be used as a black box" (idem) But: Can be inefficient and converge slowly Can fail to predict contributions due to rare events
34 Monte Carlo Flavors: macroscopic Macroscopic Monte Carlo: Instead of simulating interactions, use parametrizations of the reaction product distributions, obtained from fits to data and extrapolations Fast, especially when reactions are complex Can be more accurate than microscopic MC if the theory contains uncertainties/approximations But: The single probability distribution functions are reproduced, but the correlations among interaction products are not Cannot be extended outside the data range
35 Monte Carlo Flavors: biased Biased Monte Carlo: samples from artificial distributions, and applies a weight to the particles to correct for the bias (similar to an integration by a change of variable) predicts average quantities, but not the higher moments (on the contrary, its goal is to minimize the second moment!) same mean with smaller variance faster convergence allows sometimes to obtain acceptable statistics where an analog Monte Carlo would take years of CPU time to converge But: cannot reproduce correlations and fluctuations only privileged observables converge faster (some regions of phase space are sampled more than others)
36 The Geometry The algorithms to build a geometry and to track particles inside it differ from code to code. In general: The geometry is built from basic solids and/or surfaces: EGS: original CG, or user-defined FLUKA (extended CG): solid bodies and surfaces GEANT: solid volumes MARS: 5 different geometries, including FLUKA and MCNP geometry MCNP(X): originally surfaces, now also macrobodies MORSE, SAM-CE (original Combinatorial Geometry, CG): solid bodies PENELOPE: surfaces PHITS: MCNP geometry TRIPOLI: solid bodies and surfaces The geometry must have a way to limit the tracking, in some codes by defining an external boundary Defined by input cards (MCNP, CG) or by user-written routines (GEANT) In some codes it can allow for repetition of structures In some codes it can allow for voxel representation (CT import)
37 Discrete events Atomic interactions: Physics Photons: Compton, photoelectric, pair, coherent, Charged particles: bremsstrahlung, ray emission, large angle Coulomb scattering (all above given thresholds) Nuclear interactions: Nuclear elastic scattering Nuclear non-elastic interactions Absorption Decays Boundary crossing Escape Geometry
38 Continuous processes Physics: Charged particles lose energy and change direction as a result of thousands of discrete collisions with atomic electrons. To simulate in detail each collision would require prohibitive CPU times, except at very low particle energies. Many discrete scatterings are then replaced by a straight continuous step, and the corresponding energy losses and changes of direction are condensed into a sum of losses (de/dx) and an overall scattering angle (condensed-history technique) Ionization energy losses: all losses lower than a preset threshold are continuously distributed along the step. Any loss larger than threshold is simulated as a discrete energy imparted to an electron ( ray, transported) energy loss fluctuations can be simulated for the losses below threshold Multiple Coulomb Scattering: a deflection angle, sampled from a theoretical distribution, is applied to each particle step. Some corrections are needed: Path Length Correction (PLC) lateral displacement Transport: neutral particle displacements, particle displacements in vacuum ( steps ) When magnetic/electric field are present, substeps are necessary to follow the curvature of the trajectory, de/dx and MCS
39 Thresholds and cut-off s Transport and production thresholds are needed because of: Limits of validity of the physics models CPU time Transport thresholds: they depend on the granularity of the geometry and/or of the scoring mesh and on the interest in a given region. To reproduce correctly electronic equilibrium, neighboring regions should have the same electron energy threshold (not the same range threshold!) Photon thresholds should be lower than electron thresholds (they travel more) Production thresholds: δ-ray threshold: sets the limit between discrete and continuous ionization energy losses e/ production threshold : similar to δ-ray for other electro-magnetic processes, including bremsstrahlung
40 Results from a MC calculation Estimators It is often said that Monte Carlo is a mathematical experiment. Each aspect of a real experiment has its Monte Carlo equivalent: Experimental technique Estimator Instrument Detector Measurement Score, or Tally Result of an experiment Monte Carlo result Just as a real measurement, a score is obtained by sampling from a statistical distribution As an experimental result consists in an average of the measurements, a statistical error and a systematic error, a MC result is an average of scores, a statistical error (and a systematic error, generally unknown) There are often several different techniques to measure the same physical quantity: in the same way the same quantity can be calculated with different kinds of estimators
41 Estimator types Various types of estimators, depending on the quantity to be estimated and on the topology (phase space region over which the quantity is integrated) Boundary Crossing. Quantity: the fluence or the current of particles Phase space: a physical boundary between two space regions. Result: mono or multi-differential fluence spectra, function of energy, angle, particle type, Track length. Quantity: the fluence of particles Phase space: a region of real space Result: fluence spectra as a function of particle energy, based on the path length of the particles within the region volume
42 Estimator types Pulse-height detector (e.g. simulation of a Ge spectrometer): Quantity: energy deposited Phase space: a region of real space Result: spectrum of deposited energy within the region volume Scalar integral estimator: Quantity: deposited energy, inelastic interactions (stars), activity... Phase space: a region of real space Result: amount of given quantity (or its density) within the volume Mesh: Quantity: fluence, energy deposition, stars Phase space: regular subdivision of a portion of real space in subvolumes, generally independent from the tracking geometry Result: a 2D or 3D spatial distribution of the estimated quantity
43 Statistical Errors The variance of the mean of an estimated quantity x (e.g., fluence), calculated in N batches, is: mean of squares square of means N 1 where: n i = number of histories in the i th batch n = Σn i = total number of histories in the N batches x i = average of x in the i th batch: In the limit N = n, n i =1, the formula applies to single history statistics n x n n x n N N i i i N i x n i j i ij i n x x 1
44 Statistical Errors Can be calculated for single histories, or for batches of several histories each (not necessarily the same identical number) Distribution of scoring contributions by single histories can be very asymmetric (many histories contribute little or zero) Scoring distribution from batches tends to Gaussian for N, provided 2. (Central Limit Theorem) The standard deviation of an estimator calculated from batches or from single histories is an estimate of the standard deviation of the actual distribution ( error of the mean ) How good is such an estimate depends on the type of estimator and on the particular problem (but it converges to the true value for N )
45 Statistical Errors Practical tips: Use always at least 5-10 batches of comparable size (it is not at all mandatory that they be of equal size) Never forget that the variance itself is a stochastic variable subject to fluctuations Be careful about the way convergence is achieved: often (particularly with biasing) apparently good statistics with few isolated spikes could point to a lack of sampling of the most relevant phase-space part Plot 2D and 3D distributions! In those cases the eye is the best tool in judging the quality of the result
46 Statistical errors, systematic errors, and... mistakes Statistical errors, due to sampling (in)efficiency Relative error Quality of Tally (from the MCNP Manual) 50 to 100% Garbage 20 to 50% Factor of a few 10 to 20% Questionable < 10% Generally reliable except for point detectors < 5% Generally reliable for point detectors Why does a 30% σ mean an uncertainty of a factor of a few? Because σ in fact corresponds to the sum (in quadrature) of two uncertainties: one due to the fraction of histories which give a zero contribution, and one which reflects the spread of the non-zero contributions The MCNP guideline is empirically based on experience, not on a mathematical proof. But it has been generally confirmed also working with other codes Small penetrations and cracks are very difficult to handle by MC, because the detector is too small and too few non-zero contributions can be sampled, even by biasing
47 Statistical errors, systematic errors, and... mistakes Systematic errors, due to code weaknesses Apart from the statistical error, which other factors affect the accuracy of MC results? physics: different codes are based on different physics models. Some models are better than others. Some models are better in a certain energy range. Model quality is best shown by benchmarks at the microscopic level (e.g. thin targets) artifacts: due to imperfect algorithms, e.g., energy deposited in the middle of a step, inaccurate path length correction for multiple scattering, missing correction for cross section and de/dx change over a step, etc. Algorithm quality is best shown by benchmarks at the macroscopic level (thick targets, complex geometries) data uncertainty: an error of 1% in the absorption cross section can lead to an error of a factor 2.8 in the effectiveness of a thick shielding wall (10 attenuation lengths). Results can never be better than allowed by available experimental data!
48 Statistical errors, systematic errors, and... mistakes Systematic errors, due to user ignorance Missing information: material composition not always well known. In particular concrete/soil composition (how much water content? Can be critical) beam losses: most of the time these can only be guessed. Close interaction with engineers and designers is needed presence of additional material, not well defined (cables, supports...) Is it worth to do a very detailed simulation when some parameters are unknown or badly known? Systematic errors, due to simplification Geometries that cannot be reproduced exactly (or would require too much effort) Air contains humidity and pollutants, has a density variable with pressure
49 Statistical errors, systematic errors, and... mistakes MC codes can contain bugs: Code mistakes ( bugs ) Physics bugs: I have seen pair production cross sections fitted by a polynomial... and oscillating instead of saturating at high energies, nonuniform azimuthal scattering distributions, energy non-conservation... Programming bugs (as in any other software, of course) User mistakes mis-typing the input: some codes are more or less good at checking, but the final responsibility is the user s error in user code: use the built-in features as much as possible wrong units wrong normalization: quite common unfair biasing: energy/space cuts cannot be avoided, but must be done with much care forgetting to check that gamma production is available in the neutron cross sections
50 Some of the most popular codes EGS5: EGSnrc: (e- e + γ) FLUKA : (multiparticle) Geant4: (multiparticle) MARS: (multiparticle) MCNP: (neutrons, e- e+ γ) MCNPX: (multiparticle) PENELOPE: (e- e+ γ) PHITS: (multiparticle) TRIPOLI: (neutrons, e- e+ γ) Specialized: CORSIKA: (cosmic rays) PEREGRINE: (radiotherapy)
51 General information for various all-particle transport codes General MCNPX GEANT4 FLUKA MARS PHITS Version p Lab. Affiliation LANL CERN ESA IN2P3 PPARC INFN LIP KEK SLAC TRIUMF CERN INFN FNAL JAEA RIST GSI Chalmers Univ. Language Fortran 90/C C++ Fortran 77 Fortran 95/C Fortran 77 Cost Free Free Free Free Free Release Format Source & binary Source & binary Source & binary Binary Source & binary Availability Conditions RSICC Beta test team Open web ne User s Agreement User Manual 470 pages 280 pages 387 pages 150 pages 176 pages Users ~2000 ~2000 ~1000 ~ Web Site mcnpx.lanl.gov cern.ch/geant4 www-ap.fnal. gov/mars Under const. Workshops ~7/year ~4/year ~1/year ~2/year ~1/year Input Format Free C++ main Fixed geometry Fixed or free Free Free Input Cards ~120 N/A ~85 0 to 100 ~100 Parallel Execution
52 Geometry Capabilities Geometry MCNPX GEANT4 FLUKA MARS PHITS Description MCNP-based Solids (CSG, Boolean, some BREP/STEP) Combinatorial Fixed shapes or MCNP-based MCNP-based MORSE-based Extensions Twisted Nested Repeated Voxel (universes) Lattice (rec, hex) (logical vol.) (rec, cyl) (universes) Lattice (rec, hex) Reflections 3 types Neutron albedo Viewer Debugger Built-in: 2-D Interactive X-Windows External: Vised Moritz Built-in: 3-D Interactive OpenGL OpenInventor RayTracer External: WIRED VRML DAWN Overlap tools Built-in: ne External: Custom (X11) Debugger built in Built-in: 2-D Interactive Tcl/Tl 3-D Interactive OpenGL External: Built-in Built-in: 2,3-D Command PS via Angel External: Angel PS Setup GUI Vised Moritz GGE Tcl/Tl CAD STEP via GUI STEP via Tool Fields (E/B) 2.6.D Moving Solid Objects 2.6.D
53 Sources Source MCNPX GEANT4 FLUKA MARS PHITS Fixed General Explicit Distribution Dep. Dist. External User Sub. SSW/SSR GPS Eigenvalue Burnup (2.6.A)
54 Physics Capabilities Physics MCNPX GEANT4 FLUKA MARS PHITS Particles Charged particles Energy loss Scatter Straggling XTR/Cheren. CSDA Bethe-Bloch Rossi Vavilov CSDA Bethe-Bloch Lewis Urban CSDA Bethe-Bloch Moliere improved Custom /yes CSDA Bethe-Bloch Moliere improved Custom CSDA Bethe-Bloch Moliere Vavilov Baryons Neutron Low High Proton Low High Other Cont. (ENDF) Models Cont. (ENDF) Models Model List: Bertini ISABEL CEM INCL FLUKA89>3 GeV LAQGSM (2.6.C) Cont. (ENDF) Models Models Models Model list: Hadron-nucleous GHEISHA* INUCL(Bertini) BIC CHIPS QGS/FTF>8 GeV Multigroup(72) Models Models Models Model list: PEANUT(GINC) +DPM+Glauber Cont. (ENDF) Models Models Models Model list: Custom CEM LAQGSM DPMJET Cont. (ENDF) Models Models Models Model list: Bertini JAM>3 GeV Leptons Electrons Muon Neutrino Other ITS 3.0 CSDA/decay Production Decay Models/EEDL, EADL Models Production Decay Custom Models Models Decay Custom Models Models Models ITS 3.0 CSDA/decay Models Models
55 Physics Capabilities, cont. Mesons Models Models Models Models Models Photons Optical x-ray/g Photonuclear ITS 3.0 Libraries (IAEA) CEM Models or EPDL97, EADL CHIPS Custom+EPDL97 PEANUT VMDM Custom Custom CEM ITS 3.0 Ions ISABEL LAQGSM (2.6.C) AAM EDM BLIC RQMD-2.4 DPMJET-3 LAQGSM JQMD JAMQMD > 3 GeV/u Delayed n, (2.6.C) α,,?, n
56 Tallies/Scores/Edits (Results) Tallies MCNPX GEANT4 FLUKA MARS PHITS Standard Flux Volume Surface Point/ring Current Charge Kinetic energy Particle density Reaction rates Energy deposition Rapidity DPA Momentum Pulse-height Termination Modifiers HTAPE3X Partial 9 Limited Limited (user) User input?? 2 Star (inelastic) Some 2 (neutrons). Partial 2 2 Special Mesh Coincidence Residuals Activation Event logs rec, cyl, sph 2.5.D rec, cyl?? rec, cyl rec, cyl, sph rec,cyl
57 Tallies/Scores/Edits (Results) cont. Convergence Tests 10 Error Error Error Error Results Viewer Built-in: 1-D, 2-D Custom X-Windows External: IDL Tecplot GNUplot PAW Built-in: External: JAS PI Open Scientist PAW Built-in: ne External: Custom (X11) GNUplot PAW ROOT Built-in: Custom External: PAW Built-in: Angel External: Angel
58 Variance Reduction Variance Reduction MCNPX GEANT4 FLUKA MARS PHITS Population control Region biasing Weight cutoff Weight window mesh Energy biasing Modified sampling Source biasing Implicit capture Exp. transform Production biasing Angular bias DXTRAN RDM?? DXTRAN Viewer 2-D contour
59 Radioactive decay in Geant4 and Fluka Enjoy the code and the practical exercise
60 Conclusions? Use data to validate MC MC comparison is important and useful As always trust, but do not forget a critical view data!= information!= knowledge!= wisdom!= truth MC
Beta decay with different MC codes
Beta decay with different MC codes presented by Andrea Fontana INFN Pavia 23.03.2017 http://www.pv.infn.it/~fontana Main topics: - How to run the same setup in different MC: a lot of work! - How to evaluate
More informationExternal MC code : PHITS
External MC code : PHITS Particle and Heavy Ion Transport code System Koji. Niita 1, Tatsuhiko Sato 2, Hiroshi Iwase 3, Yosuke Iwamoto 2, Norihiro Matsuda 2, Yukio Sakamoto 2, Hiroshi Nakashima 2, Davide
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 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 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 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 informationBiasing. Beginners FLUKA Course
Biasing Beginners FLUKA Course Overview General concepts: Analog vs. biased Monte Carlo calculation Biasing options (only the most important / common options available in FLUKA) Importance biasing Weight
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 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 informationImprovements and developments of physics models in PHITS for radiotherapy and space applications
Improvements and developments of physics models in PHITS for radiotherapy and space applications L. Sihver 1-9, T. Sato 10, S. Hashimoto 10, T. Ogawa 10, K. Niita 11 1 Atominstitut, TU Wien, Austria, 2
More informationRadiation damage calculation in PHITS
Radiation Effects in Superconducting Magnet Materials (RESMM'12), 13 Feb. 15 Feb. 2012 Radiation damage calculation in PHITS Y. Iwamoto 1, K. Niita 2, T. Sawai 1, R.M. Ronningen 3, T. Baumann 3 1 JAEA,
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 informationUpcoming features in Serpent photon transport mode
Upcoming features in Serpent photon transport mode Toni Kaltiaisenaho VTT Technical Research Centre of Finland Serpent User Group Meeting 2018 1/20 Outline Current photoatomic physics in Serpent Photonuclear
More informationGeant4 Physics Lists: Status and Proposed Upgrades. Dennis Wright (SLAC) 25 February 2011
Geant4 Physics Lists: Status and Proposed Upgrades Dennis Wright (SLAC) 25 February 2011 Outline Contents of a few preferred Geant4 physics lists Updating/augmenting the physics lists Comparing Fluka and
More informationGeant4 Hadronic Physics Developments
Geant4 Hadronic Physics Developments José Manuel Quesada University of Sevilla on behalf of Geant4 Hadronic Working Group 9th Geant4 Space Users Workshop Barcelona, March 2013 Outline General matters Photo-nuclear
More informationICRP Symposium on the International System of Radiological Protection
ICRP Symposium on the International System of Radiological Protection October 24-26, 2011 Bethesda, MD, USA Akira Endo and Tatsuhiko Sato* ICRP Committee 2 & Task Group 4 (DOCAL) * Task Group 67 (Radiation
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 informationSimulating Gamma-Ray Telescopes in Space Radiation Environments with Geant4: Detector Activation
Simulating Gamma-Ray Telescopes in Space Radiation Environments with Geant4: Detector Activation Andreas Zoglauer University of California at Berkeley, Space Sciences Laboratory, Berkeley, USA Georg Weidenspointner
More informationMonte Carlo Tutorial. Part 2 Practical MC
Monte Carlo Tutorial Part 2 Practical MC IFJ-PAN, Cracow, May 2014 Niels Bassler Dept. Physics and Astronomy Aarhus University, Denmark General recipe to setup a MC run How to setup a MC run First: think
More informationFrom Russia with Love Monte Carlo Particle Transport Code SHIELD-HIT(10A)
From Russia with Love Monte Carlo Particle Transport Code (10A), PhD, Department of Physics and Astronomy Aarhus University, Denmark Department of Experimental Clinical Oncology Aarhus
More informationRadiation Transport Tools for Space Applications: A Review
Radiation Transport Tools for Space Applications: A Review Insoo Jun, Shawn Kang, Robin Evans, Michael Cherng, and Randall Swimm Mission Environments Group, February 16, 2008 5 th Geant4 Space Users Workshop
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 informationElectronuclear Interactions in FLUKA
Electronuclear Interactions in FLUKA Pavel Degtiarenko Jefferson Lab Contents What are the inelastic direct electronuclear (ea) reactions Why ea reactions important/critical for JLab s needs What is our
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 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 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 informationFeatures of PHITS2.82. PHITS development team, Dec. 25, 2015
1 Features of PHITS2.82 PHITS development team, Dec. 25, 2015 Map of Models used in PHITS2.82 Low Energy High Neutron Proton, Pion (other hadrons) 1 TeV 1 TeV/n Intra-nuclear cascade (JAM) Evaporation
More informationInteraction of Particles with Matter
Chapter 10 Interaction of Particles with Matter A scattering process at an experimental particle physics facility is called an event. Stable particles emerging from an event are identified and their momenta
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 informationPHYS 5020 Computation and Image Processing
PHYS 5020 and Image Processing : Monte Carlo Thursday 2 August 2012 Monte Carlo (MC) is a numerical method that uses random sampling of probability distributions to simulate stochastic processes in nature,
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 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 informationA MONTE CARLO SIMULATION OF COMPTON SUPPRESSION FOR NEUTRON ACTIVATION ANALYSIS. Joshua Frye Adviser Chris Grant 8/24/2012 ABSTRACT
A MONTE CARLO SIMULATION OF COMPTON SUPPRESSION FOR NEUTRON ACTIVATION ANALYSIS Joshua Frye Adviser Chris Grant 8/24/2012 ABSTRACT A Monte Carlo simulation has been developed using the Geant4 software
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 informationIAC-08-A MONTE CARLO SIMULATIONS OF ENERGY LOSSES BY SPACE PROTONS IN THE CRATER DETECTOR
IAC-08-A1.4.06 MONTE CARLO SIMULATIONS OF ENERGY LOSSES BY SPACE PROTONS IN THE CRATER DETECTOR Lawrence W. Townsend The University of Tennessee, Knoxville, Tennessee, United States of America ltownsen@tennessee.edu
More informationGeant4 version 10.0.p01. Hadronic Physics I. Geant4 Tutorial: version 10.0.p01. Michael Kelsey, Wed 5 Mar 2014
Michael Kelsey, Wed 5 Mar 2014 Hadronic Physics I Geant4 Tutorial: version 10.0.p01 Hadronic Physics I What is Hadronic Physics? The Hadronic Framework - Processes vs. Models - Cross sections and process
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 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 informationOn the energy deposition by electrons in air and the accurate determination of the air-fluorescence yield
arxiv:1207.2913v1 [astro-ph.im] 12 Jul 2012 On the energy deposition by electrons in air and the accurate determination of the air-fluorescence yield J. Rosado, P. Gallego, D. García-Pinto, F. Blanco and
More informationThe Development of Particle Physics. Dr. Vitaly Kudryavtsev E45, Tel.:
The Development of Particle Physics Dr. Vitaly Kudryavtsev E45, Tel.: 0114 2224531 v.kudryavtsev@sheffield.ac.uk Discovery of the muon and the pion Energy losses of charged particles. This is an important
More informationPhysics 663. Particle Physics Phenomenology. April 23, Physics 663, lecture 4 1
Physics 663 Particle Physics Phenomenology April 23, 2002 Physics 663, lecture 4 1 Detectors Interaction of Charged Particles and Radiation with Matter Ionization loss of charged particles Coulomb scattering
More informationIntroduction to Geant4 Physics Overview
Introduction to Geant4 Physics Overview Koi, Tatsumi SLAC SCCS Based on Presentations at SLAC Geant4 Tutorial 2007 Outline Geant4 Physics Overview Process Physics List Standard EM Low Energy EM Hadron
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 informationProblem P7. Stéphanie Ménard. Dosimetry Department Fontenay-aux FRANCE IRSN QUADOS IRSN
Problem P7 Stéphanie Ménard Dosimetry Department 92262 Fontenay-aux aux-roses FRANCE What are the applications of Gamma-Ray Spectrometry in Radiological Protection and in Safety? In the environment: after
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 informationRecent Developments in Geant4 Hadronic Physics
Recent Developments in Geant4 Hadronic Physics Dennis Wright (SLAC) 1 st Geant4 Australian School Recently Released Many new features added in Geant4 9.4 (December 2010) and patch 9.4 p01 2 Energy/Momentum
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 informationBETA-RAY SPECTROMETER
14 Sep 07 β-ray.1 BETA-RAY SPECTROMETER In this experiment, a 180, constant-radius magnetic spectrometer consisting of an electromagnet with a Geiger-Muller detector, will be used to detect and analyze
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 informationNeutron and/or photon response of a TLD-albedo personal dosemeter on an ISO slab phantom
Neutron and/or photon response of a TLD-albedo personal dosemeter on an ISO slab phantom Problem P4 Rick J Tanner National Radiological Protection Board Chilton, Didcot, Oxon OX11 0RQ, United Kingdom Intercomparison
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 informationLecture notes on Chapter 13: Electron Monte Carlo Simulation. Nuclear Engineering and Radiological Sciences NERS 544: Lecture 13, Slide # 1:13.
Lecture notes on Chapter 13: Electron Monte Carlo Simulation Nuclear Engineering and Radiological Sciences NERS 544: Lecture 13, Slide # 1:13.0 Topics covered Basic, relevant interaction processes Hard
More informationOn the energy deposition by electrons in air and the accurate determination of the air-fluorescence yield
EPJ Web of Conferences 53, 10001 (2013) DOI: 10.1051/epjconf/20135310001 C Owned by the authors, published by EDP Sciences, 2013 On the energy deposition by electrons in air and the accurate determination
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 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 informationRecent Developments of the PHITS code
Progress in NUCLEAR SCIENCE and TECHNOLOGY, Vol. 1, p.1-6 (2011) REVIEW Recent Developments of the PHITS code Koji NIITA 1*, Hiroshi IWASE 2, Tatsuhiko SATO 3, Yosuke IWAMOTO 3, Norihiro MATSUDA 3, Yukio
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 informationNew photon transport model in Serpent 2
New photon transport model in Serpent 2 Toni Kaltiaisenaho VTT Technical Research Centre of Finland Serpent User Group Meeting 1/20 Motivation On average, 8 prompt fission photons over an energy range
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 information^ He, Li. INITIAL STUDIES OF GEANT and COMPARISONS with EGS. ^/vo-ffr<-y^-/^ 7. 2 Material Generation. W. Frati. September 8, 1990.
/vo-ffr
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 informationBrief Introduction to: Transport Theory/Monte Carlo Techniques/MCNP
22:54 Neutron Interactions and Applications Brief Introduction to: Transport Theory/Monte Carlo Techniques/MCNP The behavior of individual neutrons and nuclei cannot be predicted. However, the average
More informationFLUKA Manual and Basic input. 7th FLUKA Course NEA Paris, Sept.29-Oct.3, 2008
FLUKA Manual and Basic input 7th FLUKA Course NEA Paris, Sept.29-Oct.3, 2008 The FLUKA Manual in continuous development, just as the program more a User Guide than a Reference Manual (only a short summary
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 informationParticle Energy Loss in Matter
Particle Energy Loss in Matter Charged particles loose energy when passing through material via atomic excitation and ionization These are protons, pions, muons, The energy loss can be described for moderately
More informationLast Lecture 1) Silicon tracking detectors 2) Reconstructing track momenta
Last Lecture 1) Silicon tracking detectors 2) Reconstructing track momenta Today s Lecture: 1) Electromagnetic and hadronic showers 2) Calorimeter design Absorber Incident particle Detector Reconstructing
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 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 informationICTP-IAEA Joint Workshop on Nuclear Data for Science and Technology: Medical Applications. 30 September - 4 October, 2013
2484-11 ICTP-IAEA Joint Workshop on Nuclear Data for Science and Technology: Medical Applications 30 September - 4 October, 2013 Experimental techniques (Nuclear reaction data, estimation of uncertainties)
More informationPhotonuclear Reaction Cross Sections for Gallium Isotopes. Serkan Akkoyun 1, Tuncay Bayram 2
Photonuclear Reaction Cross Sections for Gallium Isotopes Serkan Akkoyun 1, Tuncay Bayram 2 1 Cumhuriyet University, Vocational School of Healt, Sivas, Turkey 2 Sinop University, Department of Physics,
More informationCurrent status of EGS5
Workshop on Computational Medical Physics (Nara, September. 2, 2012) Current status of EGS5 Comparison with experiment on e backscattering Y. Namito and H. Hirayama (KEK) History of EGS System Period Program
More informationBasic physics Questions
Chapter1 Basic physics Questions S. Ilyas 1. Which of the following statements regarding protons are correct? a. They have a negative charge b. They are equal to the number of electrons in a non-ionized
More informationRadiation Detection and Measurement
Radiation Detection and Measurement June 2008 Tom Lewellen Tkldog@u.washington.edu Types of radiation relevant to Nuclear Medicine Particle Symbol Mass (MeV/c 2 ) Charge Electron e-,! - 0.511-1 Positron
More informationDIFFERENT H.E. INTERACTION MODELS
COMPARISON OF EAS SIMULATION RESULTS USING CORSIKA CODE FOR DIFFERENT H.E. INTERACTION MODELS Dept. of Physics, Gauhati University, Assam, INDIA 15th Sep, 2010 Out line of the talk Introduction Primary
More informationGeant4 and its validation
http://geant4.web.cern.ch/geant4/ Geant4 and its validation Luciano Pandola INFN Gran Sasso and University of L Aquila for the Geant4 Collaboration Siena, May 24 th, 2004 What is? OO Toolkit for the simulation
More informationLET! (de / dx) 1 Gy= 1 J/kG 1Gy=100 rad. m(kg) dose rate
Basics of Radiation Dosimetry for the Physicist http://en.wikipedia.org/wiki/ionizing_radiation I. Ionizing radiation consists of subatomic particles or electromagnetic waves that ionize electrons along
More informationPhysics sources of noise in ring imaging Cherenkov detectors
Nuclear Instruments and Methods in Physics Research A 433 (1999) 235}239 Physics sources of noise in ring imaging Cherenkov detectors For the ALICE HMPID Group Andreas Morsch EP Division, CERN, CH-1211
More informationCross-section Measurements of Relativistic Deuteron Reactions on Copper by Activation Method
Nuclear Physics Institute, Academy of Sciences of the Czech Republic Department of Nuclear Reactors, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University in Prague Cross-section
More informationDesign, Construction, Operation, and Simulation of a Radioactivity Assay Chamber
Design, Construction, Operation, and Simulation of a Radioactivity Assay Chamber Wesley Ketchum and Abe Reddy EWI Group, UW REU 2006 Outline Neutrino Physics Background Double Beta Decay and the Majorana
More informationIonization and Transport. FLUKA Beginners Course
Ionization and Transport FLUKA Beginners Course General settings Interactions of leptons/photons Photon interactions Photoelectric Compton Rayleigh Pair production Photonuclear Photomuon production Electron/positron
More informationInteraction of particles in matter
Interaction of particles in matter Particle lifetime : N(t) = e -t/ Particles we detect ( > 10-10 s, c > 0.03m) Charged particles e ± (stable m=0.511 MeV) μ ± (c = 659m m=0.102 GeV) ± (c = 7.8m m=0.139
More informationOutline. Absorbed Dose in Radioactive Media. Introduction. Radiation equilibrium. Charged-particle equilibrium
Absorbed Dose in Radioactive Media Chapter F.A. Attix, Introduction to Radiological Physics and Radiation Dosimetry Outline General dose calculation considerations, absorbed fraction Radioactive disintegration
More informationGeant Hadronic Physics I. Geant4 Tutorial at Lund University. 6 September 2018 Dennis Wright
Geant4 10.4 Hadronic Physics I Geant4 Tutorial at Lund University 6 September 2018 Dennis Wright Outline Overview of hadronic physics Precompoundand de-excitation models Cascade models 2 Hadronic Processes,
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 informationY2 Neutrino Physics (spring term 2017)
Y2 Neutrino Physics (spring term 2017) Lecture 5 Discoveries of the leptons Dr E Goudzovski eg@hep.ph.bham.ac.uk http://epweb2.ph.bham.ac.uk/user/goudzovski/y2neutrino Previous lecture In 1940s, nuclear
More informationThe Compact Muon Solenoid Experiment. Conference Report. Mailing address: CMS CERN, CH-1211 GENEVA 23, Switzerland
Available on CMS information server CMS CR -2018/225 The Compact Muon Solenoid Experiment Conference Report Mailing address: CMS CERN, CH-1211 GENEVA 23, Switzerland 27 September 2018 (v2, 19 November
More informationAtoms, Radiation, and Radiation Protection
James E. Turner Atoms, Radiation, and Radiation Protection Third, Completely Revised and Enlarged Edition BICENTENNIAL J 0 1 8 0 Q 71 z m z CAVILEY 2007 1 ;Z z ü ; m r B10ENTENNIAL WILEY-VCH Verlag GmbH
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 informationResearchers at the University of Missouri-Columbia have designed a triple crystal
Childress, N. L. and W. H. Miller, MCNP Analysis and Optimization of a Triple Crystal Phoswich Detector, Nuclear Instruments and Methods, Section A, 490(1-2), 263-270 (Sept 1, 2002). Abstract Researchers
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 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 informationThe Physics of Cosmic Rays
The Physics of Cosmic Rays QuarkNet summer workshop July 23-27, 2012 1 Recent History Most natural phenomena can be explained by a small number of simple rules. You can determine what these rules are by
More informationPAMELA satellite: fragmentation in the instrument
PAMELA satellite: fragmentation in the instrument Alessandro Bruno INFN, Bari (Italy) for the PAMELA collaboration Nuclear Physics for Galactic CRs in the AMS-02 era 3-4 Dec 2012 LPSC, Grenoble The PAMELA
More informationGeant4 Simulation of Very Low Energy Electromagnetic Interactions
Geant4 Simulation of Very Low Energy Electromagnetic Interactions R. Capra 1, Z. Francis 2, S. Incerti 3, G. Montarou 2, Ph. Moretto 3, P. Nieminen 4, M. G. Pia 1 1 INFN Sezione di Genova; I-16146 Genova,
More information3 Monte Carlo tools for proton dosimetry
3 Monte Carlo tools for proton dosimetry 3.1 Introduction As is clear from the previous chapter the problems in proton radiation transport are difficult to solve analytically, because of the statistical
More informationPhoton transport mode in Serpent 2
Photon transport mode in Serpent 2 Toni Kaltiaisenaho VTT Technical Research Centre of Finland, LTD Serpent User Group Meeting, Knoxville, TN October 13 16, 215 October 14, 215 1/21 Outline Photon physics
More informationIon-ion Physics in Geant4. Dennis Wright (SLAC) 5 th Geant4 Space Users' Workshop 14 February 2008
Ion-ion Physics in Geant4 Dennis Wright (SLAC) 5 th Geant4 Space Users' Workshop 14 February 2008 Outline Introduction and Motivation Cross sections Existing models New Models Interfaces to external models
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 informationModelling and simulation methods of particle transport.
Modelling and simulation methods of particle transport Henrik Sjöstrand henrik.sjostrand@phycics.uu.se 1 Why this course? Skills important in many research fields Detector development, space science, accelerator
More information