Message: Re: Optical problem with geometry subtraction Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Optical problem with geometry subtraction 

Forum: Geometry
Re: None Optical problem with geometry subtraction (Whitney Armstrong)
Re: None Re: Optical problem with geometry subtraction (Vladimir Grichine)
Re: None Re: Optical problem with geometry subtraction (Whitney Armstrong)
Re: Feedback Re: Optical problem with geometry subtraction (Gumplinger Peter)
Date: 21 Aug, 2008
From: alexandre Creusot <alexandre Creusot>

Hi all,

  I had exactly the same problem than Withney. I used a trick to solved it...
in the file G4OpBoundaryProcess.cc:
in the G4OpBoundaryProcess::PostStepDoIt method, I skip the boundary process for the specific case of the involved volumes (and in the right order). 
If HDPE and vacuum are the material of the two volumes after the reflection then: 

  G4StepPoint *pPreStepPoint  = aStep.GetPreStepPoint();
  G4StepPoint *pPostStepPoint = aStep.GetPostStepPoint();
  if ((pPostStepPoint->GetStepStatus() == fGeomBoundary) && 
      (pPreStepPoint->GetPhysicalVolume()->GetLogicalVolume()->GetMaterial()->GetName() == "HDPE") && 
      (pPostStepPoint->GetPhysicalVolume()->GetLogicalVolume()->GetMaterial()->GetName() == "Vacuum"))
        return G4VDiscreteProcess::PostStepDoIt(aTrack, aStep);

The track is not killed anymore... it is a bit dangerous and radical... but well, this is the only way I found!

best regards, Alexandre.

 Add Message Add Message
to: "Re: Optical problem with geometry subtraction"

 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 ]