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

More Re: fGeomBoundary and SensitiveDetector problem 

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

Hello Peter,

Thank you for your answer. I just realized that my problem might have to do something with the way I read the energy:

G4StepPoint* point1 = theStep->GetPreStepPoint();
    if(point1->GetStepStatus()==fGeomBoundary){
         G4Track* track = theStep->GetTrack();
         G4double KinEnergy = track->GetKineticEnergy();
         G4cout << point1->GetPosition() << G4endl;
    }

The coordinates seem to be right - I mean that they show the entry point of the particle. So that part is ok. But if the particle interacts in that volume I get only the reduced energy, not the original one. For example: I'm shooting a 1 MeV gamma, and it deposits 550 keV. In my output I get something like: gamma, 450 keV. It seems to me that the first process which limits the step also decreases the gamma's energy, so when my function is called, the current energy is already 450 keV.

At this point it seems to me that I have 2 solutions:

1) to put vacuum in my detector so the gammas won't interact with it (which I wouldn't like to do)

2) do this whole thing in Stepping action, and work with GetPostStepPoint (to extract the gamma's energy before it enters the detector)

Ideas or comments?

Gergely

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

1 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 ]