Message: AL27 isotopes from 100keV electrons? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question AL27 isotopes from 100keV electrons? 

Forum: Physics List
Date: 22 Dec, 2010
From: Maciej Sznajder <Maciej Sznajder>

Dear Geant4 Forum experts,

I have got a simple physical simulation. I have got a thin film made with AL, 12 mu m thick. I hit the foil with one 100 keV electron. The output for electron looks reasonable. But electron which goes through the detector generates 69 AL27 isotopes, which is from physically point of view impossible. Am I miss something in my physics list??

Please look:

const G4RegionStore* theRegionStore = G4RegionStore::GetInstance();
  G4Region* KobeDetector = theRegionStore->GetRegion("KobeDetector");

  theParticleIterator->reset();
  while( (*theParticleIterator)() ){
    G4ParticleDefinition* particle = theParticleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
    G4String particleName = particle->GetParticleName();

    if (particleName == "e-"){

      G4eIonisation* theIonisation = new G4eIonisation;
      G4PAIModel* pai = new G4PAIModel(particle, "PAIModel");
      theIonisation->AddEmModel(0,pai,pai,KobeDetector);

      pmanager->AddProcess(new G4CoulombScattering,  -1, -1, 1);
      pmanager->AddProcess(theIonisation,            -1, 2, 2);
      pmanager->AddProcess(new G4eBremsstrahlung,    -1, 3, 3);

    } 
}

I have specified the G4Region and Cuts in the *DetectorConstruction.cc :

G4ProductionCuts* cuts = 0;

  if ( RegKobeDet != 0 )
    delete RegKobeDet;

  if ( RegKobeDet == 0 )
  { 
    RegKobeDet = new G4Region("KobeDetector");
    RegKobeDet->AddRootLogicalVolume(targetBlock_Al1_log);

    cuts = new G4ProductionCuts();
    RegKobeDet->SetProductionCuts(cuts);
  }

From my opinion I miss some physical process, because one 100 keV electron can not produce 69! AL27 isotopes. Isotopes goes through the detector without any interaction, thats also extremely strange.

Best regards, and great thanks for any tip:) Maciej

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

1 Note: Re: AL27 isotopes from 100keV electrons?   (Maciej Sznajder - 22 Dec, 2010)
2 Question: Re: AL27 isotopes from 100keV electrons?   (Maciej Sznajder - 04 Jan, 2011)
1 None: Re: AL27 isotopes from 100keV electrons?   (Marc Verderi - 04 Jan, 2011)
(_ None: Re: AL27 isotopes from 100keV electrons?   (Maciej Sznajder - 05 Jan, 2011)
(_ Idea: Re: AL27 isotopes from 100keV electrons?   (Vladimir Ivanchenko - 06 Jan, 2011)
(_ Note: Re: AL27 isotopes from 100keV electrons?   (Maciej Sznajder - 11 Jan, 2011)
(_ Idea: Re: AL27 isotopes from 100keV electrons?   (Vladimir Ivanchenko - 11 Jan, 2011)
(_ Question: Re: AL27 isotopes from 100keV electrons?   (Maciej Sznajder - 24 Jan, 2011)
 Add Message Add Message
to: "AL27 isotopes from 100keV 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 ]