Problem 2090 - No backward scattering from G4eSingleCoulombScatteringModel
Summary: No backward scattering from G4eSingleCoulombScatteringModel
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/electromagnetic/standard (show other problems)
Version: 10.4
Hardware: All All
: P4 minor
Assignee: Vladimir.Ivantchenko
URL:
Depends on:
Blocks:
 
Reported: 2018-09-10 09:20 CEST by Takao Kitaguchi
Modified: 2018-09-26 12:11 CEST (History)
0 users

See Also:

Attachments
Add an attachment (proposed patch, testcase, etc.)

Note You need to log in before you can comment on or make changes to this problem.
Description Takao Kitaguchi 2018-09-10 09:20:19 CEST 
There is no backward scattering from the G4eSingleCoulombScatteringModel class. At the 210-th line in G4eSingleCoulombScatteringModel (see the below link), the scattering angle is calculated by

 G4double cost = sqrt(1.0 - sint*sint);

 http://www-geant4.kek.jp/lxr/source/processes/electromagnetic/standard/src/G4eSingleCoulombScatteringModel.cc#L210

Therefore, cost is always positive, corresponding to the forward scattering. I think that the correct calculation is simply

  G4double cost = cos(z1);
Comment 1 Vladimir.Ivantchenko 2018-09-26 12:11:27 CEST 
Hello,

thank you very much for the detailed analysis of the problem. It is fixed now and the fix will be available for the next public release of Geant4.

Vladimir