68 Hep3Vector f3x(1.0), f3y(0.0, 1.0), f3z(0.0, 0.0, 1.0);
69 Hep3Vector d30, d3x(1.0), d3y(0.0, 1.0), d3z(0.0, 0.0, 1.0);
78 if ( !
test(d2, 1.0, 1.0, 0.0, std::sqrt(2.0),
DEPS) )
exit(1);
80 if ( !
test(d3, 1.0, 1.0, 1.0, std::sqrt(3.0),
DEPS) )
exit(1);
121 std::cin >> d0;
if ( !
test(d0, 1.1, 2.2, 3.3, 4.4,
DEPS) )
exit(1);
122 std::cin >> f0;
if ( !
test(f0, 4.0, 3.0, 2.0, 1.0,
FEPS) )
exit(1);
123 std::cout << d0 << std::endl;
124 std::cout << f0 << std::endl;
136 if ( !
test(d11, 3.0, 3.0, 3.0, std::sqrt(27.0),
FEPS) )
exit(4);
138 if ( !
test(f11, 3.0, 3.0, 3.0, std::sqrt(27.0),
FEPS) )
exit(4);
140 if ( !
test(d11, 4.0, 4.0, 4.0, std::sqrt(48.0),
FEPS) )
exit(4);
142 if ( !
test(f11, 4.0, 4.0, 4.0, std::sqrt(48.0),
FEPS) )
exit(4);
144 if ( !
test(d11, 1.0, 1.0, 1.0, std::sqrt(3.0),
FEPS) )
exit(4);
145 if ( !
test(-d11, -1.0, -1.0, -1.0, -std::sqrt(3.0),
FEPS) )
exit(4);
147 if ( !
test(f11, 1.0, 1.0, 1.0, std::sqrt(3.0),
FEPS) )
exit(4);
148 if ( !
test(-f11, -1.0, -1.0, -1.0, -std::sqrt(3.0),
FEPS) )
exit(4);
190 double M = (d1 + d2).mag();
191 double dm1 = d1.
mag();
192 double dm2 = d2.
mag();
194 d30 = -(d1 + d2).boostVector();
196 double phi = d1.
phi();
197 double theta = d1.
theta();
206 if ( !
test(d1, 0.0, 0.0, std::sqrt(p2), std::sqrt(p2 +
sqr(dm1)),
DEPS) )
exit(7);
207 if ( !
test(d2, 0.0, 0.0, -std::sqrt(p2), std::sqrt(p2 +
sqr(dm2)),
DEPS) )
exit(7);
221 if ( !
test(d3, 0.0, 0.0, std::sqrt(p2), std::sqrt(p2 +
sqr(dm1)),
DEPS) )
exit(7);
222 if ( !
test(d4, 0.0, 0.0, -std::sqrt(p2), std::sqrt(p2 +
sqr(dm2)),
DEPS) )
exit(7);
237 if ( !
test(d3, 0.0, 0.0, std::sqrt(p2), std::sqrt(p2 +
sqr(dm1)),
DEPS) )
exit(7);
238 if ( !
test(d4, 0.0, 0.0, -std::sqrt(p2), std::sqrt(p2 +
sqr(dm2)),
DEPS) )
exit(7);