Problem 2308 - Incorrect interpolation in function G4EmCorrections::Value2
Summary: Incorrect interpolation in function G4EmCorrections::Value2
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/electromagnetic (show other problems)
Version: 10.7
Hardware: PC Linux
: P4 minor
Assignee: Vladimir.Ivantchenko
URL:
Depends on:
Blocks:
 
Reported: 2021-01-10 22:51 CET by Konstantin
Modified: 2021-01-11 06:06 CET (History)
1 user (show)

See Also:

Attachments
program and it's output (4.44 KB, text/plain)
2021-01-10 22:51 CET, Konstantin 		Details
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 Konstantin 2021-01-10 22:51:05 CET 
Created attachment 658 [details]
program and it's output

File G4EmCorrections.hh, line 284.

Result of interpolation:
f(x1,y1)=z11,
f(x2,y1)=(z12+z21)/2 (should be z21 ?),
f(x1,y2)=(z12+z21)/2 (should be z12 ?),
f(x2,y2)=z22.

The attached program and it's output show problem. KShellCorrection has discontinuities and negative values despite the condition of applicability of the model G4BetheBlochModel (K>2 MeV) is met.
Comment 1 Vladimir.Ivantchenko 2021-01-11 06:03:50 CET 
Hello Konstantin,

thank you very much for this report. The fix is provided, it will be available with next public release and the 1st patch to Geant4 10.7.

Vladimir