22.05.2014 Views

CERN Program Library Long Writeup W5013 - CERNLIB ...

CERN Program Library Long Writeup W5013 - CERNLIB ...

CERN Program Library Long Writeup W5013 - CERNLIB ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Geant 3.16 GEANT User’s Guide PHYS431<br />

Origin : L.Urbán Submitted: 12.03.82<br />

Revision : Revised: 16.12.93<br />

Documentation :<br />

1 Subroutines<br />

Ionisation processes for heavy ions<br />

For a description of the subroutines used to calculate the energy loss tables for protons, see [PHYS430].<br />

CALL GTHION<br />

GTHION transports heavy ions (A, Z > 1, particle type 8). For the moment the only discrete processes<br />

activated are δ-ray and Čerenkov photon generation.<br />

2 Method<br />

The mean energy loss of the heavy ions can be expressed as:<br />

where<br />

∆E = Q 2 S p<br />

( T<br />

A<br />

)<br />

ρt (1)<br />

Q<br />

( ) T<br />

S p<br />

A<br />

T,A<br />

ρ<br />

t<br />

charge of the ion;<br />

energy loss by the a proton with kinetic energy T p = T/A;<br />

kinetic energy and mass number of the ion;<br />

density of the medium;<br />

step length<br />

This formula is used to calculate the ionisation loss for all the charged hadrons.<br />

In the case of ions, the charge of the ion makes the problem more difficult. Electron-exchange processes<br />

with the atoms change the charge as the ion traverses the medium. The main features of this process can be<br />

summarised as follows:<br />

• Q loses the dependence on the initial ion charge after a very short initial step, usually ≪ 1µ for<br />

condensed media;<br />

• Q does depend on the velocity of the ion;<br />

• even in the case of constant velocity, Q fluctuates around a mean value;<br />

• Q features a small dependence on the medium, especially for small values of T .<br />

If we want to calculate the energy loss of the ion, we need a good parametrisation for the quantity Q = Q eff .<br />

We give here a relatively simple parametrisation [83, 93]:<br />

(<br />

)<br />

Q eff = Z I<br />

[1 − 1.034 − 0.1777e −0.08114 Z I<br />

e −v] (2)<br />

v = 121.4139 β<br />

Z 2/3<br />

I<br />

(<br />

+0.0378 sin<br />

190.7165 β<br />

Z 2/3<br />

I<br />

303 PHYS431 – 1<br />

)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!