Particles stuck in DOSXYZnrc following modification to ELECTR for magnetic fields

[snaqvi69]

I tried the latest version of EGSnrc. But every dosxyznrc simulation seems to be getting stuck after a few histories. Even the examples such as air_bone_interface_example.egsinp.
I am wondering if anyone has experienced this issue. Thanks.

[mchamberland]

Yes, I can confirm the behaviour, at least with the air and bone interface example.

[mchamberland]

@ftessier I think that dosxyznrc gets stuck in some of the electron transport routine and the only commit related to electron transport that I could think of is 677bbe5. So I reverted to the commit before that (de1717c) and sure enough, dosxyznrc runs properly.

[snaqvi69]

Thanks Frederic. I am still having the same issue after downloading from https://github.com/nrc-cnrc/EGSnrc.git today. Where can I get the correct version from? Thanks.

[rtownson]

For now if you replace line 1390 in $HEN_HOUSE/src/egsnrc.mortran:

IF(irnew ~= irold) [ $electron_region_change; ];

with:

$electron_region_change;

this will fix the bug and the transport will be correct if you are not using magnetic fields. After that change, recompile any apps you are using. As @mchamberland points out this was a recently introduced bug while working on magnetic fields in vacuum; we'll push a proper fix to the release asap.

[ftessier]

Thanks @snaqvi69 and @mchamberland for reporting on this rather critical bug! Will advise shortly.