|Message: GPS halves the bin content of user energy histogram||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
I have a problem with the user defined energy histogram of the GPS, in some cases it seems to halve the provided bin content.
For the following minimal working example I used the extended/eventgenerator/exgps example program of Geant4 10.2p01 with the attached macro files.
The file gpsUserHist1.in contains the energy histograms filled with experimental values, kineticEnergy1.png is a plot of the resulting kinetic energies of the primary particles. The kinetic energy shows at 0.171MeV clearly a dip which is not in the source spectrum. The only exceptional feature of the source spectrum at these energies are two consecutive bins which contain the exact same content.
This problem is also affected by the bin boundaries: The file gpsUserHist2.in contains a similar source spectrum with the same y-values but the bin boundaries are now changed to multiples of 1MeV. Now, the spectrum shown in kineticEnergy2.png resembles perfectly the source spectrum, having with ~70000 roughly twice the counts at 0.171Mev compared to kineticEnergy1.png.
How can I obtain the correct spectrum and still using the original bin boundaries? Iím grateful for any help on this matter.
Attachment: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2017/03/20/08.22-33955-gpsUserHist1.in http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2017/03/20/08.22-52902-gpsUserHist2.in http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2017/03/20/08.22-98098-kineticEnergy1.png http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2017/03/20/08.22-71893-kineticEnergy2.png