Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/externals/clhep/src/LorentzVectorB.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  1 // -*- C++ -*-
  2 // ---------------------------------------------------------------------------
  3 //
  4 // This file is a part of the CLHEP - a Class Library for High Energy Physics.
  5 //
  6 // This is the implementation of the HepLorentzVector class:
  7 // Those methods originating in ZOOM dealing with simple boosts and rotations.
  8 // Use of one of these methods will not force loading of the HepRotation or
  9 // HepLorentzRotation class.
 10 //
 11 
 12 #include "CLHEP/Vector/LorentzVector.h"
 13 
 14 #include <cmath>
 15 #include <iostream>
 16 
 17 namespace CLHEP  {
 18 
 19 //-*********
 20 // rotationOf()
 21 //-*********
 22 
 23 // Each of these is a shell over a rotate method.
 24 
 25 HepLorentzVector rotationXOf
 26   (const HepLorentzVector & vec, double phi){
 27   HepLorentzVector vv (vec);
 28   return vv.rotateX (phi);
 29 }
 30 
 31 HepLorentzVector rotationYOf
 32   (const HepLorentzVector & vec, double phi){
 33   HepLorentzVector vv (vec);
 34   return vv.rotateY (phi);
 35 }
 36 
 37 HepLorentzVector rotationZOf
 38   (const HepLorentzVector & vec, double phi){
 39   HepLorentzVector vv (vec);
 40   return vv.rotateZ (phi);
 41 }
 42 
 43 //-********
 44 // boost
 45 //-********
 46 
 47 HepLorentzVector & HepLorentzVector::boost 
 48       ( const Hep3Vector & aaxis,  double bbeta ) {
 49   if (bbeta==0) {
 50     return *this; // do nothing for a 0 boost
 51   }
 52   double r2 = aaxis.mag2();
 53   if ( r2 == 0 ) {
 54     std::cerr << "HepLorentzVector::boost() - "
 55       << "A zero vector used as axis defining a boost -- no boost done"
 56       << std::endl;
 57     return *this;
 58   } 
 59   double b2 = bbeta*bbeta;
 60   if (b2 >= 1) {
 61     std::cerr << "HepLorentzVector::boost() - "
 62       << "LorentzVector boosted with beta >= 1 (speed of light) -- \n"
 63       << "no boost done" << std::endl;
 64   } else {
 65     Hep3Vector u = aaxis.unit();
 66     double ggamma = std::sqrt(1./(1.-b2));
 67     double betaDotV = u.dot(pp)*bbeta;
 68     double tt = ee;
 69 
 70     ee = ggamma * (tt + betaDotV);
 71     pp += ( ((ggamma-1)/b2)*betaDotV*bbeta + ggamma*bbeta*tt ) * u;
 72     // Note:  I have verified the behavior of this even when beta is very
 73     //        small -- (gamma-1)/b2 becomes inaccurate by O(1), but it is then
 74     //        multiplied by O(beta**2) and added to an O(beta) term, so the
 75     //        inaccuracy does not affect the final result.
 76   }
 77   return *this;
 78 } /* boost (axis, beta) */
 79 
 80 }  // namespace CLHEP
 81