does example file 12cqmd.inp do nuclear fragmentation?

From: gin chen (gin_chen2005@yahoo.com)
Date: Tue Jan 10 2006 - 06:09:50 CET

  • Next message: Giuseppe Battistoni: "Announcemnt of a new FLUKA course"

    Hi,
    As a test, I tried to run the example file 12cqmd.inp
    given in fluka distribution with a change in the
    thickness of the rock. Then 20 GeV/n ( Total energy
    240 GeV ) Carbon-12 particle was sent to the rock see
    what happens. The output was written using EVENTDAT
    control card. This output was read by the program
    given in the manual and I see the out put as follows.
    Output is given below.

    My question is here.
    Since Carbon-12 hit on a rock, I expect it to do
    nuclear defragmentation and inelastic scattering etc.
    In the reading program, ENDIST array
    ENDIST(1) for Ionization, ENDIST(2) for pi0, electron
    etc. ENDIST(3) for energy deposited by nuclear recoil
    and and heavy fragments. In the output it always shows
    zero. Is it possible? or what am i doing wrong.

    Also output shows it does ionization and EM part only.
    most of the energy is outside.
    Any help is appreciated. 12cqmd.inp and rdevtd.f file
    to read EVENTDAT output is attached.

    Very happy new year to you all.
    Thank you,
    Gin

    =============== A part of the output ============

     Test of 12-C transport in Rock
     DATE: 1/ 9/ 6, TIME: 22:27:26
    Number of regions: 4 Number of scored
    quantities: 2
    The scored quantities are: 208 211
     
    NCASE, WEIPRU, ENETOT: 1 1.000 240.0
    ENDIST:
        1 4.163 2 1.495
        3 0.000 4 0.000
        5 234.3 6 0.000
        7 0.000 8 0.000
        9 0.000 10 0.000
       11 0.000 12 1.1084E-12
    Quantity n. 1 (208):
     Scoring per region:
          1 234.3
          3 5.658
    Quantity n. 2 (211):
     Scoring per region:
          1 0.2815
          3 1.495
     856510 0 0. 0. 0. 0.
                                                          
                             
    NCASE, WEIPRU, ENETOT: 2 1.000 240.0
    ENDIST:
        1 4.164 2 1.481
        3 0.000 4 0.000
        5 234.4 6 0.000
        7 0.000 8 0.000
        9 0.000 10 0.000
       11 0.000 12 3.6948E-13
    Quantity n. 1 (208):
     Scoring per region:
          1 234.4
          3 5.644
    Quantity n. 2 (211):
     Scoring per region:
          1 0.1189
          3 1.481
     1723015 0 0. 0. 0. 0.
                                                          
                             
    NCASE, WEIPRU, ENETOT: 3 1.000 240.0
    ENDIST:
        1 4.168 2 1.790
        3 0.000 4 0.000
        5 234.0 6 0.000
        7 0.000 8 0.000
        9 0.000 10 0.000
       11 0.000 12 -3.4106E-13
    Quantity n. 1 (208):
     Scoring per region:
          1 234.0
          3 5.957
    Quantity n. 2 (211):
     Scoring per region:
          1 0.4920
          3 1.790
     2661116 0 0. 0. 0. 0.
                                                          
                             
    NCASE, WEIPRU, ENETOT: 4 1.000 240.0
    ENDIST:
        1 4.166 2 1.518
        3 0.000 4 0.000
        5 234.3 6 0.000
        7 0.000 8 0.000
        9 0.000 10 0.000
       11 0.000 12 1.1369E-13
    Quantity n. 1 (208):
     Scoring per region:
          1 234.3
          3 5.685
    Quantity n. 2 (211):
     Scoring per region:
          1 0.6156
          3 1.518

    etc.

    ===================== 12cqmd.inp
    ==============================
    TITLE
     Test of 12-C transport in Rock
    DEFAULTS
                   PRECISIO
    *23456789 123456789 123456789 123456789 123456789
    123456789 123456789 123456789
    BEAM -20.0
                   HEAVYION
    HI-PROPE 6.0 12.0
    BEAMPOS 0.0 0.0 -50.0
    *EVENTYPE 2.0
    EVENTYPE 2.0
                   DPMJET
    DPMJET 0.0 0.0 0.0 0.0
          0.0 my card
    GEOBEGIN 0.01
                   COMBINAT
    *23456789 123456789 123456789 123456789 123456789
    123456789 123456789 123456789
                             standard rock slab
      SPH 1 0.0 0.0 0.0 +500.0
      RPP 2-100.0 +100.0 -100.0 +100.0
    -100.0 +200.0
      XYP 3 0.0
      XYP 4 32.0
      END
    * black hole
        1 5 1 -2
    * vacuum at the beginning
        2 5 2 3
    * rock
        3 5 2 4 -3
    * vacuum at the end
        4 5 2 -4
      END
    GEOEND
    *23456789 123456789 123456789 123456789 123456789
    123456789 123456789 123456789
    MATERIAL 1.0 1.0079 .0000899 3.0
           1.0 HYDROGEN
    MATERIAL 6.0 12.01 2.25 6.0
                   CARBON
    MATERIAL 8.0 15.999 0.001429 8.0
                   OXYGEN
    MATERIAL 12.0 24.305 1.738 9.0
                   MAGNESIU
    MATERIAL 13.0 26.982 2.70 10.0
                   ALUMINUM
    MATERIAL 19.0 39.102 0.031165 11.0
                   POTASSIU
    MATERIAL 20.0 40.08 1.54 12.0
                   CALCIUM
    MATERIAL 14.0 28.086 0.5825 13.0
                   SILICON
    *
    *23456789 123456789 123456789 123456789 123456789
    123456789 123456789 123456789
    MATERIAL 1.9 14.0
                   DETRIT
    COMPOUND -0.09800 10.0 -0.016000 11.0
    -0.11600 13.0 DETRIT
    COMPOUND -0.50500 8.0 -0.005000 3.0
    -0.06000 6.0 DETRIT
    COMPOUND -0.00300 9.0 -0.196000 12.0
                   DETRIT
    *
    MATERIAL 2.675 15.0
                   KARST
    COMPOUND -0.40000 12.0 -0.12000 6.0
    -0.47900 8.0 KARST
    *
    MATERIAL 2.625 16.0
                   LIMSTON
    COMPOUND -0.35000 12.0 -0.036000 9.0
    -0.12300 6.0 LIMSTON
    COMPOUND -0.49000 8.0
                   LIMSTON
    *
    MATERIAL 2.825 17.0
                   DOLOMIT
    COMPOUND -0.20000 12.0 -0.14400 9.0
    -0.13100 6.0 DOLOMIT
    COMPOUND -0.52400 8.0
                   DOLOMIT
    *
    *23456789 123456789 123456789 123456789 123456789
    123456789 123456789 123456789
    MATERIAL 2.675 18.0
                   MARL
    COMPOUND -0.40000 10.0 -0.08400 13.0
    -0.00600 11.0 MARL
    COMPOUND -0.49300 8.0 -0.08600 6.0
    -0.28800 12.0 MARL
    COMPOUND -0.00200 3.0
                   MARL
    *
    *LOW-NEUT 0.0
    *
    ASSIGNMAT 18.0 3.0
    * External Black Hole
    ASSIGNMAT 1.0 1.0
    * Vacuum
    ASSIGNMAT 2.0 2.0
    ASSIGNMAT 2.0 4.0
    *23456789 123456789 123456789 123456789 123456789
    123456789 123456789 123456789
    STEPSIZE 2.0 1.0 4.0
    SCORE 208.0 211.0
    EVENTDAT -23.
                   EVT.SCOR
    * **** Scoring ****
    *
    *---o---- ----o---- ----o---- ----o---- ----o----
    ----o---- ----o---- ----o----
    *EMF
                    EMF-OFF
    *
    RANDOMIZE 1.
    *---o---- ----o---- ----o---- ----o---- ----o----
    ----o---- ----o---- ----o----
    START 50.0 99999999. 0.0 1.0
    STOP

    ======================== Reader program
    ===============

    * This program reads the output from EVENTDAT *
    * HOW TO COMPILE
    * g77 -o rdevdtfluka rdevdt.f
    *
    * 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!)
    *
    * Default (option EVENTDAT not given): no event by
    event scoring
    *
    * Note: 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
    * (s 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).
    *
    * Note: 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).
    *
    * 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.
    *
    * The next is an example of a user program
    to read a binary file
    * written by EVENDAT.
    *
    *
    * Example:
    *...+....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
    ****************************************************************************
           PROGRAM RDEVDT
           CHARACTER*80 RUNTIT, FILNAM
           CHARACTER*32 RUNTIM
           DIMENSION ISCORE(4), ENDIST(12), REGSCO(5000,4)
           REAL*4 TOTALE
           INTEGER*4 BGOHIT
     
           WRITE(*,*) 'Name of the EVENTDAT binary file?'
           READ(*,'(A)') FILNAM
           IB = INDEX(FILNAM,' ')
           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:'
              TOTALE = 0.0
              BGOHIT = 0
              DO 500 IR = 1, NREGS
                 IF(REGSCO(IR,ISC).GT.0.0)THEN
                 WRITE(8,'(I7,3X,1P,G12.4)') IR,
    REGSCO(IR,ISC)
                 ENDIF
                 IF(IR.EQ.3)THEN
                 IF(REGSCO(IR,ISC).GT.0.0)THEN
                 IF(ISCORE(ISC).EQ.208)THEN
                 WRITE(56,*)NCASE,IR,REGSCO(IR,ISC)
                 ENDIF
                 ENDIF
                 ENDIF
    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

                    
    __________________________________________
    Yahoo! DSL Something to write home about.
    Just $16.99/mo. or less.
    dsl.yahoo.com


  • Next message: Giuseppe Battistoni: "Announcemnt of a new FLUKA course"

    This archive was generated by hypermail 2.1.6 : Tue Jan 10 2006 - 09:41:33 CET