Hello everyone.
I am running a simulation where I am shooting 122-keV gammas into a volume of liquid xenon. I have
this line in my UserSteppingAction:
G4cout << theTrack->GetDefinition()->GetParticleName()
<< " #" << theTrack->GetTrackID() << ", "
<< "step eKin = " << theStep->GetPreStepPoint()->GetKineticEnergy() << ", "
<< "track eKin = " << theTrack->GetKineticEnergy() << ", "
<< "eDep = " << theStep->GetTotalEnergyDeposit() << ", "
<< "deltaE = " << theStep->GetDeltaEnergy() << G4endl;
Here is the output from several events:
EVENT 1
gamma #1, step eKin = 0.122, track eKin = 0, eDep = 0.001122, deltaE = -0.122
e- #2, step eKin = 0.120878, track eKin = 0.0177671, eDep = 0.103111, deltaE = -0.103111
e- #2, step eKin = 0.0177671, track eKin = 0, eDep = 0.0280091, deltaE = -0.0177671
e- #2, step eKin = 0, track eKin = 0, eDep = 0, deltaE = 0
EVENT 2
gamma #1, step eKin = 0.122, track eKin = 0, eDep = 0.0047747, deltaE = -0.122
gamma #3, step eKin = 0.0297813, track eKin = 0, eDep = 0.0051039, deltaE = -0.0297813
e- #4, step eKin = 0.0246774, track eKin = 0, eDep = 0.0398681, deltaE = -0.0246774
e- #4, step eKin = 0, track eKin = 0, eDep = 0, deltaE = 0
e- #2, step eKin = 0.087444, track eKin = 0, eDep = 0.174888, deltaE = -0.087444
EVENT 3
gamma #1, step eKin = 0.122, track eKin = 0, eDep = 0.0047747, deltaE = -0.122
gamma #3, step eKin = 0.0297813, track eKin = 0, eDep = 0.0054169, deltaE = -0.0297813
e- #4, step eKin = 0.0243644, track eKin = 0, eDep = 0.0401791, deltaE = -0.0243644
e- #4, step eKin = 0, track eKin = 0, eDep = 0, deltaE = 0
e- #2, step eKin = 0.087444, track eKin = 0, eDep = 0.103789, deltaE = -0.087444
e- #2, step eKin = 0, track eKin = 0, eDep = 0, deltaE = 0
Note that in all three of these events, at some point the energy deposited by an electron is greater than
the kinetic energy of the electron. I get this kind of error occurring often in the output, not just the first
three events. Note also that in event 2, at one point the energy deposited is greater than the energy of
the primary particle.
How is this possible?
I am using the QGSP_HP_BIC and Livermore physics lists.
Any and all help would be greatly appreciated.
Kareem
|