Geant4 Cross Reference |
1 // -*- C++ -*- 1 // -*- C++ -*- 2 // 2 // 3 // This file is a part of the CLHEP - a Class 3 // This file is a part of the CLHEP - a Class Library for High Energy Physics. 4 // 4 // 5 // This is the definition of the HepBoostX cla 5 // This is the definition of the HepBoostX class for performing specialized 6 // Lorentz transformations which are pure boos 6 // Lorentz transformations which are pure boosts in the X direction, on 7 // objects of the HepLorentzVector class. 7 // objects of the HepLorentzVector class. 8 // 8 // 9 // HepLorentzRotation is a concrete implementa 9 // HepLorentzRotation is a concrete implementation of Hep4RotationInterface. 10 // 10 // 11 // .SS See Also 11 // .SS See Also 12 // RotationInterfaces.h 12 // RotationInterfaces.h 13 // LorentzVector.h LorentzRotation.h 13 // LorentzVector.h LorentzRotation.h 14 // Boost.h 14 // Boost.h 15 // 15 // 16 // .SS Author 16 // .SS Author 17 // Mark Fischler 17 // Mark Fischler 18 18 19 #ifndef HEP_BOOSTX_H 19 #ifndef HEP_BOOSTX_H 20 #define HEP_BOOSTX_H 20 #define HEP_BOOSTX_H 21 21 >> 22 #ifdef GNUPRAGMA >> 23 #pragma interface >> 24 #endif >> 25 22 #include "CLHEP/Vector/RotationInterfaces.h" 26 #include "CLHEP/Vector/RotationInterfaces.h" 23 #include "CLHEP/Vector/LorentzVector.h" 27 #include "CLHEP/Vector/LorentzVector.h" 24 28 25 namespace CLHEP { 29 namespace CLHEP { 26 30 27 // Declarations of classes and global methods 31 // Declarations of classes and global methods 28 class HepBoostX; 32 class HepBoostX; 29 inline HepBoostX inverseOf ( const HepBoostX & 33 inline HepBoostX inverseOf ( const HepBoostX & b ); 30 class HepBoost; 34 class HepBoost; 31 class HepRotation; 35 class HepRotation; 32 36 33 /** 37 /** 34 * @author 38 * @author 35 * @ingroup vector 39 * @ingroup vector 36 */ 40 */ 37 class HepBoostX { 41 class HepBoostX { 38 42 39 public: 43 public: 40 44 41 // ---------- Constructors and Assignment: 45 // ---------- Constructors and Assignment: 42 46 43 inline HepBoostX(); 47 inline HepBoostX(); 44 // Default constructor. Gives a boost of 0. 48 // Default constructor. Gives a boost of 0. 45 49 46 inline HepBoostX(const HepBoostX & b); 50 inline HepBoostX(const HepBoostX & b); 47 inline HepBoostX(HepBoostX && b) = default; << 51 // Copy constructor. 48 // Copy and move constructors. << 49 52 50 inline HepBoostX & operator = (const HepBoos 53 inline HepBoostX & operator = (const HepBoostX & m); 51 inline HepBoostX & operator = (HepBoostX && << 54 // Assignment. 52 // Copy and move assignment operators. << 53 55 54 HepBoostX & set (double beta); 56 HepBoostX & set (double beta); 55 inline HepBoostX (double beta); 57 inline HepBoostX (double beta); 56 // Constructor from beta 58 // Constructor from beta 57 59 58 // ---------- Accessors: 60 // ---------- Accessors: 59 61 60 inline double beta() const; 62 inline double beta() const; 61 inline double gamma() const; 63 inline double gamma() const; 62 inline Hep3Vector boostVector() const; 64 inline Hep3Vector boostVector() const; 63 inline Hep3Vector getDirection() const; 65 inline Hep3Vector getDirection() const; 64 66 65 inline double xx() const; 67 inline double xx() const; 66 inline double xy() const; 68 inline double xy() const; 67 inline double xz() const; 69 inline double xz() const; 68 inline double xt() const; 70 inline double xt() const; 69 inline double yx() const; 71 inline double yx() const; 70 inline double yy() const; 72 inline double yy() const; 71 inline double yz() const; 73 inline double yz() const; 72 inline double yt() const; 74 inline double yt() const; 73 inline double zx() const; 75 inline double zx() const; 74 inline double zy() const; 76 inline double zy() const; 75 inline double zz() const; 77 inline double zz() const; 76 inline double zt() const; 78 inline double zt() const; 77 inline double tx() const; 79 inline double tx() const; 78 inline double ty() const; 80 inline double ty() const; 79 inline double tz() const; 81 inline double tz() const; 80 inline double tt() const; 82 inline double tt() const; 81 // Elements of the matrix. 83 // Elements of the matrix. 82 84 83 inline HepLorentzVector col1() const; 85 inline HepLorentzVector col1() const; 84 inline HepLorentzVector col2() const; 86 inline HepLorentzVector col2() const; 85 inline HepLorentzVector col3() const; 87 inline HepLorentzVector col3() const; 86 inline HepLorentzVector col4() const; 88 inline HepLorentzVector col4() const; 87 // orthosymplectic column vectors 89 // orthosymplectic column vectors 88 90 89 inline HepLorentzVector row1() const; 91 inline HepLorentzVector row1() const; 90 inline HepLorentzVector row2() const; 92 inline HepLorentzVector row2() const; 91 inline HepLorentzVector row3() const; 93 inline HepLorentzVector row3() const; 92 inline HepLorentzVector row4() const; 94 inline HepLorentzVector row4() const; 93 // orthosymplectic row vectors 95 // orthosymplectic row vectors 94 96 95 HepRep4x4 rep4x4() const; 97 HepRep4x4 rep4x4() const; 96 // 4x4 representation: 98 // 4x4 representation: 97 99 98 HepRep4x4Symmetric rep4x4Symmetric() const; 100 HepRep4x4Symmetric rep4x4Symmetric() const; 99 // Symmetric 4x4 representation. 101 // Symmetric 4x4 representation. 100 102 101 // ---------- Decomposition: 103 // ---------- Decomposition: 102 104 103 void decompose (HepRotation & rotation, Hep 105 void decompose (HepRotation & rotation, HepBoost & boost) const; 104 void decompose (HepAxisAngle & rotation, Hep 106 void decompose (HepAxisAngle & rotation, Hep3Vector & boost) const; 105 // Find R and B such that L = R*B -- trivial 107 // Find R and B such that L = R*B -- trivial, since R is identity 106 108 107 void decompose ( HepBoost & boost, HepRotat 109 void decompose ( HepBoost & boost, HepRotation & rotation) const; 108 void decompose (Hep3Vector & boost, HepAxisA 110 void decompose (Hep3Vector & boost, HepAxisAngle & rotation) const; 109 // Find R and B such that L = B*R -- trivial 111 // Find R and B such that L = B*R -- trivial, since R is identity 110 112 111 // ---------- Comparisons: 113 // ---------- Comparisons: 112 114 113 inline int compare( const HepBoostX & b ) c 115 inline int compare( const HepBoostX & b ) const; 114 // Dictionary-order comparison, in order of 116 // Dictionary-order comparison, in order of beta. 115 // Used in operator<, >, <=, >= 117 // Used in operator<, >, <=, >= 116 118 117 inline bool operator == (const HepBoostX & b 119 inline bool operator == (const HepBoostX & b) const; 118 inline bool operator != (const HepBoostX & b 120 inline bool operator != (const HepBoostX & b) const; 119 inline bool operator <= (const HepBoostX & b 121 inline bool operator <= (const HepBoostX & b) const; 120 inline bool operator >= (const HepBoostX & b 122 inline bool operator >= (const HepBoostX & b) const; 121 inline bool operator < (const HepBoostX & b 123 inline bool operator < (const HepBoostX & b) const; 122 inline bool operator > (const HepBoostX & b 124 inline bool operator > (const HepBoostX & b) const; 123 // Comparisons. 125 // Comparisons. 124 126 125 inline bool isIdentity() const; 127 inline bool isIdentity() const; 126 // Returns true if a null boost. 128 // Returns true if a null boost. 127 129 128 inline double distance2( const HepBoostX & b 130 inline double distance2( const HepBoostX & b ) const; 129 double distance2( const HepBoost & b 131 double distance2( const HepBoost & b ) const; 130 // Defined as the distance2 between the vect 132 // Defined as the distance2 between the vectors (gamma*betaVector) 131 133 132 double distance2( const HepRotation & r ) c 134 double distance2( const HepRotation & r ) const; 133 double distance2( const HepLorentzRotation & 135 double distance2( const HepLorentzRotation & lt ) const; 134 // Decompose lt = B*R; add norm2 to distance 136 // Decompose lt = B*R; add norm2 to distance2 to between boosts. 135 137 136 inline double howNear( const HepBoostX & b 138 inline double howNear( const HepBoostX & b ) const; 137 inline double howNear( const HepBoost & b 139 inline double howNear( const HepBoost & b ) const; 138 inline double howNear( const HepRotation & 140 inline double howNear( const HepRotation & r ) const; 139 inline double howNear( const HepLorentzRot 141 inline double howNear( const HepLorentzRotation & lt ) const; 140 142 141 inline bool isNear( const HepBoostX & b, 143 inline bool isNear( const HepBoostX & b, 142 double epsilon=Hep4RotationInterf 144 double epsilon=Hep4RotationInterface::tolerance) const; 143 inline bool isNear( const HepBoost & b, 145 inline bool isNear( const HepBoost & b, 144 double epsilon=Hep4RotationInterf 146 double epsilon=Hep4RotationInterface::tolerance) const; 145 bool isNear( const HepRotation & r, 147 bool isNear( const HepRotation & r, 146 double epsilon=Hep4RotationInterf 148 double epsilon=Hep4RotationInterface::tolerance) const; 147 bool isNear( const HepLorentzRotation & lt 149 bool isNear( const HepLorentzRotation & lt, 148 double epsilon=Hep4RotationInterf 150 double epsilon=Hep4RotationInterface::tolerance) const; 149 151 150 // ---------- Properties: 152 // ---------- Properties: 151 153 152 inline double norm2() const; 154 inline double norm2() const; 153 // distance2 (IDENTITY), which is beta^2 * g 155 // distance2 (IDENTITY), which is beta^2 * gamma^2 154 156 155 void rectify(); 157 void rectify(); 156 // sets according to the stored beta 158 // sets according to the stored beta 157 159 158 // ---------- Application: 160 // ---------- Application: 159 161 160 inline HepLorentzVector operator()( const He 162 inline HepLorentzVector operator()( const HepLorentzVector & w ) const; 161 // Transform a Lorentz Vector. 163 // Transform a Lorentz Vector. 162 164 163 inline HepLorentzVector operator* ( const He 165 inline HepLorentzVector operator* ( const HepLorentzVector & w ) const; 164 // Multiplication with a Lorentz Vector. 166 // Multiplication with a Lorentz Vector. 165 167 166 // ---------- Operations in the group of 4-R 168 // ---------- Operations in the group of 4-Rotations 167 169 168 HepBoostX operator * (const HepBoostX & b) 170 HepBoostX operator * (const HepBoostX & b) const; 169 HepLorentzRotation operator * (const HepBoos 171 HepLorentzRotation operator * (const HepBoost & b) const; 170 HepLorentzRotation operator * (const HepRota 172 HepLorentzRotation operator * (const HepRotation & r) const; 171 HepLorentzRotation operator * (const HepLore 173 HepLorentzRotation operator * (const HepLorentzRotation & lt) const; 172 // Product of two Lorentz Rotations (this) * 174 // Product of two Lorentz Rotations (this) * lt - matrix multiplication 173 // Notice that the product of two pure boost 175 // Notice that the product of two pure boosts in different directions 174 // is no longer a pure boost. 176 // is no longer a pure boost. 175 177 176 inline HepBoostX inverse() const; 178 inline HepBoostX inverse() const; 177 // Return the inverse. 179 // Return the inverse. 178 180 179 inline friend HepBoostX inverseOf ( const He 181 inline friend HepBoostX inverseOf ( const HepBoostX & b ); 180 // global methods to invert. 182 // global methods to invert. 181 183 182 inline HepBoostX & invert(); 184 inline HepBoostX & invert(); 183 // Inverts the Boost matrix. 185 // Inverts the Boost matrix. 184 186 185 // ---------- I/O: 187 // ---------- I/O: 186 188 187 std::ostream & print( std::ostream & os ) co 189 std::ostream & print( std::ostream & os ) const; 188 // Output form is BOOSTX (beta=..., gamma=.. 190 // Output form is BOOSTX (beta=..., gamma=...); 189 191 190 // ---------- Tolerance 192 // ---------- Tolerance 191 193 192 static inline double getTolerance(); 194 static inline double getTolerance(); 193 static inline double setTolerance(double tol 195 static inline double setTolerance(double tol); 194 196 195 protected: 197 protected: 196 198 197 inline HepLorentzVector vectorMultiplication 199 inline HepLorentzVector vectorMultiplication 198 ( cons 200 ( const HepLorentzVector & w ) const; 199 // Multiplication with a Lorentz Vector. 201 // Multiplication with a Lorentz Vector. 200 202 201 HepLorentzRotation matrixMultiplication (con 203 HepLorentzRotation matrixMultiplication (const HepRep4x4 & m) const; 202 HepLorentzRotation matrixMultiplication (con 204 HepLorentzRotation matrixMultiplication (const HepRep4x4Symmetric & m) const; 203 205 204 inline HepBoostX (double beta, double gamma) 206 inline HepBoostX (double beta, double gamma); 205 207 206 double beta_; 208 double beta_; 207 double gamma_; 209 double gamma_; 208 210 209 }; // HepBoostX 211 }; // HepBoostX 210 212 211 inline 213 inline 212 std::ostream & operator << 214 std::ostream & operator << 213 ( std::ostream & os, const HepBoostX& b ) {r 215 ( std::ostream & os, const HepBoostX& b ) {return b.print(os);} 214 216 215 } // namespace CLHEP 217 } // namespace CLHEP 216 218 217 #include "CLHEP/Vector/BoostX.icc" 219 #include "CLHEP/Vector/BoostX.icc" 218 220 219 #endif /* HEP_BOOSTX_H */ 221 #endif /* HEP_BOOSTX_H */ 220 222