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

Feedback Re: GetStepStatus() == fGeomBoundary - Problem 

Forum: Event and Track Management
Re: None GetStepStatus() == fGeomBoundary - Problem (Asim)
Date: 20 Jul, 2007
From: Gumplinger Peter <Gumplinger Peter>


How thick is your "thin layer"? If you put it at 100mm and it is:

G4box: z_width/2 = 1mm(0.1mm)

the boundary crossing will be at 99mm(99.9mm) and not >99.9*mm

I suggest you write your prePoint.z() when pre->GetStepStatus() == fGeomBoundary and see what you get.

In any case, the cleaner solution is:

G4String name = pre->GetPhysicalVolume()->GetLogicalVolume()->GetName();

if(pre->GetStepStatus() == fGeomBoundary && name == "test_logical"){}

assuming you called your test volume "test_logical". This way your code will work anywhere you place it not just at 100cm.

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

1 None: Re: GetStepStatus() == fGeomBoundary - Problem   (Asim - 20 Jul, 2007)
 Add Message Add Message
to: "Re: 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 ]