JResultTransformer.hh

Go to the documentation of this file.

#ifndef __JTOOLS__JRESULTTRANSFORMER__
#define __JTOOLS__JRESULTTRANSFORMER__

#include "JTools/JMultiMapTransformer.hh"
#include "JTools/JResult.hh"


/**
 * \author mdejong
 */

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

namespace JTOOLS {

  /**
   * Auxiliary class to handle multidimensional map transformations for given result type.
   */
  template<class JResult_t>
  struct JResultTransformer {

    typedef JResult_t                    result_type;

    template<unsigned int N, class JArgument_t>
    static result_type transform(const JMultiMapTransformer<N, JArgument_t>&                     transformer,
                                 typename JMultiMapTransformer<N, JArgument_t>::const_array_type array,
                                 const result_type&                                              result) 
    {
      return result;
    }
  };


  /**
   * Auxiliary class to handle multidimensional map transformations for given result type.
   */
  template<class JResult_t>
  struct JResultTransformer< JResultHesse<JResult_t> > {

    typedef JResultHesse<JResult_t>      result_type;

    template<unsigned int N, class JArgument_t>
    static result_type transform(const JMultiMapTransformer<N, JArgument_t>&                     transformer,
                                 typename JMultiMapTransformer<N, JArgument_t>::const_array_type array,
                                 const result_type&                                              result) 
    {
      const JArgument_t z = transformer.getXn(array, 1.0) - transformer.getXn(array, 0.0);
       
      result_type value(result);

      value.fp /= z;

      return value;
    }
  };


  /**
   * Auxiliary class to handle multidimensional map transformations for given result type.
   */
  template<class JResult_t>
  struct JResultTransformer< JResultPDF<JResult_t> > {

    typedef JResultPDF<JResult_t>        result_type;

    template<unsigned int N, class JArgument_t>
    static result_type transform(const JMultiMapTransformer<N, JArgument_t>&                     transformer,
                                 typename JMultiMapTransformer<N, JArgument_t>::const_array_type array,
                                 const result_type&                                              result) 
    {
      typedef JArgument_t argument_type;

      const argument_type z = transformer.getXn(array, 1.0) - transformer.getXn(array, 0.0);

      result_type value(result);

      value.fp /= z;
      value.v  *= z;
      value.V  *= z;

      return value;
    }
  };
}

#endif