29 (
double bx,
double by,
double bz) {
30 double bp2 = bx*bx + by*by + bz*bz;
33 "Boost Vector supplied to set HepLorentzRotation represents speed >= c."));
35 double gamma = 1.0 / std::sqrt(1.0 - bp2);
36 double bgamma = gamma * gamma / (1.0 + gamma);
37 mxx = 1.0 + bgamma * bx * bx;
38 myy = 1.0 + bgamma * by * by;
39 mzz = 1.0 + bgamma * bz * bz;
40 mxy =
myx = bgamma * bx * by;
41 mxz =
mzx = bgamma * bx * bz;
42 myz =
mzy = bgamma * by * bz;
240 std::setw(11) << std::setprecision(6) <<
xx() <<
" " <<
241 std::setw(11) << std::setprecision(6) <<
xy() <<
" " <<
242 std::setw(11) << std::setprecision(6) <<
xz() <<
" " <<
243 std::setw(11) << std::setprecision(6) <<
xt() <<
")\n"
245 std::setw(11) << std::setprecision(6) <<
yx() <<
" " <<
246 std::setw(11) << std::setprecision(6) <<
yy() <<
" " <<
247 std::setw(11) << std::setprecision(6) <<
yz() <<
" " <<
248 std::setw(11) << std::setprecision(6) <<
yt() <<
")\n"
250 std::setw(11) << std::setprecision(6) <<
zx() <<
" " <<
251 std::setw(11) << std::setprecision(6) <<
zy() <<
" " <<
252 std::setw(11) << std::setprecision(6) <<
zz() <<
" " <<
253 std::setw(11) << std::setprecision(6) <<
zt() <<
")\n"
255 std::setw(11) << std::setprecision(6) <<
tx() <<
" " <<
256 std::setw(11) << std::setprecision(6) <<
ty() <<
" " <<
257 std::setw(11) << std::setprecision(6) <<
tz() <<
" " <<
258 std::setw(11) << std::setprecision(6) <<
tt() <<
") ]\n";