The Gaisser–Hillas function is used in astroparticle physics. It parameterizes the longitudinal particle density in a cosmic ray air shower. The function was proposed in 1977 by Thomas K. Gaisser and Anthony Michael Hillas. The number of particles N ( X ) {\displaystyle N(X)} as a function of traversed atmospheric depth X {\displaystyle X} is expressed as
N ( X ) = N max ( X − X 0 X max − X 0 ) X max − X 0 λ exp ( X max − X λ ) , {\displaystyle N(X)=N_{\text{max}}\left({\frac {X-X_{0}}{X_{\text{max}}-X_{0}}}\right)^{\frac {X_{\text{max}}-X_{0}}{\lambda }}\exp \left({\frac {X_{\text{max}}-X}{\lambda }}\right),}
where N max {\displaystyle N_{\text{max}}} is maximum number of particles observed at depth X max {\displaystyle X_{\text{max}}} , and X 0 {\displaystyle X_{0}} and λ {\displaystyle \lambda } are primary mass and energy dependent parameters. Using substitutions
n = N N max {\displaystyle n={\frac {N}{N_{\text{max}}}}} , x = X − X 0 λ {\displaystyle x={\frac {X-X_{0}}{\lambda }}} and m = X max − X 0 λ {\displaystyle m={\frac {X_{\text{max}}-X_{0}}{\lambda }}}
the function can be written in an alternative one-parametric (m) form as
n ( x ) = ( x m ) m exp ( m − x ) = x m e − x m m e − m = exp [ m ( ln x − ln m ) − ( x − m ) ] . {\displaystyle n(x)=\left({\frac {x}{m}}\right)^{m}\exp(m-x)={\frac {x^{m}\,e^{-x}}{m^{m}\,e^{-m}}}=\exp \left[m\,(\ln x-\ln m)-(x-m)\right]\,.}
References
