[fluka-discuss]: Reading USERDUMP

From: Helga Holmestad <helga.holmestad_at_cern.ch>
Date: Fri, 22 Aug 2014 10:02:15 +0000

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