Message: GetStepStatus() == fGeomBoundary - Problem Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None GetStepStatus() == fGeomBoundary - Problem 

Forum: Event and Track Management
Date: 20 Jul, 2007
From: Asim <herr.aghdiri@gmx.de>

Hi there,

I'm still simulating a proton-beam and I want to count all particles which pass z=100.

Therefor, I defined a thin layer of vacuum at z=100 (worldmaterial is air) and added the following lines to "steppingaction.cc":

const G4StepPoint* pre = step->GetPreStepPoint(); G4ThreeVector prePoint = step->GetPreStepPoint()->GetPosition(); G4double z = prePoint.z(); if(pre->GetStepStatus() == fGeomBoundary && z>99.9*mm) {...}

unfortunately, it doesn't work.

If I do type: if(pre->GetStepStatus() == fGeomBoundary) {...}

it seems to work. But since I have serveral volumes in front of the vacuumlayer, not only those boundary-crossings which I'm looking for will be countet (at least that's what I'm afraid of).

Can somebody tell me what's wrong with the first way ?

Asim.

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

1 Feedback: Re: GetStepStatus() == fGeomBoundary - Problem   (Peter Gumplinger - 20 Jul, 2007)
(_ None: Re: GetStepStatus() == fGeomBoundary - Problem   (Asim - 20 Jul, 2007)
 Add Message Add Message
to: "GetStepStatus() == fGeomBoundary - Problem"

 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 ]