pyrex.askaryan.AskaryanSignal.had_shower_profile

static AskaryanSignal.had_shower_profile(z, energy, density=0.92, crit_energy=0.17006, rad_length=39.562, int_length=113.03, scale_factor=0.11842)

Calculates the hadronic shower longitudinal profile.

The longitudinal charge is calculated for a given distance, density, critical energy, hadron radiation length, and interaction length in ice, plus a scale factor for the number of particles.

Parameters

zarray_like

Distance (m) along the shower at which to calculate the charge.

energyfloat

Energy (GeV) of the shower.

densityfloat, optional

Density (g/cm^3) of ice.

crit_energyfloat, optional

Critical energy (GeV) for shower formation.

rad_lengthfloat, optional

Hadron radiation length (g/cm^2) in ice.

int_lengthfloat, optional

Interaction length (g/cm^2) in ice.

scale_factorfloat, optional

Scale factor S_0 which multiplies the number of particles in the shower.

Returns

array_like

The number of particles at the given distance along the shower.

Notes

Profile calculated by the Gaisser-Hillas model [1], based on equation 1 of the Alvarez hadronic shower paper [2].

References

1

T.K. Gaisser & A.M. Hillas “Reliability of the Method of Constant Intensity Cuts for Reconstructing the Average Development of Vertical Showers.” ICRC proceedings, 353 (1977).

2

J. Alvarez-Muniz & E. Zas, “EeV Hadronic Showers in Ice: The LPM effect.” ICRC proceedings, 17-25 (1999). arXiv:astro-ph/9906347