5.2.26. Dumpall option

With the setting of icntl=12, PHITS reads the detailed information of a completed transport calculation from a dumpall file created beforehand by using the dumpall option, and reproduces all of the transport result event by event. As such, the input file except the tally sections must be identical as used in the previous calculation when the dumpall option is activated. maxcas and maxbch cannot be changed either, because they are read from the dump file.

This can be a very useful function in order to calculate different tallies which were not used in the previous calculation. However, please be careful that the calculation with the dumpall option may create a huge dump file.

The dumped data written in ascii can be used on the computers, but the data written in binary cannot. The data sequence and meaning are given in the following.

(1). NCOL

NCOL is an intrinsic variable in the program and denotes identification of processes.

NCOL


1
:
start of calculation
2
:
end of calculation
3
:
end of a batch
4
:
source
5
:
detection of geometry error
6
:
recovery of geometry error
7
:
termination by geometry error
8
:
termination by weight cut-off
9
:
termination by time cut-off
10
:
geometry boundary crossing
11
:
termination by energy cut-off
12
:
termination by escape or leakage
13
:
(n,x) reaction
14
:
(n,n’x) reaction
15
:
sequential transport only for tally
16
:
surface cross for WW of xyz mesh

When NCOL=1, 2, 3, the output is finished. The followings are for NCOL \(\ge\) 4.

(2). NOCAS, NOBCH, RCASC, RSOUIN

These four data are written only for NCOL=4 and their meaning are,

NOCAS
:
current history number in this batch
NOBCH
:
current batch number
RCASC
:
real number of NOCAS+maxcas*(NOBCH-1)
RSOUIN
:
sum of the weight of source particle

(3). NO, MAT, ITYP, KTYP, JTYP, MTYP, RTYP, OLDWT

These mean

NO
:
cascade id in this history
MAT
:
material id
ITYP
:
particle type
KTYP
:
particle kf-code
JTYP
:
charge number of the particle
MTYP
:
baryon number of the particle
RTYP
:
rest mass of the particle (MeV)
OLDWT
:
weight of the particle at (x,y,z)

(a). QS

This data is written only for ITYP=12, 13, namely electron and positron. QS is dE/dx for electrons.

QS
:
dE/dx of electron at (x,y,z)

(4). IBLZ1, IBLZ2, ILEV1, ILEV2

These mean

IBLZ1
:
cell id at (x,y,z)
IBLZ2
:
cell id after crossing
ILEV1
:
level structure id of the cell at (x,y,z)
ILEV2
:
level structure id of the cell after crossing

(a). ILAT1

This is a variable of level structure of cell. The next data are written only for ILEV1 \(>\) 0 as

write(io) ( ( ILAT1(i,j), i=1,5 ), j=1,ILEV1 )

(b). ILAT2

This is a variable of level structure of cell. The next data are written only for ILEV2 \(>\) 0 as

write(io) ( ( ILAT2(i,j), i=1,5 ), j=1,ILEV2 )

(5). COSTH, UANG(1), UANG(2), UANG(3), NSURF

These variables mean as follows. These had come to be output from ver. 2.30.

COSTH
:
cosine of an angle of incidence in a surface crossing
UANG(1,2,3)
:
x,y,z component of a normal vector of its surface, respectively
NSURF
:
internal number of the surface
Note that this is different from the surface number defined in the [surface] section

(6). NAME, NCNT(1), NCNT(2), NCNT(3)

These mean

NAME
:
collision number of the particle
NCNT(1,2,3)
:
values of counter 1, 2, and 3

(7). WT, U, V, W

These mean

WT
:
weight of the particle at (xc,yc,zc)
U, V, W
:
unit vector of momentum of the particle

(8). E, T, X, Y, Z

These mean

E
:
energy of the particle at (x,y,z) (MeV)
T
:
time of the particle at (x,y,z) (nsec)
X, Y, Z
:
position coordinate of the preceding event point (cm)

(9). EC, TC, XC, YC, ZC

These mean

EC
:
energy of the particle at (xc,yc,zc) (MeV)
TC
:
time of the particle at (xc,yc,zc) (nsec)
XC, YC, ZC
:
position coordinate of the particle (cm)

(10). SPX, SPY, SPZ

These mean

SPX, SPY, SPZ
:
unit vector of spin direction of the particle

(11). NZST

This is the charge state of particles.

(12). NCLSTS

This variable is written only for NCOL=13, 14, collision case, and means the number of produced particle and nucleus. The next data are written for NCLSTS \(>\) 0 case.

(a). MATHZ, MATHN, JCOLL, KCOLL

These mean

MATHZ
:
Z number of the mother nucleus
MATHN
:
N number of the mother nucleus
JCOLL
:
reaction type id1
KCOLL
:
reaction type id2

JCOLL and KCOLL indicate the following meaning.

JCOLL


0
:
nothing happen
1
:
Hydrogen collisions
2
:
Particle Decays
3
:
Elastic collisions
4
:
High Energy Nuclear collisions
5
:
Heavy Ion reactions
6
:
Neutron reactions by data
7
:
Photon reactions by data
8
:
Electron reactions by data
9
:
P,d,a, and photo-nuclear reactions by data
10
:
Neutron event mode
11
:
Delta Ray production
12
:
Muon atomic interaction
13
:
Photon reactions by EGS5
14
:
Electron reactions by EGS5
15
:
Photon photonuclear interaction
16
:
Negative muon captured by nucleon
17
:
Muon photonuclear interaction
18
:
Electron recoil by track structure mode
19
:
Muon pair production (photon -> mu+ mu-)
20
:
User defined interaction
KCOLL


0
:
normal
1
:
high energy fission
2
:
high energy absorption
3
:
low energy n elastic
4
:
low energy n non-elastic
5
:
low energy n fission
6
:
low energy n absorption

(b). ICLUSTS, JCLUSTS, QCLUSTS, JCOUNT

These variables have an array and denote the information on the produced particle and nucleus.

do i = 1, NCLSTS
write(io) ICLUSTS(i)
write(io) ( JCLUSTS(j,i), j=0,8)
write(io) ( QCLUSTS(j,i), j=0,12)
write(io) ( JCOUNT(j,i),  j=1,3)
end do

These mean

ICLUSTS

kind of particle
0
:
nucleus
1
:
proton
2
:
neutron
3
:
pion
4
:
photon
5
:
kaon
6
:
muon
7
:
others
JCLUSTS(i)


i = 0
:
angular momentum
= 1
:
proton number
= 2
:
neutron number
= 3
:
ityp
= 4
:
status of the particle 0: real, \(<\) 0 : dead
= 5
:
charge number
= 6
:
baryon number
= 7
:
kf code
= 8
:
isomer level (0: ground state, 1: first isomer, 2: second isomer)
QCLUSTS(i)


i = 0
:
impact parameter
= 1
:
x-component of unit vector of momentum
= 2
:
y-component of unit vector of momentum
= 3
:
z-component of unit vector of momentum
= 4
:
e \(_{tot}\) = \(\sqrt{ {p}^2 + {m}^2 }\) (GeV)
= 5
:
rest mass (GeV)
= 6
:
excitation energy (MeV)
= 7
:
kinetic energy (MeV)
= 8
:
weight
= 9
:
time (nsec)
= 10
:
x coordinate (cm)
= 11
:
y coordinate (cm)
= 12
:
z coordinate (cm)