ABSum.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_ABSum
#define ROOT_Minuit2_ABSum
 
#include "Minuit2/ABObj.h"
 
namespace ROOT {
 
   namespace Minuit2 {
 
 
template<class M1, class M2>
class ABSum {
 
private:
 
  ABSum() : fA(M1()), fB(M2()) {}
 
  ABSum& operator=(const ABSum&) {return *this;}
 
  template<class MI1, class MI2>
  ABSum& operator=(const ABSum<MI1,MI2>&) {return *this;}
 
public:
 
  ABSum(const M1& a, const M2& b): fA(a), fB(b) {}
 
  ~ABSum() {}
 
  ABSum(const ABSum& sum) : fA(sum.fA), fB(sum.fB) {}
 
  template<class MI1, class MI2>
  ABSum(const ABSum<MI1,MI2>& sum) : fA(M1(sum.A() )), fB(M2(sum.B() )) {}
 
  const M1& A() const {return fA;}
  const M2& B() const {return fB;}
 
private:
 
  M1 fA;
  M2 fB;
};
 
// ABObj + ABObj
template<class atype, class A, class btype, class B, class T>
inline ABObj<typename AlgebraicSumType<atype, btype>::Type, ABSum<ABObj<atype,A,T>, ABObj<btype,B,T> >,T> operator+(const ABObj<atype,A,T>& a, const ABObj<btype,B,T>& b) {
 
  return ABObj<typename AlgebraicSumType<atype,btype>::Type, ABSum<ABObj<atype,A,T>, ABObj<btype,B,T> >,T>(ABSum<ABObj<atype,A,T>, ABObj<btype,B,T> >(a, b));
}
 
// ABObj - ABObj
template<class atype, class A, class btype, class B, class T>
inline ABObj<typename AlgebraicSumType<atype, btype>::Type, ABSum<ABObj<atype,A,T>, ABObj<btype,B,T> >,T> operator-(const ABObj<atype,A,T>& a, const ABObj<btype,B,T>& b) {
 
  return ABObj<typename AlgebraicSumType<atype,btype>::Type, ABSum<ABObj<atype,A,T>, ABObj<btype,B,T> >,T>(ABSum<ABObj<atype,A,T>, ABObj<btype,B,T> >(a, ABObj<btype,B,T>(b.Obj(), T(-1.)*b.f())));
}
 
  }  // namespace Minuit2
 
}  // namespace ROOT
 
#endif  // ROOT_Minuit2_ABSum