Problem 2553 - New neutron inelastic thermal scattering data for: h_water contains wrong cos value
Summary: New neutron inelastic thermal scattering data for: h_water contains wrong cos...
Status: RESOLVED DUPLICATE of problem 2552
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic/models/neutron_hp (show other problems)
Version: 11.1
Hardware: All All
: P4 normal
Assignee: dennis.herbert.wright
URL:
Depends on:
Blocks:
 
Reported: 2023-07-13 14:24 CEST by Vasilis Vlachoudis
Modified: 2023-07-15 08:46 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 Vasilis Vlachoudis 2023-07-13 14:24:19 CEST 
The neutron Thermal Scattering Inelastic scattering for h_water
G4NDL4.7/ThermalScattering/Inelastic/FS/h_water.z
ThermalScattering_XS_JEFF33_TSL_mix_JEFF33_ENDFB8_LowPrecision

(maybe also the other files I didn't check them)
for some energies the last iso probable cos-theta is less than the previous entry. This leads to sample the cos-theta from the wrong iso-probable angular interval.
 
Inelastic: T=293.6
Ein=2.855 Eout=2.85501 j=7 x=0.999813 xprev=0.999828
Ein=3.12 Eout=3.11999 j=7 x=0.999813 xprev=0.999838
Ein=3.12 Eout=3.12001 j=7 x=0.999817 xprev=0.999831
Ein=3.75 Eout=3.74999 j=7 x=0.999825 xprev=0.999841
Ein=3.75 Eout=3.74999 j=7 x=0.999838 xprev=0.999839
Ein=3.75 Eout=3.75001 j=7 x=0.999825 xprev=0.999841
Ein=4.07 Eout=4.06999 j=7 x=0.999816 xprev=0.999865
Ein=4.07 Eout=4.07001 j=7 x=0.999816 xprev=0.999865
Inelastic: T=323.6
Ein=3.75 Eout=3.74999 j=7 x=0.99983 xprev=0.999846
Ein=3.75 Eout=3.75001 j=7 x=0.999819 xprev=0.999856
Ein=4.07 Eout=4.06999 j=7 x=0.99983 xprev=0.999842
Ein=4.07 Eout=4.07001 j=7 x=0.999834 xprev=0.999837
Inelastic: T=373.6
Ein=2.59 Eout=2.59001 j=7 x=0.999811 xprev=0.999817
Ein=2.855 Eout=2.85499 j=7 x=0.999816 xprev=0.999826
Ein=2.855 Eout=2.85501 j=7 x=0.999813 xprev=0.999826
Ein=3.12 Eout=3.11999 j=7 x=0.999812 xprev=0.999816
Ein=3.12 Eout=3.12001 j=7 x=0.99982 xprev=0.999833
Ein=3.42 Eout=3.41999 j=7 x=0.999818 xprev=0.999841
Ein=3.42 Eout=3.41999 j=7 x=0.999818 xprev=0.999823
Ein=3.42 Eout=3.42001 j=7 x=0.999818 xprev=0.999823
Ein=3.42 Eout=3.42001 j=7 x=0.999829 xprev=0.999836
Ein=3.42 Eout=3.42001 j=7 x=0.999821 xprev=0.999837
Ein=3.75 Eout=3.74999 j=7 x=0.999828 xprev=0.999835
Ein=3.75 Eout=3.74999 j=7 x=0.999838 xprev=0.999839
Ein=3.75 Eout=3.75001 j=7 x=0.999835 xprev=0.999842
Ein=3.75 Eout=3.75001 j=7 x=0.99983 xprev=0.999831
Ein=4.07 Eout=4.06999 j=7 x=0.999819 xprev=0.999861
Ein=4.07 Eout=4.06999 j=7 x=0.999839 xprev=0.999846
Ein=4.07 Eout=4.06999 j=7 x=0.999829 xprev=0.999834
Ein=4.07 Eout=4.07001 j=7 x=0.999829 xprev=0.999834
Ein=4.07 Eout=4.07001 j=7 x=0.999843 xprev=0.999843
Ein=4.07 Eout=4.07001 j=7 x=0.999835 xprev=0.999849
Ein=4.07 Eout=4.07001 j=7 x=0.999821 xprev=0.999857
Inelastic: T=423.6
Ein=2.855 Eout=2.85499 j=7 x=0.999812 xprev=0.999829
Ein=2.855 Eout=2.85499 j=7 x=0.999821 xprev=0.999822
Ein=2.855 Eout=2.85501 j=7 x=0.999821 xprev=0.999822
Ein=2.855 Eout=2.85501 j=7 x=0.999813 xprev=0.999826
Ein=3.12 Eout=3.11999 j=7 x=0.999816 xprev=0.999833
Ein=3.12 Eout=3.11999 j=7 x=0.999811 xprev=0.999817
Ein=3.12 Eout=3.12001 j=7 x=0.999812 xprev=0.999814
Ein=3.12 Eout=3.12001 j=7 x=0.999818 xprev=0.99983
Ein=3.42 Eout=3.41999 j=7 x=0.999821 xprev=0.999837
Ein=3.42 Eout=3.41999 j=7 x=0.999821 xprev=0.999842
Ein=3.42 Eout=3.41999 j=7 x=0.999815 xprev=0.999822
Ein=3.42 Eout=3.42001 j=7 x=0.999813 xprev=0.999826
Ein=3.42 Eout=3.42001 j=7 x=0.999815 xprev=0.999822
Ein=3.75 Eout=3.74999 j=7 x=0.999826 xprev=0.999839
Ein=3.75 Eout=3.74999 j=7 x=0.999835 xprev=0.999842
Ein=3.75 Eout=3.74999 j=7 x=0.999815 xprev=0.999836
Ein=3.75 Eout=3.74999 j=7 x=0.999821 xprev=0.999832
Ein=3.75 Eout=3.75001 j=7 x=0.999813 xprev=0.999839
Ein=3.75 Eout=3.75001 j=7 x=0.999815 xprev=0.999836
Ein=3.75 Eout=3.75001 j=7 x=0.999837 xprev=0.99984
Ein=3.75 Eout=3.75001 j=7 x=0.999831 xprev=0.999845
Ein=3.75 Eout=3.75001 j=7 x=0.999828 xprev=0.999836
Ein=4.07 Eout=4.06999 j=7 x=0.999817 xprev=0.999865
Ein=4.07 Eout=4.06999 j=7 x=0.999839 xprev=0.999846
Ein=4.07 Eout=4.06999 j=7 x=0.999842 xprev=0.999844
Ein=4.07 Eout=4.06999 j=7 x=0.999814 xprev=0.999847
Ein=4.07 Eout=4.06999 j=7 x=0.999818 xprev=0.999845
Ein=4.07 Eout=4.06999 j=7 x=0.999831 xprev=0.999834
Ein=4.07 Eout=4.07001 j=7 x=0.999826 xprev=0.999838
Ein=4.07 Eout=4.07001 j=7 x=0.999814 xprev=0.999847
Ein=4.07 Eout=4.07001 j=7 x=0.999842 xprev=0.999844
Ein=4.07 Eout=4.07001 j=7 x=0.999836 xprev=0.999848
Ein=4.07 Eout=4.07001 j=7 x=0.999818 xprev=0.999862
Inelastic: T=473.6
Ein=2.855 Eout=2.85499 j=7 x=0.999816 xprev=0.999821
Ein=2.855 Eout=2.85499 j=7 x=0.999813 xprev=0.999826
Ein=2.855 Eout=2.85501 j=7 x=0.999813 xprev=0.999826
Ein=2.855 Eout=2.85501 j=7 x=0.999815 xprev=0.999824
Ein=2.855 Eout=2.85501 j=7 x=0.999817 xprev=0.999818
Ein=3.12 Eout=3.11999 j=7 x=0.999817 xprev=0.99983
Ein=3.12 Eout=3.11999 j=7 x=0.999813 xprev=0.999838
Ein=3.12 Eout=3.11999 j=7 x=0.999818 xprev=0.999835
Ein=3.12 Eout=3.12001 j=7 x=0.999818 xprev=0.999835
Ein=3.42 Eout=3.41999 j=7 x=0.999821 xprev=0.999842
Ein=3.42 Eout=3.41999 j=7 x=0.999816 xprev=0.999821
Ein=3.42 Eout=3.42001 j=7 x=0.999821 xprev=0.999822
Ein=3.42 Eout=3.42001 j=7 x=0.999817 xprev=0.999819
Ein=3.42 Eout=3.42001 j=7 x=0.999821 xprev=0.999842
Ein=3.42 Eout=3.42001 j=7 x=0.999815 xprev=0.999849
Ein=3.75 Eout=3.74998 j=7 x=0.999829 xprev=0.999834
Ein=3.75 Eout=3.74999 j=7 x=0.999815 xprev=0.999859
Ein=3.75 Eout=3.74999 j=7 x=0.999828 xprev=0.999848
Ein=3.75 Eout=3.74999 j=7 x=0.99982 xprev=0.999826
Ein=3.75 Eout=3.74999 j=7 x=0.999818 xprev=0.999832
Ein=3.75 Eout=3.74999 j=7 x=0.999824 xprev=0.999832
Ein=3.75 Eout=3.75001 j=7 x=0.999817 xprev=0.999837
Ein=3.75 Eout=3.75001 j=7 x=0.999819 xprev=0.999829
Ein=3.75 Eout=3.75001 j=7 x=0.99982 xprev=0.999826
Ein=3.75 Eout=3.75001 j=7 x=0.999828 xprev=0.999848
Ein=3.75 Eout=3.75001 j=7 x=0.999822 xprev=0.999852
Ein=3.75 Eout=3.75001 j=7 x=0.999816 xprev=0.999857
Ein=3.75 Eout=3.75002 j=7 x=0.999829 xprev=0.999834
Ein=4.07 Eout=4.06998 j=7 x=0.99982 xprev=0.99986
Ein=4.07 Eout=4.06998 j=7 x=0.99983 xprev=0.999853
Ein=4.07 Eout=4.06998 j=7 x=0.99982 xprev=0.999838
Ein=4.07 Eout=4.06998 j=7 x=0.999819 xprev=0.99984
Ein=4.07 Eout=4.06998 j=7 x=0.999817 xprev=0.999843
Ein=4.07 Eout=4.06998 j=7 x=0.999814 xprev=0.999848
Ein=4.07 Eout=4.07001 j=7 x=0.999816 xprev=0.999846
Ein=4.07 Eout=4.07001 j=7 x=0.99982 xprev=0.999838
Ein=4.07 Eout=4.07001 j=7 x=0.99983 xprev=0.999853
Ein=4.07 Eout=4.07002 j=7 x=0.999826 xprev=0.999857
Ein=4.07 Eout=4.07002 j=7 x=0.99982 xprev=0.99986
Inelastic: T=523.6
Ein=3.12 Eout=3.11998 j=7 x=0.999818 xprev=0.999828
Ein=3.12 Eout=3.12001 j=7 x=0.999818 xprev=0.999828
Ein=3.12 Eout=3.12001 j=7 x=0.999819 xprev=0.999826
Ein=3.42 Eout=3.41998 j=7 x=0.999818 xprev=0.999842
Ein=3.42 Eout=3.41999 j=7 x=0.999816 xprev=0.999821
Ein=3.42 Eout=3.42001 j=7 x=0.999817 xprev=0.999819
Ein=3.42 Eout=3.42002 j=7 x=0.999818 xprev=0.999842
Ein=3.75 Eout=3.74998 j=7 x=0.999829 xprev=0.999834
Ein=3.75 Eout=3.74998 j=7 x=0.999823 xprev=0.999845
Ein=3.75 Eout=3.74998 j=7 x=0.999819 xprev=0.999851
Ein=3.75 Eout=3.74998 j=7 x=0.999818 xprev=0.999853
Ein=3.75 Eout=3.74998 j=7 x=0.999816 xprev=0.999835
Ein=3.75 Eout=3.75002 j=7 x=0.999816 xprev=0.999835
Ein=3.75 Eout=3.75002 j=7 x=0.999821 xprev=0.999825
Ein=3.75 Eout=3.75002 j=7 x=0.999818 xprev=0.999853
Ein=3.75 Eout=3.75002 j=7 x=0.99983 xprev=0.999832
Ein=4.07 Eout=4.06998 j=7 x=0.999819 xprev=0.99986
Ein=4.07 Eout=4.06998 j=7 x=0.999826 xprev=0.999828
Ein=4.07 Eout=4.06998 j=7 x=0.999825 xprev=0.99983
Ein=4.07 Eout=4.06998 j=7 x=0.999822 xprev=0.999834
Ein=4.07 Eout=4.06998 j=7 x=0.999816 xprev=0.999846
Ein=4.07 Eout=4.07002 j=7 x=0.999817 xprev=0.999844
Ein=4.07 Eout=4.07002 j=7 x=0.999823 xprev=0.999832
Ein=4.07 Eout=4.07002 j=7 x=0.999819 xprev=0.99986
Ein=4.07 Eout=4.07002 j=7 x=0.999821 xprev=0.999858
Inelastic: T=573.6
Ein=3.42 Eout=3.41998 j=7 x=0.999823 xprev=0.999831
Ein=3.42 Eout=3.41998 j=7 x=0.999822 xprev=0.999833
Ein=3.42 Eout=3.41998 j=7 x=0.999814 xprev=0.999826
Ein=3.42 Eout=3.42002 j=7 x=0.999814 xprev=0.999826
Ein=3.42 Eout=3.42002 j=7 x=0.999822 xprev=0.999833
Ein=3.42 Eout=3.42002 j=7 x=0.999824 xprev=0.99983
Ein=3.75 Eout=3.74998 j=7 x=0.999827 xprev=0.999836
Ein=3.75 Eout=3.74998 j=7 x=0.999825 xprev=0.999842
Ein=3.75 Eout=3.74998 j=7 x=0.999824 xprev=0.999844
Ein=3.75 Eout=3.74998 j=7 x=0.999821 xprev=0.999826
Ein=3.75 Eout=3.75002 j=7 x=0.999822 xprev=0.999825
Ein=3.75 Eout=3.75002 j=7 x=0.999824 xprev=0.999844
Ein=4.07 Eout=4.06998 j=7 x=0.999834 xprev=0.999835
Ein=4.07 Eout=4.06998 j=7 x=0.999829 xprev=0.999845
Ein=4.07 Eout=4.06998 j=7 x=0.999826 xprev=0.99985
Ein=4.07 Eout=4.06998 j=7 x=0.999825 xprev=0.999852
Ein=4.07 Eout=4.07002 j=7 x=0.999825 xprev=0.999841
Ein=4.07 Eout=4.07002 j=7 x=0.999822 xprev=0.999835
Ein=4.07 Eout=4.07002 j=7 x=0.999825 xprev=0.999852
Inelastic: T=623.6
Ein=3.75 Eout=3.74998 j=7 x=0.999828 xprev=0.999834
Ein=3.75 Eout=3.75002 j=7 x=0.999828 xprev=0.999834
Ein=3.75 Eout=3.75002 j=7 x=0.999829 xprev=0.999833
Ein=3.75 Eout=3.75002 j=7 x=0.99983 xprev=0.999831
Ein=4.07 Eout=4.06998 j=7 x=0.999831 xprev=0.99984
Ein=4.07 Eout=4.06998 j=7 x=0.99983 xprev=0.999842
Ein=4.07 Eout=4.06998 j=7 x=0.999816 xprev=0.999848
Ein=4.07 Eout=4.07002 j=7 x=0.999817 xprev=0.999846
Ein=4.07 Eout=4.07002 j=7 x=0.99983 xprev=0.999842
Ein=4.07 Eout=4.07002 j=7 x=0.999832 xprev=0.999838
Inelastic: T=647.2
Ein=3.12 Eout=3.11998 j=7 x=0.999813 xprev=0.999814
Ein=3.12 Eout=3.12002 j=7 x=0.999813 xprev=0.999814
Ein=3.75 Eout=3.74998 j=7 x=0.999819 xprev=0.99983
Ein=3.75 Eout=3.75002 j=7 x=0.99982 xprev=0.999828
Ein=4.07 Eout=4.06998 j=7 x=0.999834 xprev=0.999835
Ein=4.07 Eout=4.06998 j=7 x=0.999833 xprev=0.999837
Ein=4.07 Eout=4.06998 j=7 x=0.999826 xprev=0.999828
Ein=4.07 Eout=4.06998 j=7 x=0.999832 xprev=0.999836
Ein=4.07 Eout=4.07002 j=7 x=0.99983 xprev=0.999837
Ein=4.07 Eout=4.07002 j=7 x=0.999833 xprev=0.999837
Inelastic: T=800
Ein=3.12 Eout=3.12002 j=7 x=0.999816 xprev=0.999819
Ein=3.42 Eout=3.41998 j=7 x=0.999814 xprev=0.99983
Ein=3.42 Eout=3.41998 j=7 x=0.999819 xprev=0.999828
Ein=3.42 Eout=3.41998 j=7 x=0.999822 xprev=0.999826
Ein=3.42 Eout=3.42002 j=7 x=0.999822 xprev=0.999826
Ein=3.42 Eout=3.42002 j=7 x=0.999816 xprev=0.999831
Ein=3.42 Eout=3.42002 j=7 x=0.999815 xprev=0.999829
Ein=3.75 Eout=3.74998 j=7 x=0.999814 xprev=0.999845
Ein=3.75 Eout=3.74998 j=7 x=0.999821 xprev=0.999839
Ein=3.75 Eout=3.74998 j=7 x=0.999824 xprev=0.999837
Ein=3.75 Eout=3.75002 j=7 x=0.999824 xprev=0.999837
Ein=3.75 Eout=3.75002 j=7 x=0.999815 xprev=0.999844
Ein=4.07 Eout=4.06998 j=7 x=0.999816 xprev=0.999851
Ein=4.07 Eout=4.06998 j=7 x=0.999822 xprev=0.999847
Ein=4.07 Eout=4.06998 j=7 x=0.999811 xprev=0.999815
Ein=4.07 Eout=4.07002 j=7 x=0.999811 xprev=0.999813
Ein=4.07 Eout=4.07002 j=7 x=0.999822 xprev=0.999847
Ein=4.07 Eout=4.07002 j=7 x=0.999818 xprev=0.99985
Ein=4.07 Eout=4.07002 j=7 x=0.999816 xprev=0.999851
Comment 1 Gabriele Cosmo 2023-07-15 08:46:05 CEST 

*** This problem has been marked as a duplicate of problem 2552 ***