How to compute energy losses with MMC

when you run Amanda with "-user" option, a user line is added to each event. A
definition will also be added near the beginning of the output. You need to
consult the f2k manual to understand the user line, but it might already be
obvious once you see the definition.

You have the Eini and Efin energies: these, when positive, are the energies of
the muon at points of entry and exit of the detector (defined as a cylinder
with -radius and -length options). When negative or zero, the muon did not
reach or go through these points.

Then, you also have the muon length that is recorded in the TR line for the
muon. Knowing where it started you can calculate where it stopped.

You also have a record of all secondaries (energy losses, let's call them E_i)
and where they happened. As an example, when both Eini and Efin are positive,
you can now estimate the rate of continuous energy losses as
contLossRate=(Eini-Efin-SUM(E_i for all i))/(muon track length between ini and
fin points). If one of Eini or Efin is not positive, you need to pay attention
to where the muon stopped (and/or started, if inside the detector volume) in
the above calculation.

Now, the muon energy at point p is calculated as

Eini-SUM(E_i for i before point p)-(distance between ini and p)*contLossRate.

You can then compute the differences between muon energies at different points,
e.g., to compute energy loss at a given depth increment.

You can also do this calculation with a (relatively new) option -recc, which
could simplify the tracking of the continuous energy losses (as it prints the
continuous energy losses after each muon propagation step between pairs of
secondaries).