Message: Re: Deposited energy greater than kinetic energy for electrons? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Note Re: Deposited energy greater than kinetic energy for electrons? 

Forum: Electromagnetic Processes
Re: Warning Deposited energy greater than kinetic energy for electrons? (Kareem Kazkaz)
Date: Oct 01, 13:23
From: Kareem Kazkaz <Kareem Kazkaz>

I think I found my problem. Here is the relevant section from my physics list:

  G4PhysListFactory factory;
  G4VModularPhysicsList *phys = factory.GetReferencePhysList( "QGSP_BIC_HP" );
  for( G4int i=0; ; ++i ) {
    G4VPhysicsConstructor *elem =
        const_cast<G4VPhysicsConstructor*> (phys->GetPhysics(i));
    if( elem == NULL ) break;
    G4cout << "RegisterPhysics: " << elem->GetPhysicsName() << G4endl;
    RegisterPhysics( elem );
  }
  
  physicsVector->push_back( new G4RadioactiveDecayPhysics() );
  physicsVector->push_back( new G4EmExtraPhysics() );
  physicsVector->push_back( new LUXSimPhysicsOpticalPhysics() );
  physicsVector->push_back( new PhysListEmLivermore() );

I set the tracking verbosity to 3, and it appeared that for electrons, eIoni is being forced as the process,
rather than LowEnergyIoni. I tried moving all the "physicsVector->" lines to before the factory loop, and
that fixed the energy non-conservation. What's strange, though is that eIoni is still the process governing
electrons, rather than LowEnergyIoni. If someone could explain how THAT's happening, I'd love to know
how to fix it.

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

1 Disagree: Re: Deposited energy greater than kinetic energy for electrons?   (Gumplinger Peter - Oct 02, 14:19)
 Add Message Add Message
to: "Re: Deposited energy greater than kinetic energy for electrons?"

 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 ]