IAP GITLAB

ObservationPlanes show non-existing particles

This is a follow-up to the discussion on the slack development channel.

I've been simulating electromagnetic (1TeV initial electron) showers with the examples/corsika.cpp script. The only option I have changed is the emcut, which I changed to 20MeV.

If I look at the longitudinal profile (here for the charged particles), they look like this:

image

This means that I should expect zero particles on my observation plane. But if I look at the output from the ObservationPlane, I have 8931 entries in the parquet file! The energy distribution of these particles look like this:

image

Note that a lot of these particles have an energy very close to the cut energy, but not all of them.

For a 10 TeV shower, we would expect some charged particles according to the observation plane:

image

But the output of the observation plane shows 91676 particles!

This might, but doesn't have to be, related to issue #572 (closed).

Edited by Jean-Marco Alameddine
To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information