JMultiHistogram.hh

Go to the documentation of this file.

#ifndef __JTOOLS__JMULTIHISTOGRAM__
#define __JTOOLS__JMULTIHISTOGRAM__

#include <stdarg.h>

#include "JTools/JMultiMap.hh"
#include "JTools/JHistogram.hh"
#include "JTools/JArray.hh"


/**
 * \author mdejong
 */

namespace JTOOLS {}
namespace JPP { using namespace JTOOLS; }

namespace JTOOLS {


  /**
   * Multidimensional histogram.
   *
   * The template parameters respectively refer to:
   * -# histogram object used for the lowest dimension(s);
   * -# list of histogram maps used for the higher dimension(s); and
   * -# distance operator.
   *
   * The number of dimensions of this histogram object is equal to the length of
   * the given map list plus the extra dimensions of the first histogram object.
   */
  template<class JHistogram_t,
           class JMaplist_t,
           class JDistance_t = JDistance<typename JHistogram_t::abscissa_type> >
  class JMultiHistogram :
    public JMultiMap<typename JHistogram_t::abscissa_type, JHistogram_t, JMaplist_t, JDistance_t>
  {
  public:

    typedef JMultiMap<typename JHistogram_t::abscissa_type,
                      JHistogram_t,
                      JMaplist_t,
                      JDistance_t>                                        multimap_type;

    enum { NUMBER_OF_DIMENSIONS = JMapLength<JMaplist_t>::value + JHistogram_t::NUMBER_OF_DIMENSIONS };

    typedef typename multimap_type::abscissa_type                         abscissa_type;
    typedef typename multimap_type::ordinate_type                         ordinate_type;

    typedef typename multimap_type::const_iterator                        const_iterator;
    typedef typename multimap_type::const_reverse_iterator                const_reverse_iterator;
    typedef typename multimap_type::iterator                              iterator;
    typedef typename multimap_type::reverse_iterator                      reverse_iterator;

    typedef typename multimap_type::super_iterator                        super_iterator;
    typedef typename multimap_type::super_const_iterator                  super_const_iterator;

    typedef typename JHistogram_t::contents_type                          contents_type;


    /**
     * Default constructor.
     */
    JMultiHistogram()
    {}


    /**
     * Get multidimensional histogram.
     *
     * \return                      this multidimensional histogram
     */
    const JMultiHistogram& getMultiHistogram() const
    {
      return static_cast<const JMultiHistogram&>(*this);
    }


    /**
     * Get multidimensional histogram.
     *
     * \return                      this multidimensional histogram
     */
    JMultiHistogram& getMultiHistogram()
    {
      return static_cast<JMultiHistogram&>(*this);
    }


    /**
     * Multi-dimensional fill method call.
     *
     * \param  x      comma seperated argument list
     */
    void fill(const abscissa_type x, ...)
    {
      va_start(ap, x);

      buffer[0] = x;

      for (int i = 1; i != NUMBER_OF_DIMENSIONS; ++i) {
        buffer[i] = va_arg(ap, abscissa_type);
      }

      const contents_type w = va_arg(ap, contents_type);

      va_end(ap);

      this->evaluate(buffer.begin(), w);
    }

    
  protected:
    mutable va_list ap;
    mutable JArray<NUMBER_OF_DIMENSIONS, abscissa_type> buffer;
  };
}

#endif