Reconstruction/MdcPatRec/MdcRecoUtil/include/MdcRecoUtil/DifComplex.h

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//  $Id: DifComplex.h,v 1.2 2009/12/23 02:59:56 zhangy Exp $
//
// Description:
//  Class Header for |DifComplex|
//      Represent a complex diffential number
//
// Environment:
//  Software developed for the BaBar Detector at the SLAC B-Factory.
//
// Author List:
//  A. Snyder
//
// Copyright Information:
//  Copyright (C) 1997  SLAC
//
// History:
//  Migration for BESIII MDC
//
//------------------------------------------------------------------------Dif
 
#ifndef DifComplex_H
#define DifComplex_H
 
#include "MdcRecoUtil/DifNumber.h"
#include <assert.h>
#include <stdlib.h>
 
class ostream;
 
class DifComplex {
 
public:
  // constructors
  inline DifComplex();
  inline DifComplex( const DifNumber& r );
  inline DifComplex( double r );
  inline DifComplex( const DifNumber& r, const DifNumber& i );
  inline DifComplex( double r, double i );
  inline DifComplex( const DifComplex& c );
 
  // access
  inline const DifNumber& real() const { return _r; }
  inline DifNumber&       gimeReal() { return _r; }
  inline const DifNumber& imag() const { return _i; }
  inline DifNumber&       gimeImag() { return _i; }
 
  // set
  void setReal( const DifNumber& r ) { _r = r; }
  void setImag( const DifNumber& i ) { _i = i; }
 
  // operators and such
  inline DifComplex& operator+=( const DifComplex& a );
  inline DifComplex& operator-=( const DifComplex& a );
  inline DifComplex& operator*=( const DifComplex& a );
  inline DifComplex& operator/=( const DifComplex& a );
  inline DifComplex& operator/=( const DifNumber& a );
  inline DifComplex& operator=( const DifComplex& a );
  inline DifComplex  operator-() const;
 
  // manupulate *this guy - reverse polish style
  inline DifComplex& conjugate();                  // take complex conjugate
  inline DifComplex& inverse();                    // invert
  inline DifComplex& flipsign();                   // flip sign
  inline DifComplex& takeLog();                    // take log
  inline DifComplex& power( const DifNumber& p );  // raise to a power
  inline DifComplex& power( const DifComplex& p ); // raise to a power
 
  // complex specfic function
  inline DifNumber mag() const;                    // magnitude
  inline DifNumber magsq() const;                  // magnitude squared
  inline DifNumber arg() const;                    // phase
  inline DifNumber phase() const { return arg(); } // a rose by ...
 
  // io
  inline void print() const;
 
  // friends
  inline friend DifComplex operator+( const DifComplex& a, const DifComplex& b );
  inline friend DifComplex operator-( const DifComplex& a, const DifComplex& b );
  inline friend DifComplex operator*( const DifComplex& a, const DifComplex& b );
  inline friend DifComplex operator/( const DifComplex& a, const DifComplex& b );
  inline friend bool       operator==( const DifComplex& a, const DifComplex& b );
  inline friend bool       operator!=( const DifComplex& a, const DifComplex& b );
  inline friend DifComplex sin( const DifComplex& a );
  inline friend DifComplex cos( const DifComplex& a );
  inline friend DifComplex tan( const DifComplex& a );
  inline friend DifComplex sec( const DifComplex& a );
  inline friend DifComplex asin( const DifComplex& a );
  inline friend DifComplex acos( const DifComplex& a );
  inline friend DifComplex atan( const DifComplex& a );
  inline friend DifComplex atan2( const DifComplex& a, const DifComplex& b );
  inline friend DifComplex exp( const DifComplex& a );
  inline friend DifComplex cosh( const DifComplex& a );
  inline friend DifComplex sinh( const DifComplex& a );
  inline friend DifComplex tanh( const DifComplex& a );
  inline friend DifComplex sqrt( const DifComplex& a );
  inline friend DifComplex log( const DifComplex& a );
  inline friend DifComplex pow( const DifComplex& a, const DifComplex& p );
  inline friend DifComplex pow( const DifComplex& a, const DifNumber& p );
  inline friend DifComplex cc( const DifComplex& a );
  inline friend DifComplex fromPolar( const DifNumber& mag, const DifNumber& phase );
 
  // destructor
  inline virtual ~DifComplex(){};
 
private:
  // data
  DifNumber _r;
  DifNumber _i;
 
  // functions
  inline DifNumber&       gimeR() { return _r; }
  inline DifNumber&       gimeI() { return _i; }
  inline const DifNumber& R() const { return _r; }
  inline const DifNumber& I() const { return _i; }
};
 
#include "MdcRecoUtil/DifComplex.icc"
/*del**
inline ostream& operator<<(ostream& o,const DifComplex& c){
  c.print(o);
  return o;
}
 */
#endif