Electronuclear interactions ICODE

Advanced Features and User Routines
3

Hello Laney,

thank you for your question, I will go through it point by point.

  1. entry USDRAW in mgdraw is called as soon as any kind of interaction has taken place in your simulation. Your result shows that some of the primary electrons don’t interact in the target region, therefore USDRAW is not called and no printing in the output.dat file is done. In fact, if you would score the particles exiting the target region, you would see that some primary electrons survive and don’t interact in the target.
  2. you can find already the answer to this question in one of your previous post. Repeating what said there, the correct filtering for electro-nuclear reactions in USDRAW is ICODE = 101 .AND. JTRACK = 3. If you are then interested in the secondary particles produced in the electro-nuclear reactions, you can access their properties with the variables in common genstk.inc
  3. continuing the previous answer, in genstk.inc you can find also variables describing the laboratory kinetic energy of the secondary particle (TKI(I)), and its direction (CXR(I), CYR(I), CZR(I) with I = 1,NP). These variables are accessible from USDRAW , so you can implement the desired logic already there before printing in the output file (no need of struprf.f).

Hope this helps you,
Giuseppe