// @(#)root/minuit2:$Id: ABObj.h 34019 2010-06-21 15:19:37Z rdm $
// Authors: M. Winkler, F. James, L. Moneta, A. Zsenei   2003-2005  

/**********************************************************************
 *                                                                    *
 * Copyright (c) 2005 LCG ROOT Math team,  CERN/PH-SFT                *
 *                                                                    *
 **********************************************************************/

#ifndef ROOT_Minuit2_ABObj
#define ROOT_Minuit2_ABObj

#include "Minuit2/ABTypes.h"

namespace ROOT {

   namespace Minuit2 {


template<class mtype, class M, class T>
class ABObj {

public:

  typedef mtype Type;

private:

  ABObj() : fObject(M()), fFactor(T(0.)) {}

  ABObj& operator=(const ABObj&) {return *this;}

  template<class a, class b, class c>
  ABObj(const ABObj<a,b,c>&) : fObject(M()), fFactor(T(0.)) {}

  template<class a, class b, class c>
  ABObj& operator=(const ABObj<a,b,c>&) {return *this;}
  
public:

  ABObj(const M& obj) : fObject(obj), fFactor(T(1.)) {}

  ABObj(const M& obj, T factor) : fObject(obj), fFactor(factor) {}

  ~ABObj() {}

  ABObj(const ABObj& obj) : 
    fObject(obj.fObject), fFactor(obj.fFactor) {}

  template<class b, class c>
  ABObj(const ABObj<mtype,b,c>& obj) : 
     fObject(M(obj.Obj() )), fFactor(T(obj.f() )) {}

  const M& Obj() const {return fObject;}

  T f() const {return fFactor;}

private:

  M fObject;
  T fFactor;
};

class LAVector;
template <> class ABObj<vec, LAVector, double> {

public:

  typedef vec Type;

private:

  ABObj& operator=(const ABObj&) {return *this;}
  
public:

  ABObj(const LAVector& obj) : fObject(obj), fFactor(double(1.)) {}

  ABObj(const LAVector& obj, double factor) : fObject(obj), fFactor(factor) {}

  ~ABObj() {}

  // remove copy constructure to Fix a problem in AIX 
  // should be able to use the compiler generated one
//   ABObj(const ABObj& obj) : 
//     fObject(obj.fObject), fFactor(obj.fFactor) {}

  template<class c>
  ABObj(const ABObj<vec,LAVector,c>& obj) : 
    fObject(obj.fObject), fFactor(double(obj.fFactor)) {}

  const LAVector& Obj() const {return fObject;}

  double f() const {return fFactor;}

private:

  const LAVector& fObject;
  double fFactor;
};

class LASymMatrix;
template <> class ABObj<sym, LASymMatrix, double> {

public:

  typedef sym Type;

private:

  ABObj& operator=(const ABObj&) {return *this;}
  
public:

  ABObj(const LASymMatrix& obj) : fObject(obj), fFactor(double(1.)) {}

  ABObj(const LASymMatrix& obj, double factor) : fObject(obj), fFactor(factor) {}

  ~ABObj() {}

  ABObj(const ABObj& obj) : 
    fObject(obj.fObject), fFactor(obj.fFactor) {}

  template<class c>
  ABObj(const ABObj<vec,LASymMatrix,c>& obj) : 
    fObject(obj.fObject), fFactor(double(obj.fFactor)) {}

  const LASymMatrix& Obj() const {return fObject;}

  double f() const {return fFactor;}

private:

  const LASymMatrix& fObject;
  double fFactor;
};

// templated scaling operator *
template<class mt, class M, class T>
inline ABObj<mt, M, T> operator*(T f, const M& obj) {
  return ABObj<mt, M, T>(obj, f);
}

// templated operator /
template<class mt, class M, class T>
inline ABObj<mt, M, T> operator/(const M& obj, T f) {
  return ABObj<mt, M, T>(obj, T(1.)/f);
}

// templated unary operator -
template<class mt, class M, class T>
inline ABObj<mt,M,T> operator-(const M& obj) {
  return ABObj<mt,M,T>(obj, T(-1.));
}

/*
// specialization for LAVector

inline ABObj<vec, LAVector, double> operator*(double f, const LAVector& obj) {
  return ABObj<vec, LAVector, double>(obj, f);
}

inline ABObj<vec, LAVector, double> operator/(const LAVector& obj, double f) {
  return ABObj<vec, LAVector, double>(obj, double(1.)/f);
}

inline ABObj<vec,LAVector,double> operator-(const LAVector& obj) {
  return ABObj<vec,LAVector,double>(obj, double(-1.));
}
*/

  }  // namespace Minuit2

}  // namespace ROOT

#endif  // ROOT_Minuit2_ABObj
