Message: Error at boundary Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Error at boundary 

Forum: Fast Simulation, Transportation & Others
Date: 17 Mar, 2018
From: Mikhail <Mikhail>

Hi!

I faced to a strange thing while simulating particle propagation through water slab (z axe is perpendicular to slab). Sometimes particle passes through the boundary in a wrong way (pre and post points are different) and on the next step it is pushed back with transportation process. The kinetic energy is around 1 MeV.

As I understand particle should reach the boundary at the first step, placed to the same point in another volume on the second step, and the then propagate in new volume on next steps.

Here is an example. Particle pass the boundary and change momentum direction along Z to opposite. By the way in backward direction it is transported in correct way(

run prepoint tName e- tID 15 stepNum 4 procName msc preVol water_slab_PV pEnergy 0.660823 x -0.880212 y 0.0969188 z 0.0402352 Mx -0.723382 My 0.613168 Mz 0.461292

run postpoint tName e- tID 15 stepNum 4 procName Transportation postVol World_PV pEnergy 0.58421 x -0.880505 y 0.0971674 z 0.0404222 Mx -0.697791 My 0.610063 Mz -0.281555

run prepoint tName e- tID 15 stepNum 5 procName Transportation preVol World_PV pEnergy 0.58421 x -0.880505 y 0.0971674 z 0.0404222 Mx -0.697791 My 0.610063 Mz -0.281555

run postpoint tName e- tID 15 stepNum 5 procName Transportation postVol water_slab_PV pEnergy 0.58421 x -0.880505 y 0.0971674 z 0.0404222 Mx -0.697791 My 0.610063 Mz -0.281555

Is my situation a bug or not and what should I do?

Is it correct to set momentum direction in next volume manually like step -> GetTrack()->SetMomentumDirection(step -> GetPreStepPoint()-> GetMomentum()); ?

Or I should apply TrackingAction and AppendStep() ? Buy the way, I appreciate if one can share an example of adding point to the track))

Hope for your suggestions! Mike

*Please correct me if I posted in wrong thread

 Add Message Add Message
to: "Error at boundary"

 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 ]