Message: Re: fGeomBoundary and SensitiveDetector problem Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Disagree Re: fGeomBoundary and SensitiveDetector problem 

Forum: Event and Track Management
Re: Question fGeomBoundary and SensitiveDetector problem (Gergely Soti)
Date: 19 Jan, 2009
From: Gumplinger Peter <Gumplinger Peter>

Hi Gergely,

I cannot reproduce your problem! The first step inside a volume has the PreStepPoint StepStatus == fGeomBoundary.

I have modified LXe extended example to shoot e+:

in LXePrimaryGeneratorAction.cc -> FindParticle(particleName="e+"));

and so it will continuously track the e+:

LXeOpticalPhysics.cc -> // theScintProcess->SetTrackSecondariesFirst(true);

and similarly for theCerenkovProcess

I added this printout in ProcessHits (LXeScintSD.cc):

  if( aStep->GetPostStepPoint()->GetStepStatus() == fGeomBoundary) {
    G4cout << "PostStepPoint on Boundary" << G4endl;
  }
  if( aStep->GetPreStepPoint()->GetStepStatus() == fGeomBoundary) {
    G4cout << "PreStepPoint on Boundary" << G4endl;
  }
  G4cout << "PreStepPoint Step Status: "
         << aStep->GetPreStepPoint()->GetStepStatus()
         << G4endl;

"scintillator" is the SD - I get this printout, exactly what I'd expect:

.....
   16   -161 um  -14.9 um  -11.3 cm    337 keV  36.5 keV  30.8 um   8.71 cm      housing         msc
   17   -166 um  -14.5 um  -11.3 cm    337 keV   446 eV   4.94 um   8.71 cm      housing  Transportation
   18   -177 um  -14.2 um  -11.3 cm    334 keV   3.1 keV  11.5 um   8.71 cm   scintillator         msc
PreStepPoint on Boundary
PreStepPoint Step Status: 1
   19   -188 um  -13.2 um  -11.3 cm    328 keV  5.52 keV  11.5 um   8.71 cm   scintillator         msc
PreStepPoint Step Status: 3
......
......
   28   -251 um   45.7 um  -11.3 cm    284 keV  15.7 keV  11.5 um   8.72 cm   scintillator         msc
PreStepPoint Step Status: 3
   29   -256 um   47.9 um  -11.3 cm    284 keV   488 eV   8.41 um   8.72 cm   scintillator  Transportation
PostStepPoint on Boundary
PreStepPoint Step Status: 3
   30   -262 um   52.4 um  -11.3 cm    280 keV  3.61 keV  9.04 um   8.72 cm      housing         msc
....

Peter

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

1 More: Re: fGeomBoundary and SensitiveDetector problem   (Gergely Soti - 20 Jan, 2009)
(_ Feedback: Re: fGeomBoundary and SensitiveDetector problem   (Gumplinger Peter - 20 Jan, 2009)
(_ More: Re: fGeomBoundary and SensitiveDetector problem   (Gergely Soti - 22 Jan, 2009)
(_ Sad: Re: fGeomBoundary and SensitiveDetector problem   (Gumplinger Peter - 22 Jan, 2009)
 Add Message Add Message
to: "Re: fGeomBoundary and SensitiveDetector 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 ]