Message: Updating Direction Vector Post collision Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Updating Direction Vector Post collision 

Forum: Particles
Date: 13 Aug, 2012
From: David <David>

I am working on an analog scattering model for electrons. Essentially, I am trying to obtain the Goudsmit-Saunderson distribution using the monte carlo method. The problem doesn't model energy loss i.e. a thin foil geometry. I figured out how to implement an event by event sampling method, but I do not understand how to update the direction vector or rather the momentum vector. To outline what I am doing:

1-Generate a particle at a specific location (x,y,z) with direction (u,v,w) and energy E_0.

2-Sample a distance to collision (based on a cross sectio I provide in the physics model)

3-Transport particle to collision site

4-Sample the deflection cosine

This is where I am stuck:

5-Given deflection cosine, update particle direction and repeat

In the past, I used a the same method given in the penelope manual pg 26:

http://www.oecd-nea.org/dbprog/penelope-2003.pdf

The quantity I am sampling is cos(theta). I am looking for a way to pass cos(theta) and update the direction vector.

Thanks

David

 Add Message Add Message
to: "Updating Direction Vector Post collision"

 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 ]