Re: [fluka-discuss]: Reading USERDUMP

From: Stefan Roesler <sroesler_at_mail.cern.ch>
Date: Fri, 22 Aug 2014 16:34:13 +0200

Dear Helga

Difficult to comment on your way of reading the data.

However, they are written by the user routine mgdraw, that is available in
the usermvax directory and which you can adapt to your needs. For example,
to be sure on the data you could change the WRITE statements to formatted
output for a few test runs and inspect with your favorite editor.

Furthermore, the manual contains in Chapter 12.2.6} a short Fortran code
for converting binary into ascii format that you could use to cross check
your own postprocessing routine.

Cheers
Stefan


On Fri, 22 Aug 2014, Helga Holmestad wrote:

> 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 - 17:47:27 CEST