Message: Photoabsorption in Neon, which phyisics list is recommended Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Photoabsorption in Neon, which phyisics list is recommended 

Forum: Physics List
Date: 21 Jan, 2015
From: Robert Schaefer <Robert Schaefer>

Dear Geant4 community,

I'm currently working on my PhD thesis, which includes photoabsorption of photons (energy range 20ev < E < 70eV) inside a low pressure ionization chamber.

Which physics list would you recommend?

Currently I'm trying to implement the PAI model, however I'm not very successful. The following code is an altered version of the "/examples/extended/electromagnetic/TestEm8/src/PhysicsList.cc" file:

// ******************************************************************** #include "B7PhysicsList.hh"

PhysicsList::PhysicsList() : G4VModularPhysicsList(), fConfig(0) {

  fConfig = G4LossTableManager::Instance()->EmConfigurator();
  G4LossTableManager::Instance()->SetVerbose(0);

  AddPAIModel("pai_photon");

  defaultCutValue = 1*mm;

  SetVerboseLevel(0);

}

PhysicsList::~PhysicsList() {}

void PhysicsList::ConstructParticle()
{
  // gamma
  G4Gamma::GammaDefinition();

  // leptons
  G4Electron::ElectronDefinition();

 // barions
  G4BaryonConstructor bConstructor;
  bConstructor.ConstructParticle();

  // ions
  G4IonConstructor iConstructor;
  iConstructor.ConstructParticle();
}

void PhysicsList::AddPAIModel(const G4String& modname)
{
  theParticleIterator->reset();  
  while ((*theParticleIterator)())
  {
    G4ParticleDefinition* particle = theParticleIterator->value();
    G4String partname = particle->GetParticleName();
    if(partname == "e-") { NewPAIModel(particle, modname, "eIoni"); } 
  }
}

void PhysicsList::NewPAIModel(const G4ParticleDefinition* part, 
                              const G4String& modname,
                              const G4String& procname)
{
  G4String partname = part->GetParticleName();

  if(modname == "pai") {
    G4PAIModel* pai = new G4PAIModel(part,"PAIModel");
	G4cout << "Adding model PAI with " << partname << G4endl;
    fConfig->SetExtraEmModel(partname,procname,pai,"GasDetector",0.0,100.*keV,pai);//Error!
	G4cout << "Configurating " << procname << G4endl;
  } else if(modname == "pai_photon") {
    G4PAIPhotModel* pai = new G4PAIPhotModel(part,"PAIPhotModel");
	G4cout << "Adding model PAI_PHOTON with " << partname << G4endl;
    fConfig->SetExtraEmModel(partname,procname,pai,"GasDetector",
                              0.0,100.*keV,pai);
	G4cout << "Configurating " << procname << G4endl;
  }
}

void PhysicsList::SetCuts()
{
  SetCutsWithDefault();
  if (verboseLevel>0) DumpCutValuesTable();
}
// ********************************************************************

The PAI model is the only one I want to use, however when I run my code and start shooting photons on my target, I have no interaction and they go straight through. Can anyone help me out here? The reason I use the PAI model is because I'm very much interested in the interactions of the electrons (are the photoelectrons capable of ionizing again and so on and even count secondary ionization processes).

Best regards, Robert

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

1 Idea: Re: Photoabsorption in Neon, which phyisics list is recommended   (Vladimir Ivanchenko - 27 Jan, 2015)
2 Question: Re: Photoabsorption in Neon, which phyisics list is recommended   (Robert Schaefer - 03 Feb, 2015)
 Add Message Add Message
to: "Photoabsorption in Neon, which phyisics list is recommended"

 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 ]