[fluka-discuss]: Reading USERDUMP
Dear experts
I need to study the anti-proton annihilation in silicon in detail. I only have two primaries and dump the events to file by using USERDUMP. The binary file is read by the the rutine Data.py found in the flair folder. I need some help to understand the output I am getting.
I print the info to screen with a small python routine like this:
import Data
datafile = Data.Mgdraw("example001_TRAKFILE")
while(True):
evtType = datafile.readEvent()
if evtType == None:
break
elif evtType == 0:
print "tracking"
#tracking
print "\t", datafile.ntrack
print "\t", datafile.mtrack
print "\t", datafile.data
I have problems interpreting the output I get from datafile.data
Here are three examples of the tracking output I get:
tracking
1
0
(0.0, 0.0, -50.0, 0.0, 0.0, 0.0, 50.0)
tracking
1
1
(0.0, 0.0, 0.0, 0.0, 0.0, 3.337009547976777e-05, 1.3322263839654624e-05, 3.337009547976777e-05)
I have problems with the last line, where I print out datafile.data.
As far as I understand this is explained in the manual chapter 11, and it should be the second record. First of all I don't understand what the variable I is. I think the three first variables are the xtrack, ytrack and ztrack, but I don't understand how to interpret the 4 or 5 last entries. I hope someone is able to answer me. Below is the part of the manual that I am not able to connect with the output that I get.
Kind regards,
Helga
In Case 1, the following variables are written: First record:
NTRACK, MTRACK, JTRACK, ETRACK, WTRACK,
(three integers and two real variables)
Next record:
(XTRACK(I), YTRACK(I), ZTRACK(I), I = 0, NTRACK),
(DTRACK(J), J = 1, MTRACK), CTRACK
(NTRACK+MTRACK+1 real variables)
where:
NTRACK = number of trajectory segments
MTRACK = number of energy deposition events along the trajectory
JTRACK = particle type (see (5))
ETRACK = total energy of the particle (rest + kinetic)
WTRACK = particle weight
XTRACK(I), YTRACK(I), ZTRACK(I) = coordinates defining the upstream
end of the (I+1)th segment; for I = NTRACK, the end of the
trajectory
DTRACK(J) = energy deposition in the Jth deposition event along the
trajectory
CTRACK = total curved path
Received on Fri Aug 22 2014 - 13:34:43 CEST
This archive was generated by hypermail 2.3.0
: Fri Aug 22 2014 - 13:34:43 CEST