Magnetic Field in BEAMnrc/DOSXYZnrc

[Imim-phys]

Hello Experts,

I am currently working with Linac simulation using BEAMnrc and DOSXYZnrc, and I would like to know if there is a way to incorporate a magnetic field in BEAMnrc or DOSXYZnrc. If so, could you please provide guidance on how to implement it? and can I specify the direction of magnetic field and its intensity?

[rtownson]

First edit the Makefile in egs_home/dosxyznrc/Makefile. Update the sources to add the EEMF line:

SOURCES = $(EGS_SOURCEDIR)egsnrc.macros \
          $(EGS_UTILS)timing.macros \
          $(MACHINE_MACROS) $(RANDOM).macros \
          dosxyznrc_user_macros.mortran srcxyznrc.macros \
          $(EGS_SOURCEDIR)transportp.macros \
          $(EGS_SOURCEDIR)pegs4_macros.mortran \
          $(EGS_UTILS)phsp_macros.mortran $(IAEA_PHSP_MACROS)\
          $(EGS_SOURCEDIR)EEMF_macros.mortran\
          $(USER_CODE).mortran \
          $(EGS_SOURCEDIR)get_inputs.mortran \
          $(EGS_SOURCEDIR)get_media_inputs.mortran $(MACHINE_MORTRAN) \
          srcxyznrc.mortran $(RANDOM).mortran \
          $(EGS_UTILS)nrcaux.mortran \
          $(EGS_SOURCEDIR)egs_utilities.mortran \
          $(EGS_SOURCEDIR)egs_parallel.mortran \
          $(EGS_SOURCEDIR)pegs4_routines.mortran \
          $(EGS_SOURCEDIR)egsnrc.mortran

Then recompile dosxyznrc. In your .egsinp file, add these lines to the transport parameters section:

:Start MC Transport Parameter:
Global ECUT= 0.7 # just an example
Global PCUT= 0.01 # just an example
Magnetic Field = 1.5 0 0 # just an example, for 1.5 T in the x direction
EM ESTEPE = 0.1
:Stop MC Transport Parameter: