Message: Simulating electron emission from a radioactive gas with Geant4.9.4 p02 TestEM5 (added question mark icon) Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Simulating electron emission from a radioactive gas with Geant4.9.4 p02 TestEM5 (added question mark icon) 

Forum: Physics List
Re: None Simulating electron emission from a radioactive gas with Geant4.9.4 p02 TestEM5 (Steve Asztalos)
Date: 12 Sep, 2011
From: Steve Asztalos <Steve Asztalos>

Hi,

I am using Geant4.9.4 p02 and have run TestEM5 using an isotropically emitting electron source (particleGun). I am now ready to make my simulation more realistic.

I wish to generate conversion electrons of fixed energy due to internal capture from various radioxenon isotopes enclosed by a cube of thin Si. I diligently read the post below for tritium gas but still would like to know

a) if I embed my isotope of interest in a gas of naturally occurring Xe will it be randomly distributed?

b) will Geant simulate a Xe cascade automatically (e.g, conversion electron + x-ray + Auger electrons) or must it be wired in by hand?

c) Will conversion electron slowing down in the gas be automatically handled?

d) Does (some facsimile ) of the code below get added to ConstructProcess() or AddDecay() (both gave compile-time errors).

e) How can I record/display/tally/output those events where the electrons are 1) stopped in the gas, hit or pass through my 450 um Si detector? I have spent hours reading the fora without deep insight into this question.

Thanks - Steve Asztalos -----------------------

void ExN01PhysicsList::ConstructProcess()
{
  // Define transportation process

  AddTransportation();
  ConstructEM();
  ConstructGeneral();

// Add Decay Process
  G4RadioactiveDecay*  theRadioactiveDecay = new G4RadioactiveDecay();
  G4GenericIon* ion = G4GenericIon::GenericIon();

  G4Decay* theDecayProcess = new G4Decay();
  theParticleIterator->reset();
  while( (*theParticleIterator)() ){
    G4ParticleDefinition* particle = theParticleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();

    if (particle == ion){
      pmanager->AddProcess(theRadioactiveDecay, 0, -1, 3);
      pmanager ->SetProcessOrdering(theRadioactiveDecay, idxPostStep);
      pmanager ->SetProcessOrdering(theRadioactiveDecay, idxAtRest);

    }
    if (theDecayProcess->IsApplicable(*particle)) { 
      pmanager ->AddProcess(theDecayProcess);
      // set ordering for PostStepDoIt and AtRestDoIt
      pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep);
      pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest);
    }
  }

 Add Message Add Message
to: "Simulating electron emission from a radioactive gas with Geant4.9.4 p02 TestEM5 (added question mark icon)"

 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 ]