ARGUS distribution

From HandWiki
ARGUS
Probability density function

c = 1.
Cumulative distribution function

c = 1.
Parameters c>0 cut-off (real)
χ>0 curvature (real)
Support x(0,c)
PDF see text
CDF see text
Mean μ=cπ/8χeχ24I1(χ24)Ψ(χ)

where I1 is the Modified Bessel function of the first kind of order 1, and Ψ(x) is given in the text.
Mode c2χ(χ22)+χ4+4
Variance c2(13χ2+χϕ(χ)Ψ(χ))μ2

In physics, the ARGUS distribution, named after the particle physics experiment ARGUS,[1] is the probability distribution of the reconstructed invariant mass of a decayed particle candidate in continuum background[clarification needed].

Definition

The probability density function (pdf) of the ARGUS distribution is:

f(x;χ,c)=χ32πΨ(χ)xc21x2c2exp{12χ2(1x2c2)},

for 0x<c. Here χ and c are parameters of the distribution and

Ψ(χ)=Φ(χ)χϕ(χ)12,

where Φ(x) and ϕ(x) are the cumulative distribution and probability density functions of the standard normal distribution, respectively.

Cumulative distribution function

The cumulative distribution function (cdf) of the ARGUS distribution is

F(x)=1Ψ(χ1x2/c2)Ψ(χ).

Parameter estimation

Parameter c is assumed to be known (the kinematic limit of the invariant mass distribution), whereas χ can be estimated from the sample X1, …, Xn using the maximum likelihood approach. The estimator is a function of sample second moment, and is given as a solution to the non-linear equation

13χ2+χϕ(χ)Ψ(χ)=1ni=1nxi2c2.

The solution exists and is unique, provided that the right-hand side is greater than 0.4; the resulting estimator χ^ is consistent and asymptotically normal.

Generalized ARGUS distribution

Sometimes a more general form is used to describe a more peaking-like distribution:

f(x)=2pχ2(p+1)Γ(p+1)Γ(p+1,12χ2)xc2(1x2c2)pexp{12χ2(1x2c2)},0xc,c>0,χ>0,p>1
F(x)=Γ(p+1,12χ2(1x2c2))Γ(p+1,12χ2)Γ(p+1)Γ(p+1,12χ2),0xc,c>0,χ>0,p>1

where Γ(·) is the gamma function, and Γ(·,·) is the upper incomplete gamma function.

Here parameters c, χ, p represent the cutoff, curvature, and power respectively.

The mode is:

c2χ(χ22p1)+χ2(χ24p+2)+(1+2p)2

The mean is:

μ=cpπΓ(p)Γ(52+p)χ2p+22p+2M(p+1,52+p,χ22)Γ(p+1)Γ(p+1,12χ2)

where M(·,·,·) is the Kummer's confluent hypergeometric function.[2][circular reference]

The variance is:

σ2=c2(χ2)p+1χp+3eχ22+(χ22(p+1)){Γ(p+2)Γ(p+2,12χ2)}χ2(p+1)(Γ(p+1)Γ(p+1,12χ2))μ2

p = 0.5 gives a regular ARGUS, listed above.

References

  1. Albrecht, H. (1990). "Search for hadronic b→u decays". Physics Letters B 241 (2): 278–282. doi:10.1016/0370-2693(90)91293-K. Bibcode1990PhLB..241..278A.  (More formally by the ARGUS Collaboration, H. Albrecht et al.) In this paper, the function has been defined with parameter c representing the beam energy and parameter p set to 0.5. The normalization and the parameter χ have been obtained from data.
  2. Confluent hypergeometric function

Further reading