Hello All,
I have noticed some peculiar behaviour when using GPS point-wise with log interpolation to generate particles. I am unsure if I am missing something in my macro construction or if it is a bug in the backend code or an issue with my machine/installs.
I have created a minimal working example for use with examples/extended/eventgenerator/exgps
step_pointwise_forum.txt (530 Bytes)
I have run this example with v10.3.3 and v10.6.2 and see the same results.
The plot below compares the input pointwise data in red to the simulated results with the interpolation /gps/hist/inter Log and /gps/hist/inter Lin. The dotted black line is a log interpolation I performed in python to check what output should be occurring.
The Log data shown (in blue) is similar to that of the linear dataset however it stops between 2 and 3 MeV with asymptotic behaviour in the 1.945-1.95 MeV bin and the 2.715-2.720 MeV bin.
If anyone has any ideas as to what may be causing this, I would be very grateful.
Many thanks
Michael