[Gate-users] Dose actor NbOfHits and non-integer weights

[Lana Beck]

Hi there,

I noticed when running a small number of events from my phase space through to a piece of silicon, that the dose actor attached to the silicon collected energy deposits in the -Edep.root file but no hits in the -NbOfHits.root file.

When I ran more events there were a small number of hits in the -NbOfHits.root file which corresponded to a small fraction of deposits in the Edep file.

I see from my phase space that most events have a weight of 0.02 and a small fraction have a weight of 1.
So I checked the source code and I see: GateImageInt mNumberOfHitsImage;
https://github.com/OpenGATE/Gate/blob/develop/source/digits_hits/include/GateDoseActor.hh#L162

So it seems that when mNumberOfHitsImage.AddValue(index, weight) is called, it cannot handle non-integer weights.

I use the -NbOfHits.root output files to calculate the uncertainty on the energy deposited.

Best wishes,
Lana
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200414/ce67131f/attachment-0001.html>