Modern Methods of Data Analysis - WS 07/08
|
|
- Dwight Jordan
- 5 years ago
- Views:
Transcription
1 Modern Methods of Data Analysis Lecture V ( ) Contents: Central Limit Theorem Uncertainties: concepts, propagation and properties
2 Central Limit Theorem Consider the sum X of n independent variables, with i = 1,2,3,..., each taken from a different distribution with expectation value and variance Then the distribution for properties: has the following Its expectation value is Its variance is It becomes Gaussian distributed for n
3 Uniform distribution
4 Uneven Distribution
5 Exponential Distribution
6 But... In order to work No source contributes significant to overall variance multiple scattering Distribution with few events in large tails need many more statistic to converge. CLT unfortunately doesn't work for many physics applications!
7 Application: Many repeated Measurements
8 m(b0) = ± 0.53 (stat) ± 0.33 (sys) CDF has a mass resolution of 16 MeV the reconstructed mass of a single B meson is spread around the true B mass with σ=16 MeV The B mass can be measured with way better precision
9
10
11 Law of large numbers...
12 Weighted Mean (I) Combining measurements with different uncertainties: Twice time measurement of the seed of a car: v1 = 67 ± 4 m/s v2 = 62 ± 2 m/s Uncertainty on mean is larger than single uncertainty...???
13 Weighted Mean (II) To reach σ = 2 m/s, 4 single measurements with σ = 4 m/s are needed. Therefore this measurement should get 4 times the weight! More general:
14 Weighted Mean (III) v1 = 67 ± 4 m/s v2 = 53 ± 2 m/s????
15 New Scientist, 31 March 1988
16 Neutron Lifetime (PDG 2006)
17 Particle Data Group: World Average 1) Calculate weighted mean 2) Calculate there are 3 cases to consider < 1: all fine, simple weighted mean is OK 1: Depending on the reason, either make no average at all, or quote calculated average and make educated guess of error, taking into account known problems with data > 1: Errors on some or all measurements may have been underestimated, scale all of them with: to compute S, reject ones with larger uncertainties
18 Ideogram
19 Reminder: Covariance covariance of two variables x, y is defined as: if two variables are uncorrelated cov(x,y) = 0 If cov(x,y)=0 x,y are uncorrelated Correlation is defined as:
20 Properties Correlations Example of correlations for two random variables x and y the covariance V(x,y) or cov(x,y) can be represented by a matrix (often called error matrix ) General case of correlations for n random variables, there is covariance between each pair of variables analogous, correlation matrix is defined
21 Error Propagation (I) x= Vi,j and µi known y(x) is function of first order Taylor expansion...
22 Error Propagation (II)
23 Gaussian error propagation Error estimates for functions of several correlated variables : Normal errors for uncorrelated variables Additional terms accounting for correlations Special case, uncorrelated variables: This is called Gaussian error propagation, however has nothing to do with Gaussian distributions
24 And the same in more dimensions (A is Jacobi matrix)
25 Example: Track Parametrization (I) Typical tracking chamber measures 3D, due to symmetry uses cylindrical coordinates: (r,φ,z) Want to know what are the uncertainties in Cartesian coordinates (x,y,z): x = r cosφ; y = r sinφ
26 Example: Track Parametrization(II)
27 Exercise: x,y are two random variables with the correlation matrix V: a = 3/7 x +1/7 y; b = 1/7 x - 2/7y; Give the correlation matrix of a,b (note this time it is not an approximation...)
28 Orthogonal Transformation For n variables one can always find a linear transformation such that the in the covariance matrix of the new set of variables is diagonal.
29 Exercise: Some standard formulae
30 Example: Measure Asymmetry (I) Define A = (F-B)/(F+B) an asymmetry of events in forward and backward hemisphere of the detector LEP). Here F (B) is the number of events in forward (backward) hemisphere Case I: If the events in forward and backward hemisphere are uncorrelated then: if errors are individually Poisson distributed for F and B, then errors are dominated by smaller counting rate
31 Example: Measure Asymmetry (II) case II: events in forward and backward hemisphere are completely anti-correlated, N is thus fixed - distribution of events F and B are then given by Binomial distribution; let p be probability that event is in forward h.s.: - this is the same expression as before demonstrates relationship between Poisson & Binomial - either Poisson prob of obtaining N events altogether times binomial prob. of having F events in forward - or: two independent Poisson prob. in the number of backward and forward events
32 Repetition: Histogram Interpretation of bins A histogram can be equivalently regarded as: 1. A Poisson distribution in the overall number of events N and the corresponding multinomial distribution of obtaining events in each bin 2. An independent Poisson distribution of the number of entries in each bin of the histogram.
33 Exercise: Detector Efficiency Compute the error on measuring detector efficiency in two different ways: Binomial distribution; p: probability to detect traversing particle, N number of events Using error propagation, : number of detected events :number of not detected events, treat as independent variables
34 Be aware... The approximation using Taylor expansion breaks down if the function is significantly not linear in the region ± 1σ around the mean value. Example: momentum estimate in B field; p ~ 1/κ 10 % momentum bias!
35 st Failure of (1 order) error propagation experiment had data to look for non-zero mass of electron neutrino Quantity R was measured: Don't need to know details. Sufficient to know: a, b, c, d, e are measured quantities and K is constant. If R < 0.42, then neutrino must have mass. The experiment measured R=0.165 and found with error propagation σ(r) = > 3σ evidence for neutrino mass! however, the formula is highly non-linear... 1st order error propagation not applicable!
36 What to do instead? Use MC methods! Throw Gaussian distributed values for a,b,c,d,e Compute R; repeat toy many times and check how often you are below the measured value of R. In example this happens in 4% of the cases. This is the so-called p-value of this result. In many physics analysis, simple Gaussian error propagation not valid or too complicated... (e.g. highly non linear functions, many correlated variables) p-value MC method always works!
Modern Methods of Data Analysis - WS 07/08
Modern Methods of Data Analysis Lecture VIa (19.11.07) Contents: Uncertainties (II): Re: error propagation Correlated uncertainties Systematic uncertainties Re: Error Propagation (I) x = Vi,j and µi known
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 informationModern Methods of Data Analysis - WS 07/08
Modern Methods of Data Analysis Lecture VII (26.11.07) Contents: Maximum Likelihood (II) Exercise: Quality of Estimators Assume hight of students is Gaussian distributed. You measure the size of N students.
More informationStatistics. Lent Term 2015 Prof. Mark Thomson. 2: The Gaussian Limit
Statistics Lent Term 2015 Prof. Mark Thomson Lecture 2 : The Gaussian Limit Prof. M.A. Thomson Lent Term 2015 29 Lecture Lecture Lecture Lecture 1: Back to basics Introduction, Probability distribution
More informationStatistics for Data Analysis. Niklaus Berger. PSI Practical Course Physics Institute, University of Heidelberg
Statistics for Data Analysis PSI Practical Course 2014 Niklaus Berger Physics Institute, University of Heidelberg Overview You are going to perform a data analysis: Compare measured distributions to theoretical
More informationLectures on Statistical Data Analysis
Lectures on Statistical Data Analysis London Postgraduate Lectures on Particle Physics; University of London MSci course PH4515 Glen Cowan Physics Department Royal Holloway, University of London g.cowan@rhul.ac.uk
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 informationE. Santovetti lesson 4 Maximum likelihood Interval estimation
E. Santovetti lesson 4 Maximum likelihood Interval estimation 1 Extended Maximum Likelihood Sometimes the number of total events measurements of the experiment n is not fixed, but, for example, is a Poisson
More informationYETI IPPP Durham
YETI 07 @ IPPP Durham Young Experimentalists and Theorists Institute 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_yeti.html
More informationPhysics 6720 Introduction to Statistics April 4, 2017
Physics 6720 Introduction to Statistics April 4, 2017 1 Statistics of Counting Often an experiment yields a result that can be classified according to a set of discrete events, giving rise to an integer
More informationLecture 3. G. Cowan. Lecture 3 page 1. Lectures on Statistical Data Analysis
Lecture 3 1 Probability (90 min.) Definition, Bayes theorem, probability densities and their properties, catalogue of pdfs, Monte Carlo 2 Statistical tests (90 min.) general concepts, test statistics,
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 informationNumerical Methods Lecture 7 - Statistics, Probability and Reliability
Topics Numerical Methods Lecture 7 - Statistics, Probability and Reliability A summary of statistical analysis A summary of probability methods A summary of reliability analysis concepts Statistical Analysis
More informationPHYSICS 2150 EXPERIMENTAL MODERN PHYSICS. Lecture 3 Rejection of Data; Weighted Averages
PHYSICS 15 EXPERIMENTAL MODERN PHYSICS Lecture 3 Rejection of Data; Weighted Averages PREVIOUS LECTURE: GAUSS DISTRIBUTION 1.5 p(x µ, )= 1 e 1 ( x µ ) µ=, σ=.5 1. µ=3, σ=.5.5 µ=4, σ=1 4 6 8 WE CAN NOW
More informationLecture 2: Repetition of probability theory and statistics
Algorithms for Uncertainty Quantification SS8, IN2345 Tobias Neckel Scientific Computing in Computer Science TUM Lecture 2: Repetition of probability theory and statistics Concept of Building Block: Prerequisites:
More informationStatistics and Data Analysis
Statistics and Data Analysis The Crash Course Physics 226, Fall 2013 "There are three kinds of lies: lies, damned lies, and statistics. Mark Twain, allegedly after Benjamin Disraeli Statistics and Data
More informationUse of the likelihood principle in physics. Statistics II
Use of the likelihood principle in physics Statistics II 1 2 3 + Bayesians vs Frequentists 4 Why ML does work? hypothesis observation 5 6 7 8 9 10 11 ) 12 13 14 15 16 Fit of Histograms corresponds This
More informationExperimental Signatures for Lorentz Violation (in Neutral B Meson Mixing)
Experimental Signatures for Lorentz Violation (in Neutral B Meson Mixing) Theory recap (see Mike Berger's talk for details) Experimental situation Experimental issues: particle/antiparticle differences
More informationElectroweak Physics at the Tevatron
Electroweak Physics at the Tevatron Adam Lyon / Fermilab for the DØ and CDF collaborations 15 th Topical Conference on Hadron Collider Physics June 2004 Outline Importance Methodology Single Boson Measurements
More informationPhysics Sep Example A Spin System
Physics 30 7-Sep-004 4- Example A Spin System In the last lecture, we discussed the binomial distribution. Now, I would like to add a little physical content by considering a spin system. Actually this
More informationRare B Meson Decays at Tevatron
Rare B Meson Decays at Tevatron Walter Hopkins CDF and DØ Collaborations Cornell University FPCP 2011 Walter Hopkins (CDF and DØ Collaborations, Cornell University)Rare B Meson Decays at Tevatron FPCP
More informationLAB 2 1. Measurement of 2. Binomial Distribution
LAB 2 Gan Phys 3700 1. Measurement of π In this exercise we will determine a value for π by throwing darts: a) Determine π by throwing a dart 100 or more times. Use an 8x11 inch sheet of paper with a circle
More informationStatistical Data Analysis 2017/18
Statistical Data Analysis 2017/18 London Postgraduate Lectures on Particle Physics; University of London MSci course PH4515 Glen Cowan Physics Department Royal Holloway, University of London g.cowan@rhul.ac.uk
More informationb hadron properties and decays (ATLAS)
b hadron properties and decays (ATLAS) Milind V. Purohit Univ. of South Carolina, USA for the ATLAS collaboration M. V. Purohit SM@LHC 2012 1 Outline Detector, Di muon trigger, Performance: Mass, vertex
More informationStatistics, Probability Distributions & Error Propagation. James R. Graham
Statistics, Probability Distributions & Error Propagation James R. Graham Sample & Parent Populations Make measurements x x In general do not expect x = x But as you take more and more measurements a pattern
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 informationPHYSICS 2150 LABORATORY
PHYSICS 2150 LABORATORY Instructors: John Cumalat Jiayan Pheonix Dai Lab Coordinator: Jerry Leigh Lecture 2 September 2, 2008 PHYS2150 Lecture 2 Need to complete the Radiation Certification The Gaussian
More informationHQL Virginia Tech. Bob Hirosky for the D0 Collaboration. Bob Hirosky, UNIVERSITY of VIRGINIA. 26May, 2016
Bs CP-odd lifetime in Bs J/ψf0 and Afb for baryons at D0 2016 Virginia Tech Bob Hirosky for the D0 Collaboration 1 Tevatron Data D0 continues a rich physics program analyzing ~10fb-1 of recorded data from
More informationMeasuring CP violation in. B s φφ with LHCb. Jim Libby (University of Oxford) for the LHCb collaboration. 14/12/2006 CKM 2006 Nagoya 1
Measuring CP violation in B s φφ with LHCb Jim Libby (University of Oxford) for the LHCb collaboration 4//6 CKM 6 Nagoya Outline Motivation for the study B s φφ Brief introduction to LHCb Expected signal
More informationThe Solar Neutrino Day-Night Effect. Master of Science Thesis Mattias Blennow Division of Mathematical Physics Department of Physics KTH
The Solar Neutrino Day-Night Effect Master of Science Thesis Mattias Blennow Division of Mathematical Physics Department of Physics KTH 1 Why This Interest in Neutrinos? Massless in SM of particle physics
More informationECT Lecture 2. - Reactor Antineutrino Detection - The Discovery of Neutrinos. Thierry Lasserre (Saclay)
ECT Lecture 2 - Reactor Antineutrino Detection - The Discovery of Neutrinos Thierry Lasserre (Saclay) Reactor Neutrino Detection Inverse Beta Decay p + anti-v e à e + + n cross section @2 MeV : 5 10-43
More informationAlgorithms for Uncertainty Quantification
Algorithms for Uncertainty Quantification Tobias Neckel, Ionuț-Gabriel Farcaș Lehrstuhl Informatik V Summer Semester 2017 Lecture 2: Repetition of probability theory and statistics Example: coin flip Example
More informationRecent Heavy Flavors results from Tevatron. Aleksei Popov (Institute for High Energy Physics, Protvino) on behalf of the CDF and DØ Collaborations
Recent Heavy Flavors results from Tevatron Aleksei Popov (Institute for High Energy Physics, Protvino) on behalf of the CDF and DØ Collaborations March 27, 2017 Outline Tevatron, CDF and DØ Confirmation
More informationLecture 5. G. Cowan Lectures on Statistical Data Analysis Lecture 5 page 1
Lecture 5 1 Probability (90 min.) Definition, Bayes theorem, probability densities and their properties, catalogue of pdfs, Monte Carlo 2 Statistical tests (90 min.) general concepts, test statistics,
More informationScientific Measurement
Scientific Measurement SPA-4103 Dr Alston J Misquitta Lecture 5 - The Binomial Probability Distribution Binomial Distribution Probability of k successes in n trials. Gaussian/Normal Distribution Poisson
More informationFYST17 Lecture 8 Statistics and hypothesis testing. Thanks to T. Petersen, S. Maschiocci, G. Cowan, L. Lyons
FYST17 Lecture 8 Statistics and hypothesis testing Thanks to T. Petersen, S. Maschiocci, G. Cowan, L. Lyons 1 Plan for today: Introduction to concepts The Gaussian distribution Likelihood functions Hypothesis
More informationStat Lecture 20. Last class we introduced the covariance and correlation between two jointly distributed random variables.
Stat 260 - Lecture 20 Recap of Last Class Last class we introduced the covariance and correlation between two jointly distributed random variables. Today: We will introduce the idea of a statistic and
More informationChapter 5 continued. Chapter 5 sections
Chapter 5 sections Discrete univariate distributions: 5.2 Bernoulli and Binomial distributions Just skim 5.3 Hypergeometric distributions 5.4 Poisson distributions Just skim 5.5 Negative Binomial distributions
More informationParity violation. no left-handed ν$ are produced
Parity violation Wu experiment: b decay of polarized nuclei of Cobalt: Co (spin 5) decays to Ni (spin 4), electron and anti-neutrino (spin ½) Parity changes the helicity (H). Ø P-conservation assumes a
More informationModern Methods of Data Analysis - WS 07/08
Modern Methods of Data Analysis Lecture VIc (19.11.07) Contents: Maximum Likelihood Fit Maximum Likelihood (I) Assume N measurements of a random variable Assume them to be independent and distributed according
More informationLower Bound Techniques for Statistical Estimation. Gregory Valiant and Paul Valiant
Lower Bound Techniques for Statistical Estimation Gregory Valiant and Paul Valiant The Setting Given independent samples from a distribution (of discrete support): D Estimate # species Estimate entropy
More informationPHYSICS 2150 LABORATORY
PHYSICS 2150 LABORATORY Instructors: Noel Clark James G. Smith Eric D. Zimmerman Lab Coordinator: Jerry Leigh Lecture 2 January 22, 2008 PHYS2150 Lecture 2 Announcements/comments The Gaussian distribution
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 informationStatistical Data Analysis Stat 3: p-values, parameter estimation
Statistical Data Analysis Stat 3: p-values, parameter estimation London Postgraduate Lectures on Particle Physics; University of London MSci course PH4515 Glen Cowan Physics Department Royal Holloway,
More informationLecture 8 Hypothesis Testing
Lecture 8 Hypothesis Testing Taylor Ch. 6 and 10.6 Introduction l The goal of hypothesis testing is to set up a procedure(s) to allow us to decide if a mathematical model ("theory") is acceptable in light
More informationSTATISTICS OF OBSERVATIONS & SAMPLING THEORY. Parent Distributions
ASTR 511/O Connell Lec 6 1 STATISTICS OF OBSERVATIONS & SAMPLING THEORY References: Bevington Data Reduction & Error Analysis for the Physical Sciences LLM: Appendix B Warning: the introductory literature
More informationData Analysis and Monte Carlo Methods
Lecturer: Allen Caldwell, Max Planck Institute for Physics & TUM Recitation Instructor: Oleksander (Alex) Volynets, MPP & TUM General Information: - Lectures will be held in English, Mondays 16-18:00 -
More informationMath Review Sheet, Fall 2008
1 Descriptive Statistics Math 3070-5 Review Sheet, Fall 2008 First we need to know about the relationship among Population Samples Objects The distribution of the population can be given in one of the
More informationB Hadron lifetimes using a displaced track trigger
B Hadron lifetimes using a displaced track trigger Why B Lifetimes are interesting Overview one analysis Results and prospects Sneha Malde University of Oxford 1 B Hadrons B Mesons B Baryons Name B + B
More informationLecture 13: Covariance. Lisa Yan July 25, 2018
Lecture 13: Covariance Lisa Yan July 25, 2018 Announcements Hooray midterm Grades (hopefully) by Monday Problem Set #3 Should be graded by Monday as well (instead of Friday) Quick note about Piazza 2 Goals
More informationUncertainty Quantification and Validation Using RAVEN. A. Alfonsi, C. Rabiti. Risk-Informed Safety Margin Characterization. https://lwrs.inl.
Risk-Informed Safety Margin Characterization Uncertainty Quantification and Validation Using RAVEN https://lwrs.inl.gov A. Alfonsi, C. Rabiti North Carolina State University, Raleigh 06/28/2017 Assumptions
More informationIntroduction to Statistical Methods for High Energy Physics
Introduction to Statistical Methods for High Energy Physics 2011 CERN Summer Student Lectures Glen Cowan Physics Department Royal Holloway, University of London g.cowan@rhul.ac.uk www.pp.rhul.ac.uk/~cowan
More informationProton Decays. -- motivation, status, and future prospect -- Univ. of Tokyo, Kamioka Observatory Masato Shiozawa
Proton Decays -- motivation, status, and future prospect -- Univ. of Tokyo, Kamioka Observatory Masato Shiozawa Look for Baryon number violation B number conservation is experimental subject B number conservation
More informationPolarized muon decay asymmetry measurement: status and challenges
Polarized muon decay asymmetry measurement: status and challenges Glen Marshall, for the TWIST Collaboration Muon Physics in the LHC Era Symposium at the Institute of Nuclear Theory Seattle, October 008
More informationNormal Distributions Rejection of Data + RLC Circuits. Lecture 4 Physics 2CL Summer 2011
Normal Distributions Rejection of Data + RLC Circuits Lecture 4 Physics 2CL Summer 2011 Outline Reminder of simple uncertainty propagation formulae Hidden useful formula for estimating uncertainties More
More informationFINITE DIFFERENCES. Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations.
FINITE DIFFERENCES Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations. 1. Introduction When a function is known explicitly, it is easy
More informationarxiv: v1 [hep-ex] 30 Nov 2009
On extraction of oscillation parameters Jan Sobczyk and Jakub Zmuda Institute of Theoretical Physics, University of Wroclaw, plac Maxa Borna 9, 50-204 Wroclaw, Poland (Dated: November 28, 2017) arxiv:0912.0021v1
More informationStatistics 224 Solution key to EXAM 2 FALL 2007 Friday 11/2/07 Professor Michael Iltis (Lecture 2)
NOTE : For the purpose of review, I have added some additional parts not found on the original exam. These parts are indicated with a ** beside them Statistics 224 Solution key to EXAM 2 FALL 2007 Friday
More informationIntroduction to Error Analysis
Introduction to Error Analysis Part 1: the Basics Andrei Gritsan based on lectures by Petar Maksimović February 1, 2010 Overview Definitions Reporting results and rounding Accuracy vs precision systematic
More informationCharm Baryon Studies at BABAR
W.Mader@Physik.TU-Dresden.de Institut für Kern- und Teilchenphysik Technische Universität Dresden Institutsseminar IKTP 15. Juni 006 Outline 1 Introduction The BABAR Detector and PEP-II 3 Reconstruction
More informationSearch for a heavy gauge boson W e
Search for a heavy gauge boson W e Cornell University LEPP Journal Club Seminar April 1, 2011 The LHC Machine 2 The beginning of the LHC era First collisions at 7 TeV confirmed on March 30, 2010 There
More informationGlauber modelling in high-energy nuclear collisions. Jeremy Wilkinson
Glauber modelling in high-energy nuclear collisions Jeremy Wilkinson 16/05/2014 1 Introduction: Centrality in Pb-Pb collisions Proton-proton collisions: large multiplicities of charged particles produced
More informationThe Lund Model. and some extensions. Department of Astronomy and Theoretical Physics Lund University Sölvegatan 14A, SE Lund, Sweden
The Lund Model and some extensions Torbjörn Sjöstrand Department of Astronomy and Theoretical Physics Lund University Sölvegatan 4A, SE-223 62 Lund, Sweden Workshop on Collective effects in small collisions
More informationBrandon C. Kelly (Harvard Smithsonian Center for Astrophysics)
Brandon C. Kelly (Harvard Smithsonian Center for Astrophysics) Probability quantifies randomness and uncertainty How do I estimate the normalization and logarithmic slope of a X ray continuum, assuming
More informationIntroduction to Statistics and Error Analysis II
Introduction to Statistics and Error Analysis II Physics116C, 4/14/06 D. Pellett References: Data Reduction and Error Analysis for the Physical Sciences by Bevington and Robinson Particle Data Group notes
More informationIntroduction to least square fits for Geant 4 Simulation and ROOT Analysis of a Silicon Beam Telescope Mar , DESY Hamburg Olaf Behnke, DESY
Introduction to least square fits for Geant 4 Simulation and ROOT Analysis of a Silicon Beam Telescope Mar 4-6 2, DESY Hamburg Olaf Behnke, DESY Literature: Roger Barlow: Statistics, A Guide To The Use
More informationFall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede
Hypothesis Testing: Suppose we have two or (in general) more simple hypotheses which can describe a set of data Simple means explicitly defined, so if parameters have to be fitted, that has already been
More informationThe W-mass Measurement at CDF
2010-05 - 10 The W-mass Measurement at CDF Ilija Bizjak, University College London 1/33 Outline 1) Motivation for a W mass measurement Implications for the EW constraints on Higgs mass 2) Measurement of
More informationMeasurements of Particle Production in pp-collisions in the Forward Region at the LHC
2011 Europhysics Conference on High-Energy Physics, Grenoble Measurements of Particle Production in pp-collisions in the Forward Region at the LHC Thomas Ruf for the LHCb collaboration Introduction Particle
More informationECE 510 Lecture 7 Goodness of Fit, Maximum Likelihood. Scott Johnson Glenn Shirley
ECE 510 Lecture 7 Goodness of Fit, Maximum Likelihood Scott Johnson Glenn Shirley Confidence Limits 30 Jan 2013 ECE 510 S.C.Johnson, C.G.Shirley 2 Binomial Confidence Limits (Solution 6.2) UCL: Prob of
More informationSome Statistics. V. Lindberg. May 16, 2007
Some Statistics V. Lindberg May 16, 2007 1 Go here for full details An excellent reference written by physicists with sample programs available is Data Reduction and Error Analysis for the Physical Sciences,
More informationModern Methods of Data Analysis - WS 07/08
Modern Methods of Data Analysis Lecture III (29.10.07) Contents: Overview & Test of random number generators Random number distributions Monte Carlo The terminology Monte Carlo-methods originated around
More informationSIMULATION SEMINAR SERIES INPUT PROBABILITY DISTRIBUTIONS
SIMULATION SEMINAR SERIES INPUT PROBABILITY DISTRIBUTIONS Zeynep F. EREN DOGU PURPOSE & OVERVIEW Stochastic simulations involve random inputs, so produce random outputs too. The quality of the output is
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 informationPRECISION&MEASUREMENTS&
PRECISION&MEASUREMENTS& AT&Z&RESONANCE Z&Lineshape&and&number&of&neutrinos Lecture'2 Shahram&Rahatlou Fisica&delle&Par,celle&Elementari,&Anno&Accademico&2138214 http://www.roma1.infn.it/people/rahatlou/particelle/
More informationStatistical Tests: Discriminants
PHY310: Lecture 13 Statistical Tests: Discriminants Road Map Introduction to Statistical Tests (Take 2) Use in a Real Experiment Likelihood Ratio 1 Some Terminology for Statistical Tests Hypothesis This
More informationLecture 03 Positive Semidefinite (PSD) and Positive Definite (PD) Matrices and their Properties
Applied Optimization for Wireless, Machine Learning, Big Data Prof. Aditya K. Jagannatham Department of Electrical Engineering Indian Institute of Technology, Kanpur Lecture 03 Positive Semidefinite (PSD)
More informationb and c production in CMS
b and c production in CMS and ATLAS Francesco Fiori on behalf of CMS and ATLAS collaborations VII Meeting on B-physics (LAL) LHC luminosity evolution LHC startup in March 2010 9.67 pb -1 delivered lumi
More informationWhy is the field of statistics still an active one?
Why is the field of statistics still an active one? It s obvious that one needs statistics: to describe experimental data in a compact way, to compare datasets, to ask whether data are consistent with
More informationUNIT NUMBER PROBABILITY 6 (Statistics for the binomial distribution) A.J.Hobson
JUST THE MATHS UNIT NUMBER 19.6 PROBABILITY 6 (Statistics for the binomial distribution) by A.J.Hobson 19.6.1 Construction of histograms 19.6.2 Mean and standard deviation of a binomial distribution 19.6.3
More informationError propagation. Alexander Khanov. October 4, PHYS6260: Experimental Methods is HEP Oklahoma State University
Error propagation Alexander Khanov PHYS660: Experimental Methods is HEP Oklahoma State University October 4, 017 Why error propagation? In many cases we measure one thing and want to know something else
More informationKEKB collider and Belle detector
Belle B τν 1 KEKB collider and Belle detector Located(at(KEK(in(Tsukuba. e + e Υ(S) BB Resistive plate chamber + Fe (detection of μ and K L ) CsI(Tl) calorimeter (detection of e ± and γ) e! e! Four layers
More informationRadiative penguins at hadron machines. Kevin Stenson! University of Colorado!
Radiative penguins at hadron machines Kevin Stenson! University of Colorado! Introduction γ + b q u,c,t W s q b q u,c,t W + s q b q ν W W + u,c,t s q!! Radiative penguins involve b to s,d transitions with
More informationSearches for BSM Physics in Rare B-Decays in ATLAS
Searches for BSM Physics in Rare B-Decays in ATLAS Iskander Ibragimov on behalf of the ATLAS Collaboration SUSY 14 The nd International Conference on Supersymmetry and Unification of Fundamental Interactions
More informationLecture 2. Binomial and Poisson Probability Distributions
Durkin, Lecture 2, Page of 6 Lecture 2 Binomial and Poisson Probability Distributions ) Bernoulli Distribution or Binomial Distribution: Consider a situation where there are only two possible outcomes
More informationStatistics of Radioactive Decay
Statistics of Radioactive Decay Introduction The purpose of this experiment is to analyze a set of data that contains natural variability from sample to sample, but for which the probability distribution
More informationStandard Model at the LHC (Lecture 3: Measurement of Cross-Sections) M. Schott (CERN)
Standard Model at the LHC (Lecture 3: Measurement of Cross-Sections) M. Schott (CERN) Content 1 Measuring Cross-Sections 2 Hands-On 3 Background Estimations 4 Detector Calibration 5 Recent LHC Results
More informationDELPHI Collaboration DELPHI PHYS 656. Measurement of mass dierence between? and + and mass of -neutrino from three-prong -decays
DELPHI Collaboration DELPHI 96-167 PHYS 656 2 December, 1996 Measurement of mass dierence between? and + and mass of -neutrino from three-prong -decays M.Chapkin V.Obraztsov IHEP, Protvino Abstract The
More informationSimple Linear Regression for the Climate Data
Prediction Prediction Interval Temperature 0.2 0.0 0.2 0.4 0.6 0.8 320 340 360 380 CO 2 Simple Linear Regression for the Climate Data What do we do with the data? y i = Temperature of i th Year x i =CO
More informationThe Binomial distribution. Probability theory 2. Example. The Binomial distribution
Probability theory Tron Anders Moger September th 7 The Binomial distribution Bernoulli distribution: One experiment X i with two possible outcomes, probability of success P. If the experiment is repeated
More informationCLEO c. Anders Ryd Cornell University June 7, e e cc D D. D K,D K e
CLEOc Cornell University June 7, 25 e e cc D D D K,D K e K K e 1 Outline CLEOc experiment and the physics program Some early results D> Absolute hadronic branching fractions Semileptonic decays 2 Testing
More informationProspects for quarkonium studies at LHCb. «Quarkonium Production at the LHC» workshop
Prospects for quarkonium studies at LHCb «Quarkonium Production at the LHC» workshop Patrick Robbe, LAL Orsay, 19 Feb 2010 LHCb Detector Forward Spectrometer Geometry, with angular acceptance 15
More informationLikelihood Template Fitting
Momentum Scale Estimation Using Maximum Likelihood Template Fitting by Yu Zeng Department of Statistics Duke University Date: Approved: David L. Banks, Chair David B. Dunson Ashutosh V. Kotwal A thesis
More informationBayes Theorem (Recap) Adrian Bevan.
Bayes Theorem (Recap) Adrian Bevan email: a.j.bevan@qmul.ac.uk Overview Bayes theorem is given by The probability you want to compute: The probability of the hypothesis B given the data A. This is sometimes
More informationStandard Model of Particle Physics SS 2013
Lecture: Standard Model of Particle Physics Heidelberg SS 013 Weak Interactions II 1 Important Experiments Wu-Experiment (1957): radioactive decay of Co60 Goldhaber-Experiment (1958): radioactive decay
More informationStatistical Methods in Particle Physics
Statistical Methods in Particle Physics Lecture 10 December 17, 01 Silvia Masciocchi, GSI Darmstadt Winter Semester 01 / 13 Method of least squares The method of least squares is a standard approach to
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 informationJoint Probability Distributions, Correlations
Joint Probability Distributions, Correlations What we learned so far Events: Working with events as sets: union, intersection, etc. Some events are simple: Head vs Tails, Cancer vs Healthy Some are more
More informationPractical Statistics
Practical Statistics Lecture 1 (Nov. 9): - Correlation - Hypothesis Testing Lecture 2 (Nov. 16): - Error Estimation - Bayesian Analysis - Rejecting Outliers Lecture 3 (Nov. 18) - Monte Carlo Modeling -
More informationMeasurement And Uncertainty
Measurement And Uncertainty Based on Guidelines for Evaluating and Expressing the Uncertainty of NIST Measurement Results, NIST Technical Note 1297, 1994 Edition PHYS 407 1 Measurement approximates or
More information