| Geant4 Cross Reference |
1 // -*- C++ -*- 1 // -*- C++ -*-
2 // ------------------------------------------- 2 // ---------------------------------------------------------------------------
3 // 3 //
4 // This file is a part of the CLHEP - a Class 4 // This file is a part of the CLHEP - a Class Library for High Energy Physics.
5 // 5 //
6 // This is the implementation of methods of th 6 // This is the implementation of methods of the HepRotation class which
7 // were introduced when ZOOM PhysicsVectors wa 7 // were introduced when ZOOM PhysicsVectors was merged in, which might cause
8 // pulling in of LorentzTransformation related 8 // pulling in of LorentzTransformation related code units.
9 // 9 //
10 10
11 #include "CLHEP/Vector/Rotation.h" 11 #include "CLHEP/Vector/Rotation.h"
12 #include "CLHEP/Vector/LorentzRotation.h" 12 #include "CLHEP/Vector/LorentzRotation.h"
13 13
14 #include <cmath> 14 #include <cmath>
15 15
16 namespace CLHEP { 16 namespace CLHEP {
17 17
18 // ---------- distance2 and related member fu 18 // ---------- distance2 and related member functions:
19 // 19 //
20 // WHy do we have forms for HepLorentzRotation 20 // WHy do we have forms for HepLorentzRotation and HepBoost but not for
21 // HepBoostX, HepBoostY, HepBoostZ? Because t 21 // HepBoostX, HepBoostY, HepBoostZ? Because the latter can be gotten by
22 // implicit conversion to HepBoost; but if we 22 // implicit conversion to HepBoost; but if we just had HepLorentzRotation
23 // then this would involve double conversion w 23 // then this would involve double conversion when HepBoostX was used.
24 24
25 double HepRotation::distance2( const HepLorent 25 double HepRotation::distance2( const HepLorentzRotation & lt ) const {
26 HepAxisAngle a; 26 HepAxisAngle a;
27 Hep3Vector b; 27 Hep3Vector b;
28 lt.decompose(b, a); 28 lt.decompose(b, a);
29 double bet = b.beta(); 29 double bet = b.beta();
30 double bet2 = bet*bet; 30 double bet2 = bet*bet;
31 HepRotation r(a); 31 HepRotation r(a);
32 return bet2/(1-bet2) + distance2(r); 32 return bet2/(1-bet2) + distance2(r);
33 } 33 }
34 34
35 double HepRotation::distance2( const HepBoost 35 double HepRotation::distance2( const HepBoost & lt ) const {
36 return distance2( HepLorentzRotation(lt)); 36 return distance2( HepLorentzRotation(lt));
37 } 37 }
38 38
39 double HepRotation::howNear( const HepLorentzR 39 double HepRotation::howNear( const HepLorentzRotation & lt ) const {
40 return std::sqrt( distance2( lt ) ); 40 return std::sqrt( distance2( lt ) );
41 } 41 }
42 42
43 double HepRotation::howNear( const HepBoost & 43 double HepRotation::howNear( const HepBoost & lt ) const {
44 return std::sqrt( distance2( lt ) ); 44 return std::sqrt( distance2( lt ) );
45 } 45 }
46 46
47 bool HepRotation::isNear( const HepLorentzRo 47 bool HepRotation::isNear( const HepLorentzRotation & lt,
48 double ep 48 double epsilon) const {
49 return distance2( lt ) <= epsilon*epsilon; 49 return distance2( lt ) <= epsilon*epsilon;
50 } 50 }
51 51
52 bool HepRotation::isNear( const HepBoost & l 52 bool HepRotation::isNear( const HepBoost & lt,
53 double ep 53 double epsilon) const {
54 return distance2( lt ) <= epsilon*epsilon; 54 return distance2( lt ) <= epsilon*epsilon;
55 } 55 }
56 56
57 } // namespace CLHEP 57 } // namespace CLHEP
58 58
59 59