Problems implementing mgdraw

Advanced Features and User Routines
2

Hello,

I suggest you to have a look at the forum thread here. They explain how to obtain the values for the energy of the particle at the start and end of the step, and the same logic can be followed for the position of the particle.

Some relevant point I will highlight:

  • use the entry MGDRAW in mgdraw.f to print the information of any particle at the start and end of the step
  • have a look at trackr.inc since you are saving the particle information there described
  • as described in the thread, be aware of the fact that the photons don’t lose energy along the step. This changes the logic needed to save their energy at the start and end of the step

Cheers,
Giuseppe