From: gin chen (gin_chen2005@yahoo.com)
Date: Tue Jan 10 2006 - 06:09:50 CET
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
This archive was generated by hypermail 2.1.6 : Tue Jan 10 2006 - 09:41:33 CET