Message: Re: Particle entering volume, ProcessHits() Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Feedback Re: Particle entering volume, ProcessHits() 

Forum: Event and Track Management
Re: Question Particle entering volume, ProcessHits()
Re: Feedback Re: Particle entering volume, ProcessHits() (Gumplinger Peter)
Re: News Re: Particle entering volume, ProcessHits()
Date: 16 Feb, 2012
From: Gumplinger Peter <Gumplinger Peter>

 *********************************************************************************************************
> * G4Track Information:   Particle = gamma,   Track ID = 55,   Parent ID = 2
> *********************************************************************************************************
> 
> Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
>     0   1.16 m   -30.2 cm   2.08 m     502 keV     0 eV      0 fm      0 fm   CellWrapPhys    initStep
>     1   1.16 m   -30.2 cm   2.08 m     502 keV     0 eV   1.22 mm   1.22 mm   CellWrapPhys  Transportation
>     2   1.16 m   -30.2 cm   2.08 m     502 keV     0 eV   3.98 mm    5.2 mm   CellLucitePhys  Transportation
>     3   1.18 m   -30.9 cm   2.12 m     502 keV     0 eV   4.31 cm   4.83 cm   CellActivePhys  Transportation
> cellNumber 1
> anNumber 0
> touch1->GetVolume()->GetName() CellActivePhys
> vertpos (1162.01,-301.553,2078.33)
> just_entered 1
>     4   1.19 m     -31 cm   2.13 m     173 keV  13.6 eV   1.22 cm   6.05 cm   CellLucitePhys       compt
>     5   1.19 m   -30.9 cm   2.13 m     173 keV     0 eV   3.35 mm   6.39 cm   CellLucitePhys  Transportation
>     6   1.18 m   -30.4 cm   2.12 m     133 keV  8.99 eV   1.51 cm    7.9 cm   CellActivePhys       compt
> cellNumber 1
> anNumber 0
> touch1->GetVolume()->GetName() CellActivePhys
> vertpos (1162.01,-301.553,2078.33)
> just_entered 1

Everything looks good to me .... except, I don't know why you continue to print: vertpos (1162.01,-301.553,2078.33) in the two cases. I don't know what vertpos is but it doesn't seem to change for the two entry points: (1.16 m -30.2 cm 2.08 m) and (1.19 m -30.9 cm 2.13 m) - that's your programming.

When you look careful you see what the photon does. It starts in CellWrapPhys then exits CellWrapPhys (going into CellLucitePhys), then exits CellLucitePhys (going into CellActivePhys), then exits CellActivePhys (it goes in a single step through CellActivePhys and it goes into CellLucitePhys) - it is this step where the PreStepPoint is at the entrance and the PostStepPoint at the exit of your active volume. It then steps in CellLucitePhys but not all the way through - there is a compton scatter inside CellLucitePhys. This must change its direction. On the next step it leaves CellLucitePhys and enters CellActivePhys again. The next step is in CellActivePhys but only to another location where a compton scatter occurs - you print again because the PreStepPoint is fBoundary.

Anyway, this seems all reasonable to me.

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Question: Re: Particle entering volume, ProcessHits()   (Santos - 17 Feb, 2012)
(_ Feedback: Re: Particle entering volume, ProcessHits()   (Gumplinger Peter - 23 Feb, 2012)
 Add Message Add Message
to: "Re: Particle entering volume, ProcessHits()"

 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 ]