[fluka-discuss]: Counting secondary neutrons using usdraw

From: Felicia Sutanto <fsutanto_at_umich.edu>
Date: Tue, 27 Nov 2018 06:50:55 -0800

Hello FLUKA users,

I am a newbie and would really appreciate any help!
I am simulating fast neutrons interactions in a water detector.
My goal is to count all secondary neutrons created.

Here is my current effort:
In mgdraw.f, I add the following in the USDRAW:
      IF ( ICODE .EQ. 101 ) THEN
 WRITE( IODRAW,* ) "# Inelastic coll at ", XSCO, YSCO, ZSCO
      do 10 INP = 1, NP
      WRITE( IODRAW,* ) usIn, TKI(INP), MREG, KPART(INP),
     & CXR(INP),CYR(INP),CZR(INP),LTRACK, ICODE
 10 continue
      END IF
Here is an example of an output for one event:

(erg,region,pdgCode): (0.0614377,3,8) (Tx,Ty,Tz):
(0.391015,0.00256892,-0.920381) (generation no,interaction): (1,101)

(erg,region,pdgCode): (0.0110293,3,1) (Tx,Ty,Tz):
(0.450272,0.211508,-0.867479) (generation no,interaction): (1,101)

(erg,region,pdgCode): (0.00220192,3,8) (Tx,Ty,Tz):
(0.058387,-0.934017,0.352425) (generation no,interaction): (1,101)

(erg,region,pdgCode): (0.00677268,3,1) (Tx,Ty,Tz):
(-0.411252,0.366271,0.834696) (generation no,interaction): (1,101)

(erg,region,pdgCode): (0.00378545,3,7) (Tx,Ty,Tz):
(-0.27657,-0.292502,0.915397) (generation no,interaction): (1,101)

(erg,region,pdgCode): (0.0026211,3,7) (Tx,Ty,Tz):
(0.643253,-0.111001,0.757565) (generation no,interaction): (1,101)

(erg,region,pdgCode): (0.000197383,3,7) (Tx,Ty,Tz):
(-0.20775,0.535758,-0.818415) (generation no,interaction): (1,101)

(erg,region,pdgCode): (0.0550734,3,8) (Tx,Ty,Tz):
(0.400754,0.24182,-0.883696) (generation no,interaction): (2,101)

(erg,region,pdgCode): (0.00610732,3,7) (Tx,Ty,Tz):
(-0.751712,0.657577,-0.050218) (generation no,interaction): (2,101)

(erg,region,pdgCode): (0.0266198,3,8) (Tx,Ty,Tz):
(0.833339,-0.0743435,-0.54774) (generation no,interaction): (3,101)

(erg,region,pdgCode): (0.000992194,3,8) (Tx,Ty,Tz):
(0.534859,0.646249,-0.544324) (generation no,interaction): (3,101)

(erg,region,pdgCode): (0.00150951,3,1) (Tx,Ty,Tz):
(0.0855775,0.54351,0.835029) (generation no,interaction): (3,101)

(erg,region,pdgCode): (0.000419258,3,7) (Tx,Ty,Tz):
(-0.934926,-0.0832789,-0.344932) (generation no,interaction): (3,101)

(erg,region,pdgCode): (0.0202234,3,8) (Tx,Ty,Tz):
(0.808244,0.204056,-0.552361) (generation no,interaction): (4,101)

(erg,region,pdgCode): (0.00448876,3,7) (Tx,Ty,Tz):
(0.0615197,0.334407,0.940419) (generation no,interaction): (4,101)

(erg,region,pdgCode): (0.00177942,3,7) (Tx,Ty,Tz):
(0.943334,0.25021,-0.217983) (generation no,interaction): (4,101)
For each event, I count all neutrons that are produced (pdgCode = 8 =
neutron). I then create a histogram of the neutron multiplicity.

The problem with this is that I double count these neutrons...
For example, the first two neutrons (highlighted in red) are the same
parent neutron (generation no = 1).
In Geant4, I can avoid double counting by knowing the trackID.
Is there a way to get "trackID" in FLUKA as well?
Apologies if the question is very obvious, I am still new to FLUKA..

Thanks and best regards,
Felicia


__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Wed Nov 28 2018 - 22:18:06 CET

This archive was generated by hypermail 2.3.0 : Wed Nov 28 2018 - 22:18:10 CET