TEvtGen/EvtGenBase/EvtVector3R.hh

   1 //--------------------------------------------------------------------------
   2 //
   3 // Environment:
   4 //      This software is part of the EvtGen package developed jointly
   5 //      for the BaBar and CLEO collaborations.  If you use all or part
   6 //      of it, please give an appropriate acknowledgement.
   7 //
   8 // Copyright Information: See EvtGen/COPYRIGHT
   9 //      Copyright (C) 1998      Caltech, UCSB
  10 //
  11 // Module: EvtGen/EvtVector3R.hh
  12 //
  13 // Description: Class to describe real 3 vectors
  14 //
  15 // Modification history:
  16 //
  17 //    RYD    Sept. 5, 1997      Module created
  18 //
  19 //------------------------------------------------------------------------
  20
  21 #ifndef EVTVECTOR3R_HH
  22 #define EVTVECTOR3R_HH
  23
  24 #include <iosfwd>
  25
  26 class EvtVector3R {
  27
  28   friend EvtVector3R rotateEuler(const EvtVector3R& v,
  29                                  double phi,double theta,double ksi);
  30
  31   inline friend EvtVector3R operator*(double c,const EvtVector3R& v2);
  32   inline friend double operator*(const EvtVector3R& v1,const EvtVector3R& v2);
  33   inline friend EvtVector3R operator+(const EvtVector3R& v1,const EvtVector3R& v2);
  34   inline friend EvtVector3R operator-(const EvtVector3R& v1,const EvtVector3R& v2);
  35   inline friend EvtVector3R operator*(const EvtVector3R& v1,double c);
  36   inline friend EvtVector3R operator/(const EvtVector3R& v1,double c);
  37   friend EvtVector3R cross(const EvtVector3R& v1,const EvtVector3R& v2);
  38
  39 public:
  40   EvtVector3R();
  41   EvtVector3R(double x,double y ,double z);
  42   virtual ~EvtVector3R();
  43   inline EvtVector3R& operator*=(const double c);
  44   inline EvtVector3R& operator/=(const double c);
  45   inline EvtVector3R& operator+=(const EvtVector3R& v2);
  46   inline EvtVector3R& operator-=(const EvtVector3R& v2);
  47   inline void set(int i,double d);
  48   inline void set(double x,double y ,double z);
  49   void applyRotateEuler(double phi,double theta,double ksi);
  50   inline double get(int i) const;
  51   friend std::ostream& operator<<(std::ostream& s,const EvtVector3R& v);
  52   double dot(const EvtVector3R& v2);
  53   double d3mag() const;
  54
  55 private:
  56
  57   double v[3];
  58
  59 };
  60
  61 inline EvtVector3R& EvtVector3R::operator*=(const double c){
  62
  63   v[0]*=c;
  64   v[1]*=c;
  65   v[2]*=c;
  66   return *this;
  67 }
  68
  69 inline EvtVector3R& EvtVector3R::operator/=(const double c){
  70
  71   v[0]/=c;
  72   v[1]/=c;
  73   v[2]/=c;
  74   return *this;
  75 }
  76
  77 inline EvtVector3R& EvtVector3R::operator+=(const EvtVector3R& v2){
  78
  79   v[0]+=v2.v[0];
  80   v[1]+=v2.v[1];
  81   v[2]+=v2.v[2];
  82   return *this;
  83 }
  84
  85 inline EvtVector3R& EvtVector3R::operator-=(const EvtVector3R& v2){
  86
  87   v[0]-=v2.v[0];
  88   v[1]-=v2.v[1];
  89   v[2]-=v2.v[2];
  90   return *this;
  91 }
  92
  93 inline EvtVector3R operator*(double c,const EvtVector3R& v2){
  94
  95   return EvtVector3R(v2)*=c;
  96 }
  97
  98 inline EvtVector3R operator*(const EvtVector3R& v1,double c){
  99
 100   return EvtVector3R(v1)*=c;
 101 }
 102
 103 inline EvtVector3R operator/(const EvtVector3R& v1,double c){
 104
 105   return EvtVector3R(v1)/=c;
 106 }
 107
 108 inline double operator*(const EvtVector3R& v1,const EvtVector3R& v2){
 109
 110   return v1.v[0]*v2.v[0]+v1.v[1]*v2.v[1]+v1.v[2]*v2.v[2];
 111 }
 112
 113 inline EvtVector3R operator+(const EvtVector3R& v1,const EvtVector3R& v2) {
 114
 115   return EvtVector3R(v1)+=v2;
 116 }
 117
 118 inline EvtVector3R operator-(const EvtVector3R& v1,const EvtVector3R& v2) {
 119
 120   return EvtVector3R(v1)-=v2;
 121
 122 }
 123
 124 inline double EvtVector3R::get(int i) const {
 125   return v[i];
 126 }
 127
 128 inline void EvtVector3R::set(int i,double d){
 129
 130   v[i]=d;
 131 }
 132
 133 inline void EvtVector3R::set(double x,double y, double z){
 134
 135   v[0]=x;
 136   v[1]=y;
 137   v[2]=z;
 138 }
 139
 140 #endif
 141