[Gate-users] proton trajectory parameters actor

[Nils Krah]

Hi Dmitry,
  
​To extend Nicolas' suggestion:
  
​
  
Bear in mind that a PhaseSpaceActor by default stores the particles' kinematic properties once it enters the volume to which the actor is attached.
  
​
  
​​I see two ways in which you can use the   PhaseSpaceActor   for your purpose:
  
​
  
​1) Insert many   PhaseSpaceActors   at different depths into the volume in which you want to track the protons. You need to be a bit careful with the geometry because each actor needs a volume to attach to which in turn needs to be inserted into the volume across which the protons are tracked. Maybe very thin slices of a suitable material would work. Essentially, you are placing virtual tracking detectors into your volume. As a result you will get a list of events per actor, i.e. a list of root files, from which you need to merge information in post-processing, e.g. by matching the EventID. Per proton you will then get a set of transverse positions in discretized depths (where you placed the actors). Bear in mind that nuclear interactions can produce secondary particle with the same EventID, so that multiple positions can appear per actor. You will want some filtering strategy, based on what you try to analyze.
  
​
  
​2) Option 2: Attach one actor to the entire volume across which the protons are tracked and turn on the option "storeAllSteps". This will write an entry for each step of each proton   to a root file, and not only when the proton enters the volume. You get one root file. You get the full information about where the proton and its descendants were at each tracking step. Clearly, output data is much larger. And the depths at which information is stored is not regularly spaced because Geant4 dynamically adjusts the step length. So again, you need post-processing to merge entries with the same EventID in the root file and take care of secondary particles according to your needs. If you want to compare transverse positions of the protons among each other, you might need to interpolate to common depths.
  
​
  
​Hope my explanations are clear, please ask if not.
  
​
  
​Cheers,
  
​Nils
  
  
On Feb 11 2021, at 10:39 am, Nicolas Arbor  <nicolas.arbor at iphc.cnrs.fr>  wrote:
  
>   
> Hi Dmitry,
>   
>   
> I think the PhaseSpace actor could be suitable for that purpose, no ?
>   
>   
>   https://opengate.readthedocs.io/en/latest/tools_to_interact_with_the_simulation_actors.html#phase-space-actor
>   
>   
> Cheers,
>   
> Nicolas
>   
>   
>   
> Le 11/02/2021 à 10:12, Dmitry Hits a écrit   :
>   
>   
> >   
> >   
> > Dear Gate-Users,
> >   
> >   
> > I was wondering if someone has already written an actor that stores the parameters of proton trajectory?
> >   
> > In particular, I am interested to know the energy and direction of a proton at different points its trajectory.
> >   
> >   
> > Cheers,
> >   
> >   
> > Dmitry
> >   
> >   
> >   
> > _______________________________________________
> >   
> > Gate-users mailing list
> >   
> > Gate-users at lists.opengatecollaboration.org (mailto:Gate-users at lists.opengatecollaboration.org)http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
> >   
> >   
>   
>   
>   
> --
>   
> ***********************************************************************
>   
> Nicolas ARBOR
>   
> Maître de Conférences
>   
> Université de Strasbourg - Faculté de Physique et Ingénierie
>   
> Groupe de recherche DeSIs - Institut Pluridisciplinaire Hubert Curien
>   
> 23, rue du Loess BP28 67037 Strasbourg Cedex 2
>   
> mail :  nicolas.arbor at iphc.cnrs.fr (mailto:nicolas.arbor at iphc.cnrs.fr)  / tel : 03.88.10.64.27
>   
> ***********************************************************************
>   
>   
>   
> _______________________________________________
>   
> 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/pipermail/gate-users/attachments/20210211/37ff95b9/attachment.html>