Kinematic Fitting. Mark M. Ito. July 11, Jefferson Lab. Mark Ito (JLab) Kinematic Fitting July 11, / 25
|
|
- Sherilyn Hampton
- 5 years ago
- Views:
Transcription
1 Kinematic Fitting Mark M. Ito Jefferson Lab July 11, 2013 Mark Ito (JLab) Kinematic Fitting July 11, / 25
2 Introduction Motivation Particles measured independently Often: known relationships between their kinematics energy and momentum conservation daughters of common parent particle of known mass particles from a common vertex Measured kinematics do not respect relationships How to incoporate knowledge of relationships? Why? better measurements hypothesis testing Mark Ito (JLab) Kinematic Fitting July 11, / 25
3 PrimEx Experiment in Hall B PrimEx Detector Mark Ito (JLab) Kinematic Fitting July 11, / 25
4 PrimEx Experiment in Hall B Hybrid Calorimeter Mark Ito (JLab) Kinematic Fitting July 11, / 25
5 PrimEx Experiment in Hall B Event Display Mark Ito (JLab) Kinematic Fitting July 11, / 25
6 PrimEx Experiment in Hall B Two-photon Mass Distribution Mark Ito (JLab) Kinematic Fitting July 11, / 25
7 Example Example: Decay of Neutral Pion into Two Photons π 0 γγ both photons detected assume photon directions are known precisely energies have relative uncertainty σ E /E = 5%/ E for simplicity: look at 500 MeV/c π 0 moving in z-direction but, will not assume 500 MeV/c nor z-direction for pion Mark Ito (JLab) Kinematic Fitting July 11, / 25
8 Example two-photon invariant mass Mark Ito (JLab) Kinematic Fitting July 11, / 25
9 Example relative error on single photon energy Mark Ito (JLab) Kinematic Fitting July 11, / 25
10 Example two-photon measured momentum Mark Ito (JLab) Kinematic Fitting July 11, / 25
11 Example What is the problem? want to improve resolution assume that photons came from π 0, so mass is known can we use this information? could adjust one photon (why just one?) could scale them both (high energy γ: better E measurement) could minimize ( ) χ 2 E1,fit E 2 ( 1,meas E2,fit E 2,meas = + σ 1 but minimum is clear: E fit = E meas (something is missing!) must introduce constraint: (k 1 + k 2 ) 2 = m 2 π gives 2E 1 E 2 (1 cos θ) = m 2 π Problem: minimize χ 2 while simultaneously satisfying constraint Mark Ito (JLab) Kinematic Fitting July 11, / 25 σ 2 ) 2
12 Example Minimization Strategy multi-variable minimization with constraints: Lagrange multipliers instead of minimizing over two variables, minimize over three E 1,fit, E 2,fit, and λ ( ) χ 2 E1,fit E 2 ( 1,meas E2,fit E 2,meas = + σ 1 σ 2 +2λ [ 2E 1,fit E 2,fit (1 cos θ) mπ 2 ] ) 2 Mark Ito (JLab) Kinematic Fitting July 11, / 25
13 Results fit two-photon invariant mass Mark Ito (JLab) Kinematic Fitting July 11, / 25
14 Results fit relative error on single photon energy Mark Ito (JLab) Kinematic Fitting July 11, / 25
15 Results fit two-photon measured momentum Mark Ito (JLab) Kinematic Fitting July 11, / 25
16 General Formalism Variable Definitions measured variables N number of measured variables (input) y vector of measurements, N-dimensional (input) V covariance matrix, N N (input) η vector of fit values of measured variables, N-dimensional χ 2 = (y η) T V 1 (y η) unmeasured variables J number of unmeasured variables (input) ξ vector of unmeasured variable values, J-dimensional Mark Ito (JLab) Kinematic Fitting July 11, / 25
17 General Formalism Variable Definitions 2 constraints K number of constraints (input) f vector of constraint functions, K dimensional (input) Each constraint a function of measured and unmeasured variables. When constraint is satisfied f k (η 1,..., η N, ξ 1,..., ξ J ) = 0 for k = 1,..., K Lagrange multipliers λ vector of mulipliers, K-dimensional Extended χ 2 to be minimized: χ 2 (η, ξ, λ) = (y η) T V 1 (y η) + 2λ T f (η, ξ) Mark Ito (JLab) Kinematic Fitting July 11, / 25
18 General Formalism Minimization Condition Set all partial derivatives to zero: where χ 2 [ ] = 2V 1 (y η) + 2Fη T λ η = 0, n n χ 2 ] = [2Fξ T ξ λ = 0, j j χ 2 λ k = [2f ] k = 0, j = 1,..., J k = 1,..., K (F η ) kn = f k η n and (F ξ ) kj = f k ξ j n = 1,..., N All equations manifestly linear, except for the multiplier derivatives. Mark Ito (JLab) Kinematic Fitting July 11, / 25
19 General Formalism Linearize the Constraints; Iterate the Solutions Go to a linear form of the constraints: 1 Change variables: η = η 0 + δη, ξ = ξ 0 + δξ, and λ = λ 0 + δλ 2 Take first order in Taylor expansion: f (δη, δξ, δλ) = f (η 0, ξ 0, λ 0 ) + (F η )(δη) + (F ξ )(δξ) Now have a linear system of N + J + K equations with N + J + K unknowns, i. e. δη, δξ, δλ But an approximation was made: need to iterate Also: must initialize variables Mark Ito (JLab) Kinematic Fitting July 11, / 25
20 General Formalism Errors covariance matrices: fit variables V η = V [ ] I (G HUH T )V unmeasured variables V ξ = U covariances between fit and unmeasured variables V η,ξ = VHU where G = Fη T S 1 F η, H = Fη T S 1 F ξ, U 1 = Fξ T S 1 F ξ and S = Fη T VF η Mark Ito (JLab) Kinematic Fitting July 11, / 25
21 Diagnostics Stretch Functions or Pulls How to tell if the thing is working? look at use these N quantities: z n = y n η n σ 2 (y n ) σ 2 (η n ) n = 1,..., N Gaussian with mean at 0, σ of 1 If not there are problems: offset mean: measurements biased wrong width: errors not correct tails: non-gaussian tails in measurements, background in sample Mark Ito (JLab) Kinematic Fitting July 11, / 25
22 Diagnostics Unmeasured Variables, Number of Constraints go back to π 0 decay could have introduced unmeasured variables: p π,x, p π,y, p π,z but then would have to apply 3-momentum conservation now have 4 constraints with 3 measured variables used to have 1 constraint with 0 measured variables same problem recast: 1-C fit C = K J Mark Ito (JLab) Kinematic Fitting July 11, / 25
23 Diagnostics Check χ 2 χ 2 should have a standard probablility density distribution: f (χ 2 ) Convenient to check χ 2 probability: P(χ 2 0) = P runs from 0 to 1 χ 2 0 f (χ 2 )dχ 2 for nominal χ 2 distribution P will be uniform non-uniformity: problem with errors, check the pulls often see peaks near 0: bad χ 2, background in sample Mark Ito (JLab) Kinematic Fitting July 11, / 25
24 Diagnostics Example Probability Distribution Mark Ito (JLab) Kinematic Fitting July 11, / 25
25 Summary Summary measured variables, with or without statistical correlation, may have physical relationships kinematic fit varies values of measured quantities to satisfy relationships minimize χ 2 with constraints improved measurements diagnostics about bias and errors in measurements are generated goodness of fit a handle on correctness of physical relationships assumed Reference: A. G. Frodesen, O. Skjeggestad, H. Tøfte. Probability and Statistics in Particle Physics. Universitetsforlaget, ISBN Mark Ito (JLab) Kinematic Fitting July 11, / 25
Studies of the KS and KL lifetimes and
Studies of the KS and KL lifetimes and BR(K ) with KLOE ± ± + Simona S. Bocchetta* on behalf of the KLOE Collaboration KAON09 Tsukuba June 9th 2009 * INFN and University of Roma Tre Outline DA NE and KLOE
More informationDeeply Virtual Compton Scattering at Jefferson Lab
Deeply Virtual Compton Scattering at Jefferson Lab June 16-17, 2016 Frederic Georges (PhD student) PhD Supervisor: Carlos Muñoz Camacho Institut de Physique Nucléaire d Orsay CNRS-IN2P3 Université Paris-Sud,
More informationMIP Reconstruction Techniques and Minimum Spanning Tree Clustering
SLAC-PUB-11359 July 25 MIP Reconstruction Techniques and Minimum Spanning Tree Clustering Wolfgang F. Mader The University of Iowa, 23 Van Allen Hall, 52242 Iowa City, IA The development of a tracking
More informationStudies of e γx and of e π0x with clas and clas12
Studies of e γx and of e π0x with clas and clas12 DeepPWG meeting, JLab, June 15, 2017 Hard exclusive photon production BH propagators Possibilities with e γx Extracting PBPT from double spin asymmetry
More informationPerformance of the GlueX Detector Systems
Performance of the GlueX Detector Systems GlueX-doc-2775 Gluex Collaboration August 215 Abstract This document summarizes the status of calibration and performance of the GlueX detector as of summer 215.
More informationPhysics 736. Experimental Methods in Nuclear-, Particle-, and Astrophysics. - Statistical Methods -
Physics 736 Experimental Methods in Nuclear-, Particle-, and Astrophysics - Statistical Methods - Karsten Heeger heeger@wisc.edu Course Schedule and Reading course website http://neutrino.physics.wisc.edu/teaching/phys736/
More informationTracking and Vertex reconstruction at LHCb for Run II
Tracking and Vertex reconstruction at LHCb for Run II Hang Yin Central China Normal University On behalf of LHCb Collaboration The fifth Annual Conference on Large Hadron Collider Physics, Shanghai, China
More information05/09/07 CHEP2007 Stefano Spataro. Simulation and Event Reconstruction inside the PandaRoot Framework. Stefano Spataro. for the collaboration
for the collaboration Overview Introduction on Panda Structure of the framework Event generation Detector implementation Reconstruction The Panda experiment AntiProton Annihilations at Darmstadt Multi
More informationCharged Particle Reconstruction in HIC Detectors
Charged Particle Reconstruction in HIC Detectors Ralf-Arno Tripolt, Qiyan Li [http://de.wikipedia.org/wiki/marienburg_(mosel)] H-QM Lecture Week on Introduction to Heavy Ion Physics Kloster Marienburg/Mosel,
More informationTheoretical Concepts of Machine Learning
Theoretical Concepts of Machine Learning Part 2 Institute of Bioinformatics Johannes Kepler University, Linz, Austria Outline 1 Introduction 2 Generalization Error 3 Maximum Likelihood 4 Noise Models 5
More informationStarting a Data Analysis
03/20/07 PHY310: Statistical Data Analysis 1 PHY310: Lecture 17 Starting a Data Analysis Road Map Your Analysis Log Exploring the Data Reading the input file (and making sure it's right) Taking a first
More informationTOOLS FOR DATA ANALYSIS INVOLVING
TOOLS FOR DATA ANALYSIS INVOLVING µ-vertex DETECTORS KalmanFitter package : Primary vertex fit Secondary vertex fit Decay chain TMVA package : Multivariate analysis 1 J. Bouchet Kent State University cτ
More informationTracking and Vertexing in 3D B-field
Tracking and Vertexing in 3D B-field Norman Graf (SLAC) HPS Collaboration Meeting, JLab October 26, 2015 Track Extrapolation At the heart of both track and vertex fitting in the presence of a non-uniform
More informationMixture Models and the EM Algorithm
Mixture Models and the EM Algorithm Padhraic Smyth, Department of Computer Science University of California, Irvine c 2017 1 Finite Mixture Models Say we have a data set D = {x 1,..., x N } where x i is
More informationStudy of t Resolution Function
Belle-note 383 Study of t Resolution Function Takeo Higuchi and Hiroyasu Tajima Department of Physics, University of Tokyo (January 6, 200) Abstract t resolution function is studied in detail. It is used
More informationFirst LHCb measurement with data from the LHC Run 2
IL NUOVO CIMENTO 40 C (2017) 35 DOI 10.1393/ncc/i2017-17035-4 Colloquia: IFAE 2016 First LHCb measurement with data from the LHC Run 2 L. Anderlini( 1 )ands. Amerio( 2 ) ( 1 ) INFN, Sezione di Firenze
More informationEvent reconstruction in STAR
Chapter 4 Event reconstruction in STAR 4.1 Data aquisition and trigger The STAR data aquisition system (DAQ) [54] receives the input from multiple detectors at different readout rates. The typical recorded
More informationOPERA: A First ντ Appearance Candidate
OPERA: A First ντ Appearance Candidate Björn Wonsak On behalf of the OPERA collaboration. 1 Overview The OPERA Experiment. ντ Candidate Background & Sensitivity Outlook & Conclusions 2/42 Overview The
More informationMarkov chain Monte Carlo methods
Markov chain Monte Carlo methods (supplementary material) see also the applet http://www.lbreyer.com/classic.html February 9 6 Independent Hastings Metropolis Sampler Outline Independent Hastings Metropolis
More informationPoS(Baldin ISHEPP XXII)134
Implementation of the cellular automaton method for track reconstruction in the inner tracking system of MPD at NICA, G.A. Ososkov and A.I. Zinchenko Joint Institute of Nuclear Research, 141980 Dubna,
More informationCalorimeter Object Status. A&S week, Feb M. Chefdeville, LAPP, Annecy
Calorimeter Object Status St A&S week, Feb. 1 2017 M. Chefdeville, LAPP, Annecy Outline Status of Pi0 calibration Calibration survey with electrons (E/p) Calorimeter performance with single photons from
More informationPerformance Testing and Tuning of Kalman Track-Fitting for CLEO III
Performance Testing and Tuning of Kalman Track-Fitting for CLEO III Daniela Silva Department of Mathematics, Wayne State University, Detroit, MI, 48202 Abstract CLEO III will use a Kalman track fitter.
More informationDijet A LL Dijet cross section 2006 Dijet A LL Preparation for Tai Sakuma MIT
Dijet A LL 2005 Dijet cross section 2006 Dijet A LL Preparation for 2008 Tai Sakuma MIT Motivation for dijet A LL is to constrain ΔG G with initial event kinematics The initial state variables can be written
More informationA new inclusive secondary vertex algorithm for b-jet tagging in ATLAS
A new inclusive secondary vertex algorithm for b-jet tagging in ATLAS Giacinto Piacquadio, Christian Weiser Albert-Ludwig University of Freiburg, Hermann-Herder-Str. 3, D-79104 Freiburg E-mail: giacinto.piacquadio@physik.uni-freiburg.de
More informationπ ± Charge Exchange Cross Section on Liquid Argon
π ± Charge Exchange Cross Section on Liquid Argon Kevin Nelson REU Program, College of William and Mary Mike Kordosky College of William and Mary, Physics Dept. August 5, 2016 Abstract The observation
More informationInference and Representation
Inference and Representation Rachel Hodos New York University Lecture 5, October 6, 2015 Rachel Hodos Lecture 5: Inference and Representation Today: Learning with hidden variables Outline: Unsupervised
More informationAnalysis of Σ 0 baryon, or other particles, or detector outputs from the grid data at ALICE
Analysis of Σ 0 baryon, or other particles, or detector outputs from the grid data at ALICE Introduction Analysis Chain Current status of Σ 0 analysis Examples of root files from the data and MC Discussion
More informationGAMs semi-parametric GLMs. Simon Wood Mathematical Sciences, University of Bath, U.K.
GAMs semi-parametric GLMs Simon Wood Mathematical Sciences, University of Bath, U.K. Generalized linear models, GLM 1. A GLM models a univariate response, y i as g{e(y i )} = X i β where y i Exponential
More informationElectron and Photon Reconstruction and Identification with the ATLAS Detector
Electron and Photon Reconstruction and Identification with the ATLAS Detector IPRD10 S12 Calorimetry 7th-10th June 2010 Siena, Italy Marine Kuna (CPPM/IN2P3 Univ. de la Méditerranée) on behalf of the ATLAS
More informationThe BABAR Vertexing. Vertexing and Composition Tools Group 1. BABAR Analysis Document #102, Version 7.00 February 7, 2005
BABAR Analysis Document #102, Version 7.00 February 7, 2005 The BABAR Vertexing Vertexing and Composition Tools Group 1 Abstract This document describes the vertexing algorithms in BABAR, including the
More informationThursday, February 16, More C++ and root
More C++ and root Today s Lecture Series of topics from C++ Example of thinking though a problem Useful physics classes in ROOT Functions by this point you are used to the syntax of a function pass by
More informationAlignment of the ATLAS Inner Detector
Alignment of the ATLAS Inner Detector Heather M. Gray [1,2] on behalf of the ATLAS ID Alignment Group [1] California Institute of Technology [2] Columbia University The ATLAS Experiment tile calorimeter
More informationCMS Conference Report
Available on CMS information server CMS CR 2005/021 CMS Conference Report 29 Septemebr 2005 Track and Vertex Reconstruction with the CMS Detector at LHC S. Cucciarelli CERN, Geneva, Switzerland Abstract
More informationTime of CDF (II)
TOF detector lecture, 19. august 4 1 Time of Flight @ CDF (II) reconstruction/simulation group J. Beringer, A. Deisher, Ch. Doerr, M. Jones, E. Lipeles,, M. Shapiro, R. Snider, D. Usynin calibration group
More information3D-Triplet Tracking for LHC and Future High Rate Experiments
3D-Triplet Tracking for LHC and Future High Rate Experiments André Schöning Physikalisches Institut, Universität Heidelberg Workshop on Intelligent Trackers WIT 2014 University of Pennsylvania May 14-16,
More informationSupport Vector Machines (a brief introduction) Adrian Bevan.
Support Vector Machines (a brief introduction) Adrian Bevan email: a.j.bevan@qmul.ac.uk Outline! Overview:! Introduce the problem and review the various aspects that underpin the SVM concept.! Hard margin
More informationADAPTIVE METROPOLIS-HASTINGS SAMPLING, OR MONTE CARLO KERNEL ESTIMATION
ADAPTIVE METROPOLIS-HASTINGS SAMPLING, OR MONTE CARLO KERNEL ESTIMATION CHRISTOPHER A. SIMS Abstract. A new algorithm for sampling from an arbitrary pdf. 1. Introduction Consider the standard problem of
More informationQuick Start with CASSY Lab. Bi-05-05
Quick Start with CASSY Lab Bi-05-05 About this manual This manual helps you getting started with the CASSY system. The manual does provide you the information you need to start quickly a simple CASSY experiment
More informationSegmentation: Clustering, Graph Cut and EM
Segmentation: Clustering, Graph Cut and EM Ying Wu Electrical Engineering and Computer Science Northwestern University, Evanston, IL 60208 yingwu@northwestern.edu http://www.eecs.northwestern.edu/~yingwu
More informationAlignment of the CMS silicon tracker using Millepede II
Journal of Physics: Conference Series Alignment of the CMS silicon tracker using Millepede II To cite this article: Peter Schleper et al 2008 J. Phys.: Conf. Ser. 119 032040 Related content - CMS silicon
More informationSavitsky-Golay Filters
Savitsky-Golay Filters Robert DeSerio March 24, 2008 References 1. Abraham Savitsky and Marcel J. E. Golay, Smoothing and differentiation of data by simplified least squares procedures, Anal. Chem. 36,
More informationData Mining Chapter 3: Visualizing and Exploring Data Fall 2011 Ming Li Department of Computer Science and Technology Nanjing University
Data Mining Chapter 3: Visualizing and Exploring Data Fall 2011 Ming Li Department of Computer Science and Technology Nanjing University Exploratory data analysis tasks Examine the data, in search of structures
More informationarxiv:hep-ex/ v1 24 Jun 1994
MULTIPLE SCATTERING ERROR PROPAGATION IN PARTICLE TRACK RECONSTRUCTION M. Penţia, G. Iorgovan INSTITUTE OF ATOMIC PHYSICS RO-76900 P.O.Box MG-6, Bucharest, ROMANIA e-mail: pentia@roifa.bitnet arxiv:hep-ex/9406006v1
More informationNote Set 4: Finite Mixture Models and the EM Algorithm
Note Set 4: Finite Mixture Models and the EM Algorithm Padhraic Smyth, Department of Computer Science University of California, Irvine Finite Mixture Models A finite mixture model with K components, for
More informationPhysics 736. Experimental Methods in Nuclear-, Particle-, and Astrophysics. - Statistics and Error Analysis -
Physics 736 Experimental Methods in Nuclear-, Particle-, and Astrophysics - Statistics and Error Analysis - Karsten Heeger heeger@wisc.edu Feldman&Cousin what are the issues they deal with? what processes
More informationTHIS paper describes vertex reconstruction in the AT-
Vertex Reconstruction in the Experiment at the LHC E. Bouhova-Thacker, T. Koffas, V. Kostyukhin, W. Liebig, M. Limper, G. Piacquadio, K. Prokofiev, C. Weiser, A. Wildauer, on behalf of the Collaboration.
More informationTrack reconstruction for the Mu3e experiment based on a novel Multiple Scattering fit Alexandr Kozlinskiy (Mainz, KPH) for the Mu3e collaboration
Track reconstruction for the Mu3e experiment based on a novel Multiple Scattering fit Alexandr Kozlinskiy (Mainz, KPH) for the Mu3e collaboration CTD/WIT 2017 @ LAL-Orsay Mu3e Experiment Mu3e Experiment:
More informationPoS(EPS-HEP2017)492. Performance and recent developments of the real-time track reconstruction and alignment of the LHCb detector.
Performance and recent developments of the real-time track reconstruction and alignment of the LHCb detector. CERN E-mail: agnieszka.dziurda@cern.ch he LHCb detector is a single-arm forward spectrometer
More informationPhysics 736. Experimental Methods in Nuclear-, Particle-, and Astrophysics. Lecture 14
Physics 736 Experimental Methods in Nuclear-, Particle-, and Astrophysics Lecture 14 Karsten Heeger heeger@wisc.edu Course Schedule and Reading course website http://neutrino.physics.wisc.edu/teaching/phys736/
More informationMonte Carlo Tuning of LUARLW Mode
Monte Carlo Tuning of LUARLW Mode Lei Xia 1, Guangshun Huang 1, Ronggang Ping 2, Xi an Xiong 2 1 University of Science and Technology of China 2 Institute of High Energy Physics, Chinese Academy of Sciences
More informationMissing Data Analysis for the Employee Dataset
Missing Data Analysis for the Employee Dataset 67% of the observations have missing values! Modeling Setup For our analysis goals we would like to do: Y X N (X, 2 I) and then interpret the coefficients
More informationTPC tracking and particle identification in high-density environment
TPC tracking and particle identification in high-density environment Y.Belikov, M.Ivanov, K.Safarik CERN, Switzerland J.Bracinik Comenius University, Bratislava Track finding and fitting algorithm in the
More informationFeature Detectors and Descriptors: Corners, Lines, etc.
Feature Detectors and Descriptors: Corners, Lines, etc. Edges vs. Corners Edges = maxima in intensity gradient Edges vs. Corners Corners = lots of variation in direction of gradient in a small neighborhood
More informationfor the DESY/ ECFA study detector
The TPC Tracker for the DESY/ ECFA study detector Ties Behnke DESY 1-May-1999 the TPC tracker requirements from physics a TPC at TESLA: can this work? results from simulation technical issues conclusion
More informationA Novel Strip Energy Splitting Algorithm for the Fine Granular Readout of a Scintillator Strip Electromagnetic Calorimeter
1 3 A Novel Strip Energy Splitting Algorithm for the Fine Granular Readout of a Scintillator Strip Electromagnetic Calorimeter 4 Katsushige Kotera a, Daniel Jeans b, Akiya Miyamoto c, and Tohru Takeshita
More informationA Novel Strip Energy Splitting Algorithm for the Fine Granular Readout of a Scintillator Strip Electromagnetic Calorimeter
1 3 A Novel Strip Energy Splitting Algorithm for the Fine Granular Readout of a Scintillator Strip Electromagnetic Calorimeter 4 Katsushige Kotera a, Daniel Jeans b, Akiya Miyamoto c, and Tohru Takeshita
More informationKristoffer H. Rose. Version /01/07
Summary of qsymbols Kristoffer H. Rose Version 1.12 1997/01/07 Abstract qsymbols is a L A TEX [1] package defining systematic mnemonic abbreviations, starting with a single open quote for symbols, and
More informationUNIVERSITY OF OSLO. Please make sure that your copy of the problem set is complete before you attempt to answer anything. k=1
UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Eam in: STK4051 Computational statistics Day of eamination: Thursday November 30 2017. Eamination hours: 09.00 13.00. This problem set consists
More informationExpectation-Maximization Methods in Population Analysis. Robert J. Bauer, Ph.D. ICON plc.
Expectation-Maximization Methods in Population Analysis Robert J. Bauer, Ph.D. ICON plc. 1 Objective The objective of this tutorial is to briefly describe the statistical basis of Expectation-Maximization
More informationRadial Basis Function Networks: Algorithms
Radial Basis Function Networks: Algorithms Neural Computation : Lecture 14 John A. Bullinaria, 2015 1. The RBF Mapping 2. The RBF Network Architecture 3. Computational Power of RBF Networks 4. Training
More informationFirst results from the LHCb Vertex Locator
First results from the LHCb Vertex Locator Act 1: LHCb Intro. Act 2: Velo Design Dec. 2009 Act 3: Initial Performance Chris Parkes for LHCb VELO group Vienna Conference 2010 2 Introducing LHCb LHCb is
More informationDetector Alignment with Tracks. Wouter Hulsbergen (Nikhef, BFYS)
Detector Alignment with Tracks Wouter Hulsbergen (Nikhef, BFYS) Detector alignment LHC silicon detectors provide
More informationEstimation of Item Response Models
Estimation of Item Response Models Lecture #5 ICPSR Item Response Theory Workshop Lecture #5: 1of 39 The Big Picture of Estimation ESTIMATOR = Maximum Likelihood; Mplus Any questions? answers Lecture #5:
More informationShort-Cut MCMC: An Alternative to Adaptation
Short-Cut MCMC: An Alternative to Adaptation Radford M. Neal Dept. of Statistics and Dept. of Computer Science University of Toronto http://www.cs.utoronto.ca/ radford/ Third Workshop on Monte Carlo Methods,
More informationMODEL SELECTION AND REGULARIZATION PARAMETER CHOICE
MODEL SELECTION AND REGULARIZATION PARAMETER CHOICE REGULARIZATION METHODS FOR HIGH DIMENSIONAL LEARNING Francesca Odone and Lorenzo Rosasco odone@disi.unige.it - lrosasco@mit.edu June 6, 2011 ABOUT THIS
More informationOptimal designs for comparing curves
Optimal designs for comparing curves Holger Dette, Ruhr-Universität Bochum Maria Konstantinou, Ruhr-Universität Bochum Kirsten Schorning, Ruhr-Universität Bochum FP7 HEALTH 2013-602552 Outline 1 Motivation
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationCloth Hair. and. soft bodies
Cloth Hair Lesson 11 and soft bodies Lesson 08 Outline Problem definition and motivations Modeling deformable solids with mass-spring model Position based dynamics Modeling cloths with mass-spring model
More informationRobot Mapping. Least Squares Approach to SLAM. Cyrill Stachniss
Robot Mapping Least Squares Approach to SLAM Cyrill Stachniss 1 Three Main SLAM Paradigms Kalman filter Particle filter Graphbased least squares approach to SLAM 2 Least Squares in General Approach for
More informationGraphbased. Kalman filter. Particle filter. Three Main SLAM Paradigms. Robot Mapping. Least Squares Approach to SLAM. Least Squares in General
Robot Mapping Three Main SLAM Paradigms Least Squares Approach to SLAM Kalman filter Particle filter Graphbased Cyrill Stachniss least squares approach to SLAM 1 2 Least Squares in General! Approach for
More informationEUDET Telescope Geometry and Resolution Studies
EUDET EUDET Telescope Geometry and Resolution Studies A.F.Żarnecki, P.Nieżurawski February 2, 2007 Abstract Construction of EUDET pixel telescope will significantly improve the test beam infrastructure
More informationAnalogue, Digital and Semi-Digital Energy Reconstruction in the CALICE AHCAL
Analogue, Digital and Semi-Digital Energy Reconstruction in the AHCAL Deutsches Elektronen Synchrotron (DESY), Hamburg, Germany E-mail: coralie.neubueser@desy.de Within the collaboration different calorimeter
More informationUpdate on Energy Resolution of
Update on Energy Resolution of the EMC Using µµγ Sample David Hopkins Royal Holloway, University of London EMC Reconstruction Workshop, December 5 th, 2004 Outline Study of photon energy resolution Compare
More informationAnimation Lecture 10 Slide Fall 2003
Animation Lecture 10 Slide 1 6.837 Fall 2003 Conventional Animation Draw each frame of the animation great control tedious Reduce burden with cel animation layer keyframe inbetween cel panoramas (Disney
More information2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123
2.7 Cloth Animation 320491: Advanced Graphics - Chapter 2 123 Example: Cloth draping Image Michael Kass 320491: Advanced Graphics - Chapter 2 124 Cloth using mass-spring model Network of masses and springs
More informationMonte Carlo Tuning of LUARLW Mode
Monte Carlo Tuning of LUARLW Mode Lei Xia 1, Guangshun Huang 1, Ronggang Ping 2, Xi an Xiong 2 1 University of Science and Technology of China 2 Institute of High Energy Physics, Chinese Academy of Sciences
More informationMixture Models and EM
Table of Content Chapter 9 Mixture Models and EM -means Clustering Gaussian Mixture Models (GMM) Expectation Maximiation (EM) for Mixture Parameter Estimation Introduction Mixture models allows Complex
More informationChap.12 Kernel methods [Book, Chap.7]
Chap.12 Kernel methods [Book, Chap.7] Neural network methods became popular in the mid to late 1980s, but by the mid to late 1990s, kernel methods have also become popular in machine learning. The first
More informationAn Introduction to Markov Chain Monte Carlo
An Introduction to Markov Chain Monte Carlo Markov Chain Monte Carlo (MCMC) refers to a suite of processes for simulating a posterior distribution based on a random (ie. monte carlo) process. In other
More informationCMS Alignement and Calibration workflows: lesson learned and future plans
Available online at www.sciencedirect.com Nuclear and Particle Physics Proceedings 273 275 (2016) 923 928 www.elsevier.com/locate/nppp CMS Alignement and Calibration workflows: lesson learned and future
More informationFunction approximation using RBF network. 10 basis functions and 25 data points.
1 Function approximation using RBF network F (x j ) = m 1 w i ϕ( x j t i ) i=1 j = 1... N, m 1 = 10, N = 25 10 basis functions and 25 data points. Basis function centers are plotted with circles and data
More informationSolution Sketches Midterm Exam COSC 6342 Machine Learning March 20, 2013
Your Name: Your student id: Solution Sketches Midterm Exam COSC 6342 Machine Learning March 20, 2013 Problem 1 [5+?]: Hypothesis Classes Problem 2 [8]: Losses and Risks Problem 3 [11]: Model Generation
More informationDirect photon measurements in ALICE. Alexis Mas for the ALICE collaboration
Direct photon measurements in ALICE Alexis Mas for the ALICE collaboration 1 Outline I - Physics motivations for direct photon measurements II Direct photon measurements in ALICE i - Conversion method
More informationBoltzmann Statistics REDOR (BS-REDOR) Data Analysis Guide
Boltzmann Statistics REDOR (BS-REDOR) Data Analysis Guide BS-REDOR manual version 1.0, pertaining to BS-REDOR version 1.0 12 July 2007 Melbourne, Australia BS-REDOR has been published: John D. Gehman,
More informationImproving the way neural networks learn Srikumar Ramalingam School of Computing University of Utah
Improving the way neural networks learn Srikumar Ramalingam School of Computing University of Utah Reference Most of the slides are taken from the third chapter of the online book by Michael Nielson: neuralnetworksanddeeplearning.com
More informationPrimary Vertex Reconstruction at LHCb
LHCb-PUB-214-44 October 21, 214 Primary Vertex Reconstruction at LHCb M. Kucharczyk 1,2, P. Morawski 3, M. Witek 1. 1 Henryk Niewodniczanski Institute of Nuclear Physics PAN, Krakow, Poland 2 Sezione INFN
More informationModelling of non-gaussian tails of multiple Coulomb scattering in track fitting with a Gaussian-sum filter
Modelling of non-gaussian tails of multiple Coulomb scattering in track fitting with a Gaussian-sum filter A. Strandlie and J. Wroldsen Gjøvik University College, Norway Outline Introduction A Gaussian-sum
More informationEicRoot for tracking R&D studies
EicRoot for tracking R&D studies Alexander Kiselev EIC Software Meeting Jefferson Lab September,24 2015 Contents of the talk Tracking code implementation in EicRoot Few particular applications: Basic forward
More informationWhat the Dual-Dirac Model is and What it is Not
What the Dual-Dirac Model is and What it is Not Ransom Stephens October, 006 Abstract: The dual-dirac model is a simple tool for estimating total jitter defined at a bit error ratio, TJ(BER), for serial
More informationProgramming assignment 3 Mean-shift
Programming assignment 3 Mean-shift 1 Basic Implementation The Mean Shift algorithm clusters a d-dimensional data set by associating each point to a peak of the data set s probability density function.
More informationPoS(ACAT)049. Alignment of the ATLAS Inner Detector. Roland Haertel Max-Planck-Institut für Physik, Munich, Germany
Max-Planck-Institut für Physik, Munich, Germany E-mail: haertel@mppmu.mpg.de The ATLAS experiment at the LHC is currently under construction at CERN and will start operation in summer 2008. The Inner Detector
More informationSetting Up a C++ Analyzer Replay
Setting Up a C++ Analyzer Replay Ole Hansen Jefferson Lab Joint Hall A & Hall C Data Analysis Workshop January 14, 2015 Ole Hansen (Jefferson Lab) C++ Analyzer Replay Tutorial January 14, 2015 1 / 14 Brief
More informationModel validation T , , Heli Hiisilä
Model validation T-61.6040, 03.10.2006, Heli Hiisilä Testing Neural Models: How to Use Re-Sampling Techniques? A. Lendasse & Fast bootstrap methodology for model selection, A. Lendasse, G. Simon, V. Wertz,
More informationINLA: Integrated Nested Laplace Approximations
INLA: Integrated Nested Laplace Approximations John Paige Statistics Department University of Washington October 10, 2017 1 The problem Markov Chain Monte Carlo (MCMC) takes too long in many settings.
More informationMeasurement of Skin Dose
Measurement of Skin Dose Sources of Uncertainty Kenneth A. Fetterly, Ph.D. William Pavlicek, Ph.D. Dan Bednarek, PhD 2014 AAPM Annual Meeting, Austin Texas 2013 MFMER slide-1 Purpose 1. Present a framework
More informationEdge and local feature detection - 2. Importance of edge detection in computer vision
Edge and local feature detection Gradient based edge detection Edge detection by function fitting Second derivative edge detectors Edge linking and the construction of the chain graph Edge and local feature
More informationAlignment of the CMS Silicon Tracker
Alignment of the CMS Silicon Tracker Tapio Lampén 1 on behalf of the CMS collaboration 1 Helsinki Institute of Physics, Helsinki, Finland Tapio.Lampen @ cern.ch 16.5.2013 ACAT2013 Beijing, China page 1
More informationAlignment of the CMS silicon tracker
Journal of Physics: Conference Series Alignment of the CMS silicon tracker To cite this article: Gero Flucke and the CMS Collaboration 22 J. Phys.: Conf. Ser. 368 236 View the article online for updates
More informationIn Homework 1, you determined the inverse dynamics model of the spinbot robot to be
Robot Learning Winter Semester 22/3, Homework 2 Prof. Dr. J. Peters, M.Eng. O. Kroemer, M. Sc. H. van Hoof Due date: Wed 6 Jan. 23 Note: Please fill in the solution on this sheet but add sheets for the
More informationPassive Differential Matched-field Depth Estimation of Moving Acoustic Sources
Lincoln Laboratory ASAP-2001 Workshop Passive Differential Matched-field Depth Estimation of Moving Acoustic Sources Shawn Kraut and Jeffrey Krolik Duke University Department of Electrical and Computer
More information