JToolsSupportkit.hh

Go to the documentation of this file.

#ifndef __JROOT__JTOOLSSUPPORTKIT__
#define __JROOT__JTOOLSSUPPORTKIT__
 
#include "TH1.h"
#include "TH2.h"
#include "TH3.h"
 
#include "JTools/JElement.hh"
#include "JTools/JMapList.hh"
#include "JTools/JDistance.hh"
#include "JTools/JHistogram1D.hh"
#include "JTools/JMultiHistogram.hh"
#include "JTools/JHistogramMap_t.hh"
 
 
/**
 * \author bjung
 */
 
namespace JROOT {}
namespace JPP { using namespace JROOT; }
 
namespace JROOT {
 
  using JTOOLS::JDistance;
  using JTOOLS::JBin1D;
  using JTOOLS::JHistogram1D;
  using JTOOLS::JMapList;
  using JTOOLS::JMAPLIST;
  using JTOOLS::JMultiHistogram;
  using JTOOLS::JHistogramGridMap_t;
 
 
  /**
   * Auxiliary class for copying ROOT histograms to Jpp histograms.
   */
  class JHistogramHelper
  {
  public:
    
    /**
     * Default constructor.
     */
    JHistogramHelper()
    {}
    
 
    /**
     * Copy histogram from <tt>from</tt> to <tt>to</tt>.
     *
     * \param  from        ROOT histogram
     * \param  to          Jpp  histogram
     */
    template<class JAbscissa_t,
             class JOrdinate_t,
             template<class, class> class JContainer_t,
             class JDistance_t = JDistance<JAbscissa_t> >
    static void copy(const TH1& from, JHistogram1D<JBin1D<JAbscissa_t, JOrdinate_t>, JContainer_t, JDistance_t>& to)
    {
      to.clear();
 
      const Int_t NbinsX = from.GetXaxis()->GetNbins();
      
      to.underflow = from.GetBinContent(0);
      to.overflow  = from.GetBinContent(NbinsX+1);
      to.integral  = from.Integral(0, NbinsX+1);
 
      for (Int_t i = 1; i <= NbinsX; ++i) {
 
        const double x0 = from.GetXaxis()->GetBinLowEdge(i);
        const double xc = from.GetXaxis()->GetBinCenter (i);      
      
        const double y  = from.GetBinContent(i);
        const double ye = from.GetBinError  (i);
 
        JBin1D<JAbscissa_t, JOrdinate_t> bin(x0, y, xc*y, ye*ye);
      
        to.insert(bin);
      }
    }
 
 
    /**
     * Copy histogram from <tt>from</tt> to <tt>to</tt>.
     *
     * \param  from        ROOT histogram
     * \param  to          Jpp  histogram
     */
    template<class JAbscissa_t,
             class JOrdinate_t,
             template<class, class> class JContainer_t,
             template<class, class, class> class JMap_t = JHistogramGridMap_t,
             class JDistance_t = JDistance<JAbscissa_t> >
    static void copy(const TH2& from, JMultiHistogram<JHistogram1D<JBin1D<JAbscissa_t, JOrdinate_t>, JContainer_t, JDistance_t>, JMapList<JMap_t>, JDistance_t>& to)
    {
      typedef JHistogram1D<JBin1D<JAbscissa_t, JOrdinate_t>, JContainer_t, JDistance_t>   JHistogram1D_t;
 
      const Int_t NbinsX = from.GetXaxis()->GetNbins();
      const Int_t NbinsY = from.GetYaxis()->GetNbins();      
      
      for (Int_t i = 0; i <= NbinsX+1; ++i) {
 
        JHistogram1D_t histogram;
 
        const double xc = from.GetXaxis()->GetBinCenter(i);
 
        for (Int_t j = 0; j <= NbinsY+1; ++j) {
 
          const Int_t N = from.GetBin(i,j);
 
          const double y0 = from.GetYaxis()->GetBinLowEdge(j);
          const double yc = from.GetYaxis()->GetBinCenter (j);
          
          const double z  = from.GetBinContent(N);
          const double ze = from.GetBinError  (N);
 
          histogram.integral    += z;
          
          if        (from.IsBinUnderflow(N)) {
            histogram.underflow += z;
          } else if (from.IsBinOverflow(N)) {
            histogram.overflow  += z;
          } else {
            JBin1D<JAbscissa_t, JOrdinate_t> bin(y0, z, yc*z, ze*ze);     
            histogram.insert(bin);
          }
        }
 
        if (histogram.getSize() > 0) {
          to[xc] = histogram;
        }
      }
    }
 
 
    /**
     * Copy histogram from <tt>from</tt> to <tt>to</tt>.
     *
     * \param  from        ROOT histogram
     * \param  to          Jpp  histogram
     */
    template<class JAbscissa_t,
             class JOrdinate_t,
             template<class, class> class JContainer_t,
             template<class, class, class> class JMap_t1 = JHistogramGridMap_t,
             template<class, class, class> class JMap_t2 = JHistogramGridMap_t,
             class JDistance_t = JDistance<JAbscissa_t> >
    static void copy(const TH3& from, JMultiHistogram<JHistogram1D<JBin1D<JAbscissa_t, JOrdinate_t>, JContainer_t, JDistance_t>, typename JMAPLIST<JMap_t1, JMap_t2>::maplist, JDistance_t>& to)
    {
      typedef JHistogram1D<JBin1D<JAbscissa_t, JOrdinate_t>, JContainer_t, JDistance_t>   JHistogram1D_t;
      
      const Int_t NbinsX = from.GetXaxis()->GetNbins();
      const Int_t NbinsY = from.GetYaxis()->GetNbins();
      const Int_t NbinsZ = from.GetZaxis()->GetNbins();
 
      for (Int_t i = 0; i <= NbinsX+1; ++i) {
 
        const double xc = from.GetXaxis()->GetBinCenter(i);
 
        for (Int_t j = 0; j <= NbinsY+1; ++j) {
 
          JHistogram1D_t histogram;
 
          const double yc = from.GetYaxis()->GetBinCenter(j);
 
          for (Int_t k = 0; k <= NbinsZ+1; ++k) {
 
            const Int_t N = from.GetBin(i,j,k);
 
            const double z0 = from.GetYaxis()->GetBinLowEdge(j);
            const double zc = from.GetYaxis()->GetBinCenter (j);
          
            const double w  = from.GetBinContent(N);
            const double we = from.GetBinError  (N);
 
            histogram.integral    += w;
 
            if        (from.IsBinUnderflow(N)) {
              histogram.underflow += w;
            } else if (from.IsBinOverflow(N)) {
              histogram.overflow  += w;
            } else {
              JBin1D<JAbscissa_t, JOrdinate_t> bin(z0, w, zc*w, we*we);
              histogram.insert(bin);
            }
          }
          
          if (histogram.getSize() > 0) {
            to[xc][yc] = histogram;
          }
        }
      }
    }
 
    
    /**
     * Copy histogram from <tt>from</tt> to <tt>to</tt>.
     *
     * \param  from        ROOT histogram
     * \param  to          Jpp  histogram
     */
    template<class JAbscissa_t,
             class JOrdinate_t,
             template<class, class> class JContainer_t,
             class JDistance_t = JDistance<JAbscissa_t> >
    void operator()(const TH1& from, JHistogram1D<JBin1D<JAbscissa_t, JOrdinate_t>, JContainer_t, JDistance_t>& to) const
    {
      return copy(from, to);
    }
 
 
    /**
     * Copy histogram from <tt>from</tt> to <tt>to</tt>.
     *
     * \param  from        ROOT histogram
     * \param  to          Jpp  histogram
     */
    template<class JAbscissa_t,
             class JOrdinate_t,
             template<class, class> class JContainer_t,
             template<class, class, class> class JMap_t = JHistogramGridMap_t,
             class JDistance_t = JDistance<JAbscissa_t> >
    void operator()(const TH2& from, JMultiHistogram<JHistogram1D<JBin1D<JAbscissa_t, JOrdinate_t>, JContainer_t, JDistance_t>, JMapList<JMap_t>, JDistance_t>& to) const
    {
      return copy(from, to);
    }
    
 
    /**
     * Copy histogram from <tt>from</tt> to <tt>to</tt>.
     *
     * \param  from        ROOT histogram
     * \param  to          Jpp  histogram
     */
    template<class JAbscissa_t,
             class JOrdinate_t,
             template<class, class> class JContainer_t,
             template<class, class, class> class JMap_t1 = JHistogramGridMap_t,
             template<class, class, class> class JMap_t2 = JHistogramGridMap_t,
             class JDistance_t = JDistance<JAbscissa_t> >
    void operator()(const TH3& from, JMultiHistogram<JHistogram1D<JBin1D<JAbscissa_t, JOrdinate_t>, JContainer_t, JDistance_t>, typename JMAPLIST<JMap_t1, JMap_t2>::maplist, JDistance_t>& to) const
    {
      return copy(from, to);
    }
  };
 
 
  /**
   * Function object for copying histograms.
   */
  static const JHistogramHelper copyHistogram;
}
 
#endif