[Gate-users] DoseActor not recording any dose (empty 0's files)

[G Tom]

Thanks all. thanks David. I will start checking volumes one by one and
report back.

GT

On Wed, Mar 3, 2021 at 2:20 PM David Sarrut <
David.Sarrut at creatis.insa-lyon.fr> wrote:

> Hi,
>
> hard to say, please check any overlap in the volumes? If a volume A is
> included in a volume B, but not defined as a daughter of B, it may be
> ignored by Geant4 ...
> Maybe start with as few volume as possible and add them one by one to
> check ...
>
> HTH,
> David
>
> On Wed, Mar 3, 2021 at 6:47 PM G Tom <gregthom992 at gmail.com> wrote:
>
>> I am using GATE version 9.0 (March 2020), Geant4 version Name:
>> geant4-10-06-patch-02    (29-May-2020)
>> My phase space files that produce the source particles are from IAEAphsp
>> files and the headers look like this:
>>
>> $RECORD_CONTENTS:
>>   1   // X is stored ?
>>   1   // Y is stored ?
>>   0   // Z is stored ?
>>   1   // U is stored ?
>>   1   // V is stored ?
>>   1   // W is stored ?
>>   0   // Weight is stored ?
>>   0   // Extra floats stored ?
>>   0   // Extra longs stored ?
>>
>> Thanks
>>
>> GT
>>
>>
>>
>> On Wed, Mar 3, 2021 at 12:43 PM G Tom <gregthom992 at gmail.com> wrote:
>>
>>> Hi All,
>>>
>>> This has been puzzling me all week. I set up my simulation and I can
>>> visualize it working and interactions happening in the phantom to which I
>>> have attached a doseActor.
>>>
>>> [image: image.png]
>>>
>>> When I read the output after 10^8 histories, the dose actor recorded
>>> only 0s. What could I be missing ?
>>>
>>> Thanks
>>>
>>> GT
>>>
>>> My macro is below and placed just before the /gate/run/initialize line
>>> in my main macro and I start the simulations with
>>> #=====================================================
>>> # START
>>> #=====================================================
>>> /gate/random/setEngineName MersenneTwister
>>> /gate/random/setEngineSeed auto
>>>
>>> /gate/application/setTotalNumberOfPrimaries {NumParticles}
>>> /gate/application/start
>>>
>>>
>>> #=====================================================
>>> # ACTORS: stats, dose, fluence etc
>>> #=====================================================
>>> #<ACTORS>
>>> /control/alias BluePhantom2XLength 480.0
>>> /control/alias BluePhantom2YLength 480.0
>>> /control/alias BluePhantom2ZLength 310.0
>>> /control/alias BluePhantom2Material Water
>>>
>>>
>>> /control/divide HalfBluePhantom2XLength {BluePhantom2XLength} 2.0
>>> /control/divide HalfBluePhantom2YLength {BluePhantom2YLength} 2.0
>>> /control/divide HalfBluePhantom2ZLength {BluePhantom2ZLength} 2.0
>>>
>>> /control/alias BluePhantom2XPlacement 0.00
>>> /control/alias BluePhantom2YPlacement 0.00
>>>
>>> /control/subtract BluePhantom2ZPlacement 1000 {SSD}
>>> /control/subtract BluePhantom2ZPlacement {BluePhantom2ZPlacement}
>>> {HalfBluePhantom2ZLength}
>>>
>>> /gate/world/daughters/name BluePhantom2
>>> /gate/world/daughters/insert box
>>> /gate/BluePhantom2/setMaterial {BluePhantom2Material}
>>> /gate/BluePhantom2/geometry/setXLength {BluePhantom2XLength} mm
>>> /gate/BluePhantom2/geometry/setYLength {BluePhantom2YLength} mm
>>> /gate/BluePhantom2/geometry/setZLength {BluePhantom2ZLength} mm
>>> /gate/BluePhantom2/placement/setTranslation {BluePhantom2XPlacement}
>>> {BluePhantom2YPlacement} {BluePhantom2ZPlacement} mm
>>> /gate/BluePhantom2/vis/setColor cyan
>>> /gate/BluePhantom2/vis/forceWireframe 1
>>> /gate/BluePhantom2/vis/setVisible 1
>>>
>>> /gate/physics/Positron/SetCutInRegion BluePhantom2 0.01 mm
>>> /gate/physics/Electron/SetCutInRegion BluePhantom2 0.01 mm
>>> /gate/physics/Gamma/SetCutInRegion BluePhantom2 0.01 mm
>>>
>>>
>>>
>>> /control/alias DoseActorXVoxelSize 2 #mm
>>> /control/alias DoseActorYVoxelSize 2 #mm
>>> /control/alias DoseActorZVoxelSize 1 #mm
>>>
>>> /control/alias DoseActorXPlacement 0
>>> /control/alias DoseActorYPlacement 0
>>> /control/alias DoseActorZPlacement 0
>>>
>>> /control/alias DoseActorXSize {BluePhantom2XLength}
>>> /control/alias DoseActorYSize {BluePhantom2YLength}
>>> /control/alias DoseActorZSize {BluePhantom2ZLength}
>>>
>>> /control/divide DoseActorXResolution {DoseActorXSize}
>>> {DoseActorXVoxelSize}
>>> /control/divide DoseActorYResolution {DoseActorYSize}
>>> {DoseActorYVoxelSize}
>>> /control/divide DoseActorZResolution {DoseActorZSize}
>>> {DoseActorZVoxelSize}
>>> # Dose
>>> /gate/actor/addActor DoseActor doseMeasurement
>>> /gate/actor/doseMeasurement/attachTo BluePhantom2
>>> /gate/actor/doseMeasurement/save {OutputDirName}/3d.mhd
>>> /gate/actor/doseMeasurement/stepHitType random
>>> /gate/actor/doseMeasurement/setPosition {DoseActorXPlacement}
>>> {DoseActorYPlacement} {DoseActorZPlacement} mm
>>> /gate/actor/doseMeasurement/setVoxelSize   {DoseActorXVoxelSize}
>>> {DoseActorYVoxelSize} {DoseActorZVoxelSize} mm
>>> #/gate/actor/doseMeasurement/setResolution  {DoseActorXResolution}
>>> {DoseActorYResolution} {DoseActorZResolution}
>>> /gate/actor/doseMeasurement/setSize  {DoseActorXSize} {DoseActorYSize}
>>> {DoseActorZSize} mm
>>> /gate/actor/doseMeasurement/enableEdep  false
>>> /gate/actor/doseMeasurement/enableUncertaintyEdep  false
>>> /gate/actor/doseMeasurement/enableSquaredEdep  false
>>> /gate/actor/doseMeasurement/enableDose  true
>>> /gate/actor/doseMeasurement/enableSquaredDose  true
>>> /gate/actor/doseMeasurement/enableUncertaintyDose  false
>>> /gate/actor/doseMeasurement/enableNumberOfHits false
>>>
>>>
>>>
>>>
>>> _______________________________________________
>> Gate-users mailing list
>> Gate-users at lists.opengatecollaboration.org
>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>
>
>
> --
> David Sarrut, Phd
> Directeur de recherche CNRS
> CREATIS, UMR CNRS 5220, Inserm U1294
> Centre de lutte contre le cancer Léon Bérard
> 28 rue Laënnec, 69373 Lyon cedex 08
> Tel : 04 78 78 51 51 / 06 74 72 05 42
> http://dsarrut.github.io <http://www.creatis.insa-lyon.fr/~dsarrut>
> _________________________________
>  "2 + 2 = 5,  for extremely large values of 2"
> _________________________________
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210303/eb53682e/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 31916 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210303/eb53682e/attachment-0001.png>