EvtVector3R.hh

Go to the documentation of this file.

//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of the EvtGen package developed jointly
//      for the BaBar and CLEO collaborations.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/COPYRIGHT
//      Copyright (C) 1998      Caltech, UCSB
//
// Module: EvtGen/EvtVector3R.hh
//
// Description: Class to describe real 3 vectors
//
// Modification history:
//
//    RYD    Sept. 5, 1997      Module created
//
//------------------------------------------------------------------------
 
#ifndef EVTVECTOR3R_HH
#define EVTVECTOR3R_HH
 
// #include <iostream.h>
// #include <math.h>
#include <iosfwd>
 
class EvtVector3R {
 
  friend EvtVector3R rotateEuler( const EvtVector3R& v, double phi, double theta, double ksi );
 
  inline friend EvtVector3R operator*( double c, const EvtVector3R& v2 );
  inline friend double      operator*( const EvtVector3R& v1, const EvtVector3R& v2 );
  inline friend EvtVector3R operator+( const EvtVector3R& v1, const EvtVector3R& v2 );
  inline friend EvtVector3R operator-( const EvtVector3R& v1, const EvtVector3R& v2 );
  inline friend EvtVector3R operator*( const EvtVector3R& v1, double c );
  inline friend EvtVector3R operator/( const EvtVector3R& v1, double c );
  friend EvtVector3R        cross( const EvtVector3R& v1, const EvtVector3R& v2 );
 
public:
  EvtVector3R();
  EvtVector3R( double x, double y, double z );
  virtual ~EvtVector3R();
  inline EvtVector3R&  operator*=( const double c );
  inline EvtVector3R&  operator/=( const double c );
  inline EvtVector3R&  operator+=( const EvtVector3R& v2 );
  inline EvtVector3R&  operator-=( const EvtVector3R& v2 );
  inline void          set( int i, double d );
  inline void          set( double x, double y, double z );
  void                 applyRotateEuler( double phi, double theta, double ksi );
  inline double        get( int i ) const;
  friend std::ostream& operator<<( std::ostream& s, const EvtVector3R& v );
  double               dot( const EvtVector3R& v2 );
  double               d3mag() const;
 
private:
  double v[3];
};
 
inline EvtVector3R& EvtVector3R::operator*=( const double c ) {
 
  v[0] *= c;
  v[1] *= c;
  v[2] *= c;
  return *this;
}
 
inline EvtVector3R& EvtVector3R::operator/=( const double c ) {
 
  v[0] /= c;
  v[1] /= c;
  v[2] /= c;
  return *this;
}
 
inline EvtVector3R& EvtVector3R::operator+=( const EvtVector3R& v2 ) {
 
  v[0] += v2.v[0];
  v[1] += v2.v[1];
  v[2] += v2.v[2];
  return *this;
}
 
inline EvtVector3R& EvtVector3R::operator-=( const EvtVector3R& v2 ) {
 
  v[0] -= v2.v[0];
  v[1] -= v2.v[1];
  v[2] -= v2.v[2];
  return *this;
}
 
inline EvtVector3R operator*( double c, const EvtVector3R& v2 ) {
 
  return EvtVector3R( v2 ) *= c;
}
 
inline EvtVector3R operator*( const EvtVector3R& v1, double c ) {
 
  return EvtVector3R( v1 ) *= c;
}
 
inline EvtVector3R operator/( const EvtVector3R& v1, double c ) {
 
  return EvtVector3R( v1 ) /= c;
}
 
inline double operator*( const EvtVector3R& v1, const EvtVector3R& v2 ) {
 
  return v1.v[0] * v2.v[0] + v1.v[1] * v2.v[1] + v1.v[2] * v2.v[2];
}
 
inline EvtVector3R operator+( const EvtVector3R& v1, const EvtVector3R& v2 ) {
 
  return EvtVector3R( v1 ) += v2;
}
 
inline EvtVector3R operator-( const EvtVector3R& v1, const EvtVector3R& v2 ) {
 
  return EvtVector3R( v1 ) -= v2;
}
 
inline double EvtVector3R::get( int i ) const { return v[i]; }
 
inline void EvtVector3R::set( int i, double d ) { v[i] = d; }
 
inline void EvtVector3R::set( double x, double y, double z ) {
 
  v[0] = x;
  v[1] = y;
  v[2] = z;
}
 
#endif