[Gate-users] DoseActor not recording any dose (empty 0's files)
G Tom
gregthom992 at gmail.com
Thu Mar 4 15:18:40 CET 2021
After stripping down my simulation and testing and using a gPS source I
could score dose.
When I used the phase space files , again I could not score dose.
I finally figured out it was the missing Constant Weight line in the phase
space file header
======================
$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 ? <============================ Changed from 0
to 1
0 // Weight is stored ?
0 // Extra floats stored ?
0 // Extra longs stored ?
$RECORD_CONSTANT:
26.7 // Constant Z
1.0000 // Constant Weight <============================ ADDED this
line
I works now I can score dose !!
[image: image.png]
Thanks all
GT
On Wed, Mar 3, 2021 at 3:07 PM David Sarrut <
David.Sarrut at creatis.insa-lyon.fr> wrote:
>
> 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/20210304/b222612b/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/20210304/b222612b/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 32373 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210304/b222612b/attachment-0003.png>
More information about the Gate-users
mailing list