[Gate-users] Attach DoseActor to Cylinder

[Andreas Resch]

Hi,

just an additional complementary info for cylindrical geometry: you may 
also use the dose scoring method and simply rescale by pi/4 (if I 
calculated correctly). Gate assumes it is a rectangular region with side 
length a and area a^2, whereas a cylinder correctly would have a base of 
(a/2)^2*pi. If I remember correctly edep you would need to re-scale by 
1.6e-7/mass [g/cm^3] (but better check those numbers again).

Best,
Andreas

Am 2020-12-10 08:52, schrieb Grevillot Loic:
> Hello,
> 
> just complementary info: unfortunately the dose actor dose not manage
> (yet) cylindrical geometries, so to store the dose in the cylinder,
> you can score the energy deposited in the cylinder and divide “by
> hand” by the cylinder mass to get the dose. Your understanding was
> correct (the dose is being calculated in a box as opposed to a
> cylinder. I ran the same simulation but scoring dose in a cuboid
> geometry whose X and Y lengths were the same as the diameter (for the
> same depth resolution) of the cylindrical volume and I get the exact
> same curve)
> 
> Best,
> 
> Loïc
> 
> FROM: Gate-users
> [mailto:gate-users-bounces at lists.opengatecollaboration.org] ON BEHALF
> OF David Sarrut
> SENT: Donnerstag, 10. Dezember 2020 06:55
> TO: Botnariuc, Daniela <daniela.botnariuc.19 at ucl.ac.uk>
> CC: gate-users at lists.opengatecollaboration.org
> SUBJECT: Re: [Gate-users] Attach DoseActor to Cylinder
> 
> Hello Daniela,
> 
> the DoseActor is attached to the PMMA cylinder, so it will only store
> deposited energy of events that occur *inside* the cylinder, whatever
> it is larger or smaller than the cylinder.
> 
> You can change the size of the DoseActor with something like:
> 
> /gate/actor/MyActor/setSize 5 5 5 cm
> 
> If the size is larger than the cylinder, areas in the (parallelepiped)
> DoseActor that are outside the cylinder will store zero.
> 
> If the size is smaller than the cylinder (I think it is the case by
> default), the DoseActor will not store all deposited energy in the
> cylinder.
> 
> HTH,
> 
> best,
> 
> David
> 
> On Wed, Dec 9, 2020 at 10:07 PM Botnariuc, Daniela
> <daniela.botnariuc.19 at ucl.ac.uk> wrote:
> 
>> Dear Gate users,
>> 
>> I am simulating a proton pencil beam and I am trying to score the
>> depth-dose curve in a cylindrical geometry. When I attach the dose
>> actor to the cylindrical geometry, it appears that the dose is being
>> calculated in a box as opposed to a cylinder. I ran the same
>> simulation but scoring dose in a cuboid geometry whose X and Y
>> lengths were the same as the diameter (for the same depth
>> resolution) of the cylindrical volume and I get the exact same
>> curve. How would I define the dose actor in order to score dose in
>> the cylinder?
>> 
>> #GEOMETRY
>> 
>> /gate/world/daughters/name   PMMA
>> 
>> /gate/world/daughters/insert cylinder
>> 
>> /gate/PMMA/geometry/setRmax 4.08 cm
>> 
>> /gate/PMMA/geometry/setRmin 0. mm
>> 
>> /gate/PMMA/geometry/setHeight 20 cm
>> 
>> /gate/PMMA/placement/setTranslation 0 0 0.10 m
>> 
>> /gate/PMMA/setMaterial Water
>> 
>> /gate/PMMA/vis/setColor cyan
>> 
>> /gate/PMMA/describe
>> 
>> #DOSE ACTOR
>> 
>> /gate/actor/addActor   DoseActor   depthdose
>> 
>> /gate/actor/depthdose/attachTo   PMMA
>> 
>> /gate/actor/depthdose/stepHitType   random
>> 
>> /gate/actor/depthdose/setPosition   0 0 0 mm
>> 
>> /gate/actor/depthdose/setResolution   1  1 1000 mm
>> 
>> /gate/actor/depthdose/saveEveryNSeconds   60
>> 
>> /gate/actor/depthdose/enableDose   true
>> 
>> /gate/actor/depthdose/enableUncertaintyDose   true
>> 
>> Thank you in advance for your time.
>> 
>> Best wishes,
>> 
>> Daniela
>> 
>> _______________________________________________
>> Gate-users mailing list
>> Gate-users at lists.opengatecollaboration.org
>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>> [1]
> 
> --
> 
> David Sarrut, Phd
> Directeur de recherche CNRS
> CREATIS, UMR CNRS 5220, Inserm U1206
> 
> 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://www.creatis.insa-lyon.fr/~dsarrut [2]
> _________________________________
> 
>  "2 + 2 = 5,  for extremely large values of 2"
> _________________________________
> 
> DISCLAIMER:
> Please notice our E-Mail Disclaimer
> http://www.medaustron.at/email-disclaimer/ [3]
> 
> 
> Links:
> ------
> [1]
> https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.opengatecollaboration.org%2Fmailman%2Flistinfo%2Fgate-users&data=04%7C01%7Cloic.grevillot%40medaustron.at%7Cf40c05e9a9d742f22add08d89cd0346e%7Cadb749e06dfd4f6fbea8cd9f8e0c9b89%7C0%7C0%7C637431765343432261%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=uj8I90npp6OnKSGJAHBLklqiyJklGS2qY7Sqiocucz4%3D&reserved=0
> [2]
> https://eur03.safelinks.protection.outlook.com/?url=http:%2F%2Fwww.creatis.insa-lyon.fr%2F~dsarrut&data=04%7C01%7Cloic.grevillot%40medaustron.at%7Cf40c05e9a9d742f22add08d89cd0346e%7Cadb749e06dfd4f6fbea8cd9f8e0c9b89%7C0%7C0%7C637431765343437259%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=UcQOfMzlKy9wEgt0IofIERcHeNCj%2FBeAfE%2BriKRoqHQ%3D&reserved=0
> [3] http://www.medaustron.at/email-disclaimer/
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users