ABObj.h

Go to the documentation of this file.

// @(#)root/minuit2:$Id$
// 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 &) = delete;
 
public:
   ABObj(const LAVector &obj) : fObject(obj), fFactor(double(1.)) {}
 
   ABObj(const LAVector &obj, double factor) : fObject(obj), fFactor(factor) {}
 
   ~ABObj() {}
 
   // remove copy constructor 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