GPU Accelerated Reweighting Calculations for Neutrino Oscillation Analyses with the T2K Experiment
|
|
- Dale Simpson
- 5 years ago
- Views:
Transcription
1 GPU Accelerated Reweighting Calculations for Neutrino Oscillation Analyses with the T2K Experiment Oxford Many Core Seminar - 19th February Richard Calland rcalland@hep.ph.liv.ac.uk
2 Outline of Talk Neutrino Physics The T2K experiment Neutrino Oscillation analysis Statistical & Computing techniques How GPUs can help Feasibility of analysis Improve sensitivity of results Financially Conclusions 2
3 Neutrino Physics A branch of High Energy Physics (CERN etc) What are neutrinos? Why are they interesting? How do we study them? 3
4 What Neutrinos are the lightest particles we know of, but they are some of the most abundant They have tiny cross-sections Flooded by neutrinos from the sun, supernova explosions, the earths core and even produced in the upper atmosphere 100 billion solar neutrinos passing through your thumbnail every second! They pass straight through the earth! They have no charge (like a neutron) and only interact via the weak interaction (responsible for nuclear decays) There are 3 flavours (types), electron neutrino, muon neutrino and tau neutrino Analogues to the electron, muon and tau particles and are deeply related 4
5 5
6 Why Until recently, neutrinos were thought to have no mass (like light particles i.e. photons) The standard model of particle physics predicts massless neutrinos, and cannot explain the neutrino masses Don't let those CERN guys fool you into thinking the Higgs boson was the last piece of the puzzle Neutrinos are the ONLY particles we know that show signs of physics beyond the standard model 6
7 Neutrino Oscillations The reason we know neutrinos have mass is because we have observed a phenomenon called neutrino oscillation Simply put, neutrinos of one flavour (e,μ,τ) can transform into another Breaking a conservation law! Hint at physics beyond the standard model? This is explained in a quantum mechanical way 7
8 Neutrino Oscillations Neutrinos have 3 flavours (e,μ,τ) and 3 different masses m1,2,3 But they are more complicated than that Actually exist as a quantum superposition of all 3 Long story short, a neutrino created (e.g. in a radioactive decay inside the sun) as an electron neutrino could be observed later on as a muon neutrino 8
9 After ~400 km, non-zero probability for the electron neutrinos to be detected as muon neutrinos A macroscopic manifestation of a purely quantum mechanical phenomenon! 9
10 Oscillation Probability Accelerator Nuclear Reactor Solar 10
11 Neutrino Oscillations in Matter The interactions with ambient electrons in matter cause the neutrinos to feel an extra potential These so-called matter effects must be modelled Addition of an extra potential term to calculations creates a different oscillation probability compared to vacuum 2 neutrino approximation 11
12 Arbitrary state vector in flavour space Initial state transition amplitude Time evolution of the state Has 3 independent solutions for row vector ψj, assemble into matrix X Probability = A 2 12 Barger et al. Phys. Rev. D22(1980) 2718
13 How to study neutrino oscillations 13
14 The Tokai-to-Kamioka Experiment Measure the beam here Fire a beam of neutrinos from here 14
15 Over 500 members 11 Countries UK institutions one of the largest contributors (including Oxford) 15
16 How to measure neutrino oscillations At T2K we want to measure θ23 and θ13. We have sensitivity to both through νμ and νe samples at the far detector (SuperK) νμ = 20% νμ = 99% νe = 10% νe = 1% Can't see ντ = 0% νμ beam ντ = 70% 295 km Electron-like sample Electron neutrinos interact to produce an electron. Observe an excess of electrons compared to prediction muon-like sample Muon neutrinos interact to produce a muon. Observe a deficit of muons compared to prediction 16
17 17
18 Cerenkov Radiation 18
19 19
20 March 2011 Earthquake 20
21 Neutrino Oscillation Analysis 21
22 Basic Idea Simulate the neutrino beam, and how it looks in at Super-K, i.e. what kind of energy spectra it produces Build a probability distribution function (PDF) empirically from the detector Monte Carlo We then change the energy spectra to represent changing oscillation parameters Wobble these oscillation parameters (along with ~200 systematic uncertainty parameters) until a good fit has been found I.e. construct a likelihood function and minimize 22
23 Changing The MC to Simulate Neutrino Oscillation Apply a weight to each event which is equal to the probability of neutrino oscillation + matter effects + matter effects 23
24 Reweight method Samples of Monte Carlo events used to create probability distribution function Construct histogram, lose information but gain speed Treat every event separately, retain all information, lose speed Of course, I chose the latter 24
25 Binned Reweighting Fill once Reweight each bin using bin center 1,000,000 Monte Carlo Events Pro: Only have to do <Nbin> calculations per iteration No longer have to store all Monte Carlo in memory Con: All information for each event is lost 25
26 Event-by-event Reweighting 1,000,000 Monte Carlo Events For each event, calculate weight, fill histogram Pro: Retain all information about individual events Calculated weight for each event is more accurate Con: Have to store all Monte Carlo in memory Have to do 1,000,000 calculations per iteration 26
27 Flaws with binned reweighting Large binning to ensure adequate statistics per bin. If using binned reweighting, an event lying on the edge of a wide bin would not be represented fairly by the bin center value. We find differences ~1% 27
28 Flaws with binned reweighting (2) Monte Carlo events are broken down into sub categories If using the binned method, all the information about events is lost Can't model the effects of events migrating between different interaction modes / samples properly Proper modelling would give better sensitivity to final result! 28
29 3 Neutrino Oscillations with Matter Effects Luckily, someone already wrote a library to do this Unfortunately, rather slow ~ 2-3 seconds to reweight all 1,000,000 MC entries Need to do this ~20 million times for one analysis 2,000,000,000,000 calculations needed! (2 x 1012) ~2 years on a single machine with no GPU! Ported some functions to GPU using CUDA Vacuum calculations no speedup Full Matter Effects (constant density) speedup! 29
30 Benchmark Code compared with original Original test yielded x300, turns out wasn't compiled with -O2... Code also compared with quick OpenMP implementation In general, CUDA sees ~75-130x OpenMP sees ~2-4x speed improvement Double Precision Machine dependent 30
31 Benchmark Code Snippet clock.start(); CPU for (int i = 0; i < N; ++i) { bnu >SetMNS( nominal[0], nominal[2], nominal[1], nominal[3], nominal[4], nominal[5], 100.0, true ); bnu >propagatelinear( 2, 295, 2.6 ); sample_weights[i] = bnu >GetProb(2, 2); } clock.stop(); clock.start(); OpenMP #pragma omp parallel for num_threads(4) for (int i = 0; i < N; ++i) { bnu >SetMNS( nominal[0], nominal[2], nominal[1], nominal[3], nominal[4], nominal[5], 100.0, true ); bnu >propagatelinear( 2, 295, 2.6 ); sample_weights[i] = bnu >GetProb(2, 2); } clock.stop(); clock.start(); CUDA setmns(nominal[0], nominal[2], nominal[1], nominal[3], nominal[4], nominal[5], true); GetProb(2, 2, 295, 2.6, energy, N, sample_weights); clock.stop(); 31
32 extern "C" host void GetProb(int Alpha, int Beta, double Path, double Density, double *Energy, int n, double *oscw) { size_t dmsize = 3*3*sizeof(double); typedef double dmarray[3]; dmarray *d = (dmarray*)malloc(dmsize); memcpy(d, &dm, dmsize); dmarray *dm_device; Copy mixing matrix and mass cudamalloc((void **) &dm_device, dmsize); matrix (matter effects) to cudamemcpy(dm_device, dm, dmsize, cudamemcpyhosttodevice); size_t mixsize = 3*3*2*sizeof(double); typedef double mixarray[3][2]; mixarray *m = (mixarray*)malloc(mixsize); memcpy(m, &mix, mixsize); mixarray *mix_device; cudamalloc((void **) &mix_device,mixsize); cudamemcpy(mix_device, mix, mixsize, cudamemcpyhosttodevice); size_t size = n * sizeof(double); double *energy_device = NULL; cudamalloc((void **) &energy_device, size); cudamemcpy(energy_device, Energy, size, cudamemcpyhosttodevice); device Could copy to constant / texture memory Copy Monte Carlo event energies to device double *osc_weights; cudamalloc((void **) &osc_weights, size); dim3 block_size; block_size.x = 512; dim3 grid_size; grid_size.x = (n / block_size.x) + 1; Execute kernel propagatelinear<<<grid_size, block_size>>>(alpha, Beta, Path, Density, mix_device, dm_device, energy_device, osc_weights,n); } cudamemcpy(oscw, osc_weights, size, cudamemcpydevicetohost); clean_up(); // cudafree everything Copy oscillation weights32 back to host
33 Benchmark Strange features Xeon 2.67 GHz All threads occupied here Binned coarse reweighting around here Using all information with event33 by-event reweighting
34 Speedup Factor Xeon 2.67 GHz 34
35 Validation 10 million random comparisons between CPU and GPU calculations Agreement to precision, more than good enough for this application Difference attributed to extended ALU of CPU and different hardware implementations of non-associative calculations 35
36 Publication In the process of submission to Journal of Instrumentation (JINST) Preprint:
37 Non-linear Parameter Response For many parameters, the way they change is not gaussian Our models are too computing intensive to add into the analysis We create splines for each variation and use that instead Interaction mode C Interaction mode B Interaction mode A Reconstructed energy True energy 37
38 Splines A spline is a series of polynomials connected by points that have the same differential at the node Ensures a smooth changing function Want to sample each spline at various values 38
39 However As shown 2 pages back, there are a lot of splines (~4,000,000) A large chunk of memory (~Gbs) The next significant bottleneck of our code Off-loading work onto the GPU As we know, GPUs are compute- rather than data-intensive 39
40 Reformatting HEP experiments use a C++ library named ROOT Our splines are formatted as TSpline3 cubic spline objects Lots of bloat: 4M instantiations of a C++ class Try to reformat to perform better on GPU Instead of a class, format as an array and use a kernel function to evaluate 40
41 Convert TSpline3 Into an array 41
42 The Spline Monolith Convert TSpline3 objects into a monolithic array Eval(x) becomes Eval(i,x) i is the spline ID Already 2x faster on CPU! Smaller, takes up about 1.2 Gb RAM 42
43 GPU Implementation Copy all 1.2 Gb single array onto GPU RAM at initialization, keep it there Every iteration, evaluate all splines with a CUDA kernel and push the weight from each spline back onto CPU RAM Crude GPU implementation yields ~20x speed-up for the evaluation of 4,000,000 splines Issues with branching that are unavoidable Not using any shared memory which could be useful Would like to revisit this, but for now it helps 43
44 Bringing Everything Together Since neutrinos interact very rarely, neutrino experimental datasets are relatively small (~100 events) Need to make the most of what we have Event-by-event reweighting takes care of this Need a way to treat systematic uncertainties properly Must be able to handle large number of parameters Must be able to scale across a cluster 44
45 Markov Chain Monte Carlo To evaluate the posterior distribution, need to integrate over high-dimensions MCMC provides an efficient way to perform the ~200-dimensional integral MCMC performs a semi-random walk through parameter space, following the path of the likelihood function Can run multiple chains on a cluster and combine output 45
46 Method 1 Throw 200 proposed new values of parameters Evaluate likelihood function 2 Complete cycle ~20 million times 3 Is new state better than current? Accept or reject performed on GPU Write to disk Time per step ~0.3 seconds ~ 20x speed-up 46
47 Emerald Cluster MC samples, Splines etc MCMC chain MCMC chain job_script.py MCMC chain Combine output MCMC chain Data file (real or fake) MCMC chain MCMC chain Chains are considered independent samples and can simply be added together. 47
48 Final Result Best fit points and errors extracted from the posterior density of the combined markov chains A lot of information can be extracted from this high dimensional distribution (covariance, correlations, model comparisons etc) 48
49 Financial & Environmental Benefit Like to compare total dissipated power (TDP) and total costs for a full analysis M2090 = 225 W X5850 = 95 W kwh = 0.15 Assuming running at max TDP CPU version 95 W * 100 jobs * (11*20) hours = MW = 9.5 kwh GPU version 225 W + 95 W * 100 jobs * 11 hours = MW = 32 kwh ~ 6x less power / less expensive 49
50 Conclusions Wanted to write an analysis framework that makes full use of the information in our detector Monte Carlo Off-loaded bottlenecking algorithms to GPU to make this a possibility Combination of neutrino oscillation weight calculations and spline evaluations done on GPU ~20x faster than CPU only Full analysis takes 12 hours on Emerald Would take 10 days on an equivalent CPU cluster! Needed to produce 100 toy experiments for validation, again, without Emerald this would have taken 3 years on a CPU cluster! ~2 months with Emerald 50
51 The Future Currently writing my thesis High energy physics taking GPUs very seriously to solve big data problems they will face in the future CERN are investing in computing, they cannot handle data from the LHC upgrade presently! This analysis only off-loads a few components, with a rewrite, could make much better use of GPUs for reweighting Using shared memory CUDA + MPI 51
52 Thanks for listening! 52
Accelerated Neutrino Oscillation Probability Calculations and Reweighting on GPUs. Richard Calland
Accelerated Neutrino Oscillation Probability Calculations and Reweighting on GPUs Richard Calland University of Liverpool GPU Computing in High Energy Physics University of Pisa, 11th September 2014 Introduction
More informationLeïla Haegel University of Geneva
Picture found at: http://www.ps.uci.edu/~tomba/sk/tscan/pictures.html Leïla Haegel University of Geneva Seminar at the University of Sheffield April 27th, 2017 1 Neutrinos are: the only particle known
More informationHigh-performance processing and development with Madagascar. July 24, 2010 Madagascar development team
High-performance processing and development with Madagascar July 24, 2010 Madagascar development team Outline 1 HPC terminology and frameworks 2 Utilizing data parallelism 3 HPC development with Madagascar
More informationHow to find a Higgs boson. Jonathan Hays QMUL 12 th October 2012
How to find a Higgs boson Jonathan Hays QMUL 12 th October 2012 Outline Introducing the scalar boson Experimental overview Where and how to search Higgs properties Prospects and summary 12/10/2012 2 The
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 informationRecent Results from Alysia Marino, University of Colorado at Boulder Neutrino Flux Workshop, University of Pittsburgh, Dec 6 8,2012
Recent Results from Alysia Marino, University of Colorado at Boulder Neutrino Flux Workshop, University of Pittsburgh, Dec 6 8,2012 Outline T2K Beamline and Detectors Uncertainties from prior measurements
More informationDr. Andrea Bocci. Using GPUs to Accelerate Online Event Reconstruction. at the Large Hadron Collider. Applied Physicist
Using GPUs to Accelerate Online Event Reconstruction at the Large Hadron Collider Dr. Andrea Bocci Applied Physicist On behalf of the CMS Collaboration Discover CERN Inside the Large Hadron Collider at
More informationAtmospheric Neutrinos and Neutrino Oscillations
FEATURE Principal Investigator Takaaki Kajita Research Area Experimental Physics Atmospheric Neutrinos and Neutrino Oscillations Introduction About a hundred years ago Victor Hess aboard a balloon measured
More informationParticle Physics: Neutrinos part I
Particle Physics: Neutrinos part I José I. Crespo-Anadón Week 8: November 10, 2017 Columbia University Science Honors Program Course policies Attendance record counts Up to four absences Lateness or leaving
More informationHydra. A library for data analysis in massively parallel platforms. A. Augusto Alves Jr and Michael D. Sokoloff
Hydra A library for data analysis in massively parallel platforms A. Augusto Alves Jr and Michael D. Sokoloff University of Cincinnati aalvesju@cern.ch Presented at NVIDIA s GPU Technology Conference,
More informationNEW νe Appearance Results from the. T2K Experiment. Matthew Malek Imperial College London. University of Birmingham HEP Seminar 13 June 2012
NEW νe Appearance Results from the T2K Experiment Matthew Malek Imperial College London University of Birmingham HEP Seminar 13 June 2012 Outline Physics motivation: Neutrinos & Oscillations Overview of
More informationSystematic uncertainties in statistical data analysis for particle physics. DESY Seminar Hamburg, 31 March, 2009
Systematic uncertainties in statistical data analysis for particle physics DESY Seminar Hamburg, 31 March, 2009 Glen Cowan Physics Department Royal Holloway, University of London g.cowan@rhul.ac.uk www.pp.rhul.ac.uk/~cowan
More information27 : Distributed Monte Carlo Markov Chain. 1 Recap of MCMC and Naive Parallel Gibbs Sampling
10-708: Probabilistic Graphical Models 10-708, Spring 2014 27 : Distributed Monte Carlo Markov Chain Lecturer: Eric P. Xing Scribes: Pengtao Xie, Khoa Luu In this scribe, we are going to review the Parallel
More informationNeutrino Event Tagging Based On Nucleon Energy Spectra
Neutrino Event Tagging Based On Nucleon Energy Spectra Joshua Gevirtz Dr. Robert Svoboda UC Davis REU Program 2009 October 20, 2009 Abstract Since they were first theorized in 1930 by Wolfgang Pauli, much
More informationCS-206 Concurrency. Lecture 13. Wrap Up. Spring 2015 Prof. Babak Falsafi parsa.epfl.ch/courses/cs206/
CS-206 Concurrency Lecture 13 Wrap Up Spring 2015 Prof. Babak Falsafi parsa.epfl.ch/courses/cs206/ Created by Nooshin Mirzadeh, Georgios Psaropoulos and Babak Falsafi EPFL Copyright 2015 EPFL CS-206 Spring
More informationResults from T2K. Alex Finch Lancaster University ND280 T2K
Results from T2K Alex Finch Lancaster University ND280 T2K 1 Contents T2K Overview Tokai Kamioka Neutrino Oscillations Results Muon neutrino disappearance Electron neutrino appearance Charged Current inclusive
More informationIntro to Particle Physics and The Standard Model. Robert Clare UCR
Intro to Particle Physics and The Standard Model Robert Clare UCR Timeline of particle physics Ancient Greeks Rutherford 1911 Rutherford Chadwick Heisenberg 1930 s Hofstader Gell-Mann Ne eman 1960 s Timeline
More informationS XMP LIBRARY INTERNALS. Niall Emmart University of Massachusetts. Follow on to S6151 XMP: An NVIDIA CUDA Accelerated Big Integer Library
S6349 - XMP LIBRARY INTERNALS Niall Emmart University of Massachusetts Follow on to S6151 XMP: An NVIDIA CUDA Accelerated Big Integer Library High Performance Modular Exponentiation A^K mod P Where A,
More informationData Analysis II. PHY310: Lecture 18. Road Map
03/22/07 PHY310: Statistical Data Analysis 1 PHY310: Lecture 18 Data Analysis II Road Map Choosing the number of bins in a histogram Correcting Histograms for Systematic Uncertainty Weighting Events 03/22/07
More informationThe first one second of the early universe and physics beyond the Standard Model
The first one second of the early universe and physics beyond the Standard Model Koichi Hamaguchi (University of Tokyo) @ Colloquium at Yonsei University, November 9th, 2016. Credit: X-ray: NASA/CXC/CfA/M.Markevitch
More informationLecture 02. The Standard Model of Particle Physics. Part I The Particles
Lecture 02 The Standard Model of Particle Physics Part I The Particles The Standard Model Describes 3 of the 4 known fundamental forces Separates particles into categories Bosons (force carriers) Photon,
More informationRecent Results from T2K and Future Prospects
Recent Results from TK and Future Prospects Konosuke Iwamoto, on behalf of the TK Collaboration University of Rochester E-mail: kiwamoto@pas.rochester.edu The TK long-baseline neutrino oscillation experiment
More informationAnalyzing CMS events
Quarknet University of Rochester, March 23, 2012 Analyzing CMS events Questions in Particle Physics Introducing the Standard Model The Large Hadron Collider The CMS detector W and Z bosons: decays ispy
More informationStatistical Methods in Particle Physics Lecture 1: Bayesian methods
Statistical Methods in Particle Physics Lecture 1: Bayesian methods SUSSP65 St Andrews 16 29 August 2009 Glen Cowan Physics Department Royal Holloway, University of London g.cowan@rhul.ac.uk www.pp.rhul.ac.uk/~cowan
More informationPHYS 3446 Lecture #1. 6. Lab 7. Evaluation Policy. Syllabus Semester projects. Wednesday, Jan. 19, 2005 PHYS 3446, Spring 2005 Jae Yu
PHYS 3446 Lecture #1 Wednesday, Jan. 19, 2005 Dr. 1. Who am I? 2. Class time and location 3. Information and communication sources 4. Class specifications and style 5. Class plans Syllabus Semester projects
More information6-8 February 2017 Hotel do Mar Sesimbra. Hands on Neutrinos
6-8 February 2017 Hotel do Mar Sesimbra Hands on Neutrinos Hands on Neutrinos 1 I. BRIEF HISTORY OF NEUTRINOs The neutrinowas first postulated by Wolfgang Pauli in 1930 to explain how β particles emitted
More informationParticle + Physics at ATLAS and the Large Hadron Coillder
Particle + Physics at ATLAS and the Large Hadron Coillder Discovering the elementary particles of the Universe Kate Shaw The International Centre for Theoretical Physics + Overview Introduction to Particle
More informationNeutrinos and the Universe
Neutrinos and the Universe Susan Cartwright University of Sheffield Neutrinos and the Universe Discovering neutrinos Detecting neutrinos Neutrinos and the Sun Neutrinos and Supernovae Neutrinos and Dark
More informationPHYS 5326 Lecture #6. 1. Neutrino Oscillation Formalism 2. Neutrino Oscillation Measurements
PHYS 5326 Lecture #6 Wednesday, Feb. 14, 2007 Dr. 1. Neutrino Oscillation Formalism 2. Neutrino Oscillation Measurements 1. Solar Neutrinos 2. Atmospheric neutrinos 3. Accelerator Based Oscillation Experiments
More informationF. TASNÁDI LINKÖPING UNIVERSITY THEORETICAL PHYSICS NEUTRINO OSCILLATIONS & MASS
F. TASNÁDI LINKÖPING UNIVERSITY THEORETICAL PHYSICS NEUTRINO OSCILLATIONS & MASS the fundamental discoveries in physics con4nues 1 CONGRATULATIONS - NOBEL PRIZE IN PHYSICS 2016 the secrets of exotic matter
More informationParticles and Forces
Particles and Forces Particles Spin Before I get into the different types of particle there's a bit more back story you need. All particles can spin, like the earth on its axis, however it would be possible
More informationThe Discovery of the Higgs Boson: one step closer to understanding the beginning of the Universe
The Discovery of the Higgs Boson: one step closer to understanding the beginning of the Universe Anna Goussiou Department of Physics, UW & ATLAS Collaboration, CERN Kane Hall, University of Washington
More informationAn introduction to Bayesian reasoning in particle physics
An introduction to Bayesian reasoning in particle physics Graduiertenkolleg seminar, May 15th 2013 Overview Overview A concrete example Scientific reasoning Probability and the Bayesian interpretation
More informationNeutrinos. Invisible particles all around us. Amol Dighe. Department of Theoretical Physics Tata Institute of Fundamental Research, Mumbai
Neutrinos Invisible particles all around us Amol Dighe Department of Theoretical Physics Tata Institute of Fundamental Research, Mumbai DTP Presentation, VSRP 2013, TIFR, May 29, 2013 Why care about invisible
More informationIntroduction to the Standard Model
Introduction to the Standard Model Bill Murray, RAL, Quarks and leptons Bosons and forces The Higgs March 2002 1 Outline: An introduction to particle physics What is the Higgs Boson? Some unanswered questions
More informationNeutrino Physics: Lecture 1
Neutrino Physics: Lecture 1 Overview: discoveries, current status, future Amol Dighe Department of Theoretical Physics Tata Institute of Fundamental Research Feb 1, 2010 Plan of the course Omnipresent
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 informationStatistical Methods for Particle Physics (I)
Statistical Methods for Particle Physics (I) https://agenda.infn.it/conferencedisplay.py?confid=14407 Glen Cowan Physics Department Royal Holloway, University of London g.cowan@rhul.ac.uk www.pp.rhul.ac.uk/~cowan
More informationIntroduction to Particle Physics and the Standard Model. Robert Clare UCR
Introduction to Particle Physics and the Standard Model Robert Clare UCR Timeline of particle physics Ancient Greeks Rutherford 1911 Rutherford Chadwick Heisenberg 1930 s Hofstader Gell-Mann Ne eman 1960
More informationThe God particle at last? Astronomy Ireland, Oct 8 th, 2012
The God particle at last? Astronomy Ireland, Oct 8 th, 2012 Cormac O Raifeartaigh Waterford Institute of Technology CERN July 4 th 2012 (ATLAS and CMS ) A new particle of mass 125 GeV I The Higgs boson
More informationParallelization of the QC-lib Quantum Computer Simulator Library
Parallelization of the QC-lib Quantum Computer Simulator Library Ian Glendinning and Bernhard Ömer VCPC European Centre for Parallel Computing at Vienna Liechtensteinstraße 22, A-19 Vienna, Austria http://www.vcpc.univie.ac.at/qc/
More informationAn Introduction to Modern Particle Physics. Mark Thomson University of Cambridge
An Introduction to Modern Particle Physics Mark Thomson University of Cambridge Science Summer School: 30 th July - 1 st August 2007 1 Course Synopsis Introduction : Particles and Forces - what are the
More informationNeutrino oscillation physics potential of Hyper-Kamiokande
Neutrino oscillation physics potential of Hyper-Kamiokande on behalf of the Hyper-Kamiokande Collaboration Queen Mary University of London E-mail: l.cremonesi@qmul.ac.uk Hyper-Kamiokande (Hyper-K) is a
More informationA Quantum Chemistry Domain-Specific Language for Heterogeneous Clusters
A Quantum Chemistry Domain-Specific Language for Heterogeneous Clusters ANTONINO TUMEO, ORESTE VILLA Collaborators: Karol Kowalski, Sriram Krishnamoorthy, Wenjing Ma, Simone Secchi May 15, 2012 1 Outline!
More informationSP-CNN: A Scalable and Programmable CNN-based Accelerator. Dilan Manatunga Dr. Hyesoon Kim Dr. Saibal Mukhopadhyay
SP-CNN: A Scalable and Programmable CNN-based Accelerator Dilan Manatunga Dr. Hyesoon Kim Dr. Saibal Mukhopadhyay Motivation Power is a first-order design constraint, especially for embedded devices. Certain
More informationLecture 6: Markov Chain Monte Carlo
Lecture 6: Markov Chain Monte Carlo D. Jason Koskinen koskinen@nbi.ku.dk Photo by Howard Jackman University of Copenhagen Advanced Methods in Applied Statistics Feb - Apr 2016 Niels Bohr Institute 2 Outline
More informationSPARSE SOLVERS POISSON EQUATION. Margreet Nool. November 9, 2015 FOR THE. CWI, Multiscale Dynamics
SPARSE SOLVERS FOR THE POISSON EQUATION Margreet Nool CWI, Multiscale Dynamics November 9, 2015 OUTLINE OF THIS TALK 1 FISHPACK, LAPACK, PARDISO 2 SYSTEM OVERVIEW OF CARTESIUS 3 POISSON EQUATION 4 SOLVERS
More informationPhoton Coupling with Matter, u R
1 / 16 Photon Coupling with Matter, u R Consider the up quark. We know that the u R has electric charge 2 3 e (where e is the proton charge), and that the photon A is a linear combination of the B and
More informationStudy Sheet for Modern Physics
Study Sheet for Modern Physics Classical mechanics was meant to provide the general rules that govern the dynamics of all material bodies, such as cannon balls, planets, and pendulums, and is defined as
More informationEarly SUSY searches at the LHC
on behalf of the ATLAS and CMS Collaborations Imperial College London E-mail: a.tapper@imperial.ac.uk Supersymmetry may give rise to striking events that could be discovered early in LHC running. Search
More informationExam #3. Final Exam. Exam 3 review. How do we measure properties of a star? A detailed outline of study topics is here:
Exam #3 Exam #3 is Thursday 4/9 in this room You can bring page of notes (front and back) Bring your calculator and a # pencil Exam 3 covers material from 4/1 onward (only 8 lectures) Consequently, no
More informationRyan Stillwell Paper: /10/2014. Neutrino Astronomy. A hidden universe. Prepared by: Ryan Stillwell. Tutor: Patrick Bowman
Neutrino Astronomy A hidden universe Prepared by: Ryan Stillwell Tutor: Patrick Bowman Paper: 124.129 Date: 10 October 2014 i Table of Contents 1. Introduction pg 1 1.1 Background pg 1 2. Findings & Discussion
More informationNeutrino Oscillations and the Matter Effect
Master of Science Examination Neutrino Oscillations and the Matter Effect RAJARSHI DAS Committee Walter Toki, Robert Wilson, Carmen Menoni Overview Introduction to Neutrinos Two Generation Mixing and Oscillation
More informationNeutrinos: Yesterday, Today and Tomorrow. Stanley Wojcicki SLAC Summer Institute 2010 August 13, 2010
Neutrinos: Yesterday, Today and Tomorrow August 13, 2010 1 My Marching Orders 2 My Marching Orders...the summary talk should be visionary, rather than a dedicated summary of the SSI program. 2 My Marching
More informationObservation of a New Particle with a Mass of 125 GeV
Observation of a New Particle with a Mass of 125 GeV CMS Experiment, CERN 4 July 2012 Summary In a joint seminar today at CERN and the ICHEP 2012 conference[1] in Melbourne, researchers of the Compact
More informationSearch for heavy neutrinos in kaon decays
Search for heavy neutrinos in kaon decays L. Littenberg (work mainly done by A.T.Shaikhiev INR RAS) HQL-2016 Outline Motivation Previous heavy neutrino searches Experiment BNL-E949 Selection criteria Efficiency
More informationA Search for Point Sources of High Energy Neutrinos with AMANDA-B10
A Search for Point Sources of High Energy Neutrinos with AMANDA-B10 Scott Young, for the AMANDA collaboration UC-Irvine PhD Thesis: http://area51.berkeley.edu/manuscripts Goals! Perform an all-sky search
More informationTowards a unified treatment of systematic errors
Towards a unified treatment of systematic errors Anselmo Cervera Villanueva IFIC (Valencia) Andrea Donini IFT (Madrid) & IFIC (Valencia) EURO-ν general meeting Strasbourg, 2 June 2010 Ingreedients 2 To
More informationSpecial Contribution Observation of Neutrinos at Super-Kamiokande Observatory
Special Contribution Observation of Neutrinos at Super-Kamiokande Observatory Yoshinari Hayato Associate Professor Institute for Cosmic Ray Research The University of Tokyo 1. Introduction Neutrinos are
More informationStatistical Methods for Particle Physics Lecture 1: parameter estimation, statistical tests
Statistical Methods for Particle Physics Lecture 1: parameter estimation, statistical tests http://benasque.org/2018tae/cgi-bin/talks/allprint.pl TAE 2018 Benasque, Spain 3-15 Sept 2018 Glen Cowan Physics
More informationCharacterisation of Silicon Photomultipliers for the T2K Experiment
Characterisation of Silicon Photomultipliers for the T2K Experiment, 18th May 2010 Martin Haigh, University of Oxford Outline Brief introduction to the T2K experiment. Overall configuration and goals.
More informationDiscovery of the Neutrino Mass-I. P1X* Frontiers of Physics Lectures October 2004 Dr Paul Soler University of Glasgow
-I P1X* Frontiers of Physics Lectures 19-0 October 004 Dr Paul Soler University of Glasgow Outline 1. Introduction: the structure of matter. Neutrinos:.1 Neutrino interactions. Neutrino discovery and questions.3
More informationNew Results from the MINOS Experiment
University College London E-mail: annah@hep.ucl.ac.uk The MINOS experiment is a long-baseline neutrino experiment designed to study neutrino behaviour, in particular the phenomenon of neutrino oscillations.
More informationThe Physics of Neutrinos
AAAS Feb. 15, 2004 The Physics of Neutrinos Boris Kayser Neutrinos are Abundant We, and all the everyday objects around us here on earth, are made of electrons, protons, and neutrons. In the universe,
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 informationPHYSICS 107. Lecture 27 What s Next?
PHYSICS 107 Lecture 27 What s Next? The origin of the elements Apart from the expansion of the universe and the cosmic microwave background radiation, the Big Bang theory makes another important set of
More informationAdvanced statistical methods for data analysis Lecture 1
Advanced statistical methods for data analysis Lecture 1 RHUL Physics www.pp.rhul.ac.uk/~cowan Universität Mainz Klausurtagung des GK Eichtheorien exp. Tests... Bullay/Mosel 15 17 September, 2008 1 Outline
More informationThe 64th Compton Lecture Series Unsolved Mysteries of the Universe: Looking for Clues in Surprising Places
The 64th Compton Lecture Series Unsolved Mysteries of the Universe: Looking for Clues in Surprising Places Brian Odom Fall 2006 http://kicp.uchicago.edu/~odom/compton.htm Lecture 2: From the Big Bang to
More informationImplementing NNLO into MCFM
Implementing NNLO into MCFM Downloadable from mcfm.fnal.gov A Multi-Threaded Version of MCFM, J.M. Campbell, R.K. Ellis, W. Giele, 2015 Higgs boson production in association with a jet at NNLO using jettiness
More informationLong Baseline Neutrino Experiments
Physics in Collision, 27/6/2008 Mark Thomson 1 Long Baseline Neutrino Experiments Mark Thomson Cavendish Laboratory University of Cambridge Introduction Neutrino Beams Past Experiments: K2K Current Experiments:
More informationThe Compact Muon Solenoid Experiment. Conference Report. Mailing address: CMS CERN, CH-1211 GENEVA 23, Switzerland
Available on CMS information server CMS CR -2017/094 The Compact Muon Solenoid Experiment Conference Report Mailing address: CMS CERN, CH211 GENEVA 23, Switzerland 18 March 2017 (v3, 20 March 2017) Recent
More informationParticles in the Early Universe
Particles in the Early Universe David Morrissey Saturday Morning Physics, October 16, 2010 Using Little Stuff to Explain Big Stuff David Morrissey Saturday Morning Physics, October 16, 2010 Can we explain
More informationThe Particle World. This talk: What is our Universe made of? Where does it come from? Why does it behave the way it does?
The Particle World What is our Universe made of? Where does it come from? Why does it behave the way it does? Particle physics tries to answer these questions. This talk: particles as we understand them
More informationNew Results for ν µ ν e oscillations in MINOS
New Results for ν µ ν e oscillations in MINOS Jelena Ilic Rutherford Appleton Lab 4/28/10 RAL PPD Seminar 1 Neutrino Mixing Mass eigenstates flavour eigenstates Maki-Nakagawa-Sakata: Flavour composition
More informationEco517 Fall 2013 C. Sims MCMC. October 8, 2013
Eco517 Fall 2013 C. Sims MCMC October 8, 2013 c 2013 by Christopher A. Sims. This document may be reproduced for educational and research purposes, so long as the copies contain this notice and are retained
More informationDid something ν just break the speed limit?
Did something ν just break the speed limit? Amol Dighe Department of Theoretical Physics Tata Institute of Fundamental Research Chai-and-Why, Prithvi Theatre, Nov 6, 2011 Did something ν break the speed
More informationBeyond Standard Model Effects in Flavour Physics: p.1
Beyond Standard Model Effects in Flavour Physics: Alakabha Datta University of Mississippi Feb 13, 2006 Beyond Standard Model Effects in Flavour Physics: p.1 OUTLINE Standard Model (SM) and its Problems.
More informationCosmology and particle physics
Cosmology and particle physics Lecture notes Timm Wrase Lecture 5 The thermal universe - part I In the last lecture we have shown that our very early universe was in a very hot and dense state. During
More informationHow large is the "LSND anomaly"?
How large is the "LSND anomaly"? Andrey Elagin on behalf of the HARP-CDP group HEP lunch, UChicago, March 12, 2012 Outline 3.8 σ LSND Final Paper PRD 64, 112007 3.8 σ 2.9 σ HARP-CDP Paper I arxiv:1110.4265
More informationUNVEILING THE ULTIMATE LAWS OF NATURE: DARK MATTER, SUPERSYMMETRY, AND THE LHC. Gordon Kane, Michigan Center for Theoretical Physics Warsaw, June 2009
UNVEILING THE ULTIMATE LAWS OF NATURE: DARK MATTER, SUPERSYMMETRY, AND THE LHC Gordon Kane, Michigan Center for Theoretical Physics Warsaw, June 2009 OUTLINE! Some things we ve learned about the physical
More informationSearch for a SM-like Higgs boson. Zijun Xu Peking University ISHP, Beijing Aug 15, 2013
l n W H W Search for a SM-like Higgs boson decaying into WW lνq q in CMS Zijun Xu Peking University ISHP, Beijing Aug 15, 2013 Outline Introduction H WW Lep + MET + 2 jets [1] : 160-600 GeV H WW Lep +
More informationTHE FLORIDA STATE UNIVERSITY COLLEGE OF ARTS & SCIENCES W BOSON PRODUCTION CHARGE ASYMMETRY IN THE ELECTRON CHANNEL ASHLEY S HUFF
THE FLORIDA STATE UNIVERSITY COLLEGE OF ARTS & SCIENCES W BOSON PRODUCTION CHARGE ASYMMETRY IN THE ELECTRON CHANNEL By ASHLEY S HUFF A Thesis submitted to the Department of Physics In partial fulfillment
More informationIntroduction to Parallel Programming in OpenMP Dr. Yogish Sabharwal Department of Computer Science & Engineering Indian Institute of Technology, Delhi
Introduction to Parallel Programming in OpenMP Dr. Yogish Sabharwal Department of Computer Science & Engineering Indian Institute of Technology, Delhi Lecture - 33 Parallel LU Factorization So, now, how
More informationChapter 32 Lecture Notes
Chapter 32 Lecture Notes Physics 2424 - Strauss Formulas: mc 2 hc/2πd 1. INTRODUCTION What are the most fundamental particles and what are the most fundamental forces that make up the universe? For a brick
More informationThe future of neutrino physics (at accelerators)
Mauro Mezzetto, Istituto Nazionale Fisica Nucleare, Padova The future of neutrino physics (at accelerators) Present Status Concepts, strategies, challenges The two players: Dune and Hyper-Kamiokande Conclusions
More informationT2K and the MAINZ Photon Beam
TK and the MAINZ Photon Beam February, S. Boyd University of4 Warwick 11 Precis Long baseline neutrino oscillation experiment TK UK Collaboration is building an ECAL for detection We are thinking of using
More informationAnnouncements. Problem Set 6 due next Monday, February 25, at 12:50PM. Midterm graded, will be returned at end of lecture.
Turing Machines Hello Hello Condensed Slide Slide Readers! Readers! This This lecture lecture is is almost almost entirely entirely animations that that show show how how each each Turing Turing machine
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 informationCMS. S.Paktinat. School of Particles and Accelerators Institute for Studies in Theoretical Physics and Mathematics (IPM) May 23, 2013
IPM @ CMS S.Paktinat School of Particles and Accelerators Institute for Studies in Theoretical Physics and Mathematics (IPM) May 23, 2013 The CMS detector for LHC Each color shows a different layer This
More informationarxiv: v1 [hep-ex] 11 May 2017
LATEST RESULTS FROM TK arxiv:1705.0477v1 [hep-ex] 11 May 017 Marcela Batkiewicz a, for the TK collaboration Institute of Nuclear Physics Polish Academy of Sciences, Cracow, Poland Abstract. The TK (Tokai
More informationHydra. A. Augusto Alves Jr and M.D. Sokoloff. University of Cincinnati
Hydra A library for data analysis in massively parallel platforms A. Augusto Alves Jr and M.D. Sokoloff University of Cincinnati aalvesju@cern.ch Presented at the Workshop Perspectives of GPU computing
More informationSunrise: Patrik Jonsson. Panchromatic SED Models of Simulated Galaxies. Lecture 2: Working with Sunrise. Harvard-Smithsonian Center for Astrophysics
Sunrise: Panchromatic SED Models of Simulated Galaxies Lecture 2: Working with Sunrise Patrik Jonsson Harvard-Smithsonian Center for Astrophysics Lecture outline Lecture 1: Why Sunrise? What does it do?
More informationSTA 4273H: Sta-s-cal Machine Learning
STA 4273H: Sta-s-cal Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 2 In our
More informationAn Improved Energy Flow Diagram is Shown for an HCE8S Universe
1 An Improved Energy Flow Diagram is Shown for an HCE8S Universe George R. Briggs Abstract: An alternate forward-time, reverse-time energy cycle of the 9 th cycle of an HCE8S universe for a full loop of
More informationThe God particle at last? Science Week, Nov 15 th, 2012
The God particle at last? Science Week, Nov 15 th, 2012 Cormac O Raifeartaigh Waterford Institute of Technology CERN July 4 th 2012 (ATLAS and CMS ) A new particle of mass 125 GeV Why is the Higgs particle
More informationS Subdivide, Preprocess and Conquer: Micromagnetism FEM/BEM-Simulations on Single-Node/Multi-GPU Systems
S4283 - Subdivide, : Micromagnetism FEM/BEM-Simulations on Single-Node/Multi-GPU Systems Elmar Westphal - Forschungszentrum Jülich GmbH 1 Contents Micromagnetism TetraMag, a FEM/BEM Micromagnetism Simulator
More informationFigure 1 Overview of the T2K experiment
Figure 1 Overview of the T2K experiment The proton beams extracted from the main ring synchrotron at J-PARC are directed in a westward direction through the T2K primary beam line. The beams strike a target
More informationScalable and Power-Efficient Data Mining Kernels
Scalable and Power-Efficient Data Mining Kernels Alok Choudhary, John G. Searle Professor Dept. of Electrical Engineering and Computer Science and Professor, Kellogg School of Management Director of the
More informationFirst, a look at using OpenACC on WRF subroutine advance_w dynamics routine
First, a look at using OpenACC on WRF subroutine advance_w dynamics routine Second, an estimate of WRF multi-node performance on Cray XK6 with GPU accelerators Based on performance of WRF kernels, what
More informationAnalysis of Z ee with the ATLAS-Detector
DESY Summer Student Programme 2008 Hamburg Analysis of Z ee with the ATLAS-Detector Maximilian Schlupp ATLAS Group 11.09.2008 maximilian.schlupp@desy.de maximilian.schlupp@tu-dortmund.de Supervisor: Karsten
More information