<div dir="ltr"><div><br></div>great ! <div>Please check the weights of particles, if you use VRT (splitting for example), all particles will be associated with a w<0. It is used by the dose actor. </div><div>Not sure why the weights should be constant=1. (Note that I do not know IAEA phsp)</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Mar 4, 2021 at 3:19 PM G Tom <<a href="mailto:gregthom992@gmail.com">gregthom992@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><br></div><div>After stripping down my simulation and testing and using a gPS source I could score dose.</div><div>When I used the phase space files , again I could not score dose.</div><div><br></div><div>I finally figured out it was the missing
Constant Weight line in the phase space file header</div><div><br></div><div>======================<br></div><div>$RECORD_CONTENTS:<br> 1 // X is stored ?<br> 1 // Y is stored ?<br> 0 // Z is stored ?<br> 1 // U is stored ?<br> 1 // V is stored ?<br> 1 // W is stored ?
<============================ Changed from 0 to 1<br> 0 // Weight is stored ?<br> 0 // Extra floats stored ?<br> 0 // Extra longs stored ?<br><br>$RECORD_CONSTANT:<br> 26.7 // Constant Z<br> 1.0000 // Constant Weight <============================ ADDED this line<br></div><div><br></div><div>I works now I can score dose !! <br></div><div><img src="cid:ii_klui472e2" alt="image.png" width="514" height="260"></div><div><br></div><div>Thanks all</div><div><br></div><div>GT<br></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Mar 3, 2021 at 3:07 PM David Sarrut <<a href="mailto:David.Sarrut@creatis.insa-lyon.fr" target="_blank">David.Sarrut@creatis.insa-lyon.fr</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><br><div>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)</div><div>if it is wrong, for example outside the box, it will never score ...</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Mar 3, 2021 at 8:42 PM G Tom <<a href="mailto:gregthom992@gmail.com" target="_blank">gregthom992@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Thanks all. thanks David. I will start checking volumes one by one and report back. <br></div><div><br></div><div>GT<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Mar 3, 2021 at 2:20 PM David Sarrut <<a href="mailto:David.Sarrut@creatis.insa-lyon.fr" target="_blank">David.Sarrut@creatis.insa-lyon.fr</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi, <div><br></div><div>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 ... </div><div>Maybe start with as few volume as possible and add them one by one to check ...</div><div><br></div><div>HTH, </div><div>David</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Mar 3, 2021 at 6:47 PM G Tom <<a href="mailto:gregthom992@gmail.com" target="_blank">gregthom992@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>I am using GATE version 9.0 (March 2020), Geant4 version Name: geant4-10-06-patch-02 (29-May-2020)</div><div>My phase space files that produce the source particles are from IAEAphsp files and the headers look like this:</div><div><br></div><div>$RECORD_CONTENTS:<br> 1 // X is stored ?<br> 1 // Y is stored ?<br> 0 // Z is stored ?<br> 1 // U is stored ?<br> 1 // V is stored ?<br> 1 // W is stored ?<br> 0 // Weight is stored ?<br> 0 // Extra floats stored ?<br> 0 // Extra longs stored ?</div><div><br></div><div>Thanks</div><div><br></div><div>GT<br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Mar 3, 2021 at 12:43 PM G Tom <<a href="mailto:gregthom992@gmail.com" target="_blank">gregthom992@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi All, <br></div><div><br></div><div>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.</div><div><br></div><div>
<img src="cid:ii_kltpwvrp0" alt="image.png" width="541" height="326">
</div><div><br></div><div>When I read the output after 10^8 histories, the dose actor recorded only 0s. What could I be missing ? <br></div><div><br></div><div>Thanks <br></div><div><br></div><div>GT<br></div><div><br></div><div>My macro is below and placed just before the /gate/run/initialize line in my main macro and I start the simulations with <br></div><div>#=====================================================<br># START<br>#=====================================================<br>/gate/random/setEngineName MersenneTwister<br>/gate/random/setEngineSeed auto<br><br>/gate/application/setTotalNumberOfPrimaries {NumParticles}<br>/gate/application/start</div><div><br></div><div><br></div><div>#=====================================================<br># ACTORS: stats, dose, fluence etc<br>#=====================================================<br>#<ACTORS><br>/control/alias BluePhantom2XLength 480.0<br>/control/alias BluePhantom2YLength 480.0<br>/control/alias BluePhantom2ZLength 310.0<br>/control/alias BluePhantom2Material Water<br><br><br>/control/divide HalfBluePhantom2XLength {BluePhantom2XLength} 2.0<br>/control/divide HalfBluePhantom2YLength {BluePhantom2YLength} 2.0<br>/control/divide HalfBluePhantom2ZLength {BluePhantom2ZLength} 2.0<br><br>/control/alias BluePhantom2XPlacement 0.00<br>/control/alias BluePhantom2YPlacement 0.00<br><br>/control/subtract BluePhantom2ZPlacement 1000 {SSD}<br>/control/subtract BluePhantom2ZPlacement {BluePhantom2ZPlacement} {HalfBluePhantom2ZLength}<br><br>/gate/world/daughters/name BluePhantom2<br>/gate/world/daughters/insert box<br>/gate/BluePhantom2/setMaterial {BluePhantom2Material}<br>/gate/BluePhantom2/geometry/setXLength {BluePhantom2XLength} mm<br>/gate/BluePhantom2/geometry/setYLength {BluePhantom2YLength} mm<br>/gate/BluePhantom2/geometry/setZLength {BluePhantom2ZLength} mm<br>/gate/BluePhantom2/placement/setTranslation {BluePhantom2XPlacement} {BluePhantom2YPlacement} {BluePhantom2ZPlacement} mm<br>/gate/BluePhantom2/vis/setColor cyan<br>/gate/BluePhantom2/vis/forceWireframe 1<br>/gate/BluePhantom2/vis/setVisible 1<br><br>/gate/physics/Positron/SetCutInRegion BluePhantom2 0.01 mm<br>/gate/physics/Electron/SetCutInRegion BluePhantom2 0.01 mm<br>/gate/physics/Gamma/SetCutInRegion BluePhantom2 0.01 mm<br><br><br><br>/control/alias DoseActorXVoxelSize 2 #mm<br>/control/alias DoseActorYVoxelSize 2 #mm<br>/control/alias DoseActorZVoxelSize 1 #mm<br><br>/control/alias DoseActorXPlacement 0<br>/control/alias DoseActorYPlacement 0<br>/control/alias DoseActorZPlacement 0<br><br>/control/alias DoseActorXSize {BluePhantom2XLength}<br>/control/alias DoseActorYSize {BluePhantom2YLength}<br>/control/alias DoseActorZSize {BluePhantom2ZLength}<br><br>/control/divide DoseActorXResolution {DoseActorXSize} {DoseActorXVoxelSize}<br>/control/divide DoseActorYResolution {DoseActorYSize} {DoseActorYVoxelSize}<br>/control/divide DoseActorZResolution {DoseActorZSize} {DoseActorZVoxelSize}<br># Dose<br>/gate/actor/addActor DoseActor doseMeasurement<br>/gate/actor/doseMeasurement/attachTo BluePhantom2<br>/gate/actor/doseMeasurement/save {OutputDirName}/3d.mhd<br>/gate/actor/doseMeasurement/stepHitType random<br>/gate/actor/doseMeasurement/setPosition {DoseActorXPlacement} {DoseActorYPlacement} {DoseActorZPlacement} mm<br>/gate/actor/doseMeasurement/setVoxelSize {DoseActorXVoxelSize} {DoseActorYVoxelSize} {DoseActorZVoxelSize} mm<br>#/gate/actor/doseMeasurement/setResolution {DoseActorXResolution} {DoseActorYResolution} {DoseActorZResolution}<br>/gate/actor/doseMeasurement/setSize {DoseActorXSize} {DoseActorYSize} {DoseActorZSize} mm<br>/gate/actor/doseMeasurement/enableEdep false<br>/gate/actor/doseMeasurement/enableUncertaintyEdep false<br>/gate/actor/doseMeasurement/enableSquaredEdep false<br>/gate/actor/doseMeasurement/enableDose true<br>/gate/actor/doseMeasurement/enableSquaredDose true<br>/gate/actor/doseMeasurement/enableUncertaintyDose false<br>/gate/actor/doseMeasurement/enableNumberOfHits false</div><div><br></div><div><br></div><div><br></div><div><br></div></div>
</blockquote></div>
_______________________________________________<br>
Gate-users mailing list<br>
<a href="mailto:Gate-users@lists.opengatecollaboration.org" target="_blank">Gate-users@lists.opengatecollaboration.org</a><br>
<a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" rel="noreferrer" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">David Sarrut, Phd<br>Directeur de recherche CNRS<br>CREATIS, UMR CNRS 5220, Inserm U1294<div>Centre de lutte contre le cancer Léon Bérard<br>28 rue Laënnec, 69373 Lyon cedex 08<br>Tel : 04 78 78 51 51 / 06 74 72 05 42<br><a href="http://www.creatis.insa-lyon.fr/~dsarrut" target="_blank">http://dsarrut.github.io</a><br>_________________________________</div><div> "2 + 2 = 5, for extremely large values of 2"<br>_________________________________</div></div></div></div></div></div></div>
</blockquote></div>
_______________________________________________<br>
Gate-users mailing list<br>
<a href="mailto:Gate-users@lists.opengatecollaboration.org" target="_blank">Gate-users@lists.opengatecollaboration.org</a><br>
<a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" rel="noreferrer" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">David Sarrut, Phd<br>Directeur de recherche CNRS<br>CREATIS, UMR CNRS 5220, Inserm U1294<div>Centre de lutte contre le cancer Léon Bérard<br>28 rue Laënnec, 69373 Lyon cedex 08<br>Tel : 04 78 78 51 51 / 06 74 72 05 42<br><a href="http://www.creatis.insa-lyon.fr/~dsarrut" target="_blank">http://dsarrut.github.io</a><br>_________________________________</div><div> "2 + 2 = 5, for extremely large values of 2"<br>_________________________________</div></div></div></div></div></div></div>
</blockquote></div>
_______________________________________________<br>
Gate-users mailing list<br>
<a href="mailto:Gate-users@lists.opengatecollaboration.org" target="_blank">Gate-users@lists.opengatecollaboration.org</a><br>
<a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" rel="noreferrer" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">David Sarrut, Phd<br>Directeur de recherche CNRS<br>CREATIS, UMR CNRS 5220, Inserm U1294<div>Centre de lutte contre le cancer Léon Bérard<br>28 rue Laënnec, 69373 Lyon cedex 08<br>Tel : 04 78 78 51 51 / 06 74 72 05 42<br><a href="http://www.creatis.insa-lyon.fr/~dsarrut" target="_blank">http://dsarrut.github.io</a><br>_________________________________</div><div> "2 + 2 = 5, for extremely large values of 2"<br>_________________________________</div></div></div></div></div></div></div>