[Gate-users] Simple Particle Counting Actor
Robert Wilson
rkwilson at mail.sdsu.edu
Thu Jun 8 21:59:32 CEST 2017
Dear Akiv (and others),
Here is my few lines of code that I am using for just collecting counts of
particle interactions with a simple crystal volume. In my simulation I am
also using a source that emits a high amount of betas, but I am only
concerned with counting the gamma interactions, and this is the reason for
the filter; so that I count ONLY the gamma particles.
/gate/actor/addActor KillActor
counts
/gate/actor/counts/save
output/Detector_Counts.txt
/gate/actor/counts/attachTo crystal
/gate/actor/counts/addFilter
particleFilter
/gate/actor/counts/particleFilter/addParticle gamma
I hope that this helps. I know I wish I would have known about this several
months ago, but such is the process. Best of luck to you.
Cheers,
*R. Kevin Wilson*
San Diego State University | *2017*
Department of Physics | *M.S Medical Physics*
Radiation Detection Measuring & Modeling Laboratory
*PET/SPECT Modeling - Gate & GEANT4*
Clemson University | *2015*
College of Engineering & Science
Department of Physics & Astronomy | *B.S. Applied Physics*
CU EBIT & Atomic Physics Laboratory *|* *Intern Recipient and Student
Researcher*
rkwilson at mail.sdsu.edu
rkwilso at g.clemson.edu
(828) 216-5310
On Thu, Jun 8, 2017 at 9:45 AM, Robert Wilson <rkwilson at mail.sdsu.edu>
wrote:
> Dear Akiv,
>
> I know exactly what you're going through, and I'm glad to know that I'm
> actually able to help someone else on here. So the simplest solution that
> you're looking for is the Kill Track actor. It's the actor described in the
> manual just before the Dose actor. It does exactly what you're wanting. The
> output is a .txt file that just spits out the number of particles that were
> stopped by the volume that the actor is attached to.
> It should be noted that this actor is not energy or particle dependent, so
> if you have a lot of different particles flying around your simulation,
> you're going to want to add a particle filter to the actor.
> I'm currently not at my work computer that has GATE and all of my code,
> but if you like I can send you a few lines as an example for how I'm doing
> the exact same thing when I do get to it. I hope this helps, and if you
> need any other help, please don't hesitate to ask.
>
> Best regards,
>
> *R. Kevin Wilson*
>
> San Diego State University | *2017*
> Department of Physics | *M.S Medical Physics*
> Radiation Detection Measuring & Modeling Laboratory
> *PET/SPECT Modeling - Gate & GEANT4*
>
>
>
> Clemson University | *2015*
> College of Engineering & Science
> Department of Physics & Astronomy | *B.S. Applied Physics*
> CU EBIT & Atomic Physics Laboratory *|* *Intern Recipient and Student
> Researcher*
>
>
> rkwilson at mail.sdsu.edu
> rkwilso at g.clemson.edu
> (828) 216-5310
>
> On Thu, Jun 8, 2017 at 1:33 AM, Andreas Resch <
> andreas.resch at meduniwien.ac.at> wrote:
>
>> Am 2017-06-07 18:46, schrieb Akiv Jhirad:
>>
>>> Hello,
>>>
>>> I have currently been using the DoseActor attached to a small geometry
>>> (representing my detector) with a resolution of 1x1x1 along with the
>>> NumberOfHits parameter, to find the total number of particles
>>> interacting with my detector. Is there a better alternative to achieve
>>> the same task?
>>>
>>> I wanted to try the FluenceActor but don't want my output being of
>>> .mhd format. Will I simply get the number of particles detected in my
>>> detector if I just change the output format to .txt instead of .mhd?
>>>
>>> --
>>> Best Regards,
>>> Akiv Jhirad
>>> _______________________________________________
>>> Gate-users mailing list
>>> Gate-users at lists.opengatecollaboration.org
>>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>>>
>>
>> Hi,
>>
>> I didn't know that there is the number of Hits option in the dose actor.
>> However, I checked in the source code and this is the line of interest
>> (version 7.2):
>> if (mIsNumberOfHitsImageEnabled) mNumberOfHitsImage.AddValue(index,
>> weight);
>>
>> so, this is definitely not the fluence, but rather gives you the number
>> of energy depositions in a voxel. You could try to half or double your step
>> limiter, then you should see an increase or decrease of number of hits.
>>
>> The second option, the fluence actor implementation is more complex. It
>> seems, that number of particles is scored if a geometry boundary is
>> crossed. So, this actor using the option "/enableNumberOfHits true" would
>> probably give you the total number of particles entering your detector.
>>
>> Note that this number of hits is not the fluence, as fluence is according
>> to ICRU report 85 defined as:
>> "The fluence, Phi, is the quotient of dN by da, where dN is the number of
>> particles incident on a sphere of cross-sectional area da"
>>
>> Yes, if you want txt output instead of mhd, just change the extension in
>> the actor.
>>
>>
>> Cheers,
>> Andreas
>>
>> _______________________________________________
>> Gate-users mailing list
>> Gate-users at lists.opengatecollaboration.org
>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20170608/bc01b74a/attachment-0001.html>
More information about the Gate-users
mailing list