transform detector to surface coordinates in MSIS00IceCubeCentered

mceq-project / MCEq

Matrix cascade equation core code

BSD 3-Clause "New" or "Revised" License

30 stars 32 forks source link

transform detector to surface coordinates in MSIS00IceCubeCentered #8

Closed HansN87 closed 5 years ago

HansN87 commented 5 years ago

A discussion with Joeran Stettner at RWTH Aachen reminded me of this little patch from about year ago. This is a minor correction to the zenith angle as probed by IceCube. Relevant for very horizontal trajectories, where the exactly horizontal ones can not be probed.

afedynitch commented 5 years ago

Hi Hans,

thanks for this patch! I'm still a bit confused about what this transformation does. For me IceCube is not under the surface but rather above it and surface is the bedrock. The SouthPole ice suface above IceCube is at 3000m asl. It is important to make this difference for the muon fluxes and the atmosphere models that work in altitudes asl. Can you please clarify?

afedynitch commented 5 years ago

As we discussed, the geometry in MCEq does not require this transformation. IceCube would sit a few hundred meters above the bedrock = "the surface". The glacier surface above IceCube is at 3km altitude asl. in the coordinate system of MCEq.