- Timestamp:
- Dec 16, 2009, 12:14:47 PM (15 years ago)
- Location:
- trunk/documents/UserDoc/DocBookUsersGuides/PhysicsReferenceManual/latex
- Files:
-
- 20 added
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/documents/UserDoc/DocBookUsersGuides/PhysicsReferenceManual/latex/electromagnetic/lowenergy/hadrons.tex
r1211 r1222 785 785 786 786 787 \subsection{ICRU 73-based energy loss model} 788 The ICRU 73 \cite{hlei.ICRU73} report contains stopping power tables 789 for ions with atomic numbers 3--18 and 26, covering a range of different 790 elemental and compound target materials. The stopping powers derive from 791 calculations with the PASS code \cite{hlei.sigm02}, which implements the 792 binary stopping theory described in \cite{hlei.sigm02,hlei.sigm00}. Tables 793 in ICRU 73 extend over an energy range up to 1 GeV/nucleon. All stopping 794 powers were incorporated into Geant4 and are available through a 795 parameterisation model ({\tt G4IonParametrisedLossModel}). For a few 796 materials revised stopping powers were included (water, water vapor, nylon type 797 6 and 6/6 from P. Sigmund et al \cite{hlei.sigm09a} and copper from P. Sigmund 798 \cite{hlei.sigm09b}), which replace the corresponding tables of the original 799 ICRU 73 report. 800 801 To account for secondary electron production above $T_{c}$, the continuous 802 energy loss per unit path length is calculated according to 803 \begin{equation} 804 \label{hlei.rstp} 805 \frac{dE}{dx}\bigg|_{T<T_C} = \bigg(\frac{dE}{dx}\bigg)_{ICRU73} - 806 \bigg(\frac{dE}{dx}\bigg)_{\delta} 807 \end{equation} 808 where $(dE/dx)_{ICRU73}$ refers to stopping powers obtained by interpolating 809 ICRU 73 tables and $(dE/dx)_{\delta}$ is the mean energy transferred to 810 $\delta$-electrons per path length given by 811 \begin{equation} 812 \bigg(\frac{dE}{dx}\bigg)_{\delta} = \sum_{i} n_{at,i} \int_{T_c}^{T_{max}} 813 \frac{d\sigma_i(T)}{dT} T dT 814 \label{} 815 \end{equation} 816 where the index $i$ runs over all elements composing the material, $n_{at,i}$ 817 is the number of atoms of the element $i$ per volume, $T_{max}$ is the maximum 818 energy transferable to an electron according to formula (\ref{hlei.a1}) and 819 $d\sigma_i/dT$ specifies the differential cross section per atom for producing 820 an $\delta$-electron following equation (\ref{hlei.bbb}). 821 822 For compound targets not considered in the ICRU 73 report, the first term on 823 the rightern side in equation (\ref{hlei.rstp}) is computed by applying Bragg's 824 additivity rule \cite{hlei.ICRU49} if tables for all elemental components are 825 available in ICRU 73. 826 827 828 787 829 \subsection{Status of this document} 788 830 … … 793 835 19.01.2002 Minor corrections (mma) \\ 794 836 13.05.2002 Minor corrections (V.Ivanchenko) \\ 795 28.08.2002 Minor corrections (V.Ivanchenko) 837 28.08.2002 Minor corrections (V.Ivanchenko) \\ 838 11.12.2009 Modified by A. Lechner to add ICRU 73 section 796 839 797 840 \begin{latexonly} … … 862 905 \bibitem{hlei.Gryzinski1} M. Gryzinski, Phys. Rev. A 135 (1965) 305. 863 906 \bibitem{hlei.Gryzinski2} M. Gryzinski, Phys. Rev. A 138 (1965) 322. 907 \bibitem{hlei.ICRU73} 908 Stopping of Ions Heavier Than Helium, 909 ICRU Report 73, Oxford University Press (2005). 910 \bibitem{hlei.sigm02} 911 P.~Sigmund and A.~Schinner, 912 Nucl. Instr. Meth. in Phys. Res. B 195 (2002) 64. 913 \bibitem{hlei.sigm00} 914 P.~Sigmund and A.~Schinner, 915 Eur. Phys. J. D 12 (2000) 425. 916 \bibitem{hlei.sigm09a} 917 P.~Sigmund, A.~Schinner and H.~Paul, 918 Errata and Addenda for ICRU Report 73, Stopping of Ions Heavier 919 than Helium (2009). 920 \bibitem{hlei.sigm09b} 921 Personal communication with P.~Sigmund (2009). 864 922 \end{thebibliography} 865 923 … … 932 990 \item M. Gryzinski, Phys. Rev. A 135 (1965) 305. 933 991 \item M. Gryzinski, Phys. Rev. A 138 (1965) 322. 992 \item 993 Stopping of Ions Heavier Than Helium, 994 ICRU Report 73, Oxford University Press (2005). 995 \item 996 P.~Sigmund and A.~Schinner, 997 Nucl. Instr. Meth. in Phys. Res. B 195 (2002) 64. 998 \item 999 P.~Sigmund and A.~Schinner, 1000 Eur. Phys. J. D 12 (2000) 425. 1001 \item 1002 P.~Sigmund, A.~Schinner and H.~Paul, 1003 Errata and Addenda for ICRU Report 73, Stopping of Ions Heavier than Helium (2009). 1004 \item 1005 Personal communication with P.~Sigmund (2009). 1006 934 1007 \end{enumerate} 935 1008 -
trunk/documents/UserDoc/DocBookUsersGuides/PhysicsReferenceManual/latex/electromagnetic/miscellaneous/setcuts.tex
r1211 r1222 18 18 one geometry is valid also for the other geometries. 19 19 20 The value of cut is defined only for electrons, positrons and gamma.20 The value of cut is defined for electrons, positrons, gamma and protons. 21 21 At the beginning of initialization of Geant4 physics the conversion from unique 22 22 {\it cut in range} to cuts in kinetic energy for each {\it G4MaterialCutsCouple} … … 40 40 the corresponding kinetic energy can be found out. If obtained $cut in energy$ 41 41 cannot be below 42 the parameter $lowlimit$ (default $1 \; keV$) .42 the parameter $lowlimit$ (default $1 \; keV$) and above $highlimit$ (default $10 \; GeV$). 43 43 If in specific application lower cut is required, 44 44 then the allowed energy cut needs to be extended:\\ 45 45 \\ 46 46 {\it \footnotesize G4ProductionCutsTable::GetProductionCutsTable()$\to$SetEnergyRange(lowlimit,highlimit);} 47 47 or via UI commands 48 $$/cuts/setMinCutEnergy\;\; 100\;\; eV$$ 49 $$/cuts/setMaxCutEnergy\;\; 100\;\; TeV$$ 48 50 In contrary to electrons, gammas has no range, so some approximation should 49 51 be used for range to energy conversion. … … 73 75 {\it highlimit}. 74 76 77 The cut for proton is introduced with Geant4 v9.3. The main goal 78 of this cut is to limit production of all recoil ions including protons 79 in elastic scattering 80 processes. A simple linear conversion formula is used to compute energy threshold from the value 81 of cut in range, in particular, the cut in range $1~mm$ corresponds 82 to the production threshold $100 keV$. 83 75 84 The conversion from range to energy can be studied using {\it G4EmCalculator} 76 85 class. This class allows access or recalculation of energy loss, ranges and … … 90 99 04.12.04 minor re-wording by D.H. Wright \\ 91 100 18.05.07 rewritten by V. Ivanchenko \\ 92 11.12.08 minor revision by V. Ivanchenko \\ 101 11.12.08 minor revision by V. Ivanchenko, Geant4 v9.2 \\ 102 11.12.09 minor revision by V. Ivanchenko, Geant4 v9.3 \\ 93 103 94 104 \begin{latexonly} -
trunk/documents/UserDoc/DocBookUsersGuides/PhysicsReferenceManual/latex/electromagnetic/standard/hion.tex
r1211 r1222 244 244 \subsubsection{Nuclear Stopping} 245 245 246 For scaled energies below $T_{lim} = 2 MeV$ the contribution of non-ionizing energy loss 247 needs to be taken into account. The additional energy loss due to {\it nuclear stopping power} 248 $\Delta T_N \Delta s$ is added the the energy loss. The process {\it G4ionIonisation} 249 has a flag, which allows to switch on or off this correction. For that the method\\ 250 \\ 251 {\it G4ionIonisation::ActivateNuclearStopping(G4bool)} 252 \\ 253 \\ 254 can be used. By default 255 this correction is active and the ICRU'49 parameterisation \cite{hion.ICRU49} is used. 246 Nuclear stopping due to elastic ion-ion scattering since Geant4 v9.3 247 can be simulated with the continuous process 248 {\it G4NuclearStopping}. By default 249 this correction is active and the ICRU'49 parameterisation \cite{hion.ICRU49} is used, 250 which is implemented in the model class {\it G4ICRU49NuclearStoppingModel}. 256 251 257 252 -
trunk/documents/UserDoc/DocBookUsersGuides/PhysicsReferenceManual/latex/electromagnetic/standard/msc.tex
r1211 r1222 15 15 \item 16 16 G4MuMultipleScattering. 17 \item18 G4MultipleScattering;19 17 \end{itemize} 20 18 For concrete simulation the {\it G4VMscModel} interface is used, which is an extension … … 22 20 \begin{itemize} 23 21 \item 24 G4UrbanMscModel ;22 G4UrbanMscModel90 - Geant4 v9.0 applied to muons, hadrons and ions; 25 23 \item 26 G4UrbanMscModel 2;24 G4UrbanMscModel92 - Geant4 v9.2 (current default) applied for electron and positron; 27 25 \item 28 G4UrbanMscModel9 0;26 G4UrbanMscModel93 - Geant4 v9.3 applied for electron and positron for Option2, Option3 and other EM Physics Lists; 29 27 \item 30 G4WentzelVIModel. 28 G4GoudsmitSaundersonModel - for electrons and positrons (beta-version); 29 \item 30 G4WentzelVIModel - for muons and hadrons, should be included in Physics Lists together with G4CoulombScattering process; 31 31 \end{itemize} 32 The last model is not yet in the production mode, so below we will concentrate 33 on models developed by L.~Urban \cite{msc.urban}. Today default model for 34 electrons and positrons is {\it G4UrbanMscModel2}, for hadrons - {\it G4UrbanMscModel90}. 32 The last models are not yet in the production mode, so below we will concentrate 33 on models developed by L.~Urban \cite{msc.urban}. 35 34 36 35 \subsection{Introduction} … … 811 810 08.12.08 revised by L. Urb\'an, for Geant4 V9.2 \\ 812 811 11.12.08 minor revision by V. Ivanchenko \\ 812 11.12.09 minor revision by V. Ivanchenko, for Geant4 v9.3 \\ 813 813 814 814 \begin{latexonly}
Note: See TracChangeset
for help on using the changeset viewer.