diff --git a/corsikaio/longitudinal.py b/corsikaio/longitudinal.py index 63f23bc..c610254 100644 --- a/corsikaio/longitudinal.py +++ b/corsikaio/longitudinal.py @@ -4,8 +4,8 @@ import re import numpy as np -PARTICLE_HEADER_RE = re.compile(r"LONGITUDINAL DISTRIBUTION IN\s+(\d+)\s+(SLANT\s+)?STEPS OF\s+(\d+(?:.\d*)?) G\/CM\*\*2 FOR SHOWER\s+(\d+)") -ENERGY_HEADER_RE = re.compile(r"LONGITUDINAL ENERGY DEPOSIT IN\s+(\d+)\s+(SLANT\s+)?STEPS OF\s+(\d+(?:.\d*)?) G\/CM\*\*2 FOR SHOWER\s+(\d+)") +PARTICLE_HEADER_RE = re.compile(r"LONGITUDINAL DISTRIBUTION IN\s+(\d+)\s+(SLANT|VERTICAL)\s+STEPS OF\s+(\d+(?:.\d*)?) G\/CM\*\*2 FOR SHOWER\s+(\d+)") +ENERGY_HEADER_RE = re.compile(r"LONGITUDINAL ENERGY DEPOSIT IN\s+(\d+)\s+(SLANT|VERTICAL)\s+STEPS OF\s+(\d+(?:.\d*)?) G\/CM\*\*2 FOR SHOWER\s+(\d+)") ENERGY_COLUMNS = [ "depth", @@ -72,7 +72,8 @@ def read_longitudinal_distributions(path): n_steps = int(match.group(1)) longi = dict( shower=int(match.group(4)), - slant=match.group(2) is not None, + n_steps=n_steps, + slant=match.group(2) == "SLANT", step_width=float(match.group(3)), ) diff --git a/tests/resources/corsika_77500_particle.long b/tests/resources/corsika_77500_particle_slant.long similarity index 100% rename from tests/resources/corsika_77500_particle.long rename to tests/resources/corsika_77500_particle_slant.long diff --git a/tests/resources/corsika_77500_particle_vertical.long b/tests/resources/corsika_77500_particle_vertical.long new file mode 100644 index 0000000..220a149 --- /dev/null +++ b/tests/resources/corsika_77500_particle_vertical.long @@ -0,0 +1,426 @@ + LONGITUDINAL DISTRIBUTION IN 208 VERTICAL STEPS OF 5. G/CM**2 FOR SHOWER 1 + DEPTH GAMMAS POSITRONS ELECTRONS MU+ MU- HADRONS CHARGED NUCLEI CHERENKOV + 5.0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 + 10.0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 + 15.0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 + 20.0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 + 25.0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 + 30.0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 + 35.0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 + 40.0 4.60000E+01 1.00000E+00 1.00000E+00 4.00000E+00 5.00000E+00 2.90000E+01 3.50000E+01 0.00000E+00 0.00000E+00 + 45.0 5.20000E+01 7.00000E+00 7.00000E+00 1.10000E+01 9.00000E+00 1.80000E+01 4.80000E+01 0.00000E+00 0.00000E+00 + 50.0 1.08000E+02 1.30000E+01 1.60000E+01 2.20000E+01 1.60000E+01 1.80000E+01 8.10000E+01 0.00000E+00 0.00000E+00 + 55.0 1.56000E+02 2.40000E+01 3.00000E+01 2.30000E+01 1.80000E+01 1.50000E+01 1.05000E+02 0.00000E+00 0.00000E+00 + 60.0 2.57000E+02 3.40000E+01 4.20000E+01 2.30000E+01 1.80000E+01 3.60000E+01 1.45000E+02 0.00000E+00 0.00000E+00 + 65.0 3.60000E+02 6.40000E+01 7.70000E+01 2.90000E+01 2.10000E+01 2.70000E+01 2.11000E+02 0.00000E+00 0.00000E+00 + 70.0 5.26000E+02 8.20000E+01 1.01000E+02 3.00000E+01 2.30000E+01 2.40000E+01 2.53000E+02 0.00000E+00 0.00000E+00 + 75.0 7.39000E+02 1.20000E+02 1.40000E+02 3.10000E+01 2.60000E+01 2.50000E+01 3.35000E+02 0.00000E+00 0.00000E+00 + 80.0 1.00400E+03 1.57000E+02 1.97000E+02 3.20000E+01 2.80000E+01 2.20000E+01 4.30000E+02 0.00000E+00 0.00000E+00 + 85.0 1.33400E+03 1.89000E+02 2.50000E+02 3.10000E+01 2.90000E+01 2.00000E+01 5.13000E+02 0.00000E+00 0.00000E+00 + 90.0 1.68800E+03 2.41000E+02 3.17000E+02 3.30000E+01 3.10000E+01 2.40000E+01 6.39000E+02 0.00000E+00 0.00000E+00 + 95.0 2.10800E+03 2.96000E+02 3.92000E+02 3.40000E+01 3.00000E+01 2.50000E+01 7.70000E+02 0.00000E+00 0.00000E+00 + 100.0 2.65600E+03 3.83000E+02 5.18000E+02 4.10000E+01 3.70000E+01 2.80000E+01 9.98000E+02 0.00000E+00 0.00000E+00 + 105.0 3.23300E+03 4.86000E+02 6.34000E+02 4.40000E+01 4.00000E+01 3.20000E+01 1.22700E+03 0.00000E+00 0.00000E+00 + 110.0 3.96300E+03 5.95000E+02 7.79000E+02 4.30000E+01 4.00000E+01 2.90000E+01 1.47800E+03 0.00000E+00 0.00000E+00 + 115.0 4.84900E+03 6.90000E+02 9.03000E+02 4.30000E+01 3.90000E+01 2.80000E+01 1.69500E+03 0.00000E+00 0.00000E+00 + 120.0 5.83000E+03 8.01000E+02 1.07900E+03 4.30000E+01 3.70000E+01 3.00000E+01 1.98100E+03 0.00000E+00 0.00000E+00 + 125.0 6.98300E+03 9.26000E+02 1.28500E+03 4.40000E+01 3.90000E+01 3.20000E+01 2.31500E+03 0.00000E+00 0.00000E+00 + 130.0 8.22600E+03 1.10000E+03 1.50000E+03 4.90000E+01 4.10000E+01 5.30000E+01 2.72500E+03 0.00000E+00 0.00000E+00 + 135.0 9.66600E+03 1.24500E+03 1.69500E+03 5.40000E+01 5.30000E+01 5.00000E+01 3.07900E+03 0.00000E+00 0.00000E+00 + 140.0 1.13250E+04 1.50800E+03 2.06200E+03 5.80000E+01 5.90000E+01 7.70000E+01 3.74000E+03 0.00000E+00 0.00000E+00 + 145.0 1.31000E+04 1.69900E+03 2.30900E+03 6.20000E+01 6.30000E+01 6.80000E+01 4.17700E+03 0.00000E+00 0.00000E+00 + 150.0 1.52720E+04 1.91900E+03 2.66100E+03 6.60000E+01 6.70000E+01 6.30000E+01 4.75200E+03 0.00000E+00 0.00000E+00 + 155.0 1.74100E+04 2.21900E+03 3.06800E+03 6.50000E+01 7.00000E+01 6.00000E+01 5.45800E+03 0.00000E+00 0.00000E+00 + 160.0 1.99510E+04 2.47800E+03 3.41100E+03 6.90000E+01 7.00000E+01 5.40000E+01 6.05800E+03 0.00000E+00 0.00000E+00 + 165.0 2.25920E+04 2.77400E+03 3.86000E+03 7.30000E+01 7.10000E+01 4.90000E+01 6.80300E+03 0.00000E+00 0.00000E+00 + 170.0 2.53870E+04 3.15500E+03 4.35000E+03 7.50000E+01 6.80000E+01 5.20000E+01 7.67600E+03 0.00000E+00 0.00000E+00 + 175.0 2.84800E+04 3.43300E+03 4.76900E+03 7.40000E+01 6.80000E+01 5.10000E+01 8.37100E+03 0.00000E+00 0.00000E+00 + 180.0 3.17090E+04 3.76800E+03 5.20800E+03 7.80000E+01 7.30000E+01 6.10000E+01 9.16400E+03 0.00000E+00 0.00000E+00 + 185.0 3.52390E+04 4.11600E+03 5.69400E+03 8.00000E+01 7.40000E+01 6.00000E+01 9.99800E+03 0.00000E+00 0.00000E+00 + 190.0 3.88850E+04 4.49000E+03 6.35400E+03 8.10000E+01 7.80000E+01 7.10000E+01 1.10480E+04 0.00000E+00 0.00000E+00 + 195.0 4.28330E+04 4.80800E+03 6.80200E+03 8.70000E+01 8.20000E+01 6.00000E+01 1.18130E+04 0.00000E+00 0.00000E+00 + 200.0 4.69050E+04 5.25100E+03 7.44000E+03 8.60000E+01 8.40000E+01 7.50000E+01 1.29070E+04 0.00000E+00 0.00000E+00 + 205.0 5.11800E+04 5.61200E+03 8.08200E+03 8.80000E+01 8.60000E+01 7.10000E+01 1.39090E+04 0.00000E+00 0.00000E+00 + 210.0 5.56580E+04 6.05200E+03 8.70200E+03 8.90000E+01 9.00000E+01 7.30000E+01 1.49770E+04 0.00000E+00 0.00000E+00 + 215.0 6.02570E+04 6.51900E+03 9.41600E+03 9.10000E+01 9.30000E+01 7.30000E+01 1.61620E+04 0.00000E+00 0.00000E+00 + 220.0 6.50710E+04 6.89000E+03 1.00350E+04 9.70000E+01 9.70000E+01 8.40000E+01 1.71710E+04 0.00000E+00 0.00000E+00 + 225.0 6.98930E+04 7.34800E+03 1.07120E+04 1.02000E+02 1.02000E+02 9.00000E+01 1.83230E+04 0.00000E+00 0.00000E+00 + 230.0 7.52060E+04 7.84100E+03 1.15410E+04 1.04000E+02 1.06000E+02 7.90000E+01 1.96410E+04 0.00000E+00 0.00000E+00 + 235.0 8.05080E+04 8.26100E+03 1.22720E+04 1.07000E+02 1.08000E+02 1.04000E+02 2.08230E+04 0.00000E+00 0.00000E+00 + 240.0 8.61330E+04 8.77700E+03 1.27960E+04 1.11000E+02 1.14000E+02 9.50000E+01 2.18650E+04 0.00000E+00 0.00000E+00 + 245.0 9.18960E+04 9.15800E+03 1.35100E+04 1.14000E+02 1.18000E+02 8.70000E+01 2.29610E+04 0.00000E+00 0.00000E+00 + 250.0 9.75440E+04 9.67700E+03 1.42800E+04 1.19000E+02 1.17000E+02 1.22000E+02 2.42800E+04 0.00000E+00 0.00000E+00 + 255.0 1.03479E+05 1.01650E+04 1.49610E+04 1.20000E+02 1.22000E+02 1.51000E+02 2.54850E+04 0.00000E+00 0.00000E+00 + 260.0 1.09275E+05 1.07210E+04 1.58490E+04 1.34000E+02 1.29000E+02 1.54000E+02 2.69490E+04 0.00000E+00 0.00000E+00 + 265.0 1.15398E+05 1.11950E+04 1.66500E+04 1.39000E+02 1.38000E+02 1.44000E+02 2.82270E+04 0.00000E+00 0.00000E+00 + 270.0 1.21569E+05 1.17320E+04 1.73930E+04 1.42000E+02 1.42000E+02 1.47000E+02 2.95170E+04 0.00000E+00 0.00000E+00 + 275.0 1.27908E+05 1.23030E+04 1.82880E+04 1.48000E+02 1.47000E+02 1.35000E+02 3.09860E+04 0.00000E+00 0.00000E+00 + 280.0 1.34197E+05 1.27260E+04 1.90330E+04 1.52000E+02 1.52000E+02 1.28000E+02 3.21550E+04 0.00000E+00 0.00000E+00 + 285.0 1.40711E+05 1.33070E+04 1.97870E+04 1.57000E+02 1.56000E+02 1.43000E+02 3.35120E+04 0.00000E+00 0.00000E+00 + 290.0 1.47403E+05 1.38770E+04 2.06500E+04 1.68000E+02 1.67000E+02 2.35000E+02 3.50510E+04 0.00000E+00 0.00000E+00 + 295.0 1.53661E+05 1.42110E+04 2.12600E+04 1.77000E+02 1.84000E+02 2.32000E+02 3.60120E+04 0.00000E+00 0.00000E+00 + 300.0 1.60274E+05 1.46760E+04 2.19830E+04 1.89000E+02 1.97000E+02 2.06000E+02 3.72010E+04 0.00000E+00 0.00000E+00 + 305.0 1.66686E+05 1.51710E+04 2.27370E+04 1.95000E+02 2.06000E+02 2.08000E+02 3.84650E+04 0.00000E+00 0.00000E+00 + 310.0 1.73172E+05 1.56820E+04 2.35500E+04 2.03000E+02 2.15000E+02 2.13000E+02 3.98080E+04 0.00000E+00 0.00000E+00 + 315.0 1.79356E+05 1.61260E+04 2.43570E+04 2.10000E+02 2.19000E+02 2.00000E+02 4.10590E+04 0.00000E+00 0.00000E+00 + 320.0 1.85779E+05 1.65700E+04 2.51320E+04 2.12000E+02 2.21000E+02 1.90000E+02 4.22720E+04 0.00000E+00 0.00000E+00 + 325.0 1.92332E+05 1.69370E+04 2.57080E+04 2.13000E+02 2.22000E+02 2.06000E+02 4.32300E+04 0.00000E+00 0.00000E+00 + 330.0 1.98342E+05 1.75540E+04 2.63220E+04 2.22000E+02 2.33000E+02 2.48000E+02 4.45090E+04 0.00000E+00 0.00000E+00 + 335.0 2.04747E+05 1.77840E+04 2.71320E+04 2.34000E+02 2.53000E+02 2.59000E+02 4.55940E+04 0.00000E+00 0.00000E+00 + 340.0 2.11223E+05 1.81940E+04 2.75840E+04 2.49000E+02 2.68000E+02 2.41000E+02 4.64700E+04 0.00000E+00 0.00000E+00 + 345.0 2.17967E+05 1.85530E+04 2.83460E+04 2.66000E+02 2.73000E+02 2.29000E+02 4.75960E+04 0.00000E+00 0.00000E+00 + 350.0 2.24201E+05 1.89210E+04 2.89770E+04 2.73000E+02 2.79000E+02 2.19000E+02 4.85970E+04 0.00000E+00 0.00000E+00 + 355.0 2.30261E+05 1.93080E+04 2.94580E+04 2.84000E+02 2.86000E+02 2.15000E+02 4.94820E+04 0.00000E+00 0.00000E+00 + 360.0 2.36068E+05 1.96640E+04 3.02500E+04 2.90000E+02 2.94000E+02 2.51000E+02 5.06730E+04 0.00000E+00 0.00000E+00 + 365.0 2.41831E+05 1.99850E+04 3.08850E+04 2.97000E+02 3.01000E+02 2.49000E+02 5.16390E+04 0.00000E+00 0.00000E+00 + 370.0 2.47649E+05 2.03840E+04 3.14520E+04 3.14000E+02 3.10000E+02 2.59000E+02 5.26320E+04 0.00000E+00 0.00000E+00 + 375.0 2.53107E+05 2.08120E+04 3.21820E+04 3.23000E+02 3.14000E+02 2.64000E+02 5.38060E+04 0.00000E+00 0.00000E+00 + 380.0 2.58436E+05 2.11110E+04 3.26990E+04 3.27000E+02 3.23000E+02 2.62000E+02 5.46340E+04 0.00000E+00 0.00000E+00 + 385.0 2.63239E+05 2.15060E+04 3.33440E+04 3.37000E+02 3.41000E+02 2.90000E+02 5.57270E+04 0.00000E+00 0.00000E+00 + 390.0 2.68437E+05 2.18630E+04 3.39470E+04 3.44000E+02 3.50000E+02 2.79000E+02 5.66990E+04 0.00000E+00 0.00000E+00 + 395.0 2.73656E+05 2.20540E+04 3.42830E+04 3.56000E+02 3.52000E+02 2.67000E+02 5.72260E+04 0.00000E+00 0.00000E+00 + 400.0 2.78515E+05 2.21640E+04 3.46460E+04 3.63000E+02 3.58000E+02 2.52000E+02 5.76990E+04 0.00000E+00 0.00000E+00 + 405.0 2.83314E+05 2.24170E+04 3.49130E+04 3.69000E+02 3.66000E+02 2.46000E+02 5.82270E+04 0.00000E+00 0.00000E+00 + 410.0 2.87783E+05 2.25050E+04 3.51970E+04 3.70000E+02 3.71000E+02 2.33000E+02 5.85960E+04 0.00000E+00 0.00000E+00 + 415.0 2.91653E+05 2.27750E+04 3.55290E+04 3.74000E+02 3.74000E+02 2.26000E+02 5.92030E+04 0.00000E+00 0.00000E+00 + 420.0 2.94852E+05 2.29060E+04 3.54290E+04 3.77000E+02 3.75000E+02 2.12000E+02 5.92270E+04 0.00000E+00 0.00000E+00 + 425.0 2.97303E+05 2.30120E+04 3.56740E+04 3.79000E+02 3.82000E+02 2.08000E+02 5.95830E+04 0.00000E+00 0.00000E+00 + 430.0 2.99999E+05 2.30050E+04 3.57820E+04 3.84000E+02 3.87000E+02 2.08000E+02 5.96910E+04 0.00000E+00 0.00000E+00 + 435.0 3.02419E+05 2.29730E+04 3.57070E+04 3.89000E+02 3.90000E+02 1.97000E+02 5.95800E+04 0.00000E+00 0.00000E+00 + 440.0 3.04464E+05 2.28520E+04 3.57730E+04 3.95000E+02 3.89000E+02 2.04000E+02 5.95430E+04 0.00000E+00 0.00000E+00 + 445.0 3.06027E+05 2.26200E+04 3.56180E+04 3.97000E+02 3.91000E+02 2.01000E+02 5.91560E+04 0.00000E+00 0.00000E+00 + 450.0 3.07179E+05 2.25840E+04 3.59000E+04 4.01000E+02 3.87000E+02 1.95000E+02 5.93960E+04 0.00000E+00 0.00000E+00 + 455.0 3.08669E+05 2.26110E+04 3.59680E+04 4.00000E+02 3.95000E+02 2.17000E+02 5.95180E+04 0.00000E+00 0.00000E+00 + 460.0 3.09864E+05 2.26210E+04 3.58720E+04 4.02000E+02 4.01000E+02 2.28000E+02 5.94470E+04 0.00000E+00 0.00000E+00 + 465.0 3.10209E+05 2.23560E+04 3.54470E+04 4.08000E+02 4.05000E+02 2.57000E+02 5.87900E+04 0.00000E+00 0.00000E+00 + 470.0 3.10710E+05 2.22750E+04 3.57350E+04 4.15000E+02 4.12000E+02 2.40000E+02 5.89920E+04 0.00000E+00 0.00000E+00 + 475.0 3.10501E+05 2.21160E+04 3.52190E+04 4.21000E+02 4.18000E+02 2.25000E+02 5.83140E+04 0.00000E+00 0.00000E+00 + 480.0 3.09932E+05 2.18500E+04 3.51730E+04 4.23000E+02 4.21000E+02 2.21000E+02 5.80040E+04 0.00000E+00 0.00000E+00 + 485.0 3.09250E+05 2.16260E+04 3.48100E+04 4.22000E+02 4.20000E+02 2.17000E+02 5.74150E+04 0.00000E+00 0.00000E+00 + 490.0 3.08324E+05 2.14800E+04 3.46380E+04 4.19000E+02 4.19000E+02 2.16000E+02 5.70940E+04 0.00000E+00 0.00000E+00 + 495.0 3.07262E+05 2.12450E+04 3.46130E+04 4.22000E+02 4.22000E+02 2.11000E+02 5.68360E+04 0.00000E+00 0.00000E+00 + 500.0 3.06015E+05 2.07670E+04 3.39540E+04 4.23000E+02 4.29000E+02 2.43000E+02 5.57390E+04 0.00000E+00 0.00000E+00 + 505.0 3.04626E+05 2.04850E+04 3.37170E+04 4.31000E+02 4.35000E+02 2.36000E+02 5.52260E+04 0.00000E+00 0.00000E+00 + 510.0 3.02471E+05 2.01670E+04 3.32330E+04 4.31000E+02 4.42000E+02 2.63000E+02 5.44530E+04 0.00000E+00 0.00000E+00 + 515.0 3.00299E+05 2.01080E+04 3.27130E+04 4.39000E+02 4.48000E+02 2.70000E+02 5.38920E+04 0.00000E+00 0.00000E+00 + 520.0 2.98051E+05 2.00070E+04 3.25410E+04 4.47000E+02 4.54000E+02 2.44000E+02 5.36100E+04 0.00000E+00 0.00000E+00 + 525.0 2.95700E+05 1.95290E+04 3.20240E+04 4.48000E+02 4.60000E+02 2.50000E+02 5.26220E+04 0.00000E+00 0.00000E+00 + 530.0 2.93573E+05 1.93170E+04 3.16670E+04 4.53000E+02 4.62000E+02 2.56000E+02 5.20600E+04 0.00000E+00 0.00000E+00 + 535.0 2.90808E+05 1.92130E+04 3.12990E+04 4.59000E+02 4.69000E+02 2.58000E+02 5.16000E+04 0.00000E+00 0.00000E+00 + 540.0 2.88287E+05 1.88900E+04 3.09060E+04 4.62000E+02 4.74000E+02 2.52000E+02 5.08850E+04 0.00000E+00 0.00000E+00 + 545.0 2.85470E+05 1.86030E+04 3.07160E+04 4.63000E+02 4.74000E+02 2.38000E+02 5.03960E+04 0.00000E+00 0.00000E+00 + 550.0 2.82306E+05 1.83650E+04 3.03000E+04 4.67000E+02 4.80000E+02 2.38000E+02 4.97530E+04 0.00000E+00 0.00000E+00 + 555.0 2.79826E+05 1.81980E+04 3.00980E+04 4.74000E+02 4.82000E+02 2.35000E+02 4.93930E+04 0.00000E+00 0.00000E+00 + 560.0 2.76972E+05 1.77980E+04 2.95720E+04 4.77000E+02 4.88000E+02 2.23000E+02 4.84690E+04 0.00000E+00 0.00000E+00 + 565.0 2.73777E+05 1.75900E+04 2.91140E+04 4.84000E+02 4.81000E+02 2.20000E+02 4.78020E+04 0.00000E+00 0.00000E+00 + 570.0 2.70802E+05 1.73820E+04 2.86080E+04 4.81000E+02 4.78000E+02 2.28000E+02 4.70900E+04 0.00000E+00 0.00000E+00 + 575.0 2.67547E+05 1.71720E+04 2.82580E+04 4.83000E+02 4.82000E+02 2.21000E+02 4.65280E+04 0.00000E+00 0.00000E+00 + 580.0 2.63990E+05 1.69720E+04 2.79080E+04 4.90000E+02 4.84000E+02 2.29000E+02 4.59970E+04 0.00000E+00 0.00000E+00 + 585.0 2.60988E+05 1.66870E+04 2.74160E+04 4.93000E+02 4.82000E+02 2.25000E+02 4.52170E+04 0.00000E+00 0.00000E+00 + 590.0 2.57894E+05 1.63960E+04 2.68560E+04 4.96000E+02 4.84000E+02 2.24000E+02 4.43740E+04 0.00000E+00 0.00000E+00 + 595.0 2.54317E+05 1.60590E+04 2.64700E+04 4.98000E+02 4.85000E+02 2.24000E+02 4.36520E+04 0.00000E+00 0.00000E+00 + 600.0 2.51089E+05 1.59610E+04 2.63290E+04 5.02000E+02 4.84000E+02 2.26000E+02 4.34140E+04 0.00000E+00 0.00000E+00 + 605.0 2.47517E+05 1.55260E+04 2.56440E+04 5.00000E+02 4.85000E+02 2.21000E+02 4.22900E+04 0.00000E+00 0.00000E+00 + 610.0 2.44014E+05 1.52410E+04 2.52280E+04 5.01000E+02 4.90000E+02 2.11000E+02 4.15900E+04 0.00000E+00 0.00000E+00 + 615.0 2.40635E+05 1.50290E+04 2.47740E+04 5.01000E+02 4.87000E+02 2.05000E+02 4.09200E+04 0.00000E+00 0.00000E+00 + 620.0 2.36993E+05 1.48030E+04 2.45690E+04 5.05000E+02 4.93000E+02 2.38000E+02 4.05220E+04 0.00000E+00 0.00000E+00 + 625.0 2.33338E+05 1.44340E+04 2.39840E+04 5.06000E+02 5.03000E+02 2.20000E+02 3.95650E+04 0.00000E+00 0.00000E+00 + 630.0 2.29405E+05 1.41480E+04 2.37770E+04 5.07000E+02 5.05000E+02 2.06000E+02 3.90640E+04 0.00000E+00 0.00000E+00 + 635.0 2.25980E+05 1.38240E+04 2.32040E+04 5.11000E+02 5.04000E+02 2.17000E+02 3.81770E+04 0.00000E+00 0.00000E+00 + 640.0 2.22371E+05 1.36110E+04 2.27860E+04 5.11000E+02 5.06000E+02 2.05000E+02 3.75390E+04 0.00000E+00 0.00000E+00 + 645.0 2.18844E+05 1.34310E+04 2.24500E+04 5.12000E+02 5.05000E+02 2.05000E+02 3.70240E+04 0.00000E+00 0.00000E+00 + 650.0 2.15482E+05 1.32880E+04 2.20800E+04 5.11000E+02 5.07000E+02 1.93000E+02 3.65040E+04 0.00000E+00 0.00000E+00 + 655.0 2.11862E+05 1.30210E+04 2.16490E+04 5.12000E+02 5.12000E+02 2.01000E+02 3.58200E+04 0.00000E+00 0.00000E+00 + 660.0 2.08098E+05 1.27330E+04 2.15090E+04 5.15000E+02 5.13000E+02 1.84000E+02 3.53820E+04 0.00000E+00 0.00000E+00 + 665.0 2.04701E+05 1.24090E+04 2.09950E+04 5.17000E+02 5.13000E+02 1.96000E+02 3.45610E+04 0.00000E+00 0.00000E+00 + 670.0 2.01462E+05 1.22460E+04 2.04620E+04 5.19000E+02 5.15000E+02 1.83000E+02 3.38560E+04 0.00000E+00 0.00000E+00 + 675.0 1.97724E+05 1.21360E+04 2.03500E+04 5.25000E+02 5.13000E+02 1.82000E+02 3.36370E+04 0.00000E+00 0.00000E+00 + 680.0 1.94304E+05 1.19900E+04 2.00050E+04 5.27000E+02 5.15000E+02 1.78000E+02 3.31460E+04 0.00000E+00 0.00000E+00 + 685.0 1.91024E+05 1.17770E+04 1.96490E+04 5.29000E+02 5.18000E+02 1.69000E+02 3.25740E+04 0.00000E+00 0.00000E+00 + 690.0 1.87700E+05 1.15590E+04 1.94480E+04 5.29000E+02 5.21000E+02 1.75000E+02 3.21660E+04 0.00000E+00 0.00000E+00 + 695.0 1.84667E+05 1.13460E+04 1.91260E+04 5.37000E+02 5.20000E+02 1.78000E+02 3.16410E+04 0.00000E+00 0.00000E+00 + 700.0 1.81521E+05 1.10580E+04 1.86030E+04 5.40000E+02 5.25000E+02 1.77000E+02 3.08420E+04 0.00000E+00 0.00000E+00 + 705.0 1.78412E+05 1.08890E+04 1.82390E+04 5.37000E+02 5.25000E+02 1.75000E+02 3.03070E+04 0.00000E+00 0.00000E+00 + 710.0 1.75352E+05 1.06040E+04 1.78510E+04 5.39000E+02 5.26000E+02 1.69000E+02 2.96330E+04 0.00000E+00 0.00000E+00 + 715.0 1.72377E+05 1.04500E+04 1.76650E+04 5.39000E+02 5.33000E+02 1.59000E+02 2.92910E+04 0.00000E+00 0.00000E+00 + 720.0 1.69508E+05 1.00680E+04 1.72840E+04 5.36000E+02 5.37000E+02 1.57000E+02 2.85250E+04 0.00000E+00 0.00000E+00 + 725.0 1.66726E+05 1.00330E+04 1.69970E+04 5.40000E+02 5.39000E+02 1.61000E+02 2.82110E+04 0.00000E+00 0.00000E+00 + 730.0 1.63815E+05 9.87800E+03 1.66710E+04 5.48000E+02 5.40000E+02 1.60000E+02 2.77350E+04 0.00000E+00 0.00000E+00 + 735.0 1.61017E+05 9.67800E+03 1.61400E+04 5.55000E+02 5.41000E+02 1.52000E+02 2.70040E+04 0.00000E+00 0.00000E+00 + 740.0 1.58120E+05 9.60700E+03 1.59080E+04 5.53000E+02 5.39000E+02 1.52000E+02 2.66970E+04 0.00000E+00 0.00000E+00 + 745.0 1.55271E+05 9.38100E+03 1.57370E+04 5.53000E+02 5.36000E+02 1.48000E+02 2.62990E+04 0.00000E+00 0.00000E+00 + 750.0 1.52668E+05 9.16000E+03 1.51860E+04 5.49000E+02 5.33000E+02 1.39000E+02 2.55140E+04 0.00000E+00 0.00000E+00 + 755.0 1.49994E+05 9.05400E+03 1.48710E+04 5.48000E+02 5.29000E+02 1.31000E+02 2.50790E+04 0.00000E+00 0.00000E+00 + 760.0 1.47311E+05 8.80300E+03 1.46570E+04 5.46000E+02 5.31000E+02 1.29000E+02 2.46110E+04 0.00000E+00 0.00000E+00 + 765.0 1.44318E+05 8.58200E+03 1.45280E+04 5.47000E+02 5.29000E+02 1.22000E+02 2.42520E+04 0.00000E+00 0.00000E+00 + 770.0 1.41718E+05 8.46300E+03 1.42200E+04 5.43000E+02 5.26000E+02 1.25000E+02 2.38200E+04 0.00000E+00 0.00000E+00 + 775.0 1.38912E+05 8.12400E+03 1.38130E+04 5.42000E+02 5.25000E+02 1.23000E+02 2.30740E+04 0.00000E+00 0.00000E+00 + 780.0 1.36174E+05 7.95200E+03 1.35530E+04 5.41000E+02 5.26000E+02 1.19000E+02 2.26430E+04 0.00000E+00 0.00000E+00 + 785.0 1.33508E+05 7.77300E+03 1.33320E+04 5.41000E+02 5.23000E+02 1.14000E+02 2.22370E+04 0.00000E+00 0.00000E+00 + 790.0 1.30879E+05 7.68000E+03 1.29810E+04 5.39000E+02 5.21000E+02 1.10000E+02 2.17850E+04 0.00000E+00 0.00000E+00 + 795.0 1.28101E+05 7.53400E+03 1.29710E+04 5.33000E+02 5.19000E+02 1.07000E+02 2.16200E+04 0.00000E+00 0.00000E+00 + 800.0 1.25435E+05 7.40100E+03 1.26510E+04 5.29000E+02 5.17000E+02 1.10000E+02 2.11630E+04 0.00000E+00 0.00000E+00 + 805.0 1.22848E+05 7.18900E+03 1.22870E+04 5.34000E+02 5.14000E+02 1.06000E+02 2.05840E+04 0.00000E+00 0.00000E+00 + 810.0 1.20216E+05 7.09300E+03 1.19650E+04 5.31000E+02 5.09000E+02 1.08000E+02 2.01600E+04 0.00000E+00 0.00000E+00 + 815.0 1.17590E+05 6.92300E+03 1.17630E+04 5.31000E+02 5.10000E+02 1.06000E+02 1.97860E+04 0.00000E+00 0.00000E+00 + 820.0 1.15076E+05 6.69100E+03 1.13950E+04 5.29000E+02 5.11000E+02 1.02000E+02 1.91800E+04 0.00000E+00 0.00000E+00 + 825.0 1.12473E+05 6.54400E+03 1.11860E+04 5.28000E+02 5.06000E+02 1.00000E+02 1.88190E+04 0.00000E+00 0.00000E+00 + 830.0 1.10316E+05 6.29300E+03 1.09380E+04 5.27000E+02 4.99000E+02 9.60000E+01 1.83100E+04 0.00000E+00 0.00000E+00 + 835.0 1.07842E+05 6.11300E+03 1.07820E+04 5.27000E+02 4.94000E+02 9.40000E+01 1.79680E+04 0.00000E+00 0.00000E+00 + 840.0 1.05392E+05 5.87200E+03 1.03640E+04 5.25000E+02 4.93000E+02 1.03000E+02 1.73180E+04 0.00000E+00 0.00000E+00 + 845.0 1.02634E+05 5.84800E+03 1.02530E+04 5.22000E+02 4.93000E+02 1.05000E+02 1.71840E+04 0.00000E+00 0.00000E+00 + 850.0 1.00375E+05 5.68200E+03 9.92900E+03 5.19000E+02 4.91000E+02 1.03000E+02 1.66860E+04 0.00000E+00 0.00000E+00 + 855.0 9.79030E+04 5.60000E+03 9.57200E+03 5.16000E+02 4.91000E+02 1.02000E+02 1.62440E+04 0.00000E+00 0.00000E+00 + 860.0 9.56000E+04 5.45900E+03 9.25100E+03 5.13000E+02 4.89000E+02 1.00000E+02 1.57750E+04 0.00000E+00 0.00000E+00 + 865.0 9.30910E+04 5.33100E+03 8.92700E+03 5.11000E+02 4.86000E+02 1.06000E+02 1.53230E+04 0.00000E+00 0.00000E+00 + 870.0 9.08270E+04 5.10300E+03 8.65200E+03 5.12000E+02 4.89000E+02 1.00000E+02 1.48200E+04 0.00000E+00 0.00000E+00 + 875.0 8.84800E+04 4.98700E+03 8.43300E+03 5.11000E+02 4.90000E+02 8.80000E+01 1.44750E+04 0.00000E+00 0.00000E+00 + 880.0 8.61550E+04 4.90200E+03 8.29200E+03 5.09000E+02 4.88000E+02 9.50000E+01 1.42520E+04 0.00000E+00 0.00000E+00 + 885.0 8.38760E+04 4.72800E+03 8.05900E+03 5.09000E+02 4.84000E+02 9.40000E+01 1.38390E+04 0.00000E+00 0.00000E+00 + 890.0 8.17100E+04 4.52900E+03 7.86100E+03 5.07000E+02 4.81000E+02 9.00000E+01 1.34330E+04 0.00000E+00 0.00000E+00 + 895.0 7.94490E+04 4.33800E+03 7.66500E+03 5.07000E+02 4.78000E+02 8.70000E+01 1.30410E+04 0.00000E+00 0.00000E+00 + 900.0 7.73500E+04 4.20400E+03 7.40900E+03 5.08000E+02 4.76000E+02 8.40000E+01 1.26500E+04 0.00000E+00 0.00000E+00 + 905.0 7.52340E+04 4.05400E+03 7.18400E+03 5.07000E+02 4.70000E+02 8.20000E+01 1.22660E+04 0.00000E+00 0.00000E+00 + 910.0 7.30940E+04 3.89700E+03 6.92700E+03 5.06000E+02 4.65000E+02 8.10000E+01 1.18440E+04 0.00000E+00 0.00000E+00 + 915.0 7.12510E+04 3.82700E+03 6.70600E+03 5.05000E+02 4.66000E+02 8.80000E+01 1.15570E+04 0.00000E+00 0.00000E+00 + 920.0 6.92680E+04 3.75000E+03 6.58000E+03 5.06000E+02 4.67000E+02 7.90000E+01 1.13500E+04 0.00000E+00 0.00000E+00 + 925.0 6.72960E+04 3.64900E+03 6.35800E+03 5.05000E+02 4.69000E+02 6.90000E+01 1.10180E+04 0.00000E+00 0.00000E+00 + 930.0 6.54260E+04 3.60200E+03 6.34700E+03 5.00000E+02 4.66000E+02 7.10000E+01 1.09540E+04 0.00000E+00 0.00000E+00 + 935.0 6.38460E+04 3.46600E+03 6.09300E+03 4.97000E+02 4.65000E+02 7.40000E+01 1.05640E+04 0.00000E+00 0.00000E+00 + 940.0 6.20550E+04 3.34800E+03 5.90900E+03 4.95000E+02 4.65000E+02 7.10000E+01 1.02580E+04 0.00000E+00 0.00000E+00 + 945.0 6.02920E+04 3.27500E+03 5.73700E+03 4.93000E+02 4.63000E+02 6.70000E+01 1.00040E+04 0.00000E+00 0.00000E+00 + 950.0 5.86480E+04 3.13500E+03 5.49900E+03 4.93000E+02 4.59000E+02 6.50000E+01 9.62200E+03 0.00000E+00 0.00000E+00 + 955.0 5.70220E+04 3.08200E+03 5.30700E+03 4.90000E+02 4.58000E+02 6.40000E+01 9.37400E+03 0.00000E+00 0.00000E+00 + 960.0 5.52600E+04 3.02800E+03 5.15000E+03 4.87000E+02 4.57000E+02 6.10000E+01 9.15800E+03 0.00000E+00 0.00000E+00 + 965.0 5.36410E+04 2.90900E+03 4.93900E+03 4.85000E+02 4.55000E+02 6.10000E+01 8.82400E+03 0.00000E+00 0.00000E+00 + 970.0 5.19730E+04 2.78900E+03 4.83600E+03 4.84000E+02 4.51000E+02 6.50000E+01 8.59800E+03 0.00000E+00 0.00000E+00 + 975.0 5.05800E+04 2.61700E+03 4.59900E+03 4.87000E+02 4.49000E+02 6.20000E+01 8.18800E+03 0.00000E+00 0.00000E+00 + 980.0 4.90150E+04 2.52500E+03 4.48600E+03 4.88000E+02 4.49000E+02 5.90000E+01 7.98100E+03 0.00000E+00 0.00000E+00 + 985.0 4.75600E+04 2.46200E+03 4.41600E+03 4.86000E+02 4.47000E+02 5.50000E+01 7.84300E+03 0.00000E+00 0.00000E+00 + 990.0 4.60230E+04 2.42400E+03 4.29600E+03 4.86000E+02 4.45000E+02 5.30000E+01 7.68200E+03 0.00000E+00 0.00000E+00 + 995.0 4.46310E+04 2.32900E+03 4.05200E+03 4.87000E+02 4.44000E+02 5.40000E+01 7.34300E+03 0.00000E+00 0.00000E+00 + 1000.0 4.32730E+04 2.28600E+03 4.15500E+03 4.82000E+02 4.42000E+02 5.20000E+01 7.39400E+03 0.00000E+00 0.00000E+00 + 1005.0 4.19280E+04 2.21900E+03 3.95800E+03 4.80000E+02 4.40000E+02 5.00000E+01 7.12400E+03 0.00000E+00 0.00000E+00 + 1010.0 4.05990E+04 2.15700E+03 3.81500E+03 4.77000E+02 4.41000E+02 5.50000E+01 6.92100E+03 0.00000E+00 0.00000E+00 + 1015.0 3.93020E+04 2.09800E+03 3.62300E+03 4.74000E+02 4.39000E+02 5.40000E+01 6.66100E+03 0.00000E+00 0.00000E+00 + 1020.0 3.81790E+04 1.99600E+03 3.52400E+03 4.74000E+02 4.38000E+02 5.20000E+01 6.45900E+03 0.00000E+00 0.00000E+00 + 1025.0 3.67740E+04 1.95700E+03 3.39900E+03 4.73000E+02 4.38000E+02 4.80000E+01 6.29100E+03 0.00000E+00 0.00000E+00 + 1030.0 3.56540E+04 1.89100E+03 3.23200E+03 4.74000E+02 4.36000E+02 4.60000E+01 6.05500E+03 0.00000E+00 0.00000E+00 + 1035.0 3.45410E+04 1.78600E+03 3.06500E+03 4.73000E+02 4.36000E+02 4.30000E+01 5.78000E+03 0.00000E+00 0.00000E+00 + 1040.0 3.42840E+04 1.78400E+03 3.10200E+03 4.73000E+02 4.36000E+02 4.20000E+01 5.81500E+03 0.00000E+00 0.00000E+00 + LONGITUDINAL ENERGY DEPOSIT IN 208 VERTICAL STEPS OF 5. G/CM**2 FOR SHOWER 1 + DEPTH GAMMA EM IONIZ EM CUT MU IONIZ MU CUT HADR IONIZ HADR CUT NEUTRINO SUM + 2.5 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.23472E-02 0.00000E+00 0.00000E+00 2.23472E-02 + 7.5 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.23472E-02 0.00000E+00 0.00000E+00 2.23472E-02 + 12.5 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.23472E-02 0.00000E+00 0.00000E+00 2.23472E-02 + 17.5 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.23472E-02 0.00000E+00 0.00000E+00 2.23472E-02 + 22.5 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.23472E-02 0.00000E+00 0.00000E+00 2.23472E-02 + 27.5 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.23472E-02 0.00000E+00 0.00000E+00 2.23472E-02 + 32.5 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.23472E-02 0.00000E+00 0.00000E+00 2.23472E-02 + 37.5 0.00000E+00 4.69888E-03 0.00000E+00 6.31845E-03 1.69859E+00 8.18001E-02 2.29999E+00 1.17927E+01 1.58841E+01 + 42.5 0.00000E+00 7.28118E-02 2.00000E-03 1.90825E-01 0.00000E+00 2.69490E-01 3.92653E-01 2.71046E+01 2.80324E+01 + 47.5 5.95902E-02 2.70995E-01 5.43391E-03 3.45897E-01 1.70730E+00 2.57208E-01 3.14402E+00 6.39000E+01 6.96905E+01 + 52.5 1.43780E-03 5.01321E-01 1.46203E-02 4.88120E-01 2.65042E-01 1.73211E-01 6.32843E-01 1.45310E+01 1.66076E+01 + 57.5 2.86638E-03 8.68491E-01 6.05813E-02 4.92974E-01 0.00000E+00 1.65116E-01 1.65560E+00 1.81741E+01 2.14197E+01 + 62.5 5.75094E-03 1.27711E+00 4.58980E-02 5.64120E-01 0.00000E+00 3.30626E-01 0.00000E+00 2.19518E+01 2.41753E+01 + 67.5 1.08812E-02 1.97062E+00 1.18463E-01 6.16085E-01 0.00000E+00 2.81098E-01 0.00000E+00 2.82975E+01 3.12946E+01 + 72.5 6.11371E-02 2.71340E+00 1.14692E-01 6.75530E-01 0.00000E+00 2.55117E-01 1.19626E+00 1.50014E+01 2.00175E+01 + 77.5 2.76835E-02 3.72955E+00 1.70860E-01 7.07966E-01 0.00000E+00 2.67394E-01 2.49220E-01 1.59784E+01 2.11311E+01 + 82.5 4.34701E-02 4.82405E+00 2.85762E-01 7.34988E-01 0.00000E+00 2.24625E-01 2.76739E+00 2.73529E+01 3.62332E+01 + 87.5 1.16390E-01 5.99789E+00 3.94139E-01 7.58851E-01 2.55025E+00 2.60689E-01 3.85082E+00 4.56107E+02 4.70036E+02 + 92.5 7.64882E-02 7.45080E+00 4.64205E-01 7.97665E-01 2.05658E-01 2.53237E-01 1.13659E+00 2.15134E+01 3.18981E+01 + 97.5 9.73889E-02 9.41742E+00 5.71892E-01 9.00254E-01 0.00000E+00 3.38890E-01 7.25912E-01 2.17862E+01 3.38380E+01 + 102.5 3.68821E-01 1.17983E+01 9.22857E-01 9.90539E-01 2.05658E-01 2.73332E-01 3.24630E+00 1.90768E+01 3.68827E+01 + 107.5 1.40866E-01 1.49763E+01 9.77552E-01 1.03381E+00 2.05658E-01 3.02800E-01 9.29930E-01 5.04390E+00 2.36108E+01 + 112.5 1.92456E-01 1.78278E+01 1.28612E+00 1.03320E+00 2.05658E-01 2.86201E-01 0.00000E+00 2.24646E+00 2.30779E+01 + 117.5 2.54597E-01 2.06815E+01 1.32965E+00 9.94290E-01 0.00000E+00 3.03672E-01 1.03238E-01 1.47686E+01 3.84355E+01 + 122.5 3.16846E-01 2.41082E+01 1.65376E+00 1.00318E+00 0.00000E+00 2.89085E-01 2.13376E+00 1.36871E+02 1.66376E+02 + 127.5 5.05056E-01 2.79339E+01 1.95268E+00 1.05111E+00 0.00000E+00 3.46045E-01 2.49979E+00 1.69874E+01 5.12760E+01 + 132.5 4.80308E-01 3.21630E+01 2.22178E+00 1.20921E+00 0.00000E+00 4.69702E-01 3.37245E+00 3.81901E+01 7.81066E+01 + 137.5 4.71791E-01 3.81005E+01 2.56140E+00 1.34421E+00 8.12028E-01 6.21499E-01 2.46233E+00 2.80344E+01 7.44081E+01 + 142.5 6.80561E-01 4.39911E+01 3.28659E+00 1.46366E+00 9.50409E-01 6.92449E-01 1.20701E+00 2.68285E+01 7.91003E+01 + 147.5 6.56937E-01 5.00238E+01 3.58919E+00 1.58529E+00 0.00000E+00 5.68462E-01 1.63930E+00 4.35650E+01 1.01628E+02 + 152.5 8.03675E-01 5.69622E+01 4.16398E+00 1.64509E+00 1.25627E+00 5.53992E-01 1.75478E+00 2.22948E+01 8.94347E+01 + 157.5 1.02913E+00 6.41523E+01 4.79832E+00 1.68259E+00 0.00000E+00 4.78287E-01 6.11008E-02 1.71154E+01 8.93171E+01 + 162.5 1.01346E+00 7.24096E+01 5.40756E+00 1.73576E+00 4.11317E-01 4.05437E-01 3.61370E-01 1.56466E+01 9.73911E+01 + 167.5 1.54845E+00 8.15811E+01 6.19343E+00 1.77824E+00 2.05658E-01 4.07876E-01 5.26654E+00 1.87334E+01 1.15715E+02 + 172.5 1.32405E+00 9.08761E+01 7.26631E+00 1.76703E+00 0.00000E+00 3.89630E-01 1.23387E+00 2.23269E+00 1.05090E+02 + 177.5 1.53823E+00 9.87941E+01 7.97271E+00 1.82106E+00 1.07784E+00 5.49313E-01 3.19406E+00 2.24375E+02 3.39322E+02 + 182.5 1.73275E+00 1.07697E+02 8.71580E+00 1.89504E+00 2.16504E+00 4.92700E-01 1.99691E+00 4.19643E+00 1.28892E+02 + 187.5 1.89887E+00 1.18843E+02 9.44743E+00 1.91419E+00 0.00000E+00 4.92825E-01 1.52137E+00 1.25557E+01 1.46674E+02 + 192.5 2.14077E+00 1.29143E+02 1.06527E+01 2.01942E+00 2.05658E-01 5.32759E-01 1.00000E-01 1.83069E+01 1.63101E+02 + 197.5 2.28090E+00 1.40043E+02 1.16612E+01 2.09703E+00 7.33634E-01 4.71836E-01 1.93913E+00 1.42947E+01 1.73522E+02 + 202.5 2.64685E+00 1.51099E+02 1.26715E+01 2.10999E+00 4.11317E-01 5.87528E-01 4.38079E-01 1.56988E+01 1.85663E+02 + 207.5 2.84005E+00 1.61843E+02 1.34867E+01 2.19835E+00 3.75305E+00 6.25004E-01 5.52737E+00 8.37937E+01 2.74067E+02 + 212.5 3.25595E+00 1.74900E+02 1.50308E+01 2.22864E+00 7.49839E-01 5.81734E-01 2.81961E+00 3.45840E+01 2.34150E+02 + 217.5 3.46113E+00 1.88150E+02 1.62797E+01 2.31385E+00 4.82384E+00 6.20967E-01 6.32525E+00 3.64941E+01 2.58469E+02 + 222.5 3.71532E+00 1.99485E+02 1.72400E+01 2.43326E+00 1.71116E+00 6.32218E-01 4.82299E+00 3.22340E+01 2.62274E+02 + 227.5 3.93657E+00 2.12862E+02 1.86215E+01 2.51943E+00 2.05658E-01 6.81550E-01 2.98408E+00 5.10532E+01 2.92864E+02 + 232.5 4.24983E+00 2.25983E+02 2.01806E+01 2.60538E+00 8.16883E-01 6.87444E-01 7.21676E+00 2.28266E+01 2.84567E+02 + 237.5 4.59570E+00 2.39176E+02 2.16407E+01 2.69698E+00 1.67826E+00 9.31593E-01 4.99406E+00 2.10503E+01 2.96764E+02 + 242.5 4.94950E+00 2.52290E+02 2.31471E+01 2.79643E+00 4.11317E-01 8.41274E-01 3.11881E+00 1.85433E+01 3.06098E+02 + 247.5 5.26446E+00 2.65745E+02 2.43967E+01 2.85884E+00 1.19951E+00 8.48640E-01 3.39328E+00 9.27690E+00 3.12983E+02 + 252.5 5.51137E+00 2.79382E+02 2.56902E+01 2.90319E+00 2.05030E+00 1.35524E+00 6.54964E+00 2.37694E+01 3.47212E+02 + 257.5 5.99515E+00 2.93242E+02 2.73615E+01 3.06371E+00 3.86368E+00 1.60997E+00 1.00977E+01 5.58533E+01 4.01087E+02 + 262.5 6.21434E+00 3.08338E+02 2.89447E+01 3.29544E+00 3.13735E+00 1.46229E+00 6.01357E+00 3.02641E+01 3.87670E+02 + 267.5 6.65424E+00 3.22781E+02 3.03692E+01 3.41100E+00 6.16975E-01 1.32594E+00 6.28516E+00 1.28466E+01 3.84290E+02 + 272.5 7.23326E+00 3.38755E+02 3.18652E+01 3.54230E+00 9.46604E-01 1.38334E+00 8.22455E+00 6.72038E+01 4.59155E+02 + 277.5 7.89736E+00 3.53513E+02 3.40571E+01 3.61569E+00 2.90561E+00 1.32980E+00 1.13664E+01 4.15706E+01 4.56256E+02 + 282.5 8.34614E+00 3.66640E+02 3.53419E+01 3.73099E+00 1.98306E+00 1.28057E+00 6.55080E+00 5.01229E+01 4.73996E+02 + 287.5 8.10537E+00 3.83148E+02 3.64438E+01 3.88807E+00 3.74016E+00 1.94887E+00 1.15334E+01 5.92274E+01 5.08035E+02 + 292.5 8.81514E+00 3.96686E+02 3.90005E+01 4.25272E+00 4.00315E+00 2.50940E+00 8.69602E+00 5.45696E+01 5.18532E+02 + 297.5 8.96015E+00 4.08280E+02 4.02917E+01 4.52252E+00 3.57998E+00 2.24959E+00 3.31821E+00 1.58504E+02 6.29707E+02 + 302.5 9.71633E+00 4.22618E+02 4.17213E+01 4.77674E+00 4.49329E+00 2.15379E+00 1.10035E+01 4.84946E+01 5.44978E+02 + 307.5 1.02386E+01 4.35918E+02 4.33268E+01 4.93324E+00 5.13629E+00 2.22004E+00 1.32672E+01 5.01161E+01 5.65156E+02 + 312.5 1.04223E+01 4.50463E+02 4.54609E+01 5.15847E+00 1.94942E+00 2.04243E+00 8.52077E+00 8.19247E+01 6.05942E+02 + 317.5 1.07891E+01 4.62915E+02 4.64979E+01 5.23454E+00 4.27902E+00 1.90572E+00 6.47562E+00 3.53354E+01 5.73433E+02 + 322.5 1.09427E+01 4.77330E+02 4.91935E+01 5.26494E+00 2.19482E+00 1.88124E+00 6.41893E+00 8.55603E+01 6.38786E+02 + 327.5 1.14908E+01 4.88261E+02 5.03717E+01 5.36671E+00 5.91026E+00 2.20568E+00 1.50136E+01 7.29965E+01 6.51616E+02 + 332.5 1.19361E+01 5.00139E+02 5.13089E+01 5.67699E+00 3.82652E+00 2.37009E+00 1.23984E+01 1.09240E+03 1.68006E+03 + 337.5 1.20681E+01 5.14010E+02 5.28514E+01 6.05405E+00 1.70398E+00 2.33798E+00 1.05093E+01 9.24720E+01 6.92007E+02 + 342.5 1.23490E+01 5.21510E+02 5.42620E+01 6.37794E+00 4.11317E-01 2.19005E+00 1.11503E+01 6.79179E+01 6.76168E+02 + 347.5 1.27509E+01 5.35571E+02 5.55184E+01 6.58173E+00 2.01311E+00 1.96299E+00 1.12862E+01 6.07252E+01 6.86409E+02 + 352.5 1.33077E+01 5.46371E+02 5.76060E+01 6.78823E+00 2.89323E+00 1.96006E+00 1.31853E+01 5.43796E+01 6.96492E+02 + 357.5 1.36759E+01 5.56190E+02 5.84729E+01 6.94197E+00 4.42343E+00 1.98151E+00 1.76543E+01 3.19560E+01 6.91296E+02 + 362.5 1.42781E+01 5.69658E+02 6.05033E+01 7.14192E+00 5.86906E+00 2.17718E+00 9.96818E+00 3.41134E+02 1.01073E+03 + 367.5 1.45381E+01 5.80766E+02 6.12914E+01 7.39646E+00 5.47934E+00 2.17521E+00 2.03257E+01 1.82855E+02 8.74827E+02 + 372.5 1.49540E+01 5.89510E+02 6.24846E+01 7.58861E+00 5.97579E+00 2.15511E+00 1.72258E+01 1.13243E+02 8.13136E+02 + 377.5 1.53281E+01 6.02908E+02 6.45094E+01 7.79686E+00 1.84986E+00 2.19418E+00 1.06875E+01 1.39893E+02 8.45167E+02 + 382.5 1.54423E+01 6.11828E+02 6.61635E+01 7.96229E+00 3.11386E+00 2.21687E+00 1.84188E+01 1.86638E+02 9.11784E+02 + 387.5 1.58073E+01 6.22309E+02 6.69250E+01 8.24960E+00 5.54442E+00 2.43159E+00 1.61610E+01 3.86101E+01 7.76038E+02 + 392.5 1.60969E+01 6.31163E+02 6.86152E+01 8.47054E+00 3.54836E+00 2.29998E+00 1.57432E+01 4.20241E+01 7.87962E+02 + 397.5 1.65813E+01 6.36218E+02 6.94217E+01 8.62327E+00 1.13086E+00 2.11498E+00 6.68408E+00 4.21503E+01 7.82925E+02 + 402.5 1.65367E+01 6.42913E+02 7.10363E+01 8.75430E+00 1.75166E+00 2.03373E+00 1.90242E+01 4.81706E+01 8.10221E+02 + 407.5 1.67342E+01 6.48664E+02 7.14222E+01 8.90130E+00 4.11317E-01 1.92641E+00 1.10653E+01 6.01248E+01 8.19250E+02 + 412.5 1.71698E+01 6.52550E+02 7.17759E+01 8.99922E+00 8.22633E-01 1.79239E+00 1.01201E+01 4.42900E+01 8.07520E+02 + 417.5 1.76193E+01 6.57890E+02 7.39355E+01 9.04688E+00 2.61737E+00 1.78737E+00 1.15231E+01 4.35126E+01 8.17932E+02 + 422.5 1.78971E+01 6.58783E+02 7.45428E+01 9.11849E+00 0.00000E+00 1.70407E+00 7.30709E+00 4.99614E+01 8.19314E+02 + 427.5 1.80201E+01 6.60273E+02 7.49222E+01 9.21605E+00 2.91893E+00 1.63279E+00 8.43587E+00 5.13630E+01 8.26782E+02 + 432.5 1.81311E+01 6.59898E+02 7.50286E+01 9.32577E+00 3.58524E+00 1.57313E+00 7.71508E+00 2.56450E+01 8.00902E+02 + 437.5 1.84232E+01 6.59360E+02 7.50116E+01 9.41469E+00 2.28202E+00 1.62042E+00 1.54203E+01 4.89755E+01 8.30507E+02 + 442.5 1.84195E+01 6.56836E+02 7.58088E+01 9.46982E+00 3.38542E+00 1.61951E+00 6.53277E+00 3.46057E+01 8.06677E+02 + 447.5 1.83569E+01 6.56214E+02 7.57237E+01 9.54516E+00 3.13290E+00 1.58031E+00 6.85186E+00 1.91769E+01 7.90582E+02 + 452.5 1.82706E+01 6.58492E+02 7.61998E+01 9.52308E+00 2.01469E+00 1.70665E+00 6.03783E+00 5.44051E+01 8.26650E+02 + 457.5 1.91778E+01 6.58453E+02 7.61906E+01 9.64309E+00 3.16596E+00 1.89007E+00 1.51820E+01 1.27697E+02 9.11400E+02 + 462.5 1.87779E+01 6.52226E+02 7.72842E+01 9.73954E+00 2.59408E+00 2.10854E+00 1.93226E+01 7.41268E+01 8.56179E+02 + 467.5 1.86804E+01 6.48703E+02 7.60035E+01 9.90114E+00 5.32531E+00 1.99694E+00 1.10320E+01 4.07741E+01 8.12416E+02 + 472.5 1.86511E+01 6.47465E+02 7.69976E+01 1.00372E+01 4.14584E+00 1.83096E+00 9.60202E+00 6.54727E+01 8.34202E+02 + 477.5 1.94192E+01 6.42556E+02 7.67936E+01 1.00990E+01 4.05719E+00 1.76307E+00 1.68761E+01 7.31050E+01 8.44669E+02 + 482.5 1.85120E+01 6.38182E+02 7.65831E+01 1.01837E+01 5.73622E+00 1.72238E+00 1.61719E+01 2.67816E+01 7.93873E+02 + 487.5 1.87245E+01 6.31091E+02 7.57856E+01 1.01494E+01 2.47112E+00 1.73226E+00 1.10861E+01 1.84643E+01 7.69505E+02 + 492.5 1.86552E+01 6.27893E+02 7.58075E+01 1.01077E+01 6.16975E-01 1.65405E+00 1.85135E+01 4.59820E+01 7.99230E+02 + 497.5 1.87774E+01 6.20884E+02 7.54106E+01 1.02234E+01 3.12423E+00 1.83964E+00 1.11110E+01 3.66094E+01 7.77980E+02 + 502.5 1.85128E+01 6.11207E+02 7.47867E+01 1.03467E+01 2.17399E+00 2.03159E+00 8.68880E+00 2.13229E+01 7.49070E+02 + 507.5 1.85188E+01 6.04683E+02 7.42417E+01 1.04700E+01 2.04328E+00 2.09356E+00 1.36051E+01 5.45491E+01 7.80204E+02 + 512.5 1.82513E+01 5.95617E+02 7.41759E+01 1.05902E+01 4.90812E+00 2.39934E+00 1.34593E+01 5.50473E+01 7.74448E+02 + 517.5 1.79895E+01 5.91412E+02 7.26130E+01 1.07966E+01 3.58067E+00 2.19824E+00 1.38811E+01 7.03155E+01 7.82787E+02 + 522.5 1.81884E+01 5.86370E+02 7.26538E+01 1.08876E+01 5.47750E+00 2.07366E+00 2.11933E+01 4.20529E+01 7.58897E+02 + 527.5 1.80445E+01 5.75399E+02 7.11525E+01 1.09857E+01 2.05658E-01 2.01075E+00 1.36926E+01 4.52190E+01 7.36710E+02 + 532.5 1.80037E+01 5.71147E+02 7.06867E+01 1.10910E+01 4.75520E+00 2.00312E+00 1.38574E+01 9.50279E+01 7.86572E+02 + 537.5 1.74675E+01 5.63508E+02 7.04883E+01 1.12028E+01 1.87713E+00 1.96527E+00 1.62523E+01 4.65142E+01 7.29275E+02 + 542.5 1.76318E+01 5.56665E+02 7.00343E+01 1.12976E+01 3.06699E+00 1.86572E+00 1.71287E+01 2.59445E+01 7.03634E+02 + 547.5 1.75820E+01 5.50320E+02 6.97971E+01 1.13400E+01 5.98333E+00 1.79075E+00 1.65174E+01 4.43756E+01 7.17706E+02 + 552.5 1.69980E+01 5.43670E+02 6.83412E+01 1.14467E+01 3.38229E+00 1.75224E+00 9.82047E+00 8.47597E+01 7.40171E+02 + 557.5 1.77408E+01 5.37062E+02 6.78747E+01 1.15815E+01 2.05658E+00 1.69933E+00 1.25363E+01 1.05334E+02 7.55885E+02 + 562.5 1.67779E+01 5.27634E+02 6.66589E+01 1.16853E+01 2.05658E+00 1.69919E+00 1.01797E+01 1.17431E+02 7.54122E+02 + 567.5 1.65009E+01 5.19049E+02 6.56601E+01 1.16734E+01 2.05658E-01 1.62483E+00 1.40535E+01 3.50797E+01 6.63847E+02 + 572.5 1.66474E+01 5.13638E+02 6.54957E+01 1.16120E+01 4.45750E+00 1.68120E+00 1.29364E+01 2.03936E+01 6.46862E+02 + 577.5 1.61709E+01 5.06697E+02 6.45733E+01 1.16869E+01 4.21224E+00 1.72088E+00 1.63802E+01 3.41220E+01 6.55564E+02 + 582.5 1.70010E+01 4.98513E+02 6.42327E+01 1.17903E+01 4.49004E+00 1.79659E+00 1.84821E+01 6.47650E+01 6.81071E+02 + 587.5 1.58106E+01 4.90397E+02 6.26877E+01 1.18146E+01 1.91449E+00 1.70309E+00 9.53193E+00 4.86846E+01 6.42544E+02 + 592.5 1.56098E+01 4.82456E+02 6.21518E+01 1.18720E+01 5.32210E+00 1.78033E+00 8.85252E+00 1.72308E+01 6.05276E+02 + 597.5 1.51798E+01 4.75370E+02 6.10967E+01 1.19059E+01 2.73383E+00 1.71041E+00 9.35274E+00 2.13497E+01 5.98699E+02 + 602.5 1.52845E+01 4.67688E+02 6.00378E+01 1.19019E+01 1.75731E+00 1.69423E+00 1.26198E+01 2.79371E+01 5.98920E+02 + 607.5 1.50077E+01 4.58955E+02 5.95365E+01 1.19498E+01 2.11563E+00 1.63413E+00 1.62137E+01 2.51185E+01 5.90531E+02 + 612.5 1.49107E+01 4.50488E+02 5.82921E+01 1.19951E+01 2.77916E+00 1.55292E+00 1.28430E+01 7.64498E+01 6.29311E+02 + 617.5 1.46426E+01 4.44373E+02 5.77405E+01 1.19959E+01 2.84667E+00 1.64546E+00 1.34922E+01 4.52778E+01 5.92014E+02 + 622.5 1.42687E+01 4.35792E+02 5.65965E+01 1.21639E+01 5.37358E+00 1.72146E+00 1.19911E+01 2.47447E+01 5.62652E+02 + 627.5 1.41990E+01 4.27824E+02 5.57890E+01 1.22265E+01 1.08962E+00 1.65390E+00 8.81019E+00 2.42487E+01 5.45841E+02 + 632.5 1.37383E+01 4.19243E+02 5.50295E+01 1.22403E+01 1.83835E+00 1.59657E+00 9.22957E+00 3.61944E+01 5.49110E+02 + 637.5 1.36281E+01 4.12136E+02 5.35837E+01 1.23066E+01 7.15870E+00 1.63356E+00 2.06304E+01 4.42392E+01 5.65316E+02 + 642.5 1.32267E+01 4.04381E+02 5.29472E+01 1.23324E+01 1.18319E+00 1.51609E+00 1.72054E+01 5.21876E+01 5.54980E+02 + 647.5 1.33490E+01 3.99371E+02 5.18436E+01 1.23369E+01 2.86638E+00 1.49006E+00 7.71399E+00 4.87304E+01 5.37702E+02 + 652.5 1.32825E+01 3.93253E+02 5.15025E+01 1.24064E+01 2.91023E+00 1.53461E+00 1.02694E+01 4.84323E+01 5.33591E+02 + 657.5 1.27242E+01 3.84676E+02 4.99399E+01 1.24545E+01 2.81156E+00 1.43674E+00 7.29351E+00 2.76858E+01 4.99022E+02 + 662.5 1.25564E+01 3.79267E+02 5.02025E+01 1.24533E+01 8.81537E-01 1.43309E+00 1.37031E+01 2.17055E+01 4.92203E+02 + 667.5 1.25272E+01 3.70921E+02 4.85681E+01 1.25127E+01 1.60102E+00 1.43617E+00 5.50175E+00 1.90565E+01 4.72124E+02 + 672.5 1.22774E+01 3.65880E+02 4.78624E+01 1.25474E+01 1.48356E+00 1.39963E+00 1.42438E+01 2.13385E+01 4.77033E+02 + 677.5 1.23601E+01 3.61097E+02 4.76481E+01 1.25843E+01 9.56924E-01 1.38682E+00 8.56386E+00 2.49866E+01 4.69583E+02 + 682.5 1.17525E+01 3.54583E+02 4.63387E+01 1.26520E+01 8.22633E-01 1.26333E+00 2.34468E+00 1.01581E+01 4.39915E+02 + 687.5 1.15617E+01 3.48989E+02 4.54004E+01 1.26458E+01 3.66096E+00 1.32373E+00 1.16655E+01 3.93507E+01 4.74598E+02 + 692.5 1.16167E+01 3.43981E+02 4.48284E+01 1.27456E+01 1.86714E+00 1.32601E+00 8.73279E+00 5.56405E+01 4.80738E+02 + 697.5 1.12462E+01 3.36633E+02 4.44527E+01 1.28182E+01 1.73812E+00 1.34326E+00 1.10964E+01 3.69912E+01 4.56319E+02 + 702.5 1.09370E+01 3.28031E+02 4.33183E+01 1.28798E+01 3.92491E+00 1.38835E+00 5.88493E+00 1.16183E+01 4.17982E+02 + 707.5 1.08008E+01 3.20977E+02 4.21680E+01 1.28815E+01 8.22633E-01 1.35587E+00 8.21049E+00 1.68206E+01 4.14037E+02 + 712.5 1.04037E+01 3.17452E+02 4.14850E+01 1.29183E+01 2.34281E+00 1.32175E+00 5.20540E+00 2.07675E+01 4.11896E+02 + 717.5 1.04719E+01 3.10538E+02 4.14020E+01 1.29743E+01 1.24569E+00 1.27138E+00 5.78744E+00 3.71620E+01 4.20853E+02 + 722.5 1.00289E+01 3.04375E+02 3.99868E+01 1.30113E+01 0.00000E+00 1.25784E+00 6.57942E+00 4.30035E+01 4.18243E+02 + 727.5 9.96763E+00 3.00594E+02 3.99395E+01 1.30823E+01 1.49392E+00 1.24204E+00 6.57705E+00 2.20991E+01 3.94996E+02 + 732.5 1.01246E+01 2.93028E+02 3.93089E+01 1.31933E+01 1.90664E+00 1.12506E+00 7.09142E+00 1.76734E+01 3.83451E+02 + 737.5 9.76342E+00 2.86240E+02 3.81188E+01 1.32266E+01 1.02829E+00 1.11305E+00 1.16594E+01 1.98411E+01 3.80991E+02 + 742.5 9.44027E+00 2.83298E+02 3.75235E+01 1.31902E+01 1.85268E+00 1.09382E+00 1.24093E+01 1.93104E+01 3.78118E+02 + 747.5 9.42312E+00 2.76179E+02 3.69466E+01 1.31384E+01 3.03756E+00 1.07565E+00 7.88019E+00 1.59411E+01 3.63621E+02 + 752.5 9.12935E+00 2.69915E+02 3.58135E+01 1.30860E+01 1.13585E+00 9.76427E-01 4.74791E+00 1.22165E+01 3.47021E+02 + 757.5 9.12138E+00 2.65738E+02 3.54129E+01 1.30239E+01 1.02829E+00 9.10901E-01 4.97497E+00 1.95671E+01 3.49778E+02 + 762.5 8.92617E+00 2.60265E+02 3.46279E+01 1.30299E+01 6.16975E-01 8.28400E-01 3.28244E+00 2.03845E+01 3.41961E+02 + 767.5 9.02910E+00 2.55900E+02 3.40263E+01 1.30109E+01 1.68013E+00 8.25708E-01 1.19779E+01 1.96813E+01 3.46131E+02 + 772.5 8.59373E+00 2.50869E+02 3.38139E+01 1.29368E+01 2.05658E-01 8.75975E-01 1.15032E+01 5.21066E+00 3.24009E+02 + 777.5 8.47211E+00 2.44259E+02 3.30716E+01 1.29050E+01 1.54291E+00 9.74810E-01 5.48368E+00 1.80399E+01 3.24750E+02 + 782.5 8.04423E+00 2.37871E+02 3.15979E+01 1.28922E+01 1.43961E+00 8.60998E-01 6.73828E+00 1.92420E+01 3.18686E+02 + 787.5 7.84425E+00 2.33191E+02 3.11706E+01 1.28935E+01 8.22633E-01 8.06932E-01 5.39425E+00 2.57440E+01 3.17868E+02 + 792.5 8.00745E+00 2.30632E+02 3.06578E+01 1.28111E+01 1.57286E+00 7.55721E-01 2.21518E+00 1.10303E+01 2.97682E+02 + 797.5 7.65707E+00 2.27379E+02 3.04077E+01 1.27005E+01 8.22633E-01 7.78043E-01 4.34110E+00 1.44198E+01 2.98506E+02 + 802.5 7.55632E+00 2.21465E+02 3.01760E+01 1.27059E+01 3.52177E+00 7.53319E-01 4.05708E+00 1.60606E+01 2.96296E+02 + 807.5 7.33773E+00 2.15829E+02 2.90355E+01 1.26835E+01 6.16975E-01 7.00094E-01 3.27212E+00 7.58684E+00 2.77062E+02 + 812.5 7.58748E+00 2.11178E+02 2.83137E+01 1.26170E+01 2.93510E+00 7.09702E-01 2.23541E+00 3.47184E+00 2.69048E+02 + 817.5 7.13388E+00 2.05619E+02 2.77154E+01 1.26119E+01 2.10215E+00 7.20056E-01 1.50308E+00 4.78851E+00 2.62194E+02 + 822.5 7.03108E+00 2.01687E+02 2.72731E+01 1.25963E+01 1.02829E+00 6.36013E-01 4.34315E+00 1.13325E+01 2.65927E+02 + 827.5 6.90882E+00 1.94633E+02 2.62117E+01 1.25318E+01 1.68201E+00 6.39769E-01 6.60553E+00 9.69656E+00 2.58909E+02 + 832.5 6.79061E+00 1.90838E+02 2.60359E+01 1.24388E+01 4.11317E-01 6.10251E-01 5.08183E+00 1.24584E+01 2.54666E+02 + 837.5 7.04327E+00 1.85429E+02 2.54800E+01 1.23966E+01 1.23395E+00 7.39280E-01 3.13809E+00 8.78160E+00 2.44242E+02 + 842.5 6.35597E+00 1.80474E+02 2.46156E+01 1.23850E+01 8.22633E-01 7.84348E-01 5.72595E+00 1.52183E+01 2.46382E+02 + 847.5 6.24776E+00 1.77189E+02 2.41897E+01 1.23321E+01 4.11317E-01 7.98322E-01 1.05427E+00 1.84405E+01 2.40663E+02 + 852.5 5.92837E+00 1.72342E+02 2.37005E+01 1.22703E+01 1.64527E+00 7.82842E-01 4.11603E+00 6.91941E+00 2.27705E+02 + 857.5 5.81758E+00 1.68152E+02 2.29753E+01 1.22808E+01 2.45537E+00 7.48350E-01 4.39854E+00 1.49683E+01 2.31796E+02 + 862.5 5.87099E+00 1.62153E+02 2.25716E+01 1.21912E+01 1.84705E+00 7.56121E-01 4.50077E+00 1.06242E+01 2.20515E+02 + 867.5 5.67463E+00 1.55558E+02 2.14867E+01 1.21582E+01 1.41731E+00 7.83431E-01 4.49886E+00 3.88465E+01 2.40424E+02 + 872.5 5.40841E+00 1.51616E+02 2.11584E+01 1.21849E+01 1.69394E+00 7.21190E-01 7.84411E+00 1.75074E+01 2.18134E+02 + 877.5 5.33795E+00 1.48213E+02 2.07444E+01 1.21714E+01 2.01315E+00 6.59766E-01 3.01183E+00 7.76581E+00 1.99917E+02 + 882.5 5.21301E+00 1.44838E+02 2.02883E+01 1.21323E+01 4.11317E-01 6.96657E-01 2.13598E+00 2.31380E+01 2.08854E+02 + 887.5 5.07468E+00 1.41052E+02 1.95552E+01 1.20833E+01 1.35277E+00 6.45121E-01 2.64440E+00 2.94650E+01 2.11873E+02 + 892.5 4.97497E+00 1.35872E+02 1.92867E+01 1.20214E+01 6.16975E-01 6.27911E-01 2.62215E+00 1.74470E+01 1.93469E+02 + 897.5 4.77209E+00 1.32002E+02 1.87412E+01 1.20255E+01 1.85093E+00 6.05299E-01 5.49457E+00 5.98372E+00 1.81475E+02 + 902.5 4.59793E+00 1.27346E+02 1.78058E+01 1.19957E+01 1.53463E+00 5.92630E-01 4.45948E+00 1.14090E+01 1.79741E+02 + 907.5 5.18042E+00 1.22722E+02 1.72241E+01 1.18983E+01 2.05658E-01 5.82700E-01 2.70246E+00 5.84254E+00 1.66359E+02 + 912.5 4.26527E+00 1.19282E+02 1.71189E+01 1.18590E+01 0.00000E+00 5.94699E-01 4.21328E+00 3.09199E+01 1.88253E+02 + 917.5 4.25837E+00 1.16959E+02 1.66237E+01 1.18565E+01 2.22220E+00 5.85033E-01 3.75394E+00 5.33899E+00 1.61597E+02 + 922.5 4.20854E+00 1.14379E+02 1.63301E+01 1.18960E+01 2.72871E+00 4.69495E-01 5.87566E+00 9.07508E+00 1.64963E+02 + 927.5 4.02279E+00 1.11427E+02 1.54898E+01 1.18437E+01 8.22633E-01 4.22265E-01 2.91595E+00 1.19855E+01 1.58930E+02 + 932.5 3.95224E+00 1.08972E+02 1.51313E+01 1.17710E+01 4.11317E-01 4.22079E-01 5.60698E+00 6.77351E+00 1.53040E+02 + 937.5 3.83145E+00 1.04917E+02 1.48773E+01 1.17300E+01 1.02829E+00 5.05313E-01 4.16801E+00 5.09910E+00 1.46157E+02 + 942.5 3.71627E+00 1.01858E+02 1.44962E+01 1.16937E+01 4.11317E-01 4.63286E-01 2.03016E+00 6.05241E+00 1.40721E+02 + 947.5 4.17745E+00 9.87293E+01 1.41005E+01 1.16418E+01 3.39734E+00 4.20559E-01 3.52910E+00 1.36640E+01 1.49660E+02 + 952.5 3.43835E+00 9.56864E+01 1.37363E+01 1.16098E+01 6.16975E-01 4.08199E-01 4.86770E+00 9.03454E-01 1.31267E+02 + 957.5 3.39878E+00 9.21055E+01 1.32503E+01 1.15578E+01 4.11317E-01 4.15283E-01 2.73436E+00 9.85862E+00 1.33732E+02 + 962.5 3.41625E+00 8.95771E+01 1.30236E+01 1.15057E+01 6.16975E-01 4.54100E-01 8.59265E-01 4.78435E+00 1.24237E+02 + 967.5 3.26916E+00 8.63205E+01 1.26232E+01 1.14523E+01 0.00000E+00 4.01167E-01 2.96301E-01 1.78596E+00 1.16148E+02 + 972.5 3.23102E+00 8.28524E+01 1.20088E+01 1.14270E+01 2.05658E-01 4.07884E-01 4.97574E+00 6.41128E+00 1.21520E+02 + 977.5 2.99712E+00 7.98489E+01 1.15471E+01 1.14298E+01 8.22633E-01 4.16749E-01 8.75861E-01 4.77784E+00 1.12716E+02 + 982.5 2.90789E+00 7.78006E+01 1.10307E+01 1.14413E+01 4.11317E-01 4.35002E-01 1.79759E+00 1.18712E+01 1.17696E+02 + 987.5 2.83483E+00 7.57515E+01 1.09321E+01 1.13753E+01 8.22633E-01 3.36511E-01 1.53167E-01 2.05925E+00 1.04265E+02 + 992.5 2.77574E+00 7.30690E+01 1.07884E+01 1.13722E+01 1.65287E+00 3.65523E-01 3.00760E+00 4.54579E+00 1.07577E+02 + 997.5 2.68584E+00 7.15856E+01 1.01075E+01 1.13469E+01 6.16975E-01 3.26858E-01 3.76981E+00 9.42287E+00 1.09862E+02 + 1002.5 2.56046E+00 7.04657E+01 1.01629E+01 1.12815E+01 4.11317E-01 3.04909E-01 1.20059E+00 3.84790E+00 1.00235E+02 + 1007.5 2.57656E+00 6.84069E+01 9.87757E+00 1.12417E+01 4.11317E-01 3.09770E-01 1.98070E+00 6.46648E+00 1.01271E+02 + 1012.5 2.46914E+00 6.53124E+01 9.50496E+00 1.11926E+01 8.06593E-01 3.15747E-01 6.20207E+00 7.95389E+00 1.03757E+02 + 1017.5 2.37216E+00 6.30719E+01 9.30901E+00 1.11502E+01 4.11317E-01 2.92220E-01 2.16995E+00 4.04391E+00 9.28206E+01 + 1022.5 2.47082E+00 6.05777E+01 8.99484E+00 1.11477E+01 4.11317E-01 2.66996E-01 1.28472E+00 6.78486E+00 9.19390E+01 + 1027.5 2.20163E+00 5.87250E+01 8.69179E+00 1.11170E+01 1.34876E+00 2.46473E-01 0.00000E+00 1.58832E+00 8.39189E+01 + 1032.5 2.13705E+00 5.58224E+01 8.28565E+00 1.11036E+01 1.73322E+00 2.25485E-01 3.71159E+00 1.99541E+01 1.02973E+02 + 1037.5 1.09515E+03 1.18845E+01 4.03263E+02 1.10962E+01 1.48077E+04 2.09360E-01 5.25335E+02 1.16406E+02 1.69710E+04 + FIT OF THE HILLAS CURVE N(T) = P1*((T-P2)/(P3-P2))**((P3-P2)/(P4+P5*T+P6*T**2)) * EXP((P3-T)/(P4+P5*T+P6*T**2)) + TO LONGITUDINAL DISTRIBUTION OF ALL CHARGED PARTICLES + PARAMETERS = 5.7722E+04 -3.9655E+01 4.4669E+02 5.9967E+01 2.3036E-02 7.4710E-06 + CHI**2/DOF = 2.1752E+01 + AV. DEVIATION IN % = 9.0539E+00 + diff --git a/tests/test_longitudinal.py b/tests/test_longitudinal.py index 6b5dc6f..944ae74 100644 --- a/tests/test_longitudinal.py +++ b/tests/test_longitudinal.py @@ -1,20 +1,34 @@ +import pytest import numpy as np -def test_longitudinal(): +@pytest.mark.parametrize("slant", (True, False)) +def test_longitudinal(slant): from corsikaio import read_longitudinal_distributions from corsikaio import longitudinal_fit_function - path = "./tests/resources/corsika_77500_particle.long" + + if slant: + path = "./tests/resources/corsika_77500_particle_slant.long" + n_steps = 105 + width = 10 + else: + path = "./tests/resources/corsika_77500_particle_vertical.long" + n_steps = 208 + width = 5 n_showers = 0 for longi in read_longitudinal_distributions(path): n_showers += 1 assert longi["shower"] == n_showers - assert len(longi["particles"]) == 105 - assert len(longi["energy_deposition"]) == 105 - np.testing.assert_array_equal(np.diff(longi["particles"]["depth"]), 10) - np.testing.assert_array_equal(np.diff(longi["energy_deposition"]["depth"]), 10) + assert longi["slant"] is slant + + assert longi["n_steps"] == n_steps + assert len(longi["particles"]) == n_steps + assert len(longi["energy_deposition"]) == n_steps + + np.testing.assert_array_equal(np.diff(longi["particles"]["depth"]), width) + np.testing.assert_array_equal(np.diff(longi["energy_deposition"]["depth"]), width) particles = longi["particles"] @@ -25,5 +39,7 @@ def test_longitudinal(): np.testing.assert_allclose(particles["charged"][mask], fit, rtol=0.2) - assert n_showers == 5 - + if slant: + assert n_showers == 5 + else: + assert n_showers == 1