This is very raw code. As of writing, i haven't even updated the repos it depends on yet!
This will be an approximation, using only Brehmstrahlung and pair production. Not entirely sure where exactly the approximation starts.
Currently it is probably pretty much wrong due to the soft photon spectrum not being correct. Nevertheless it does seem like for instance the muon does radiate less and loses less energy.(though i think it isn't 'less enough')
The scale of spatial sizes is not quite right yet either.
Soft photon spectrum thing.. Cutting off the low-energy portion because at some point the process is not possible, and figuring what to do with the high end.
Get spatial scales right
'Laziness' when tracks dont have to be stored.
Consistency check, particularly the energy one. Also figure stuff about rounding errors.