02.05.2014 Views

PENELOPE 2003 - OECD Nuclear Energy Agency

PENELOPE 2003 - OECD Nuclear Energy Agency

PENELOPE 2003 - OECD Nuclear Energy Agency

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

5.1. Rotations and translations 155<br />

the geometry on the computer screen. These programs use specific computer graphics<br />

software and, therefore, they are not portable. The executable files included in the<br />

penelope distribution package run on personal computers under Microsoft Windows;<br />

they are simple and effective tools for debugging the geometry definition file.<br />

5.1 Rotations and translations<br />

The definition of various parts of the material system (quadric surfaces in reduced form<br />

and modules) involves rotations and translations. To describe these transformations,<br />

we shall adopt the active point of view: the reference frame remains fixed and only the<br />

space points (vectors) are translated or rotated.<br />

In what follows, and in the computer programs, all lengths are in cm. The position<br />

and direction of movement of a particle are referred to the laboratory coordinate system,<br />

a Cartesian reference frame which is defined by the position of its origin of coordinates<br />

and the unit vectors ˆx = (1, 0, 0), ŷ = (0, 1, 0) and ẑ = (0, 0, 1) along the directions of<br />

its axes.<br />

A translation T (t), defined by the displacement vector t = (t x , t y , t z ), transforms<br />

the vector r = (x, y, z) into<br />

T (t) r = r + t = (x + t x , y + t y , z + t z ). (5.1)<br />

Evidently, the inverse translation T −1 (t) corresponds to the displacement vector −t,<br />

i.e. T −1 (t) = T (−t).<br />

A rotation R is defined through the Euler angles ω, θ and φ, which specify a sequence<br />

of rotations about the coordinate axes 1 : first a rotation of angle ω about the z-axis,<br />

followed by a rotation of angle θ about the y-axis and, finally, a rotation of angle φ<br />

about the z-axis. A positive rotation about a given axis would carry a right-handed<br />

screw in the positive direction along that axis. Positive (negative) angles define positive<br />

(negative) rotations.<br />

The rotation R(ω, θ, φ) transforms the vector r = (x, y, z) into a vector<br />

whose coordinates are given by<br />

r ′ = R(ω, θ, φ) r = (x ′ , y ′ , z ′ ), (5.2)<br />

⎛<br />

x ′ ⎞<br />

⎛<br />

⎜<br />

⎝ y ′<br />

⎟<br />

⎠ = R(ω, θ, φ) ⎜<br />

⎝<br />

z ′<br />

⎞<br />

x<br />

y<br />

z<br />

⎟<br />

⎠ , (5.3)<br />

1 This definition of the Euler angles is the one usually adopted in Quantum Mechanics (see e.g.<br />

Edmonds, 1960).

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

Saved successfully!

Ooh no, something went wrong!