Dear experts,
I am currently simulating neutron bombardment on iron and need to collect event information where the neutron energy deposition in iron exceeds 10 keV.
I have implemented the ENDRAW routine in MGDRAW, but file 99 shows no output.
When I filter only for neutron, all regions recorded in file 90 show as “1” (black hole).
sorry for the delay, I thought there was something wrong with the way you write the file (no need to open units, and also there was a weird closing of the output unit in a separate entry of MGDRAW). Regardless, this is not the issue.
The main issue is that ENDRAW scores energy losses to a material, but neutrons do not interact electromagnetically, and the interactions they undergo result in nuclear recoils or other secondaries. These particles are then transported, and thus are not counted as losses, just energy transfer. A user had a similar question last year in this post. The ENDRAW is then invoked as charged secondaries slow down in the material, depositing energy in a form that is not transported (below delta-ray production, e.g.)
sorry for the long delay, I missed this message. I hope you had found a way to do this, but in case you are still working on this, here is my suggestion.
The quickest way to do this is to modify your current USERDUMP card to have Dump: User Defined # WHAT(4) = 1)
and to invoke USDRAW instead of ENDRAW. In USDRAW, you can filter for ICODE = 300, and print:
the stack of secondaries (indexed by NP, for TKI, KPART, etc, look at genstk.inc)
the stack of heavies, since it could produce a recoil nucleus (indexed by NPHEAV, and KHEAVY(NPHEAV) for some properties, look at fheavy.inc)
these stacks will be populated with the output from the neutron interaction. Look at the include files for these stack to see what information is relevant for you.