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.

• the mean number of tries is not too large.<br />

Thus the different possible decompositions of the distribution f(x) are not equivalent from the practical<br />

point of view (e.g. they can be very different in computational speed) and it can be very useful to optimise the<br />

decomposition. A remark of practical importance is that if our distribution is not normalised ( ∫ f(x)dx =<br />

C>0; C ≈ 1), the method can be used in the same manner, the only difference is that the mean number of<br />

tries in this case is given by ∑ i α i/C .<br />

2.2 Differential cross-section for pair production<br />

The Bethe-Heitler differential cross-section with the Coulomb correction for a photon of energy E to produce<br />

a e − e + -pair with one of the particles having energy ɛE (ɛ is the fraction of the photon energy carried<br />

by one particle of the pair) is given as in [12]:<br />

dσ(Z, E, ɛ)<br />

dɛ<br />

{<br />

= r2 0αZ[Z + ξ(Z)]<br />

E 2<br />

[ɛ 2 +(1− ɛ) 2 ][Φ 1 (δ) − F (Z)<br />

2<br />

]+ 2 3 ɛ(1 − ɛ)[Φ 2(δ) − F (Z) }<br />

] (2)<br />

2<br />

where Φ i (δ) are the screening functions depending on the screening variable δ<br />

δ = 136m 1<br />

Z 1/3 E ɛ(1 − ɛ)<br />

}<br />

Φ 1 (δ) =20.867 − 3.242δ +0.625δ 2<br />

Φ 2 (δ) =20.209 − 1.930δ − 0.086δ 2 δ ≤ 1<br />

Φ 1 (δ) =Φ 2 (δ) =21.12 − 4.184 ln(δ +0.952) δ>1<br />

m = electron mass<br />

F (Z) =<br />

{<br />

8/3lnZ<br />

8/3lnZ +8f c (Z)<br />

E

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

Saved successfully!

Ooh no, something went wrong!