Dear fluka users,
I’m simulating 10 MeV neutron interactions in liquid 40Ar. The goal is to count how many 40Ar(n, g)41Ar and how many 40Ar(n, n+g)40Ar.
I try to separate these two types of interactions by looking at the secondaries produced in each interaction. Below are two examples:
Example 1:
mgdraw.f:
IF (NCASE .EQ. 56 .AND. JTRACK .EQ. 8) THEN
ETNEW = ETRACK-Am(JTRACK)
do ip = 1, NP
WRITE(82,600) NCASE, JTRACK, ETNEW,
& ip, Kpart(ip), Tki(ip), NP,
& ICHTAR, IBTAR, ICRES, IBRES
end do
END IF
Output:
56 8 0.10000E-01 1 7 0.60471E-03 3 18 40 18 40
56 8 0.10000E-01 2 7 0.29167E-02 3 18 40 18 40
56 8 0.10000E-01 3 7 0.57545E-02 3 18 40 18 40
Example 2:
mgdraw.f:
IF (NCASE .EQ. 12060 .AND. JTRACK .EQ. 8) THEN
ETNEW = ETRACK-Am(JTRACK)
do ip = 1, NP
WRITE(83,600) NCASE, JTRACK, ETNEW,
& ip, Kpart(ip), Tki(ip), NP,
& ICHTAR, IBTAR, ICRES, IBRES
end do
END IF
Output:
12060 8 0.10000E-01 1 7 0.27595E-02 5 18 40 18 40
12060 8 0.10000E-01 2 7 0.14334E-02 5 18 40 18 40
12060 8 0.10000E-01 3 7 0.64536E-02 5 18 40 18 40
12060 8 0.10000E-01 4 8 0.10126E-01 5 18 40 18 40
12060 8 0.10000E-01 5 8 0.26580E-02 5 18 40 18 40
So the output shows two interactions: n + 40Ar -> 40Ar + gamma ; n + 40Ar -> 2n + 40Ar + gamma. However, they are not physical.
Does anyone know how to interpret the outputs shown here?
Thank you!
Guanying Zhu
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Fri Oct 20 2017 - 06:08:31 CEST