ComplexProduct.h

Go to the documentation of this file.

#ifndef MINT_COMPLEX_PRODUCT_HH
#define MINT_COMPLEX_PRODUCT_HH
// author: Jonas Rademacker (Jonas.Rademacker@bristol.ac.uk)
// status:  Mon 9 Feb 2009 19:17:55 GMT

#include "Mint/IComplexFitParDependent.h"
#include "Mint/BasicComplex.h"
#include <vector>
#include <complex>
#include "Mint/counted_ptr.h"
#include "Mint/FitParDependent.h"

namespace MINT{
  class ComplexProduct : virtual public IComplexFitParDependent, public FitParDependent{
  protected:
    std::vector< counted_ptr<IReturnComplex> > _facVec; 
    std::vector< double > _fixedDoubleVec; 
    std::vector< std::complex<double> > _fixedComplexVec; 
    
    inline double fixedRealProduct(){
      double prod(1); 
      for(unsigned i = 0; i < _fixedDoubleVec.size(); i++){
    prod *= _fixedDoubleVec[i];
      }      
      return prod;
    }
    inline std::complex<double> fixedComplexProduct(){
      std::complex<double> prod(1, 0); 
      for(unsigned i = 0; i < _fixedComplexVec.size(); i++){
    prod *= _fixedComplexVec[i];
      }      
      return prod;
    }
    
    inline std::complex<double> complexProduct(){
      std::complex<double> prod(1,0); 
      for(unsigned i = 0; i < _facVec.size(); i++){
    prod *= (_facVec[i])->ComplexVal();
      }      
      return prod;
    }
     
    void setToUnity();

  public:
    ComplexProduct(IFitParRegister* daddy=0);
    ComplexProduct(double initVal, IFitParRegister* daddy=0);
    ComplexProduct(const std::complex<double>& z, IFitParRegister* daddy=0);
    ComplexProduct(const ComplexProduct& other, IFitParRegister* newDaddy=0);

    void multiply(const ComplexProduct& multiplyWith);

    void addTerm(double val); // by value
    void addTerm(const std::complex<double>& z); // by value
    void addTerm(const counted_ptr<IReturnComplex>& irc); // by reference
    void addTerm(const counted_ptr<IComplexFitParDependent>& irc); // by reference

    ComplexProduct& operator*=(double val);
    ComplexProduct& operator*=(const std::complex<double>& z);
    ComplexProduct& operator*=(const counted_ptr<IComplexFitParDependent>& irc);
    ComplexProduct& operator*=(const counted_ptr<IReturnComplex>& irc);
    ComplexProduct& operator*=(const ComplexProduct& cp);

   inline std::complex<double> ComplexVal(){
     return complexProduct() *  fixedComplexProduct() * fixedRealProduct();
   }
    
  };
}
#endif
//