Message: GPS halves the bin content of user energy histogram Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None GPS halves the bin content of user energy histogram 

Forum: Particles
Date: 20 Mar, 2017
From: Holger Kluck <Holger Kluck>

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

 Add Message Add Message
to: "GPS halves the bin content of user energy histogram"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]