Tuesday, April 4, 2017

statistics - Standard Deviation in Particle Physics


I'm familiar with sigma, and how its usually calculated and used, but would like to know how it's applied to particle physics. I recall reading that the discovery of the Higgs would only be credible unless it meets a 5 sigma criterion. What are the particle physics factors/numbers that would go into a sigma formula for a credible Higgs discovery?



Answer



Errors in particle physics are of two kinds. Statistical, and systematic. Statistical is the usual standard deviation of gausian distributions, sqrt(n)/N for 1 sigma. It is the systematics that take a lot of effort, and often are not taken well into account.


Systematic errors come from


1) the background to the signal expected. The background is calculated theoretically and entered into a monte carlo program that simulates the production of the data. There are also errors in the parameters the theory has used that enter as systematics.


2)Effects from the limitations of the measuring apparatus, which is also inserted into the monte carlo simulation of events



3) from the method of analysis, i.e. the cuts made in order to isolate the signal in the data and in the simulated events.


Ideally, the monte carlo simulation will have many more events than the data and so the statistical error from the MC can be ignored. One should estimate the errors from uncertainty in theory and from defects of detectors by varying the parameters in the MC program to the 1 sigma level of the important parameters and cuts and observing the change in the distributions.


Of course there are so many parameters and cuts that this procedure is not adhered to strictly, as happened with the new 3 sigma announcement of the "fifth force" by CDF, where the errors are just statistical from the events. That is why for discovery one asks 5 sigma. It is very hard to have an effect of 5 sigma due to the systematics enumerated above, whereas 3 sigma announcements have often disappeared . The ALEPH Higgs for example was a 3 sigma one that disappeared when the other 3 LEP experiments looked. I have known a 4 sigma resonance that was the effect of not estimating the systematics of the cuts. 5 sigma is for playing it safe.


I should add here that the different systematics in different experiments are the main reason why at least two expensive experimental set ups are approved and built in the collider setups: independent confirmation in parallel .


Edit: I want to correct this entry due to the discussions on statistical significance for the very much in the physics news OPERA neutrino speed measurement. In blog discussions I learned that the current way of dealing with systematic errors is of assuming randomness and adding them in quadrature. One interested in the details of how particle physics is treating errors should go to the CERN Yellow Report 99-03, July 1999, that has an html version.


No comments:

Post a Comment

classical mechanics - Moment of a force about a given axis (Torque) - Scalar or vectorial?

I am studying Statics and saw that: The moment of a force about a given axis (or Torque) is defined by the equation: $M_X = (\vec r \times \...