Dear Felicia,
note that the generation number you are printing (LTRACK) is the one of
the parent particle inducing the reaction, not the one of the products you
are listing (which is increased by one).
At each inelastic interaction, the parent disappers and the produced
particles are all 'new' and belonging to the next generation. So I'm not
sure I fully understand your double counting concerns.
Note also that in the output file you may already find what you are
searching for, since - towards its end - there is a table giving the
Number of secondaries generated in inelastic interactions [=stars] per
beam particle
Kind regards
Francesco
**************************************************
Francesco Cerutti
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
tel. +41 22 7678962
On Wed, 28 Nov 2018, Felicia Sutanto wrote:
> 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 Fri Nov 30 2018 - 18:53:56 CET