Problem 1761 - pixe mechanism in Pelelope
Summary: pixe mechanism in Pelelope
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/electromagnetic/lowenergy (show other problems)
Version: other
Hardware: All All
: P5 normal
Assignee: Luciano Pandola
URL:
Depends on:
Blocks:
 
Reported: 2015-06-23 15:04 CEST by michel.maire
Modified: 2015-06-24 09:00 CEST (History)
0 users

See Also:


Attachments
x-rays spectrum with penelope/ionization (7.71 KB, image/gif)
2015-06-23 15:04 CEST, michel.maire
Details

Note You need to log in before you can comment on or make changes to this problem.
Description michel.maire 2015-06-23 15:04:43 CEST
Created attachment 350 [details]
x-rays spectrum with penelope/ionization

in Penelope, Ionization process triggers automatically the atomic deexcitation cascade (pixe mechanism).
If, in addition users ask explicitly for /process/em/pixe true, there is a double generation of deexcitation cascade.
In attachement, a spectrum of generated x-rays, without and with PIXE activated.
 
 # macro file for TestEm5.cc    
 #     		    
 /control/verbose 2
 /run/verbose 2
 #
 /testem/det/setAbsMat Gold      
 /testem/det/setAbsThick 1 um
 /testem/det/setAbsYZ   10 um
 #
 /testem/phys/addPhysics empenelope
 
 /run/setCut 0.1 um
 #
 /run/initialize
 #
 /process/inactivate eBrem
 /process/inactivate compt
 /process/inactivate phot
 #
 /process/em/pixe true
 #
 /testem/gun/setDefault
 /gun/particle e- 
 /gun/energy 50 keV
 # 
 /analysis/setFileName penelope-nobrem-pixe
 /analysis/h1/set  3 100 0.01 50 keV	       #gamma: energy at vertex
 /analysis/h1/set  5 100 0.1 50 keV log10      #gamma: energy at vertex (log10) 	
 #
 /run/printProgress 100000
 # 
 /run/beamOn 1000000
Comment 1 Sebastien Incerti 2015-06-23 15:21:18 CEST
Reassigned.
Comment 2 Luciano Pandola 2015-06-23 15:37:07 CEST
Thanks for the report. The bug is confirmed. In the present configuration the PenelopeIonisation model produces the correct results *only when PIXE=off*

Unlike the Livermore and Standard ionization models, the PenelopeIonisation model has a built-in management of the atomic de-excitation, since it is able to sample the shell which is ionized and to follows the de-excitation cascade. This is the same approach which is followed for the gamma processes (photoelectric, Compton). So, x-rays are always produced (provided that the Fluo flag is on), even if the PIXE flag is off.

The Livermore and standard ionization models do not sample individually the target shell which is ionized. They sample the de-excitation statistically by calculating the shell ionization cross section. This is done using the PIXE common interface and it is in place only if the PIXE flag is on. Therefore, if PIXE=false, the Penelope model will show the x-rays, while Livermore/Standard will not. 

Unfortunately, due to a bug, the Penelope model produces the de-excitation twice if the PIXE flag is on (one through the built-in mechanism and one through the PIXE interface).
Comment 3 Luciano Pandola 2015-06-24 09:00:01 CEST
Fixed.

The Penelope model will use its own internal mechanism for atomic de-excitation when PIXE=off (thus producing the x-rays), while will use the PIXE interface when PIXE=on.