[Gate-users] Attach DoseActor to Cylinder
Botnariuc, Daniela
daniela.botnariuc.19 at ucl.ac.uk
Tue Dec 15 09:45:56 CET 2020
Dear all,
Thank you for your replies, I really appreciate it!
Best wishes,
Daniela
-----Original Message-----
From: Andreas Resch <andreas.resch at meduniwien.ac.at>
Sent: 10 December 2020 08:46
To: Grevillot Loic <loic.grevillot at medaustron.at>
Cc: David.Sarrut at creatis.insa-lyon.fr; Botnariuc, Daniela <daniela.botnariuc.19 at ucl.ac.uk>; gate-users at lists.opengatecollaboration.org
Subject: Re: [Gate-users] Attach DoseActor to Cylinder
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
>> https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Flist
>> s.opengatecollaboration.org%2Fmailman%2Flistinfo%2Fgate-users&dat
>> a=04%7C01%7C%7C5155079d4ecf49cf764708d89ce7f3e9%7C1faf88fea9984c5b93c
>> 9210a11d9a5c2%7C0%7C0%7C637431867327954359%7CUnknown%7CTWFpbGZsb3d8ey
>> JWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C10
>> 00&sdata=kmxWlLeN7t8YCKO%2FzzINHXDUwbowwq%2BryoVByxdUnEU%3D&r
>> eserved=0
>> [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
> https://eur01.safelinks.protection.outlook.com/?url=http:%2F%2Fwww.cre
> atis.insa-lyon.fr%2F~dsarrut&data=04%7C01%7C%7C5155079d4ecf49cf764
> 708d89ce7f3e9%7C1faf88fea9984c5b93c9210a11d9a5c2%7C0%7C0%7C63743186732
> 7954359%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLC
> JBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=vcteHc2EKVSmLZtILk9JYX7c
> GhySOAkVPoozSmz5Nbg%3D&reserved=0 [2]
> _________________________________
>
> "2 + 2 = 5, for extremely large values of 2"
> _________________________________
>
> DISCLAIMER:
> Please notice our E-Mail Disclaimer
> https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.m
> edaustron.at%2Femail-disclaimer%2F&data=04%7C01%7C%7C5155079d4ecf4
> 9cf764708d89ce7f3e9%7C1faf88fea9984c5b93c9210a11d9a5c2%7C0%7C0%7C63743
> 1867327954359%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2lu
> MzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=GOsWq2cDPx6twOmobz
> 3DqOsLU5pZlg0oLnnjZkLsHZk%3D&reserved=0 [3]
>
>
> Links:
> ------
> [1]
> https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists
> .opengatecollaboration.org%2Fmailman%2Flistinfo%2Fgate-users&data=
> 04%7C01%7C%7C5155079d4ecf49cf764708d89ce7f3e9%7C1faf88fea9984c5b93c921
> 0a11d9a5c2%7C0%7C0%7C637431867327954359%7CUnknown%7CTWFpbGZsb3d8eyJWIj
> oiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&am
> p;sdata=kmxWlLeN7t8YCKO%2FzzINHXDUwbowwq%2BryoVByxdUnEU%3D&reserve
> d=0
> [2]
> https://eur01.safelinks.protection.outlook.com/?url=http:%2F%2Fwww.cre
> atis.insa-lyon.fr%2F~dsarrut&data=04%7C01%7C%7C5155079d4ecf49cf764
> 708d89ce7f3e9%7C1faf88fea9984c5b93c9210a11d9a5c2%7C0%7C0%7C63743186732
> 7954359%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLC
> JBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=vcteHc2EKVSmLZtILk9JYX7c
> GhySOAkVPoozSmz5Nbg%3D&reserved=0
> [3]
> https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.m
> edaustron.at%2Femail-disclaimer%2F&data=04%7C01%7C%7C5155079d4ecf4
> 9cf764708d89ce7f3e9%7C1faf88fea9984c5b93c9210a11d9a5c2%7C0%7C0%7C63743
> 1867327954359%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2lu
> MzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=GOsWq2cDPx6twOmobz
> 3DqOsLU5pZlg0oLnnjZkLsHZk%3D&reserved=0
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org
> https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists
> .opengatecollaboration.org%2Fmailman%2Flistinfo%2Fgate-users&data=
> 04%7C01%7C%7C5155079d4ecf49cf764708d89ce7f3e9%7C1faf88fea9984c5b93c921
> 0a11d9a5c2%7C0%7C0%7C637431867327964358%7CUnknown%7CTWFpbGZsb3d8eyJWIj
> oiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&am
> p;sdata=%2FMTxD6p23sSnKMpSlzKYXy3Gy2GsKWp%2FqjXWcBSghxE%3D&reserve
> d=0
More information about the Gate-users
mailing list