FiveLayerAtmosphereDensityProfile

class ctapipe.atmosphere.FiveLayerAtmosphereDensityProfile(table: astropy.table.table.Table)

Bases: ctapipe.atmosphere.AtmosphereDensityProfile

CORSIKA 5-layer atmosphere model

Layers 1-4 are modeled with:

\[T(h) = a_i + b_i \exp{-h/c_i}\]

Layer 5 is modeled with:

..math:: T(h) = a_5 - b_5 frac{h}{c_5}

References

[corsika-user] D. Heck and T. Pierog, “Extensive Air Shower Simulation with CORSIKA:

A User’s Guide”, 2021, Appendix F

Methods Summary

__call__(height)	Returns
from_array(array[, metadata])	construct from a 5x5 array as provided by eventio
integral(height)	Integral of the profile along the height axis, i.e. the atmospheric depth \(X\).

Methods Documentation

__call__(height) → astropy.units.quantity.Quantity

Returns

u.Quantity[“g cm-3”]

the density at height h

classmethod from_array(array: numpy.ndarray, metadata: Optional[Dict] = None)

construct from a 5x5 array as provided by eventio

integral(height) → astropy.units.quantity.Quantity

Integral of the profile along the height axis, i.e. the atmospheric depth \(X\).

\[X(h) = \int_{h}^{\infty} \rho(h') dh'\]

Returns

u.Quantity[“g/cm2”]:

Integral of the density from height h to infinity