Calculating error bars for neutrino mixing parameters
  One goal of contemporary particle physics is to determine the mixing angles and mass-squared differences that constitute the phenomenological constants that describe neutrino oscillations. Of great interest are not only the best fit values of these constants but also their errors. Some of the neutrino oscillation data is statistically poor and cannot be treated by normal (Gaussian) statistics. To extract confidence intervals when the statistics are not normal, one should not utilize the value for chisquare versus confidence level taken from normal statistics. Instead, we propose that one should use the normalized likelihood function as a probability distribution; the relationship between the correct chisquare and a given confidence level can be computed by integrating over the likelihood function. This allows for a definition of confidence level independent of the functional form of the !2 function; it is particularly useful for cases in which the minimum of the !2 function is near a boundary. We present two pedagogic examples and find that the proposed method yields confidence intervals that can differ significantly from those obtained by using the value of chisquare from normal statistics. For example, we find that for the first data release of the T2K experiment the probability that chisquare is not zero, as defined by the maximum confidence level at which the value of zero is not allowed, is 92%. Using the value of chisquare at zero and assigning a confidence level from normal statistics, a common practice, gives the over estimation of 99.5%. 
H. R. Burroughs, B. K. Cogswell, J. Escamilla-Roa, D. C. Latimer, and   D. J. Ernst
nucl-th hep-ph