JMath.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <iomanip>
 
#include "TRandom3.h"
 
#include "JLang/JException.hh"
#include "JMath/JMath.hh"
#include "JMath/JMathToolkit.hh"
 
#include "Jeep/JPrint.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
 
 
namespace {
 
  using namespace JPP;
 
  /**
   * A random object with some arithmetic capabilities.
   */
  struct JObject :
    public JMath<JObject>
  {
  public:
    /**
     * Default constructor.
     */
    JObject() :
      x(0.0),
      y(0.0),
      z(0.0)
    {}
 
 
    /**
     * Constructor.
     *
     * \param  __x      x
     * \param  __y      y
     * \param  __z      z
     */
    JObject(const double __x,
            const double __y,
            const double __z) :
      x(__x),
      y(__y),
      z(__z)
    {}
 
 
    /**
     * Add object.
     *
     * \param  object   object
     * \return          this object
     */
    JObject& add(const JObject& object)
    {
      x += object.x;
      y += object.y;
      z += object.z;
 
      return *this;
    }
 
 
    /**
     * Scale object.
     *
     * \param  factor   multiplication factor
     * \return          this object
     */
    JObject& mul(const double factor)
    {
      x *= factor;
      y *= factor;
      z *= factor;
 
      return *this;
    }
 
 
    /**
     * Write object to output.
     *
     * \param  out           output stream
     * \param  object        object
     * \return               output stream
     */
    friend inline std::ostream& operator<<(std::ostream& out, const JObject& object)
    {
      using namespace std;
 
      out << showpos << FIXED(5,2) << object.x << ' ' 
          << showpos << FIXED(5,2) << object.y << ' ' 
          << showpos << FIXED(5,2) << object.z;
 
      return out;
    }
 
 
    double x;
    double y;
    double z;
  };
}
 
 
/**
 * \file
 *
 * Example program to test user class with arithmetic capabilities (based on class JMATH::JMath).
 * \author mdejong
 */
int main(int argc, char**argv)
{
  using namespace std;
 
  double alpha;
  int    debug;
 
  try {
 
    JParser<> zap("Example program to test user class with arithmetic capabilities.");
 
    zap['a'] = make_field(alpha);
    zap['d'] = make_field(debug)          = 0;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
  using namespace JPP;
 
 
  gRandom->SetSeed(0);
 
  const Double_t xmin = -1.0;
  const Double_t xmax = +1.0;
 
 
  const JObject A(gRandom->Uniform(xmin, xmax),
                  gRandom->Uniform(xmin, xmax),
                  gRandom->Uniform(xmin, xmax));
 
  const JObject B(gRandom->Uniform(xmin, xmax),
                  gRandom->Uniform(xmin, xmax),
                  gRandom->Uniform(xmin, xmax));
    
  cout << "A           " << A << endl;
  cout << "B           " << B << endl;
  cout << "interpolate " << interpolate(A, B, alpha) << endl;
 
  cout << "A + B     = " << A + B << endl;
 
  cout << "A * " << noshowpos << FIXED(5,2) << alpha << " = " << A * alpha << endl;
}