Feb. 18: Low-E Iron LDF with 5MeV cutoff

First pass at parameterizing this LDF. Seems not to have worked very well -
now changing the order in which I fit the LDF parameters (starting with eta instead of alpha).

Second pass at parameterizing the LDF.

Code used:

• fit.C (fits LDF to AIRES output)
• NKG.C (holder for LDF function)
• ntup_e.C (plots electron parameter trends)
• ntup_m.C (plots muon parameter trends)

Final LDF: Plots of Iron LDF superimposed on AIRES data

The LDF has the same general form as the AGASA LDF:

The parameters, as functions of zenith angle and energy, are:

• MR_e = 82.0
• MR_m = 102.5
• delta_e = 0.4
• delta_m = -0.9
• eta_e = 0.307 + 3.656 cos(ZA)
• eta_m = 1.247 + 0.8214 cos(ZA)
• alpha_e = 1.429 + 0.6220 (log(E) - 17)
• alpha_m = 0.5647 + 0.06972 (log(E) - 17)
• log(C_e) = 1.88 + 1.0 (log(E)-17) + 5.0 (cos(ZA) - .85)
• log(C_m) = 0.78 + 0.9 (log(E)-17) + 1.2 (cos(ZA) - .85)

For comparison: The KASCADE LDF.
The electron LDF now agrees very closely with KASCADE.

• At 10^16 eV = 10^7 GeV, Ne = 4e5 = 10e5.6 -> Open diamond on KASCADE plot.
  • I(50m) ~ 10/m^2
  • I(100m) ~ 2/m^2
• At 10^16.5 eV = 10^7.5 GeV, Ne = 2e6 = 10e6.3 -> Between shaded diamond and open triangle on KASCADE plot.
  • I(50m) ~ 50/m^2
  • I(100m) ~ 10/m^2