I’m looking for a way to get a count per second on a gamma ray simulation. Here is what I have done:
- the geometry, including the sensitive detector (a NaI crystal)
- the source (a concrete formation with a custom energy spectrum)
- setup a scorer, the same size as the sensitive detector, that records energy deposit
- setup a root file output (I’m using gears.cc btw) that records everything happening in the sensitive detector
The score file gives me a number of “entry” but so does the root file with the “deposited energy” and “total energy deposited” histograms.
The problem I’m facing is that these different scoring methods give me different mesures and I don’t know which one to pick.
Bonus side question : which physics list would you recommend for gamma rays < 3MeV. I’m using QGSP_BERT via the macro /physics_lists/select QGSP_BERT