CERN homepage INFN homepage
FLUKA: EVENTDAT Previous Index Next

EVENTDAT


    For calorimetry only.
    Prints event by event the scored star production and/or energy
    deposition in each region, and the total energy balance.

    See also EVENTBIN, SCORE

    EVENTDAT requests separate scoring by region of energy and/or
    star density for each event (primary history).
    The quantities to be scored are defined via a SCORE command (see
    SCORE for details).
    As for SCORE, a maximum per run of 4 different star densities is allowed.
    The EVENTDAT output includes also a detailed energy balance event by event.

     
WHAT(1)
= output unit If < 0, the output is unformatted. Values of |
WHAT(1)
| < 21 should be avoided (with the exception of +11).
Default
= 11 (standard output)
WHAT(2)
to
WHAT(6)
: not used
SDUM
= output file name (no default!) Max. 10 characters.
Default
(option EVENTDAT not given): no event by event scoring
Notes:
1) Unformatted data are written as follows. Once, at the beginning of the run: - RUNTIT, RUNTIM, NREGS, NSCO, (ISCORE(IS), IS = 1, NSCO) Then, for each primary particle: - NCASE, WEIPRU, ENETOT - (ENDIST(IE), IE = 1, 12) Then, NSCO times: - ISC, ISCORE(ISC) - (REGSCO(IR,ISC), IR = 1, NREGS) Then one dummy record (for historical reasons): - NDUM, DUM1, DUM2 Then: - ISEED1, ISEED2, SEED1, SEED2, SOPP1, SOPP2 where: RUNTIT = title of the run (CHARACTER*80 variable), which appears also at the beginning of the standard output RUNTIM = time of the run (CHARACTER*32 variable), which appears also at the beginning of the standard output NREGS = number of regions NSCO = number of scoring distributions requested by SCORE ISCORE(I) = Ith requested (generalised) particle distribution (see 5}) NCASE = number of primaries handled so far (current one included) WEIPRU = primary weight ENETOT = primary particle total energy (GeV) ENDIST(I) are 12 energy contributions to the total energy balance, some of which appear at the end of the standard output. Here they are given separately for each primary history (in GeV) and NOT normalised to the weight of the primary. Note that some of the contributions are meaningful only in specific contexts (e.g. if low-energy neutron transport has been requested). ENDIST(1) = energy deposited by ionisation ENDIST(2) = en. depos. by pi0, electrons, positrons and photons ENDIST(3) = en. depos. by nuclear recoils and heavy fragments ENDIST(4) = energy deposited by particles below threshold ENDIST(5) = energy leaving the system ENDIST(6) = energy carried by discarded particles ENDIST(7) = residual excitation energy after evaporation ENDIST(8) = energy deposited by low-energy neutrons (kerma, proton recoil energy not included) ENDIST(9) = energy of particles out of the time limit ENDIST(10) = energy lost in endothermic nuclear reactions (gained in exothermic reactions if < 0) above 20 MeV (not implemented yet) ENDIST(11) = energy lost in endothermic low-energy neutron reactions (gained in exothermic reactions if < 0) (not implemented yet) ENDIST(12) = missing energy NDUM, DUM1, DUM2 = three dummy variables, with no meaning REGSCO(IR,ISC) = energy or stars (corresponding to the ISCth generalised particle distribution) deposited or produced in the IRth region during the current primary history. NOT normalised, neither to the the primary weight nor to the region volume ISEED1, ISEED2, SEED1, SEED2, SOPP1, SOPP2 = random number generator information to be read in order to reproduce the current sequence (skipping calls, see RANDOMIZE). 2) All the above quantities are REAL*4, except RUNTIT and RUNTIM (which are of type CHARACTER) and those with a name beginning with I,J,K,L,M,N (which are integer). 3) The different items appearing in the EVENTDAT energy balance may sometimes give overlapping information and are not all meaningful in every circumstance (for instance residual excitation energy is meaningful only if gamma de-excitation has not been requested). Unlike the balance which is printed at the end of standard output, these terms are not additive. 4) An example on how to read EVENTDAT unformatted output is given below. PROGRAM RDEVDT CHARACTER*80 RUNTIT, FILNAM CHARACTER*32 RUNTIM DIMENSION ISCORE(4), ENDIST(12), REGSCO(5000,4) WRITE(*,*) 'Name of the EVENTDAT binary file?' READ(*,'(A)') FILNAM IB = INDEX(FILNAM,' ')i OPEN(UNIT = 7, FORM = 'UNFORMATTED', FILE = FILNAM(1:IB-1), & STATUS = 'OLD') OPEN(UNIT = 8, FORM = 'FORMATTED', FILE = FILNAM(1:IB-1)//'.txt', & STATUS = 'NEW') * Once, at the beginning of the run: READ(7) RUNTIT, RUNTIM, NREGS, NSCO, (ISCORE(IS), IS = 1, NSCO) WRITE(8,'(A80)') RUNTIT WRITE(8,'(A32)') RUNTIM WRITE(8,'(A,I6,5X,A,I4)') 'Number of regions: ', NREGS, & ' Number of scored quantities: ', NSCO WRITE(8,'(A,4I6)') 'The scored quantities are: ', & (ISCORE(IS), IS = 1, NSCO) * Loop on each primary particle: 100 CONTINUE WRITE(8,*) READ(7,END=300) NCASE, WEIPRU, ENETOT WRITE(8,'(A,I10,1P,2G12.4)') 'NCASE, WEIPRU, ENETOT: ', & NCASE, WEIPRU, ENETOT READ(7) (ENDIST(IE), IE = 1, 12) WRITE(8,'(A)') 'ENDIST: ' DO 400 IE = 1, 12, 2 WRITE(8,'(2(I5,5X,1P,G12.4))') IE,ENDIST(IE),IE+1,ENDIST(IE+1) 400 CONTINUE DO 200 ISC = 1, NSCO READ(7) IISC, ISCORE(ISC) * IISC is redundant, must be equal to ISC IF(IISC .NE. ISC) STOP 'Wrong sequence' WRITE(8,'(A,I2,A,I3,A)') & 'Quantity n. ',ISC, ' (',ISCORE(ISC),'):' READ(7) (REGSCO(IR,ISC), IR = 1, NREGS) WRITE(8,*) 'Scoring per region:' DO 500 IR = 1, NREGS WRITE(8,'(I7,3X,1P,G12.4)') IR, REGSCO(IR,ISC) 500 CONTINUE 200 CONTINUE READ(7) NDUM, DUM1, DUM2 IF (DUM1 .LT. 0.) THEN * DUM1 < 0 is used to signal that seeds follow READ(7) ISEED1, ISEED2, SEED1, SEED2, SOPP1, SOPP2 WRITE(8,*) ISEED1, ISEED2, SEED1, SEED2, SOPP1, SOPP2 ELSE BACKSPACE 7 END IF * This event is finished, start again with the next one GO TO 100 300 CONTINUE WRITE(8,*) "End of a run of ", NCASE, " particles" CLOSE (UNIT = 7) CLOSE (UNIT = 8) END Example (number-based):
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+...
SCORE 208. 211. 201. 8. 0. 0. EVENTDAT -23. 0. 0. 0. 0. 0. EVT.SCOR
* In this example, the user requests (with option SCORE) scoring of
* total and electromagnetic energy deposition, total stars and
* neutron-produced stars. The average scores for each region will be
* printed on standard output (as an effect of SCORE command), and
* corresponding scores, as well as the energy balance, will be written
* separately for each primary particle on an unformatted file EVT.SCOR
The same example, name-based: SCORE ENERGY EM-ENRGY ALL-PART NEUTRON 0. 0. EVENTDAT -23. 0. 0. 0. 0. 0. EVT.SCOR * Start_Devel_seq

Previous Index Next