Modern Methods of Data Analysis - WS 07/08
|
|
- Blake McDowell
- 6 years ago
- Views:
Transcription
1 Modern Methods of Data Analysis Lecture III ( ) Contents: Overview & Test of random number generators Random number distributions
2 Monte Carlo The terminology Monte Carlo-methods originated around 1940, related to the construction of the atomic bomb. Random processes were for the first time extensively used to theoretical predict the interaction of neutrons with matter.
3 History: Buffon's Needles (1777) Determine the value π by throwing needles (length l) on a plate width (l) 2 crossing needles k p = = --- = --all needles n π Numerical Problems can be solved approximately by statistics!
4 Monte Carlo: A method of statistical experiments Often analy. treatment of physical problem difficult/impossible Must do approximation <- or -> statistical description (MC) Application: High energy physics Numerical calculations (Integration,Differentiation coding/encoding e.g. ssh-connection Reliability tests Investment banking Earth sciences (earthquakes) Method: Find statistic model Find sequence of random numbers Calculate estimators from random quantities
5 Random Numbers First look at random numbers distributed uniform in [0,1] : U(0,1) U(0,1) random numbers shall have the correct size of fluctuations for a uniform distribution real random numbers can be extracted from physics processes radioactive decays thermal noise over a resistor alternatively on can use pseudo random numbers generation via numerical algorithms not truly random, but hopefully indistinguishable from true random numbers available any time and with a reproducible sequence important for debugging purposes
6 Historical Example Von Neumann's mid-squares generator Start with m-digit number seed. Calculate seed*seed and take the central m digits as new seed. Normalization of seed to the range [0,1] yields a pseudo-random sequence
7 Linear Congruential Generators (LCG) x(n+1) = (a * x(n) + b) mod m in [0, m-1] The properties/quality of the generator is determined by the parameters a,b and m Maximal period m is reached only for good choice of a,b and m -> distribution of single values: Example: a = 1601; b = 3456; m = 10000;
8 Properties of LCG: LCG random numbers are located on hyperplanes # of hyperplanes if function of dimension d of representation minimal requirements for good generator number of hyperplanes as large as possible period as long as possible both characteristics grow with the number of bits per integer & choice of parameters a,b,m important Mathematical max. nb. of hyperplanes:
9 RANDU Generator 1960's: RANDU (IBM, 32bit): m =, a=65539, c=0; seed: n = 1 RANDU (10000 numbers) high periodicity, however very low rate of hyperplanes many scientific studies suffered from this!
10 RANLUX Generators based on Marsaglia-Zaman algorithm mathematically equivalent to a LCG completely different implementation [hep-lat/930902] RANLUX implementation: with m = z[n] = (a * z[n-1]) mod m (prime) and a = Effectively the generator is based on 576-bit integer variables. Period: ; # of hyperplanes for d=100: However: a has only very few bits set, subsequent 576-bit states are still correlated. RANLUX-luxury-level 4, therefore discards 8760 bits of the bit stream produced by the generator, before the next 576 bits are accepted and convert them to 24 single precision floating point numbers.
11 Combining Several Generators Two standard procedures to combine two generators: combination: two random numbers are generate with two separate generators. New random number obtained by combination, e.g. via +, - or bit wise exclusive OR mixing: an array is filled with random numbers from one generator. Random number from second generator determine order in which random numbers are used. Naively one could think more complex functions yield more randomness, this is not true! [D.E. Knuth: The Art of Computer Programming]
12 Generator by L'Ecuyer combination of two LCG: x[n+1] = (a1 * x[n]) mod m1 y[n+1] = (a2 * y[n]) mod m2 z[n+1] = (x[n+1]+y[n+1]) mod m1 with following choice of variables: a1 = 40014; a2 = 40692; m1 = = * a m2 = = * a with period m1, m2: m1 = 2 * 3 * 7 * 631 * m2 = 2 * 19 * 31 * 1019 * 1789 z has period of and significantly larger nb. of hyperplanes than individual generators x, y.
13 State of the Art 1998 by Makoto Matsumoto and Takuji Nishimura Marsenne Twister algorithm based on bit-shifting bit long seed value, stored in a 624 long array with 32 bit values (31 bits are free) is a Marsenne prime number if a shift register as a size which is a Marsenne prime the testing procedure for periodicity is significantly simplified periode: in decimals: up to 623 dimension maximal # of hyperplanes it is fast! Recommended generator!
14 Marsenne Twister (MT)
15 Random Generators in ROOT TRandom (LCG,, known to have correlations!) 34 ns/call TRandom1 (RANLUX, math prooven random properties, 242 ns/call TRandom2 (based on l'ecuyer, ) 37 ns/call TRandom3 ) recommended one! (Mersenne Twister, good random properties 45 ns/call )
16 Test of Random Generators (I) Test for distance of hyperplanes of random numbers Test for distance of hyperplanes of last digits of random nb. Sequence-(up-down-) test: compare and. Set bit 0/1 according to or, count how often 0 and 1 are in a row of k numbers. For N truly uniformly distributed numbers we expect:
17 Test of Random Generators (II) Test for uniform distribution: [0,1] is split in k subintervals of same size. N random numbers are generated. Compare expected number of events (N/k) to number of events in interval i. Sum the deviation divided by the expected fluctuation (σ² =N/k) in quadrature, compare to expected distribution. Random Walk-test: choose 0 < α «1, generate large set of r.nb., count numbers below α, compare to expected number many more tests: Most try to compare generated distribution to true random distribution.
18 Warning Every Random generator failing those tests is a bad one Truly comprehensive tests are impossible. Use professional random generators! Use nontrivial seeds! However, keep a little suspect in mind.
19 Random numbers for any distribution given distribution f(x), F(x) = F(x) = u is analytically invertible: x = if u[i] uniform in [0,1], then follows pdf f(x) Method is applicable, if F(x) and analytical solvable
20 Example: Exponential Distribution Simulate random numbers according to: x>0 f(x,λ) = 0 x 0 Computation of inverse pdf: thus random numbers x[i] = - ln (u[i])/λ with u[i] uniform distribution in [0,1] follow f(x,λ)!
21 Hit - & Miss Method (von Neumann) If is not computable, then use Hit & Miss: f(x) pdf in [a,b] with maximum fmax generate uniform distributed random numbers x[i] in [a,b] generate second uniform distributed random numbers y[i] in [0, fmax] accept x[i] if y[i] < f(x[i]) else reject x[i] fmax f(x) a b
22 Modified Hit & Miss Especially with extreme functions, the method is inefficient. Improve by selecting invertible function f'(x) with according F(u) and f'(x) f(x). generate x[i] according f'[x] generate y[i] in [0,1] accept x[i] if y[j]*f'(x[j]) < f(x[j]) else reject For often used pdf (e.g. Gaussian distribution), there are optimized algorithms
23 Box-Muller-Transformation (I) not analytical solvable... x,y two independent Gaussian distributions
24 Box-Muller-Transformation (II) simulate uniform distribution u1 in [0,1] and u2 [0,2π]] x,y independent Gaussian distributed random variables
25 MC Methods: Integration (I) f(x) = 1/(a-b) a < x < b 0 elsewhere xi: uniformly distributed random numbers in [a,b]
26 Monte Carlo Methods: Integration (II)
27 Reduction of V[g(x)] (I) Stratification: c = 0.5 (a+b)
28 Reduction of V[g(x)] (II) Importance Sampling: f(x) normalized in [a,b], integratable & invertible for probability density f(x) x random numbers in [a,b] according to f(x) small if g(x)/f(x) almost constant alternative g(x) = g'(x) + f(x); F(x) analytical solvable;
29 Example in two dimensions Easy to get caught in a trap x uniform in [0,1] y uniform in [0,x] Accept only if y<x Alternative efficient algorithm: u uniform in [0,1], v uniform in [0,1] x = max(u,v); y = min(u,v)
30 Comparison to numerical Integration Trapezoidal rule: Error: 1/n² Simpson's rule: Error: 1/n4 In 1-D, conventional methods are always better than MC ( In d-d, trapezoidal rule: Simpson's rule: )
31 Quasi Random Numbers smaller variance, better equally distributed exploit stratification convergence [however not yet fully understood]
32 MC Methods: Example - Ks Mesons A source of Ks mesons, with mass M=0.498 GeV, produces a narrow monoenergetic beam of particles with energy E = GeV. Most Ks decay with an average lifetime s, or cm, into pairs. In the rest system of the Ks the decay is isotropic. Located at a distance of 14 cm behind the source is a silicon detector to register the decay pions. The detector is a circular disk with radius R=7cm, centered on the Ks beam. Determine the probability that both pions from the decay hit the detector.
33 Analytical Calculation parametrize the decay in the Ks rest frame perform Lorentz boost of decay pions into the lab-system boost parameter is γ = E/M describe the flight distance z by an exponential distribution uniform in φ the beam direction uniform cosine of polar angle of the beam uniform 0 < cosθ < 1 probability density ρ(z) = 1/(<z>) exp(z/<z>) in the lab system one finds <z> = <ct> 4.07 cm determine as a function of the decay configuration the allowed z range, from which both pions reach the detector
34 Direct Simulation (I) generate decay position according to ρ(z) generate isotropic decays in the Ks rest frame perform Lorentz transformation to the lab system analyze the events, check weather pions hit detector count as success cases where both pions hit the detector count as failure if at least one particle misses Repeat N times (p = n/n): e.g. given N= generated events & n= successes: A = ± 405 versus analytical result A =
35 Direct Simulation (II): easy to adapt to more complex cases... other Ks energies an energy spectrum rectangular detectors introducing Φ dependence finite efficiency of detector direct simulation allows detailed studies albeit with finite statistical precision
36 More complex in real life (I) modern high energy experiments consist of several detectors-system CMS (LHC): silicon strip detectors with about read out channels Simulation of physics processes: [Pythia, Sherpa, Herwig, EvtGen, Bgen] all produced particles, expected energy, momentum, angular distribution. Interaction with material: [Geant] energy deposit, multiple scattering, secondaries detection efficiency for each detector component [experiment specific] position and energy resolution of each component [experiment specific]
37 More complex in real life (II) End of simulation (1-20 minutes per event!), list of digitalized signals from the various detector components. No difference to real data! Simulation used to: - develop reconstruction algorithms. (particle trajectories in the tracking system, shower in the calorimeters) - optimize trigger selection - optimize reconstruction of signal signature (combination of all detector components) - determine detector acceptance Simulation crutial for RND (research and development planing phase) of experiments, for preparation of data taking, for optimizing analysis and for evaluating significance of result.
38 More complex in real life (III) CDF: proton-antiproton accelerator at the Tevatron (Fermilab/USA) In average 100 observed particles. 3-Jet event
Uniform Random Number Generators
JHU 553.633/433: Monte Carlo Methods J. C. Spall 25 September 2017 CHAPTER 2 RANDOM NUMBER GENERATION Motivation and criteria for generators Linear generators (e.g., linear congruential generators) Multiple
More information2008 Winton. Review of Statistical Terminology
1 Review of Statistical Terminology 2 Formal Terminology An experiment is a process whose outcome is not known with certainty The experiment s sample space S is the set of all possible outcomes. A random
More informationReview of Statistical Terminology
Review of Statistical Terminology An experiment is a process whose outcome is not known with certainty. The experiment s sample space S is the set of all possible outcomes. A random variable is a function
More informationStatistische Methoden der Datenanalyse. Kapitel 3: Die Monte-Carlo-Methode
1 Statistische Methoden der Datenanalyse Kapitel 3: Die Monte-Carlo-Methode Professor Markus Schumacher Freiburg / Sommersemester 2009 Basiert auf Vorlesungen und Folien von Glen Cowan und Abbildungen
More informationPhysics 403 Monte Carlo Techniques
Physics 403 Monte Carlo Techniques Segev BenZvi Department of Physics and Astronomy University of Rochester Table of Contents 1 Simulation and Random Number Generation Simulation of Physical Systems Creating
More informationRandom numbers and generators
Chapter 2 Random numbers and generators Random numbers can be generated experimentally, like throwing dice or from radioactive decay measurements. In numerical calculations one needs, however, huge set
More informationRandom Number Generation. Stephen Booth David Henty
Random Number Generation Stephen Booth David Henty Introduction Random numbers are frequently used in many types of computer simulation Frequently as part of a sampling process: Generate a representative
More informationMonte Carlo Techniques
Physics 75.502 Part III: Monte Carlo Methods 40 Monte Carlo Techniques Monte Carlo refers to any procedure that makes use of random numbers. Monte Carlo methods are used in: Simulation of natural phenomena
More informationPhysics 403. Segev BenZvi. Monte Carlo Techniques. Department of Physics and Astronomy University of Rochester
Physics 403 Monte Carlo Techniques Segev BenZvi Department of Physics and Astronomy University of Rochester Table of Contents 1 Simulation and Random Number Generation Simulation of Physical Systems Creating
More informationMonte Carlo Methods in High Energy Physics I
Helmholtz International Workshop -- CALC 2009, July 10--20, Dubna Monte Carlo Methods in High Energy Physics CALC2009 - July 20 10, Dubna 2 Contents 3 Introduction Simple definition: A Monte Carlo technique
More informationLehmer Random Number Generators: Introduction
Lehmer Random Number Generators: Introduction Revised version of the slides based on the book Discrete-Event Simulation: a first course LL Leemis & SK Park Section(s) 21, 22 c 2006 Pearson Ed, Inc 0-13-142917-5
More information( ) ( ) Monte Carlo Methods Interested in. E f X = f x d x. Examples:
Monte Carlo Methods Interested in Examples: µ E f X = f x d x Type I error rate of a hypothesis test Mean width of a confidence interval procedure Evaluating a likelihood Finding posterior mean and variance
More informationCDF top quark " $ )(! # % & '
$% CDF quark 7 3 5 ( "#! Tevatron Run II Started Spring 1. proton-antiproton collider with (Run I :. antiproton recycler commissioning electron cooling operational by Summer 5. increase in luminosity.
More informationA Repetition Test for Pseudo-Random Number Generators
Monte Carlo Methods and Appl., Vol. 12, No. 5-6, pp. 385 393 (2006) c VSP 2006 A Repetition Test for Pseudo-Random Number Generators Manuel Gil, Gaston H. Gonnet, Wesley P. Petersen SAM, Mathematik, ETHZ,
More informationEfficiency of the CMS Level-1 Trigger to. Selected Physics Channels
Efficiency of the CMS Level-1 Trigger to Selected Physics Channels C. Sulkko University of Colorado, Boulder, CO D. Acosta University of Florida, Gainesville, FL Abstract The efficiency of the Compact
More informationGenerating Uniform Random Numbers
1 / 43 Generating Uniform Random Numbers Christos Alexopoulos and Dave Goldsman Georgia Institute of Technology, Atlanta, GA, USA March 1, 2016 2 / 43 Outline 1 Introduction 2 Some Generators We Won t
More informationUniform random numbers generators
Uniform random numbers generators Lecturer: Dmitri A. Moltchanov E-mail: moltchan@cs.tut.fi http://www.cs.tut.fi/kurssit/tlt-2707/ OUTLINE: The need for random numbers; Basic steps in generation; Uniformly
More informationContents. 1 Probability review Introduction Random variables and distributions Convergence of random variables...
Contents Probability review. Introduction..............................2 Random variables and distributions................ 3.3 Convergence of random variables................. 6 2 Monte Carlo methods
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 informationPhysics at Hadron Colliders
Physics at Hadron Colliders Part 2 Standard Model Physics Test of Quantum Chromodynamics - Jet production - W/Z production - Production of Top quarks Precision measurements -W mass - Top-quark mass QCD
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 informationCPSC 531: Random Numbers. Jonathan Hudson Department of Computer Science University of Calgary
CPSC 531: Random Numbers Jonathan Hudson Department of Computer Science University of Calgary http://www.ucalgary.ca/~hudsonj/531f17 Introduction In simulations, we generate random values for variables
More informationarxiv: v1 [hep-ex] 21 Aug 2011
arxiv:18.155v1 [hep-ex] 1 Aug 011 Early Searches with Jets with the ATLAS Detector at the LHC University of Chicago, Enrico Fermi Institute E-mail: georgios.choudalakis@cern.ch We summarize the analysis
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 informationGenerating pseudo- random numbers
Generating pseudo- random numbers What are pseudo-random numbers? Numbers exhibiting statistical randomness while being generated by a deterministic process. Easier to generate than true random numbers?
More informationClass 12. Random Numbers
Class 12. Random Numbers NRiC 7. Frequently needed to generate initial conditions. Often used to solve problems statistically. How can a computer generate a random number? It can t! Generators are pseudo-random.
More informationHow does the computer generate observations from various distributions specified after input analysis?
1 How does the computer generate observations from various distributions specified after input analysis? There are two main components to the generation of observations from probability distributions.
More informationMeasurement of the associated production of direct photons and jets with the Atlas experiment at LHC. Michele Cascella
Measurement of the associated production of direct photons and jets with the Atlas experiment at LHC Michele Cascella Graduate Course in Physics University of Pisa The School of Graduate Studies in Basic
More informationStatistical Methods in Particle Physics
Statistical Methods in Particle Physics 4. Monte Carlo Methods Prof. Dr. Klaus Reygers (lectures) Dr. Sebastian Neubert (tutorials) Heidelberg University WS 2017/18 Monte Carlo Method Any method which
More informationS.A. Teukolsky, Computers in Physics, Vol. 6, No. 5,
Physics 75.502 Part III: Monte Carlo Methods 139 Part III: Monte Carlo Methods Topics: Introduction Random Number generators Special distributions General Techniques Multidimensional simulation References:
More informationJet Energy Calibration. Beate Heinemann University of Liverpool
Jet Energy Calibration Beate Heinemann University of Liverpool Fermilab, August 14th 2006 1 Outline Introduction CDF and D0 calorimeters Response corrections Multiple interactions η-dependent corrections
More informationFuture prospects for the measurement of direct photons at the LHC
Future prospects for the measurement of direct photons at the LHC David Joffe on behalf of the and CMS Collaborations Southern Methodist University Department of Physics, 75275 Dallas, Texas, USA DOI:
More informationGenerating Uniform Random Numbers
1 / 41 Generating Uniform Random Numbers Christos Alexopoulos and Dave Goldsman Georgia Institute of Technology, Atlanta, GA, USA 10/13/16 2 / 41 Outline 1 Introduction 2 Some Lousy Generators We Won t
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 informationGenerating Uniform Random Numbers
1 / 44 Generating Uniform Random Numbers Christos Alexopoulos and Dave Goldsman Georgia Institute of Technology, Atlanta, GA, USA 10/29/17 2 / 44 Outline 1 Introduction 2 Some Lousy Generators We Won t
More informationRadiation Transport Calculations by a Monte Carlo Method. Hideo Hirayama and Yoshihito Namito KEK, High Energy Accelerator Research Organization
Radiation Transport Calculations by a Monte Carlo Method Hideo Hirayama and Yoshihito Namito KEK, High Energy Accelerator Research Organization Monte Carlo Method A method used to solve a problem with
More informationProperties of Proton-proton Collision and. Comparing Event Generators by Multi-jet Events
Properties of Proton-proton Collision and Comparing Event Generators by Multi-jet Events Author: W. H. TANG 1 (Department of Physics, The Chinese University of Hong Kong) Supervisors: Z. L. MARSHALL 2,
More informationRecent Results of + c + X and + b + X Production Cross Sections at DØ
Recent Results of + c + X and + b + X Production Cross Sections at DØ Florida State University Wednesday March 18th Virginia HEP Seminar 1 The Standard Model (SM) The Standard Model (SM) describes the
More informationMuon reconstruction performance in ATLAS at Run-2
2 Muon reconstruction performance in ATLAS at Run-2 Hannah Herde on behalf of the ATLAS Collaboration Brandeis University (US) E-mail: hannah.herde@cern.ch ATL-PHYS-PROC-205-2 5 October 205 The ATLAS muon
More informationA brief history of accelerators, detectors and experiments: (See Chapter 14 and Appendix H in Rolnick.)
Physics 557 Lecture 7 A brief history of accelerators, detectors and experiments: (See Chapter 14 and Appendix H in Rolnick.) First came the study of the debris from cosmic rays (the God-given particle
More informationThe ATLAS Experiment and the CERN Large Hadron Collider
The ATLAS Experiment and the CERN Large Hadron Collider HEP101-4 February 20, 2012 Al Goshaw 1 HEP 101 Today Introduction to HEP units Particles created in high energy collisions What can be measured in
More informationLecture 20. Randomness and Monte Carlo. J. Chaudhry. Department of Mathematics and Statistics University of New Mexico
Lecture 20 Randomness and Monte Carlo J. Chaudhry Department of Mathematics and Statistics University of New Mexico J. Chaudhry (UNM) CS 357 1 / 40 What we ll do: Random number generators Monte-Carlo integration
More informationThe LHCf data hadronic interactions and UHECR showers. Paolo Lipari LHCf meeting Catania, 6th july 2011
The LHCf data hadronic interactions and UHECR showers Paolo Lipari LHCf meeting Catania, 6th july 2011 ~50 years of UHECR Problems of determination of: Energy Mass A Hadronic interaction Modeling Measure
More information2 The ATLAS Experiment at the Large Hadron Collider at CERN
Studies of the Di-Muons Resonances at the ATLAS Experiment at CERN PhD Detailed Research Project PhD Candidate: Camilla Maiani Supervisor: Prof. Carlo Dionisi, Dott. Stefano Giagu Università di Roma La
More informationQCD and jets physics at the LHC with CMS during the first year of data taking. Pavel Demin UCL/FYNU Louvain-la-Neuve
QCD and jets physics at the LHC with CMS during the first year of data taking Pavel Demin UCL/FYNU Louvain-la-Neuve February 8, 2006 Bon appétit! February 8, 2006 Pavel Demin UCL/FYNU 1 Why this seminar?
More information( x) ( ) F ( ) ( ) ( ) Prob( ) ( ) ( ) X x F x f s ds
Applied Numerical Analysis Pseudo Random Number Generator Lecturer: Emad Fatemizadeh What is random number: A sequence in which each term is unpredictable 29, 95, 11, 60, 22 Application: Monte Carlo Simulations
More informationAlgorithms and Networking for Computer Games
Algorithms and Networking for Computer Games Chapter 2: Random Numbers http://www.wiley.com/go/smed What are random numbers good for (according to D.E. Knuth) simulation sampling numerical analysis computer
More informationDiscovery potential of the SM Higgs with ATLAS
Discovery potential of the SM Higgs with P. Fleischmann On behalf of the Collaboration st October Abstract The discovery potential of the Standard Model Higgs boson with the experiment at the Large Hadron
More informationATLAS NOTE. August 25, Electron Identification Studies for the Level 1 Trigger Upgrade. Abstract
Draft version 1.0 ATLAS NOTE August 25, 2012 1 Electron Identification Studies for the Level 1 Trigger Upgrade 2 3 4 L. Feremenga a, M.-A. Pleier b, F. Lanni b a University of Texas at Arlington b Brookhaven
More informationSources of randomness
Random Number Generator Chapter 7 In simulations, we generate random values for variables with a specified distribution Ex., model service times using the exponential distribution Generation of random
More informationV0 cross-section measurement at LHCb. RIVET analysis module for Z boson decay to di-electron
V0 cross-section measurement at LHCb. RIVET analysis module for Z boson decay to di-electron Outline of the presentation: 1. Introduction to LHCb physics and LHCb detector 2. RIVET plug-in for Z e+e- channel
More informationStudy of Multimuon Events at CDF
Study of Multimuon Events at CDF Kevin Pitts University of Illinois Outline History: b production and decay puzzles from the 1990s Recent results Inclusive B cross sections bb cross section New study of
More informationLa ricerca dell Higgs Standard Model a CDF
La ricerca dell Higgs Standard Model a CDF Melisa Rossi INFN-TS Giornata di seminari INFN Trieste - 7 Luglio 2009 FNAL: Fermi National Accelerator Lab Tevatron currently provides the highest energy proton-antiproton
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 informationMeasurement of the Inclusive Isolated Prompt Photon Cross Section at CDF
of the Inclusive Isolated Cross at IFAE Barcelona HEP Seminar University of Virginia Outline Theoretical introduction Prompt photon production The The Photon detection prediction The pqcd NLO prediction
More informationMATH5835 Statistical Computing. Jochen Voss
MATH5835 Statistical Computing Jochen Voss September 27, 2011 Copyright c 2011 Jochen Voss J.Voss@leeds.ac.uk This text is work in progress and may still contain typographical and factual mistakes. Reports
More informationModern Methods of Data Analysis - WS 07/08
Modern Methods of Data Analysis Lecture V (12.11.07) Contents: Central Limit Theorem Uncertainties: concepts, propagation and properties Central Limit Theorem Consider the sum X of n independent variables,
More informationExclusive Central π+π- Production in Proton Antiproton Collisions at the CDF
Exclusive Central π+π- Production in Proton Antiproton Collisions at the CDF Maria Żurek Research Center Jülich University of Cologne on behalf of the CDF Collaboration XLV International Symposium on Multiparticle
More informationLecture 5: Random numbers and Monte Carlo (Numerical Recipes, Chapter 7) Motivations for generating random numbers
Lecture 5: Random numbers and Monte Carlo (Numerical Recipes, Chapter 7) Motivations for generating random numbers To sample a function in a statistically controlled manner (i.e. for Monte Carlo integration)
More informationarxiv: v1 [hep-ex] 6 Jul 2007
Muon Identification at ALAS and Oliver Kortner Max-Planck-Institut für Physik, Föhringer Ring, D-005 München, Germany arxiv:0707.0905v1 [hep-ex] Jul 007 Abstract. Muonic final states will provide clean
More informationWhat detectors measure
What detectors measure As a particle goes through matter, it releases energy Detectors collect the released energy and convert it to electric signals recorded by DAQ Raw event record is a collection of
More informationTopics in Computer Mathematics
Random Number Generation (Uniform random numbers) Introduction We frequently need some way to generate numbers that are random (by some criteria), especially in computer science. Simulations of natural
More informationATLAS EXPERIMENT : HOW THE DATA FLOWS. (Trigger, Computing, and Data Analysis)
ATLAS EXPERIMENT : HOW THE DATA FLOWS (Trigger, Computing, and Data Analysis) In order to process large volumes of data within nanosecond timescales, the trigger system is designed to select interesting
More informationThe reaction p(e,e'p)π 0 to calibrate the Forward and the Large Angle Electromagnetic Shower Calorimeters
The reaction p(e,e'p)π 0 to calibrate the Forward and the Large Angle Electromagnetic Shower Calorimeters M.Battaglieri, M.Anghinolfi, P.Corvisiero, A.Longhi, M.Ripani, M.Taiuti Istituto Nazionale di Fisica
More informationExperimental Particle Physics Informal Lecture & Seminar Series Lecture 1 Detectors Overview
Experimental Particle Physics Informal Lecture & Seminar Series 2013 Lecture 1 Detectors Overview Detectors in Particle Physics Let s talk about detectors for a bit. Let s do this with Atlas and CMS in
More informationMonte Carlo integration (naive Monte Carlo)
Monte Carlo integration (naive Monte Carlo) Example: Calculate π by MC integration [xpi] 1/21 INTEGER n total # of points INTEGER i INTEGER nu # of points in a circle REAL x,y coordinates of a point in
More informationReconstruction in Collider Experiments (Part IX)
Introduction to Hadronic Final State Reconstruction in Collider Experiments Introduction to Hadronic Final State Reconstruction in Collider Experiments (Part IX) Peter Loch University of Arizona Tucson,
More informationJet Reconstruction and Energy Scale Determination in ATLAS
Jet Reconstruction and Energy Scale Determination in ATLAS Ariel Schwartzman 3 rd Top Workshop: from the Tevatron to ATLAS Grenoble, 23-Oct-2008 1 Outline ATLAS Calorimeters Calorimeter signal reconstruction:
More informationEvidence for Single Top Quark Production. Reinhard Schwienhorst
Evidence for Single Top Quark Production Reinhard Schwienhorst MSU high energy seminar, 1/9/2007 Outline Motivation Preparation Optimized Event Analysis Sensitivity Cross section measurement Vtb Conclusions/Outlook
More informationMeasurement of jet production in association with a Z boson at the LHC & Jet energy correction & calibration at HLT in CMS
Measurement of jet production in association with a Z boson at the LHC & Jet energy correction & calibration at HLT in CMS Fengwangdong Zhang Peking University (PKU) & Université Libre de Bruxelles (ULB)
More informationSusanna Costanza. (Università degli Studi di Pavia & INFN Pavia) on behalf of the ALICE Collaboration
(Università degli Studi di Pavia & INFN Pavia) on behalf of the ALICE Collaboration 102 Congresso della Società Italiana di Fisica Padova, 26-30 settembre 2016 Outline Heavy flavour physics in ALICE The
More informationRWTH Aachen Graduiertenkolleg
RWTH Aachen Graduiertenkolleg 9-13 February, 2009 Glen Cowan Physics Department Royal Holloway, University of London g.cowan@rhul.ac.uk www.pp.rhul.ac.uk/~cowan Course web page: www.pp.rhul.ac.uk/~cowan/stat_aachen.html
More informationFeasibility of a cross-section measurement for J/ψ->ee with the ATLAS detector
Feasibility of a cross-section measurement for J/ψ->ee with the ATLAS detector ATLAS Geneva physics meeting Andrée Robichaud-Véronneau Outline Motivation Theoretical background for J/ψ hadroproduction
More informationValidation of Geant4 Physics Models Using Collision Data from the LHC
Journal of Physics: Conference Series Validation of Geant4 Physics Models Using Collision from the LHC To cite this article: S Banerjee and CMS Experiment 20 J. Phys.: Conf. Ser. 33 032003 Related content
More informationLecture 2 & 3. Particles going through matter. Collider Detectors. PDG chapter 27 Kleinknecht chapters: PDG chapter 28 Kleinknecht chapters:
Lecture 2 & 3 Particles going through matter PDG chapter 27 Kleinknecht chapters: 1.2.1 for charged particles 1.2.2 for photons 1.2.3 bremsstrahlung for electrons Collider Detectors PDG chapter 28 Kleinknecht
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 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 informationConference Report Mailing address: CMS CERN, CH-1211 GENEVA 23, Switzerland
Available on the CMS information server CMS CR 212/178 The Compact Muon Solenoid Experiment Conference Report Mailing address: CMS CERN, CH211 GENEVA 23, Switzerland 212//9 Measurement of isolated photon
More informationNA62: Ultra-Rare Kaon Decays
NA62: Ultra-Rare Kaon Decays Phil Rubin George Mason University For the NA62 Collaboration November 10, 2011 The primary goal of the experiment is to reconstruct more than 100 K + π + ν ν events, over
More informationCALICE scintillator HCAL
CALICE scintillator HCAL Erika Garutti DESY (on behalf of the CALICE collaboration) OUTLINE: electromagnetic and hadronic shower analysis shower separation The test beam prototypes 10 GeV pion shower @
More informationLHC MPI and underlying event results
LHC MPI and underlying event results Marianna Testa (ATLAS) for the ALICE, ATLAS, CMS & LHCb collaborations SM@LHC 2012, Copenhagen, 10-13 April 2012 The Underlying Event Everything in hadron collisions
More informationarxiv: v1 [nucl-ex] 7 Sep 2009
CCπ Event Reconstruction at MiniBooNE R.H. Nelson University of Colorado, Dept. of Physics, 39 UCB, Boulder, CO 839, USA arxiv:99.238v [nucl-ex] 7 Sep 29 Abstract. We describe the development of a fitter
More informationSpectroscopy and Decay results from CDF
Quarkonium Spectroscopy and Decay results from CDF KIT Quarkonium Workshop December 3, 2008 Outline Tevatron and CDF Bc Mass Lifetime bc X(3872) Mass splitting and mass X page 2 Tevatron p s = 1.96 TeV
More informationMonte Carlo Event Generators
YETI 06 SM IPPP, Durham, UK 27 29 March 2006 Monte Carlo Event Generators Torbjörn Sjöstrand Lund University 1. (today) Introduction and Overview; Monte Carlo Techniques 2. (today) Matrix Elements; Parton
More informationGEANT4 simulation of the testbeam set-up for the ALFA detector
GEANT4 simulation of the testbeam set-up for the detector V. Vorobel a and H. Stenzel b a Faculty of Mathematics and Physics, Charles University in Prague, Czech Republic b II. Physikalisches Institut,
More informationTwo Early Exotic searches with dijet events at ATLAS
ATL-PHYS-PROC-2011-022 01/06/2011 Two Early Exotic searches with dijet events at ATLAS University of Toronto, Department of Physics E-mail: rrezvani@physics.utoronto.ca This document summarises two exotic
More informationAccelerators and Colliders
Accelerators and Colliders References Robert Mann: An introduction to particle physics and the standard model Tao Han, Collider Phenomenology, http://arxiv.org/abs/hep-ph/0508097 Particle Data Group, (J.
More informationQCD Measurements at DØ
QCD Measurements at DØ University of Texas -Arlington- Seminar, University of Virginia January 24th 27 TeVatron Collider at Fermilab Proton-Antiproton Collisions at Center-of-Mass s = 1.96 TeV Two Multi-Purpose
More informationLecture 3 Cross Section Measurements. Ingredients to a Cross Section
Lecture 3 Cross Section Measurements Ingredients to a Cross Section Prerequisites and Reminders... Natural Units Four-Vector Kinematics Lorentz Transformation Lorentz Boost Lorentz Invariance Rapidity
More informationIE 581 Introduction to Stochastic Simulation. One page of notes, front and back. Closed book. 50 minutes. Score
One page of notes, front and back. Closed book. 50 minutes. Score Schmeiser Page 1 of 4 Test #1, Spring 2001 1. True or false. (If you wish, write an explanation of your thinking.) (a) T Data are "binary"
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 informationQCD Studies at LHC with the Atlas detector
QCD Studies at LHC with the Atlas detector Introduction Sebastian Eckweiler - University of Mainz (on behalf of the ATLAS Collaboration) Examples of QCD studies Minimum bias & underlying event Jet-physics
More informationMeasurement of multijets and the internal structure of jets at ATLAS
Measurement of multijets and the internal structure of jets at ALAS Bilge M. Demirköz, on behalf of the ALAS Collaboration Institut de Fisica d Altes Energies, Barcelona, SPAIN also at Middle East echnical
More informationPart II: Detectors. Peter Schleper Universität Hamburg
Part II: Detectors Peter Schleper Universität Hamburg 30.05.2018 Outline of the lecture: 1. Overview on detectors 2. Particle interactions with matter 3. Scintillators and photon detectors 4. Semiconductor
More informationQCD Studies with CMS at LHC. Gunther Roland for the Collaboration
QCD Studies with CMS at LHC Gunther Roland for the Collaboration INT Seattle 5/24/2010 First 7 TeV Collisions: March 30th 2010 2 QCD Studies with CMS pp integrated luminosity ~ 10nb -1 n.b. rates for LHC
More informationHard Scattering in Hadron-Hadron Collisions: Physics and Anatomy
Hard Scattering in Hadron-Hadron Collisions: Physics and Anatomy Section 1: Introduction, Colliders and Detectors 1. Basic anatomy of a collision 2. Collider considerations 3. Detector Implications 4.
More informationPerformance of muon and tau identification at ATLAS
ATL-PHYS-PROC-22-3 22/2/22 Performance of muon and tau identification at ATLAS On behalf of the ATLAS Collaboration University of Oregon E-mail: mansoora.shamim@cern.ch Charged leptons play an important
More informationGenerating Random Variables
Generating Random Variables These slides are created by Dr. Yih Huang of George Mason University. Students registered in Dr. Huang's courses at GMU can make a single machine-readable copy and print a single
More informationLHC Detectors and their Physics Potential. Nick Ellis PH Department, CERN, Geneva
LHC Detectors and their Physics Potential Nick Ellis PH Department, CERN, Geneva 1 Part 1 Introduction to the LHC Detector Requirements & Design Concepts 2 What is the Large Hadron Collider? Circular proton-proton
More informationRandom Number Generation. CS1538: Introduction to simulations
Random Number Generation CS1538: Introduction to simulations Random Numbers Stochastic simulations require random data True random data cannot come from an algorithm We must obtain it from some process
More information