Problem 1517 - low energy electron ignores magnetic filed
Summary: low energy electron ignores magnetic filed
Status: RESOLVED FIXED
Alias: None
Product: Examples/Extended
Classification: Unclassified
Component: electromagnetic/TestEM0-4. 11-18 (show other problems)
Version: 9.6
Hardware: All All
: P5 normal
Assignee: michel.maire
URL:
Depends on:
Blocks:
 
Reported: 2013-08-11 04:44 CEST by janstar1122
Modified: 2017-04-20 11:16 CEST (History)
1 user (show)

See Also:


Attachments
full code exhibiting the problem (210.00 KB, application/x-tar)
2013-08-11 04:44 CEST, janstar1122
Details
fig 1 - for 30 MeV electrons all works well (14.61 KB, image/png)
2013-08-11 04:45 CEST, janstar1122
Details
fig 2 - 30 keV electrons ignore B-filed during inozation losses (14.42 KB, image/png)
2013-08-11 04:47 CEST, janstar1122
Details
fig 3 - zoom in of 30 keV e- (21.89 KB, image/png)
2013-08-11 04:47 CEST, janstar1122
Details

Note You need to log in before you can comment on or make changes to this problem.
Description janstar1122 2013-08-11 04:44:20 CEST
Created attachment 236 [details]
full code exhibiting the problem

Hi,
I study propagation of low energy (20 keV) electron in the 0.5T magnetic filed in the volume filled with low pressure (1/25 atm)  H2 gas.

To illustrate the issue I took TestEm4 example and made only minor modifications:
- set volume to be larger box
- change material to be low density H2 gas
- added uniform B-filed of 0.5T parallel to the Z-axis
- Particle guns fires e- at large angle vs. Z axis

To test all works well for high energy electrons I set the gun to fire 3 30 MeV e-,
and got reasonably curved tracks (see fig 1).

Next, I change only energy of electrons to 30 keV and got this nonsense, shown in fig 2. That low PT tracks should spiral tightly around the z-axis (blue).
However, the want to jump transversely by many cm, many times. They seems to spiral only  between jumps.

Fig 3 shows zoom in of 30 keV e- transport. You can compare the size of jumps with the 10 cm long axis.

The whole code is also attached. I noticed that if I disable:
 ph->RegisterProcess(new G4eIonisation,         particle);
the problem is gone. But I do want to have ionization losses in my model

Please advice how to force Geant to propagate those low energy electrons correctly.
Comment 1 janstar1122 2013-08-11 04:45:59 CEST
Created attachment 237 [details]
fig 1 - for 30 MeV electrons all works well
Comment 2 janstar1122 2013-08-11 04:47:07 CEST
Created attachment 238 [details]
fig 2 -  30 keV electrons ignore B-filed during inozation losses
Comment 3 janstar1122 2013-08-11 04:47:43 CEST
Created attachment 239 [details]
fig 3 - zoom in of 30 keV e-
Comment 4 michel.maire 2013-09-08 11:56:25 CEST
solved in 9.6 patch 2