Crystal Ball function

From HandWiki
Examples of the Crystal Ball function.

The Crystal Ball function, named after the Crystal Ball Collaboration (hence the capitalized initial letters), is a probability density function (PDF) commonly used to model various lossy processes in high-energy physics such as Bremsstrahlung by electrons. It consists of a Gaussian core portion and a power-law low-end tail, below a certain threshold. The function itself and its first derivative are both continuous.

The Crystal Ball function is given by:

f(x;α,n,x¯,σ)=N⋅{exp⁡(−(x−x¯)22σ2),for x−x¯σ>−αA⋅(B−x−x¯σ)−n,for x−x¯σ⩽−α,

where

A=(n|α|)n⋅exp⁡(−|α|22),
B=n|α|−|α|,
N=1σ(C+D),
C=n|α|⋅1n−1⋅exp⁡(−|α|22),
D=π2(1+erf⁡(|α|2)),

with the error function erf.

The parameters of the function (that are usually determined by a fit) are:

  • N is a normalization factor (Skwarnicki 1986)
  • α>0 defines the point where the PDF changes from a power-law to a Gaussian distribution
  • n>1 is the power of the power-law tail
  • x¯ and σ are the mean and the standard deviation of the Gaussian