Dear Fluka users and experts,
I am trying to store information on secondaries from inelastic interactions. I need to assign this information to the secondaries and propagate the information further to their respective daughter particles and so on. I was hoping to do this using mgdraw.f or stuprf.f, however, with the available user variables (e.g. SPAUSR) I am so far only able to assign/propagate the same information for all secondaries. An example of this is given in the code below. Here I am interested in storing (and also propagating) position and direction of neutrons created in an interaction. The information entered in SPAUSR seems to be inherited by all secondaries. This means, as far as I understand, that if there is more than one neutron created in the interaction, only the information on the last neutron (in the loop up to NP) is stored. Of course I could write all the information to a separate file, but my goal is however to separately store, and also propagate to daughter particles, information on secondaries.
More specifically, secondary neutrons (daughters) produced in inelastic proton collisions can further produce protons (granddaughters) through elastic or inelastic scatters in my setup. I need to be able to identify exactly the neutron (and its initial properties), i.e. the daughter, that produced a proton, i.e. a granddaughter, through (n,p) reactions.
Any suggestions on how to achieve this would be much apprciated.
Best regards,
Kristian Smeland Ytre-Hauge
from usdraw:
IF(ICODE .EQ. 101 .AND. MREG .EQ. 10) THEN
DO IP = 1, NP
IF(KPART(IP) .EQ. 8) THEN
* Store origin and direction of neutron (x,y,z)
SPAUSR(1)=XSCO
SPAUSR(2)=YSCO
SPAUSR(3)=ZSCO
SPAUSR(4)=Cxr(IP)
SPAUSR(5)=Cyr(IP)
SPAUSR(6)=Czr(IP)
ENDIF
END DO
END IF
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Tue May 08 2018 - 15:23:23 CEST