/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/MdcPatRec/MdcRecoUtil/MdcRecoUtil-00-01-08/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 <assert.h>
#include <stdlib.h>
#include "MdcRecoUtil/DifNumber.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

---