Hi John,
you could do something like this (of course it would require a loop for
j=NTRACK to 1 by -1),
but you will not get the desired incoming direction, but the one actually
performed due to the MCS.
Which I would assume would be almost-equal with first sub-step direction
DX = Xtrack(1) - Xtrack(0)
DY = Ytrack(1) - Ytrack(0)
DZ = Ztrack(1) - Ztrack(0)
and then get dx/sqrt(dx^2 + dy^2 + dz^2), ...
The correct incoming direction you will have it from the previous C[XYZ]TRCK variables,
which you have to save somewhere.
Vasilis Vlachoudis
Dep EN, CERN
CH-1211 GENEVA 23
SWITZERLAND
Phone: +41-22 767 9851
GSM: +41-76 487 4378
Fax: +41-22 766 9644
Received on Thu Jan 12 2012 - 17:22:53 CET
This archive was generated by hypermail 2.2.0 : Thu Jan 12 2012 - 17:23:24 CET