Hi,
In my simulation I'm shooting CsI targets with electrons which can have energies between 29-120 keV.These electrons will emit a cascade of delta electrons producing clusters of eV secondary electrons because of the extraordinarily low work function (0.1-0.2 eV) of CsI.
The thickness of CsI is about 100 nm.
I think I will get many electrons leaving the CsI bulk,but I got a few.
So what's the problem?
Here is part of my PhysicsList:
defaultCutValue = 0.000001*mm;
cutForGamma = defaultCutValue;
cutForElectron = defaultCutValue;
cutForProton = 0.001*mm;
SetVerboseLevel(1);
G4ProductionCutsTable::GetProductionCutsTable()->SetEnergyRange(100.*eV, 1*GeV);
if (particleName == "gamma") {
// gamma
pmanager->AddDiscreteProcess(new G4PenelopePhotoElectric);
pmanager->AddDiscreteProcess(new G4PenelopeCompton);
pmanager->AddDiscreteProcess(new G4PenelopeGammaConversion);
pmanager->AddDiscreteProcess(new G4PenelopeRayleigh);
} else if (particleName == "e-") {
//electron
pmanager->AddProcess(new G4MultipleScattering, -1, 1, 1);
pmanager->AddProcess(new G4PenelopeIonisation, -1, 2, 2);
pmanager->AddProcess(new G4PenelopeBremsstrahlung, -1,-1, 3);
} else if (particleName == "e+") {
//positron
pmanager->AddProcess(new G4MultipleScattering, -1, 1, 1);
pmanager->AddProcess(new G4PenelopeIonisation, -1, 2, 2);
pmanager->AddProcess(new G4PenelopeBremsstrahlung, -1,-1, 3);
pmanager->AddProcess(new G4PenelopeAnnihilation, 0,-1, 4);
}
It seems when the low energy electrons are created in the CsI bulk most of them will not leave CsI bulk because of a very short range in CsI.
I also have used LowEnergy PhysicsList, but the result is still wrong.
what's the problem?
|