analyses

A collection of functions and parametrizations of physic analyses.

radiation energy

Implementation of all parametrizations published in Glaser et al., JCAP 09(2016)024.

analyses.radiationenergy.get_clipping(dxmax)

get clipping correction

Parameters:dxmax (float) – distance to shower maximum in g/cm^2
Returns:fraction of radiation energy that is radiated in the atmosphere
Return type:float
analyses.radiationenergy.get_a(rho)

get relative charge excess strength

Parameters:rho (float) – density at shower maximum in kg/m^3
Returns:relative charge excess strength a
Return type:float
analyses.radiationenergy.get_a_zenith(zenith)

get relative charge excess strength wo Xmax information

Parameters:zentith (float) – zenith angle in rad according to radiotools default coordinate system
Returns:relative charge excess strength a
Return type:float
analyses.radiationenergy.get_S(Erad, sinalpha, density, p0=0.250524463912, p1=-2.95290494, b_scale=1.0, b=1.8)

get corrected radiation energy (S_RD)

Parameters:
  • Erad (float) – radiation energy (in eV)
  • sinalpha (float) – sine of angle between shower axis and geomagnetic field
  • density (float) – density at shower maximum in kg/m^3
Returns:

corrected radiation energy (in eV)
Return type:

float
analyses.radiationenergy.get_S_zenith(erad, sinalpha, zeniths, b_scale=1.0, p0=0.239, p1=-3.13)

get corrected radiation energy (S_RD) wo xmax information

Parameters:
  • Erad (float) – radiation energy (in eV)
  • sinalpha (float) – sine of angle between shower axis and geomagnetic field
  • density (float) – density at shower maximum in kg/m^3
Returns:

corrected radiation energy wo Xmax information (in eV)
Return type:

float
analyses.radiationenergy.get_radiation_energy(Srd, sinalpha, density, p0=0.250524463912, p1=-2.95290494, b_scale=1.0, b=1.8)

get radiation energy (S_RD)

Parameters:
  • Srd (float) – corrected radiation energy (in eV)
  • sinalpha (float) – sine of angle between shower axis and geomagnetic field
  • density (float) – density at shower maximum in kg/m^3
Returns:

radiation energy (in eV)
Return type:

float