Message: Re: Positron annihilation coordinates Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Feedback Re: Positron annihilation coordinates 

Forum: Event and Track Management
Re: Question Positron annihilation coordinates (Ottavia)
Re: Feedback Re: Positron annihilation coordinates (Gumplinger Peter)
Re: Question Re: Positron annihilation coordinates (Ottavia)
Date: 12 Mar, 2013
From: Gumplinger Peter <Gumplinger Peter>

Hi Ottavia,

Ok, so the 'feature' in your histogram is real. You should be curious as to why some events do not have an e+ annihilation. You can do this by writing out the random number at the start of one of those events and you then re-launch the event with that seed with step-wise debugging set on. If you don't know how to do this, there should be ample instructions in the G4 documentation/forum/FAQ etc.

> I even tried to use GetVertexPosition instead of GetPosition but I
> always get (0,0,0). Shouldn't it actually give me the position of the
> vertex where the track was originated?

It could well be that the vertex position is not yet filled for the track in StackingAction. It should certainly be filled as you expect in SteppingAction. Have a look.

In any event, it looks to me now that your 'feature' has a physical explanation and there is one for the missing entries as well. The only other possibility for it is the logic with which you arrive in your code to fill the histogram.

Peter

 Add Message Add Message
to: "Re: Positron annihilation coordinates"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]