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

Wed Mar 3 21:07:27 CET 2021

another idea: the position of the actor, within the BluePhantom2 volume:
setPosition is according to the coordinate system of the box BluePhantom2
(its center). (/gate/actor/doseMeasurement/setPosition XXX)
if it is wrong, for example outside the box, it will never score ...

On Wed, Mar 3, 2021 at 8:42 PM G Tom <gregthom992 at gmail.com> wrote:

> 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"
>> _________________________________
>>
> _______________________________________________
> 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/8e42d48d/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/8e42d48d/attachment-0001.png>